##################
#Oneway ANOVA in CRD
########
Germinated=c(56,48,66,62,83,78,94,93,80,72,83,85)
Insecticide=factor(c(rep(1,4),rep(2,4),rep(3,4)))
par(mfrow=c(1,3))
lm1=lm(Germinated~Insecticide)
plot(Germinated~as.numeric(Insecticide),main="Seedlings",xlab="Insecticide", xaxt = "n",ylab="# Germinated per 100")
axis(1, at = c(1,2,3), labels = c('A', 'B', 'C'), tick = TRUE)
anova(lm1)
qqnorm(lm1$resid,main="Residuals")
qqline(lm1$resid)
r=rnorm(length(Germinated),0,7.18)
qqnorm(r,main="N(0,7.18)")
qqline(r)

##################
# Randomized Complete Block Design
#####
Productivity=c(33,38,39,42,62,35,37,43,47,71,40,42,45,52,74,54,50,55,62,84)
Workshop=factor(c(rep('A',5),rep('B',5),rep('C',5),rep('D',5)))
Attitude=factor(c(rep(1:5,4)))
lm1=lm(Productivity ~ Workshop + Attitude)
anova(lm1)
interaction.plot(Workshop,Attitude,Productivity,ylab="Productivity")