How can this be achieved in 2-way or 3 way experiments? On Monday, January 18, 2016 7:31 AM, Jianling Fan <fanjianl...@gmail.com> wrote:
Thanks David, That's exactly what I need! Thanks very much for your example. Best regards, Julian On 18 January 2016 at 13:00, David L Carlson <dcarl...@tamu.edu> wrote: > Providing us with a reproducible example, makes it much easier to answer your > question. The object returned by LSD.test() is completely different from the > object described in the "Value" section of its manual page which makes things > confusing. If you look at the first example on the manual page: > >> library(agricolae) >> data(sweetpotato) >> model<-aov(yield~virus, data=sweetpotato) >> out <- LSD.test(model,"virus", p.adj="bonferroni") >> out > $statistics > Mean CV MSerror LSD > 27.625 17.1666 22.48917 13.4704 > > $parameters > Df ntr bonferroni alpha test name.t > 8 4 3.478879 0.05 bonferroni virus > > $means > yield std r LCL UCL Min Max > cc 24.40000 3.609709 3 18.086268 30.71373 21.7 28.5 > fc 12.86667 2.159475 3 6.552935 19.18040 10.6 14.9 > ff 36.33333 7.333030 3 30.019601 42.64707 28.0 41.8 > oo 36.90000 4.300000 3 30.586268 43.21373 32.1 40.4 > > $comparison > NULL > > $groups > trt means M > 1 oo 36.90000 a > 2 ff 36.33333 a > 3 cc 24.40000 ab > 4 fc 12.86667 b > > So you can get the LSD statistic with > >> out$statistics$LSD > [1] 13.4704 >> out$statistics[4] > LSD > 13.4704 >> out[[1]][4] > LSD > 13.4704 > > ------------------------------------- > David L Carlson > Department of Anthropology > Texas A&M University > College Station, TX 77840-4352 > > > -----Original Message----- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Jianling Fan > Sent: Monday, January 18, 2016 10:47 AM > To: r-help > Subject: [R] LSD value > > Hello, everyone, > > I am using LSD.test() from package "agricolae" to do my anova > analysis. I know I can calculate LSD value by its equation > t*sqrt(MSE*2/n), but I am wondering if there is a code or something > that can give the LSD value more directly in R? > > Thanks! > > Julian > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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 -- To UNSUBSCRIBE and more, see 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. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.