coag <- scan()
62 60 63 59
63 67 71 64 65 66
68 66 71 67 68 68
56 62 60 61 63 64 63 59

coag
diet <- factor(rep(LETTERS[1:4],c(4,6,6,8)))  #create a factor
diet
coag.df <- data.frame(diet,coag)
coag.df

par(mfrow=c(1,2))
plot.design(coag.df)
plot.design(coag.df, fun= median)
par(mfrow=c(1,1))
plot.factor(coag.df)

aov.coag <- aov(coag ~ diet, coag.df)
aov.coag
summary(aov.coag)

fitted.values(aov.coag)
par(mfrow=c(3.1))
hist(resid(aov.coag))
qqnorm(resid(aov.coag))
qqline(resid(aov.coag))
plot(fitted(aov.coag), resid(aov.coag))
abline(h=0)

mca.coag <- TukeyHSD(aov.coag,"diet")
mca.coag