Hi, Tried with your example dataset: Are you able to get results with the example data?
dat1<-read.table(text=" subspecies WMF 1 rowleyi 2.50 2 rowleyi 2.30 3 rowleyi 2.35 49 beatae 2.20 50 beatae 2.35 51 beatae 2.45 91 levipes 2.45 92 levipes 2.35 93 levipes 2.50 122 carletoniA 2.20 123 carletoniA 2.10 124 carletoniA 2.30 135 carletoniB 2.60 136 carletoniB 2.40 137 carletoniB 2.60 145 schmidlyi 2.50 146 schmidlyi 2.70 147 schmidlyi 2.70 161 simulus 2.40 162 simulus 2.20 163 simulus 2.30 ",sep="",header=TRUE,stringsAsFactors=TRUE) library(agricolae) fm1<-aov(WMF~subspecies,data=dat1) duncan.test(fm1,"subspecies",14,0.01131,0.05) #Study: #Duncan's new multiple range test #for WMF # #Mean Square Error: 0.01130952 #subspecies, means # # WMF std.err replication #beatae 2.333333 0.07264832 3 #carletoniA 2.200000 0.05773503 3 #carletoniB 2.533333 0.06666667 3 #levipes 2.433333 0.04409586 3 #rowleyi 2.383333 0.06009252 3 #schmidlyi 2.633333 0.06666667 3 #simulus 2.300000 0.05773503 3 #alpha: 0.05 ; Df Error: 14 #Critical Range # 2 3 4 5 6 7 #0.1862347 0.1951445 0.2006433 0.2043603 0.2070022 0.2089363 # #Means with the same letter are not significantly different. #Groups, Treatments and means #a schmidlyi 2.63333333333333 #ab carletoniB 2.53333333333333 # bc levipes 2.43333333333333 # bcd rowleyi 2.38333333333333 # bcd beatae 2.33333333333333 # cd simulus 2.3 d carletoniA 2.2 A.K. ----- Original Message ----- From: Amanda Jones <akjone...@gmail.com> To: r-help@r-project.org Cc: Sent: Sunday, November 25, 2012 12:41 AM Subject: [R] Multiple Range Means Test Hello, My boss wants me to do a Duncan's test, which is under the agricolae package. Unfortunately I am not versed enough in R to run my data. I have 7 subspecies of deer mouse for which I have 23 measurements which are my variables of interest. I have run an ANOVA for each of the set of subspecies and variables, my data look like this: subspecies WMF 1 rowleyi 2.50 2 rowleyi 2.30 3 rowleyi 2.35 49 beatae 2.20 50 beatae 2.35 51 beatae 2.45 91 levipes 2.45 92 levipes 2.35 93 levipes 2.50 122 carletoniA 2.20 123 carletoniA 2.10 124 carletoniA 2.30 135 carletoniB 2.60 136 carletoniB 2.40 137 carletoniB 2.60 145 schmidlyi 2.50 146 schmidlyi 2.70 147 schmidlyi 2.70 161 simulus 2.40 162 simulus 2.20 163 simulus 2.30 This is condensed. My data is called dat. The help for duncan.test shows the argument like this: duncan.test(y, trt, DFerror, MSerror, alpha = 0.05, group=TRUE, main = NULL) When I try to run it as duncan.test(aov, dat, 6, .297, alpha=.05, group=TRUE, main=NULL) it gives me an error message that I have differing rows, 0 and 163. 163 is the correct number, where is it getting the 0 from? Thank you, Amanda Jones ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
