coag_c(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)))

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)

coag.aov_aov(coag~diet,coag.df)

summary(coag.aov)

fitted.values(coag.aov)

resid(coag.aov)

par(mfrow=c(1,2))

hist(resid(coag.aov))

qqnorm(resid(coag.aov))

par(mfrow=c(1,1))

plot(fitted(coag.aov),resid(coag.aov))

summary(coag.aov)[2,3]

ks.gof(resid(coag.aov),dist="norm",mean=0,sd=sqrt(summary(coag.aov)[2,3]))

model.tables(coag.aov)

model.tables(coag.aov,type="means")

multicomp(coag.aov,plot=T)