[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