<div><span style="font-size: medium; color: #0000ff;">Tem este arquivo 'dados' para rodar?</span></div>
<div> </div>
<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</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="../../../undefined//compose?to=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"> </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.ufpr.br/Rpira/Rpira/node17.html</a>>) </div>
<div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"> </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"> </div>
<div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Desde já agradeço</div>
<div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr"> </div>
<div id="m_1795547057813537629yui_3_16_0_ym19_1_1478277557220_5220" dir="ltr">Luiz</div>
</div>
</div>
<br />_______________________________________________<br /> R-br mailing list<br /> <a href="../../../undefined//compose?to=R-br@listas.c3sl.ufpr.br" target="_blank">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/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-guia</a>) e forneça código mínimo reproduzível.</blockquote>
</div>
<br /><br /><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>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br />Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>