e esqueci de dar um exemplo mais tangivel das mudancas que esperam...
set.seed(1) N <- 100 x <- rnorm(N) e <- rnorm(100, 0, (x^2+1)^2) y <- 10-2*x+e plot(y~x) mod0 <- lm(y~x) mod1 <- lm(y~x, weights=1/(x^2+1)^2) summary(mod0) summary(mod1) b _______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
