De um modo bem menos elegante do que os nossos amigos da programação, pode-se ter: length(levels(factor(dados$parcela))) results <- list() for(i in 1:length(levels(factor(dados$parcela)))){   results[[i]] <- lm(y ~ x, subset(dados, parcela==i)) } results lapply(results,summary) lapply(results,anova)   Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman@yahoo.com.br