library(agricolae)
data(sweetpotato)
attach(sweetpotato)
model<-aov(yield~virus)
df<-df.residual(model)
MSerror<-deviance(model)/df
Fc<-anova(model)[1,4]
comparison <- waller.test(yield,virus,df, MSerror, Fc,
main="Yield of sweetpotato
dealt with different virus")
# std = F (default) is standard error
#startgraph
par(mfrow=c(2,2))
par(cex=1)
bar.err(comparison,horiz=TRUE,xlim=c(0,45),angle=125,density=6,
main="Standard deviation")
bar.err(comparison,std=FALSE,horiz=TRUE,xlim=c(0,45),density=8,
col="brown",main="Standard error")
bar.err(comparison,ylim=c(0,45),density=4,angle=125,col="green",
main="Standard deviation")
bar.err(comparison,std=FALSE,ylim=c(0,45),density=10,col="blue",
main="Standard error")
#endgraph
Run the code above in your browser using DataLab