[R-br] Extrair um vetor de diversos ajustes

Gustavo Marcatti vgp.gustavo em yahoo.com.br
Quarta Maio 4 10:55:38 BRT 2011


Bom dia a todos!
Quando desejo extrair os coeficientes de diversos ajustes, por exemplo, procedo 
da seguinte maneira:
x<-1:60
y<-rep(sort(abs(rnorm(20))),3)
caso <- rep(1:3,each = 20)
dados <- data.frame(caso,x,y)
exp1 <- lapply(split(dados,dados["caso"]),FUN=function(data){
            coefs <- coef(lm(y ~ x, data=data))
            unlist(coefs)
        })
do.call("rbind",exp1)

Minha dúvida é como faço para extrair um vetor de varios casos, assim como no 
exemplo abaixo:
exp2 <- lapply(split(dados,dados["caso"]),FUN=function(data){
            preditos <- fitted(lm(y ~ x, data=data))
            list(preditos)
        })
do.call("rbind",exp2)

O resultado esperado seria um único vetor de tamanho 60 x 1 contendo todos os 
valores preditos para cada ajuste.

Obs1: Esse exemplo é apenas para ilustrar, não é o que realmente estou 
precisando.
Obs2: Os casos tem tamanhos variados, não exatamente 20 observações.

Gustavo Marcatti
Eng. Florestal 
UFV
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110504/dc279ff1/attachment.html>


Mais detalhes sobre a lista de discussão R-br