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