<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:10pt"><DIV>Bom dia a todos!</DIV>
<DIV>Quando desejo extrair os coeficientes de diversos ajustes, por exemplo, procedo da seguinte maneira:</DIV>
<DIV>x<-1:60<BR>y<-rep(sort(abs(rnorm(20))),3)<BR>caso <- rep(1:3,each = 20)<BR>dados <- data.frame(caso,x,y)</DIV>
<DIV>exp1 <- lapply(split(dados,dados["caso"]),FUN=function(data){<BR> coefs <- coef(lm(y ~ x, data=data))<BR> unlist(coefs)<BR> })</DIV>
<DIV>do.call("rbind",exp1)</DIV>
<DIV> </DIV>
<DIV>Minha dúvida é como faço para extrair um vetor de varios casos, assim como no exemplo abaixo:</DIV>
<DIV>exp2 <- lapply(split(dados,dados["caso"]),FUN=function(data){<BR> preditos <- fitted(lm(y ~ x, data=data))<BR> list(preditos)<BR> })</DIV>
<DIV>do.call("rbind",exp2)</DIV>
<DIV> </DIV>
<DIV>O resultado esperado seria um único vetor de tamanho 60 x 1 contendo todos os valores preditos para cada ajuste.</DIV>
<DIV> </DIV>
<DIV>Obs1: Esse exemplo é apenas para ilustrar, não é o que realmente estou precisando.</DIV>
<DIV>Obs2: Os casos tem tamanhos variados, não exatamente 20 observações.</DIV>
<DIV> </DIV>
<DIV>Gustavo Marcatti</DIV>
<DIV>Eng. Florestal </DIV>
<DIV>UFV</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV></div></body></html>