[R] LSD comparison
arun
smartpink111 at yahoo.com
Thu Jun 27 20:05:05 CEST 2013
Hi,Without any reproducible example or code, it is difficult to understand the problem. Did you load the package?
Using the example in ?LSD.test()
library(agricolae)
data(sweetpotato)
model<-aov(yield~virus, data=sweetpotato)
out <- LSD.test(model,"virus", p.adj="bonferroni",
main="Yield of sweetpotato\ndealt with different virus")
Study: Yield of sweetpotato
dealt with different virus
LSD t Test for yield
P value adjustment method: bonferroni
Mean Square Error: 22.48917
virus, means and individual ( 95 %) CI
yield std.err r LCL UCL Min. Max.
cc 24.40000 2.084067 3 19.594134 29.20587 21.7 28.5
fc 12.86667 1.246774 3 9.991602 15.74173 10.6 14.9
ff 36.33333 4.233727 3 26.570341 46.09633 28.0 41.8
oo 36.90000 2.482606 3 31.175100 42.62490 32.1 40.4
alpha: 0.05 ; Df Error: 8
Critical Value of t: 3.478879
Least Significant Difference 13.4704
Means with the same letter are not significantly different.
Groups, Treatments and means
a oo 36.9
a ff 36.33
ab cc 24.4
b fc 12.87
sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8
[5] LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_CA.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] agricolae_1.1-4 stringr_0.6.2 reshape2_1.2.2
loaded via a namespace (and not attached):
[1] plyr_1.8 tcltk_3.0.1 tools_3.0.1
A.K.
Dear Sir
I'm a student in the university and just start to use R for agricultural research analysis.
I know my question is an old question, but I don't find any clear explanation about my problem via internet.
I'm sorry if my email is disturbing all of you.
I install R with version 3.0.1 and no problem with anova analysis,
my problem is I can't perfprm LSD comparison, it is appearing "Error:
could not find function "LSD.test"
I think I'm already installed "agricolae package" on my system. I attach some picture as shown on my screen.
if possible please guide me how to solve the problem. (how to check if agricolae package is completely installed?)
Thanks you very much for reading
More information about the R-help
mailing list