<div dir="ltr">Maurício,<div><br></div><div>Toda vez que você tem "comparação de médias" e "pressupostos" é necessário pensar no objetivo da comparação:</div><div><br></div><div><ul><li>Você está usando um teste de médias para testar parametricamente se as distribuições podem ser consideradas iguais ou não, ou;<br></li><li>você está querendo testar as médias e fazer uso o CLT para que o teste tenha robustez estatística?<br></li></ul></div><div>Se você pensa em usar uma transformação (você citou Box-Cox) então a nova pergunta é você quer testar as distribuições, e aí os testes teriam que ser nos parâmetros transformada B-C (que além de uma medida de tendência central e dispersão ter um parâmetro a mais) e quais os pressupostos para testes dessa(s) distribuição(ões), ou se você continua interessado somente nas médias: por que transformar?</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-10-26 22:59 GMT-02:00 Maurício Lordêlo via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Caros,<div>Preciso fazer um teste de comparação de médias mas estou tendo problemas com os pressupostos. Tentei uma transformação Box Cox mas não consegui êxito. Alguma sugestão?</div><div><br></div><div><div>Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681)</div><div>Cor = factor(c(rep("ambar",6),rep("<wbr>ambar_claro",3),rep("ambar",6)<wbr>,rep("ambar_claro",6),rep("<wbr>branco",6),</div><div> rep("extra_ambar_claro",3),<wbr>rep("branco",3))) </div><div>Fenois; Cor</div><div>boxplot(Fenois~Cor,ylab=<wbr>expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro"))</div><div>tapply(Fenois, Cor, mean)</div><div>tapply(Fenois, Cor, sd)</div><div>tapply(Fenois, Cor, var)<br></div><div><br></div><div>#Teste de homogeneidade de variâncias</div><div>bartlett.test(Fenois ~ Cor)</div><div>#Gráficos para avaliar a variável resposta</div><div>qqnorm(Fenois)</div><div>qqline(Fenois)</div><div>hist(Fenois)</div><div>boxplot(Fenois)</div><div><br></div><div>#Transformação Box-Cox</div><div>require(MASS)</div><div>boxcox(Fenois~Cor, lambda = seq(-2, 2, 1/10)) </div><div>boxcox(Fenois~Cor, lambda = seq(0.5, 1.5, 1/100))<br></div><div>#Optando por lambda = 0.6</div><div>Fenois2 = ((Fenois^0.6) - 1)/0.6</div><div>par(mfrow=c(1,2))<br></div><div>boxplot(Fenois~Cor,ylab=<wbr>expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro"))</div><div>boxplot(Fenois2~Cor,ylab=<wbr>expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro"))</div><div><br></div><div>bartlett.test(Fenois2 ~ Cor)</div></div><div><br></div></div>
<br>______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>