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.
