Bom dia,este é um CMR hipotético.
dados <- expand.grid(A=factor(paste("A", 1:3, sep="")),B=factor(paste("B", 1:4, sep="")), rep=1:4)dados$y <- rgamma(n=48,shape=2, scale = 2)
summary(lm(y~rep, data = subset(dados,B=="B1")))plot(y ~resp) summary(lm(y~rep, data = subset(dados,B=="B2")))plot(y ~resp)
summary(lm(y~rep, data = subset(dados,B=="B3")))plot(y ~resp)
summary(lm(y~rep, data = subset(dados,B=="B4")))
Como posso automatizar estas regressões com um for? Alguma sugestão? Sem plot já está ótimo!
Obrigado
André Oliveira Souza