<div dir="ltr"><div>Caros,</div><div>Peço auxilio para um problema de um fatorial em que há um problema em atender aos pressupostos. <br></div><div><br></div><div>Agradeço alguma sugestão.</div><div><br></div><div>Maurício</div><div><br></div><div>m_seca_raiz = c(0.000000000, 0.000000000, 0.000460000, 0.000000000, 0.000000000, 0.000140000,<br>                0.000000000, 0.000300000, 0.000150000, 0.000000000, 0.000120000, 0.000500000,<br>                0.000560000, 0.000300000, 0.000000000, 0.000160000, 0.000160000, 0.000600000,<br>                0.000200000, 0.000200000, 0.000440000, 0.000000000, 0.000080000, 0.000240000,<br>                0.000380000, 0.000080000, 0.000000000, 0.000250000, 0.000960000, 0.000800000,<br>                0.000320000, 0.000760000, 0.000920000, 0.000366667, 0.000720000, 0.000540000,<br>                0.000440000, 0.000620000, 0.000100000, 0.000600000, 0.000200000, 0.000280000,<br>                0.000000000, 0.000520000, 0.000480000, 0.000000000, 0.000280000, 0.000280000,<br>                0.000800000, 0.000280000, 0.000000000, 0.000680000, 0.000700000, 0.000380000,<br>                0.000000000, 0.000360000, 0.000300000, 0.000340000, 0.000280000, 0.000360000,<br>                0.000340000, 0.001780000, 0.000250000, 0.000260000, 0.000300000, 0.000280000,<br>                0.000160000, 0.000120000, 0.000125000, 0.000080000, 0.000220000, 0.000320000)<br>meio = factor((rep(c("m2","m3"),each=12, times=3)))<br>sacarose = factor(rep(c("7,5","15,0","30,0"), each=24))<br>carvão = factor((rep(c("c1","c2"),each=6, times=6)))<br>tratamento = meio:sacarose:carvão<br>head(data.frame(m_seca_raiz, meio, sacarose, carvão, tratamento),18)<br><br><br>hist(m_seca_raiz)<br>boxplot(m_seca_raiz ~ tratamento,ylab="Massa seca da raiz", xlab="Tratamento", <br>        cex.axis=0.85 )<br><br>max(m_seca_raiz)<br><br>modelo_m_seca_raiz = aov(m_seca_raiz ~ tratamento)<br>summary(modelo_m_seca_raiz)<br><br>hist(modelo_m_seca_raiz$residuals)<br>shapiro.test(modelo_m_seca_raiz$residuals)<br>require(car)<br>leveneTest(modelo_m_seca_raiz$residuals ~ tratamento)<br><br>anova(modelo_m_seca_raiz)$'Mean Sq'[2]<br>res_pad_m_seca_raiz = modelo_m_seca_raiz$residuals/sqrt(anova(modelo_m_seca_raiz)$'Mean Sq'[2])<br>plot(modelo_m_seca_raiz$fitted.values,  res_pad_m_seca_raiz, ylim = c(-5,5))<br>abline(h=c(-3,3),lty=3)<br><br>#Elimando a observação 62<br>shapiro.test(modelo_m_seca_raiz$residuals[-62])<br>leveneTest(modelo_m_seca_raiz$residuals[-62] ~ tratamento[-62])<br><br>#Transformação Box-Cox <br>require(forecast)<br>#require(fpp)<br>BoxCox.lambda(m_seca_raiz+0.5, method = "guerrero", lower=-3,upper = 3)<br>BoxCox.lambda(m_seca_raiz+0.5, method = "loglik", lower=-3,upper = 3)<br>#lambda = BoxCox.lambda(m_seca_raiz,lower=-3)<br>#lambda<br>m_seca_raizt = BoxCox(m_seca_raiz+0.5, lambda=-3)<br><br>hist(m_seca_raizt)<br>boxplot(m_seca_raizt ~ tratamento, ylab="m_seca_raiz (transformada)", xlab="Tratamentos")<br>modelo_m_seca_raizt = aov(m_seca_raizt ~ tratamento)<br>summary(modelo_m_seca_raizt)<br>hist(modelo_m_seca_raizt$residuals)<br>plot(modelo_m_seca_raizt$residuals)<br>plot(modelo_m_seca_raizt$fitted.values,  modelo_m_seca_raizt$residuals)<br>shapiro.test(modelo_m_seca_raizt$residuals)<br>leveneTest(modelo_m_seca_raizt$residuals ~tratamento)<br><br><br>anova(modelo_m_seca_raizt)$'Mean Sq'[2]<br>res_pad_m_seca_raizt = modelo_m_seca_raizt$residuals/sqrt(anova(modelo_m_seca_raizt)$'Mean Sq'[2])<br>plot(modelo_m_seca_raizt$fitted.values,  res_pad_m_seca_raizt, ylim = c(-6,6))<br>abline(h=c(-3,3),lty=3)<br></div></div>