<div><span style="font-size: medium; color: #0000ff;">Ele não reconhece estes comandos 'mcglm','gof', tem alguma bibliotéca para isto?</span></div>
<div> </div>
<div dir="ltr">
<div>Caros, </div>
<div><br />Alguém postou esse conjunto de dados com problema de pressupostos, principalmente heterocedasticidade. Agora a pouco veio outro e-mail com um problema similar. Fiz um exemplo um pouco mais detalhado de como isso pode</div>
ser facilmente resolvido e mostrando o efeito disso no modelo.<br /><br /># Example 2 ------------------------------------------------------------<br />Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, <br />           328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, <br />           342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, <br />           88.571, 88.571,  90.252,  41.513,  52.437,  49.076,  88.571,  <br />           88.571,  90.252,  64.202,  60.000,  61.681)<br />Cor = factor(c(rep("ambar",6), rep("ambar_claro",3), rep("ambar",6), <br />               rep("ambar_claro",6),rep("branco",6), <br />               rep("extra_ambar_claro",3),rep("branco",3)))<br /><br /># Exploratory analysis<br />boxplot(Fenois ~ Cor)<br />tapply(Fenois, Cor, sd)<br />dados <- data.frame(Fenois, Cor)<br />dados$id <- 1<br /><br /># Fitting <br />fit1 <- mcglm(c(Fenois ~ Cor), list(mc_id(dados)), data = dados)<br />fit2 <- mcglm(c(Fenois ~ Cor), list(mc_dglm(~ Cor, id = "id", data = dados)),<br />              covariance = "expm", data = dados)<br /># Goodness-of-fit <br />gof(fit1)<br />gof(fit2)<br /><br /># Comparing estimates and standard errors<br />coef(fit1, type = "beta", std.error = TRUE)<br />coef(fit2, type = "beta", std.error = TRUE)<br /><br />O interessante é que a estimativa pontual é exatamente a mesma, porém<br />olha a enorme diferença nos erros padrões dos betas.<br />
<div>
<div><br />-- <br />
<div class="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>
</div>
<div>_______________________________________________<br /><br /></div>