<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span>Bom dia senhores!</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; ">Gostaria de entender porque os resíduos de um modelo não-ponderado estão sendo os mesmos de um modelo ponderado pelo inverso da variância. Apenas os plots dos modelos estão refletindo as reais mudanças da ponderação. Com certeza, os resíduos padronizados pelo plot do modelo estão diferentes do resíduo que eu estou simulando. O que está errado? Segue um CMR.</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div><div><div>set.seed(1)</div><div>da <- expand.grid(rep=1:20,Geracao=1:5)</div><div>y1 <- rnorm(20,5,1)</div><div>y2
<- rnorm(20,10,4)</div><div>y3 <- rnorm(20,20,7)</div><div>y4 <- rnorm(20,28,12)</div><div>y5 <- rnorm(20,36,28)</div><div>da <- data.frame(da,y=matrix(c(y1,y2,y3,y4,y5),ncol=1)) # </div><div><br></div><div>boxplot(da$y ~ da$Geracao)</div><div><br></div><div>mod <- lm(y ~ factor(Geracao), data=da)</div><div>summary(mod)</div><div>plot(mod)#problema nos resíduos</div><div><br></div><div>#obtendo os resíduos padronizados</div><div>resid_mod <- resid(mod)</div><div>s2_mod <- sum(resid_mod^2)/mod$df.res </div><div>respad_mod <- resid_mod/sqrt(s2_mod)</div><div><br></div><div>library(nortest)</div><div>lillie.test(respad_mod)</div><div><br></div><div>library(car)</div><div>leveneTest(respad_mod ~ factor(Geracao),data=da,center=mean)</div><div><br></div><div>#Fazendo a ponderação</div><div>Var_ger <- rep(tapply(da$y,da$Geracao,var),each=20)</div><div>Var_ger</div><div><br></div><div>da$Var <-
Var_ger</div><div>da</div><div>#Novo ajuste como o modelo ponderado</div><div>mod1 <- lm(y ~ factor(Geracao), weights=1/Var, data=da)</div><div>summary(mod1)</div><div><br></div><div>#comparando os dois ajustes</div><div>par(mfrow=c(2,4))</div><div>plot(mod)</div><div>plot(mod1)</div><div><br></div><div>#obtendo os resíduos padronizados</div><div>resid_mod1 <- resid(mod1)</div><div>s2_mod1 <- sum(resid_mod1^2)/mod1$df.res </div><div>respad_mod1 <- resid_mod1/sqrt(s2_mod1)#são os mesmos do modelo anterior???????????? </div><div><br></div><div>lillie.test(respad_mod1)</div><div><br></div><div>leveneTest(respad_mod1 ~ factor(Geracao),data=da,center=mean)</div><div><br></div><div>qqnorm(resid_mod1)#diferente do plot???????????</div><div><br></div><div>boxplot(resid_mod ~ Geracao,data=da)</div><div>boxplot(resid_mod1 ~ Geracao,data=da)#deveria estar diferente não?????</div><div><br></div><div>Desde já, agradeço a
atenção de todos.</div><div><br></div><div>(S,f,P)</div><div>Allaman</div></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "> </div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\begin{signature}</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2"><<>>=</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Prof. Dr. Ivan
Bezerra Allaman</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Universidade Estadual de Santa Cruz</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Departamento de Ciências Exatas e Tecnológicas</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Ilhéus/BA - Brasil</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; "
align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Fone: +55 73 3680-5076</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com</font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">@</font></font></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; text-align: left; background-color: transparent; " align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif"
size="2">\end{signature}</font></font></div></div></body></html>