Dear Tu, Use library(agricolae) LSD.test(), It is necessary to aov()
comparison <- LSD.test(yield,virus,df,MSerror,group=F) Comparison between treatments means tr.i tr.j diff pvalue 1 1 2 11.5333333 0.0176 2 1 3 11.9333333 0.0151 3 1 4 12.5000000 0.0121 4 2 3 23.4666667 0.0003 5 2 4 24.0333333 0.0003 6 3 4 0.5666667 0.8873 or comparison <- LSD.test(yield,virus,df,MSerror, group=T) Least Significant Difference 8.928965 Means with the same letter are not significantly different. Groups, Treatments and means a oo 36.9 a ff 36.33333 b cc 24.4 c fc 12.86667 Regards, Felipe. -----Original Message----- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]on Behalf Of Chunhao Tu Sent: Tuesday, June 23, 2009 10:41 PM To: r-help@r-project.org Subject: [R] Help LSD multiple comparison test Hi R users, I know this is a old question but I did not see any clear answer. Is there any way or package to perform LSD test in R such as H0: A=B. H0, A=C and H0 B=C. Thank you in advance Tu -- View this message in context: http://www.nabble.com/Help-LSD-multiple-comparison-test-tp24178161p24178161.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.
