Re: [R] LSD-test

2017-12-25 Thread Kenneth Roy Cabrera Torres
Check the agricolae package. library(agricolae) El 25/12/17 a las 17:33, David Winsemius escribió: On Dec 25, 2017, at 2:09 PM, Ahmed Attia wrote: The model should be class aov or lm and my model class is aovlist. tried tidy from broom library but did not work. To make it class aov, I had t

Re: [R] LSD-test

2017-12-25 Thread David Winsemius
> On Dec 25, 2017, at 2:09 PM, Ahmed Attia wrote: > > The model should be class aov or lm and my model class is aovlist. > tried tidy from broom library but did not work. To make it class aov, > I had to remove the error term; > > model <- > aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Ni

Re: [R] LSD-test

2017-12-25 Thread Ahmed Attia
The model should be class aov or lm and my model class is aovlist. tried tidy from broom library but did not work. To make it class aov, I had to remove the error term; model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen),data=Rotationdata_R) Ahmed Attia, Ph.D. Agronomist & Soil S

Re: [R] LSD-test

2017-12-25 Thread David Winsemius
> On Dec 25, 2017, at 10:21 AM, Ahmed Attia wrote: > > LSD.test ?LSD.test No documentation for ‘LSD.test’ in specified packages and libraries: you could try ‘??LSD.test’ -- David Winsemius Alameda, CA, USA 'Any technology distinguishable from magic is insufficiently advanced.' -Gehm's Cor

[R] LSD-test

2017-12-25 Thread Ahmed Attia
LSD-test produces error for this code; code <- as.factor(Rotationdata_R$`Rot/code`) #factor in the main Nitrogen <- as.factor(Rotationdata_R$Nitrogen) #factor in the sub Rep <- as.factor(Rotationdata_R$REP) #blocks Year <- as.factor(Rotationdata_R$YEAR) #years model <- aov(Rotationdata_R$`GY(M