similar to: Kruskal-Wallis, Friedman tests and Tukey HSD

Dear all I have a dataset that looks like: fr.wt site 1 4400 glen 2 235 glen 3 225 glen ' ' ' ' ' ' ' ' ' 82 550 glen 83 550 kom 84 550 kom ' ' ' ' ' ' ' ' ' 191 820 kom 192 2000 soet ' ' ' ' ' ' I need to do a series of histograms for each of the codes, levels or factors in

Dear all I have made some lattice graphs (bwplot(dead ~ treat | group,...). 'group' is a factor with three levels (Artemia larvae, Abalone larvae and Abalone spat), and the result is a graph with three panels corresponding to the levels in 'group', ordered in alphabetical order from bottom to top, as expected. How does one re-order the order in which the levels of

I followed the threads that enquired about doing post-hoc tests after doing Kruskal testing. It took me to npmc. But npmc is giving an output I do not understand. I noticed a thread entitled "npmc function: 'x' must be atomic" but there never appeared to be a resolution. > npmc(npmcinput) Error in sort(unique.default(x), na.last = TRUE) : 'x' must be atomic

Dear all, I run a kruskal wallis test and found significant results. Then, I conducted all pairwise comparisons and found no significant results. Could anyone please give me a hint as to why this happens or redirect me towards a specific web page where I can find more info? (I used alpha=5% and made no bonferroni or other correction for the pairwise comparisons) Thank you   Dr. Iasonas

Well fellow R users, I throw myself on your mercy. Help me, the unworthy, satisfy my employer, the ungrateful. My feeble ramblings follow... I've searched R-Help, the R Website and done a GOOGLE without success for a one way ANOVA procedure to analyse data that are both non-normal in nature and which exhibit unequal variances and unequal sample sizes across the 4 treatment levels. My

Dear all I need to obtain an estimate of the 50% lethal dose (LD50) from a logistic regression model obtained by applying the glm procedure to some binomial data. The model appears to fit the data very well. I used dope.p from MASS to try and find LD50. The following output appears: > dose.p(iso.glm.logit, cf = c(1,3), p = 1:3/4) Error in dose.p(iso.glm.logit, cf = c(1, 3), p = 1:3/4) :

Dear all, I would be glad if someone could tell me if Tukey HSD (function of stats library) accept the NA values and could correct the unbalanced design (different number of sample in each group) Thanks for help Alice Le Bars

Hi everyone, I am working on analyzing spectra which requires hundreds of kruskal wallis tests to be performed. Once significant results are found post tests are performed ect. My questions are: how do I perform all these kw tests and save just the pvalues in a table? Is there anyway to have R automatically perform post tests if a result of p<0.05 is found? Currently when I run the test

Hi, I have got two groups of samples; and for every row, I wish to calculate Kruskal-Wallis' p-value. In the example below, and the stars (****) show where I am struggling to design and put things together. Any help would be appreciated. myFile <- data.frame(Sample_1a = 1:10, Sample_1b = 2:11, Sample_1c = 3:12, Sample_2a=4:13, Sample_2b=7:16, row.names=LETTERS[1:10]) groups <- rep

Dear all, I am doing a one way between subjects anova in an unbalanced data set. Suppose we have "a" levels of the one factor. I want to merge the not so significantly different levels into the same cluster. Can I do a Tukey Kramer HSD and then use the following algorithm: For i in 2 : "a" For j in 1 : i-1 if mean of level i is not significantly different

Hello, Is there like the kruskal wallis test in relation to ANOVA (no restrictions on normallity and variance homogenity) something (in R) for MANOVA? thanks -- Dr.Nicolaas Busscher Universit?t GH Kassel Nordbahnhofstrasse: 1a, D-37213 Witzenhausen Phone: 0049-(0)5542-98-1715, Fax: 0049-(0)5542-98-1713

Hi all, The Kruskal-Wallis test is a generalization of the two-sample Mann-Whitney test to *k* samples. That being the case, the Kruskal-Wallis test with *k*=2 should give an identical p-value to the Mann-Whitney test, should it not? x1<-c(1:5) x2<-c(6,8,9,11) a<-wilcox.test(x1,x2,paired=FALSE) b<-kruskal.test(list(x1,x2),paired=FALSE) a$p.value [1] 0.01587302 b$p.value [1]

Hello, I need nemenyi-test or any other post-hoc test for kruskal-wallis, but I just can't find out how to implement this in R. My data set is nitrite concentrations in four different groups which I intend to compare. The kruskal-wallis-test showed significance, but I still don't know between which groups these significances are. Can anybody help me with this? Thanks

Hi, I need to determinate HSD value from a matrix like that Thesis Days A1 Cx 0 66.07 Cx 0 60.24 Cx 0 42.86 Tw 0 66.07 Tw 0 60.24 Tw 0 42.86 Aa 0 66.07 Aa 0 60.24 Aa 0 42.86 Qe 0 66.07 Qe 0 60.24 Qe 0 42.86 Cx 56 310.86 Cx 56 223.17 Cx 56 186.77 Tw 56 149.42 Tw 56 127.75 Tw 56 138.59 Aa 56 130.24 Aa 56 214.83 Aa 56 137.95 Qe 56 186.64 Qe 56 189.09 Qe 56 187.87 this is my script

Hello All, I want "Ranks' Table in R as like in SPSS ouput in the given link. http://www.statisticssolutions.com/methods-chapter/statistical-tests/kruskal-wallis-test/ Is the code is already available? Please let me know. Thanks, Lawrence

Hi, I have searched and found a response to a question similar to mine but when I tried the code, R says it's not an actual function so I thought I'd ask here. http://r.789695.n4.nabble.com/file/n4633035/Cookies.csv Cookies.csv I have attached the data I am using. I am trying to look at two things: how moisture content changes over time, and how it changes along the length of a log

Hi, I was wondering if you could help me: The following are the first few lines of my data set: subject group condition depvar s1 c ver 114.87 s1 c feet 114.87 s1 c body 114.87 s2 c ver 73.54 s2 c feet 64.32 s2 c body 61.39 s3 a ver 114.87 s3 a feet 97.21 s3 a body 103.31 etc. I entered the following ANOVA command: >dat <- read.table("mydata.txt", header=T)

dear all effect of a, b and c on d, total 48 comparisons, got one anova result in model1=aov(d~A*B*C). can we get all the result in one command? or can we interpret the whole comparisons from this result? how it work in Tukey HSD? jose ____________________________________________________________________________________ Access over 1 million songs.

Drug US1 US2 Aptecha Celebrex 235.54 269.99 121.02 Detrol LA 157.99 190.99 55.3 Flomax 166.00 190.99 93.45 Lipitor 174.99 200.99 137.7 Novaldex 108.6 129.99 22.48 Norvasc 186.66 203.99 161.93 Plavix 107.99 106.99 64.53 Prevacid 117.39 134.99 59.83 Prilosec 115.99 126.99 57.75 Zyrtec 181.1 200.99 58.79 US1=c(235.54,157.99,166,174.99,108.60,186.66,107.99,117.39,115.99,181.10)

Hello. We have some questions concerning the statistical analysis of a dataset. We aim to compare the sample means of more than 2 independent samples; the sample sizes are unbalanced. The requirements of normality distribution and variance homogeneity were not met even after transforming the data. Thus we applied a nonparametric test: the Kruskal-Wallis-test (H-Test). The null hypothesis was