Pessoal como eu poderia modelar os Ys com um for ou de alguma outra forma? (Tenho 20 Ys)
Exemplo.
library(mvtnorm)
Sigma <- matrix(c(1,1,-.4,-.9,1,4,1.6,-1.2,-.4,1.6,4, 1.8,-.9,-1.2,1.8,9),ncol=4)
Y <- rmvnorm(10, c(2,4,6,8), Sigma)
X1=as.factor(c(rep("A",5),rep("B",5)))
X2=as.factor(c(rep("1",2),rep("3",3),rep("2",2),rep("0",3)))
dados = data.frame(Y,X1,X2)
for (i in 1:4) lm(dados[,i] ~X1*X2))
No verdade desejo usar ajustar com a função!
require(ExpDes.pt)
fat2.dic(X1, X2, dados[i], quali = c(TRUE, TRUE), mcomp = "sk", fac.names = c("Doses", "Substrato"), sigT = 0.05, sigF = 0.05)
obg
André Oliveira Souza.
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES