The Author Online Book Forums are Moving

The Author Online Book Forums will soon redirect to Manning's liveBook and liveVideo. All book forum content will migrate to liveBook's discussion forum and all video forum content will migrate to liveVideo. Log in to liveBook or liveVideo with your Manning credentials to join the discussion!

Thank you for your engagement in the AoF over the years! We look forward to offering you a more enhanced forum experience.

dksamuel (6) [Avatar] Offline
#1
Sir, I wanted to compute Fisher's least significant difference (LSD) method
I used the following code to read the table

library(RODBC)
> library(agricolae)
> channel <- odbcConnectExcel("G:/oneway.xls"smilie
> mydata <- sqlFetch(channel,"Sheet1"smilie
> library(reshape)
> mydata.m <- melt(mydata, id=c("Treatment"smilie)
> fit <- aov(value~Treatment, data=mydata.m)
> summary(fit)
Df Sum Sq Mean Sq F value Pr(>F)
Treatment 9 23494.8 2610.54 1589.5 < 2.2e-16 ***
Residuals 23 37.8 1.64
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

then I used the following code as given in the following sites

http://cran.r-project.org/web/packages/agricolae/index.html
http://tarwi.lamolina.edu.pe/~fmendiburu/


library(agricolae)
LSD.test()
comparison <- LSD.test(Treatment,Residuals,df,MSerror,group=F)


but I could not proceed further,
Sir can you please include the LSD and information on comparing means after ANOVA

Sir, if you have included the needed information in chapter 8 please excuse me, as I have not got my copy yet. Samuel
robert.kabacoff (170) [Avatar] Offline
#2
Re: Fisher's least significant difference (LSD) method in R
I am not familiar with the agricolae package, but will certainly take a look.

From a cursory look at your code, I am wondering if you need to specify the dependent (outcome) variable in the line

comparison <- LSD.test(value, Treatment,Residuals,df,MSerror,group=F)

Sincerely,

Rob