Ajustando as médias

Como o efeito da Covariável foi significativo, isso implica na correção das médias dos tratamentos.

Com o comando lm pode-se obter o valor de $\beta$ da covariável diâmetro para correção das médias.

beta<-lm(resp~cov+maq)
summary(beta)

> beta$coef[2]
      cov 
0.9539877

Depois, pode-se obter o valor da média corrigida, por exemplo da máquina 1:

media.resp<-tapply(resp,maq,mean)
media.cov<-tapply(cov,maq,mean)
media.cor1<-media.resp[1]-beta$coef[2]*(media.cov[1]-mean(cov))