[R-br] Ajuste de regressão linear

[Paulo Nogueira Starzynski]

André,
rapidamente, sem o plot e também sem testar:

dados <- expand.grid(A=factor(paste("A", 1:3, sep="")),B=factor(paste("B",
1:4, sep="")), rep=1:4)
dados$y  <- rgamma(n=48,shape=2, scale = 2)
lapply(paste("B", 1:4, sep=""), function(x){summary(lm(y~rep, data =
subset(dados,B==x)))})


Se funcionar, é só colocar as linhas do plot dentro da função.

abs
Paulo

Atenciosamente,
Paulo Nogueira Starzynski

2015-07-02 8:44 GMT-03:00 Diego Miro <d.miro1089 em gmail.com>:

> André, mão está otimizado mas já te dá um norte.
>
> Estou assumindo que o resp do plot é o rep do objeto dados
>
> b <- sort(unique(dados$B)) # ou levels(dados$B)
> modelos.summary <- vector(mode = 'list', length = length(b))
> for (i in 1:length(b)) {
>   dados.modelo <- subset(x = dados, subset = B == b[i])
>   modelos.summary[[i]] <- summary(lm(y ~ rep, data = dados.modelo))
>   png(paste0("Gráfico ", b[i], ".png"))
>   plot(dados.modelo$y ~ dados.modelo$rep)
>   dev.off()
> }
> names(modelos.summary) <- b
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150702/f317debc/attachment.html>