box<-rep(1:4,c(2,3,3,2))
box<-as.factor(box)
box
sales<-c(12,18,14,12,13,19,17,21,24,30)
BS.df_data.frame(sales,box)
BS.df


BS.aov_aov(sales~box,BS.df)
summary(BS.aov)
model.tables(BS.aov)
model.tables(BS.aov,type="means")

par(mfrow=c(2,2))
plot.design(BS.df)
hist(resid(BS.aov))
qqnorm(resid(BS.aov))
qqline(resid(BS.aov))
plot(fitted(BS.aov),resid(BS.aov))

multicomp.lm(BS.aov, method="tukey",plot=T)
multicomp.lm(BS.aov, method="scheffe",plot=T)
multicomp.lm(BS.aov, method="bon",plot=T)
multicomp.lm(BS.aov, method="best",plot=T)  
par(mfrow=c(1,1))

food<-rep(c("A","B","C","D","E"),c(7,9,10,12,11))
food
weig<-c(28,24.8,27.8,24.8,28.7,34.8,30.8,
24.7,30.7,30.5,22,26.7,28.7,28.7,22.7,34.7,
30.5,24.8,30.7,22.7,37.7,28.5,28.7,28.7,33.5,24,
30.7,24.7,24.8,30.7,20.7,30.8,30.5,22.5,22.8,21.7,28.7,38.7,
28.7,30.8,22.7,27.8,23.8,18.7,24.7,26.7,25,31,25)
WF.df_data.frame(weig,food)
WF.df


WF.aov_aov(weig~food,WF.df)
summary(BS.aov)
model.tables(WF.aov)
model.tables(WF.aov,type="means")

par(mfrow=c(2,2))
plot.design(WF.df)
hist(resid(WF.aov))
qqnorm(resid(WF.aov))
qqline(resid(WF.aov))
plot(fitted(WF.aov),resid(WF.aov))

multicomp.lm(WF.aov, method="tukey",plot=T)
multicomp.lm(WF.aov, method="scheffe",plot=T)
multicomp.lm(WF.aov, method="bon",plot=T)
multicomp.lm(WF.aov, method="best",plot=T)  
par(mfrow=c(1,1))