<div dir="ltr"><div>Complementando minha resposta.<br><br>dados <- read.table("dados.csv", header = TRUE, sep = ";", dec = ",")<br>plot(dados$y ~ dados$x)<br><br># Veja que o A33 nao tem variabilidade vou remover<br>dados <- dados[which(dados$x != "A33"),]<br>dados$x <- droplevels(dados$x)<br>tapply(dados$y, dados$x, sd)<br><br># Aparentemente o A21, A69 e A95 tem variancias diferentes<br><br># Modelo de regressao<br>fit1 <- lm(y ~ x, data = dados)<br>plot(residuals(fit1)~ fitted(fit1)) # Parece heterocedastico<br><br># Modelo de regressao media e variancia<br></div># A ligação pra média é a identidade e pra variance é o log.<br><div>require(mcglm)<br>dados$id <- 1<br>fit2 <- mcglm(c(y ~ x), list(mc_dglm(~ x, id = "id", data = dados)), <br>              covariance = "expm", data = dados)<br>summary(fit2)<br><br># Agora o modelo é heterocedastico e não precisa mais preocupar com isso.<br><br><br><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 4 de novembro de 2016 18:06, Luiz Leal 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> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220">Prezados</div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Nos casos em que a transformação box-cox não é capaz remover a heterocedasticidade há alguma outra solução?</div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"><br></div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Eu utilizei a sugestão do Leg/ UFPR nos meus dados, mas não funcionou.</div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">(exemplo: <<a href="http://www.leg.ufpr.br/Rpira/Rpira/node17.html" target="_blank">http://www.leg.<wbr>ufpr.br/Rpira/Rpira/node17.<wbr>html</a>>) </div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"><br></div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Segue em anexo o script e os dados.</div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"><br></div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Desde já agradeço</div><span class="HOEnZb"><font color="#888888"><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"><br></div><div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Luiz</div></font></span></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><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Wagner Hugo Bonat<br>----------------------------------------------------------------------------------------------<br>Department of Mathematics and Computer Science (IMADA)<br>University of Southern Denmark (SDU) and<br>Laboratório de Estatística e Geoinformação (LEG)<br>Universidade Federal do Paraná (UFPR)<br><br></div></div></div></div></div>
</div>