CheFact_list(ConReac=c("L","H"),Catal=c("L","H"))

Che.des_fac.design(c(2,2),CheFact,rep=3)

time_scan()
 28 36 18 31 
 25 32 19 30
 27 32 23 29
 
replic_rep(1:3,rep(4,3))

replic_as.factor(replic)

FacDes22.df_data.frame(time,Che.des,replic)

FacDes22.df

par(mfrow=c(1,2))

plot.factor(FacDes22.df[,1:3])

attach(FacDes22.df)

interaction.plot(ConReac,Catal,time)

interaction.plot(Catal,ConReac,time)

detach()

FacDes22.aov_aov(time~(ConReac+Catal)^2,FacDes22.df)

summary(FacDes22.aov)

model.tables(FacDes22.aov,type="feffects",se=T)

summary.lm(FacDes22.aov)

qqnorm(resid(FacDes22.aov))

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