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.

Reply via email to