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