<div dir="ltr">Se é que eu entendi sua dúvida...<div><br></div><div><div>vars=list()</div><div>for(i in 1:30) vars[[i]] <- rnorm(48)</div><div>names(vars) <- paste('var',1:30,sep='')</div><div>vars</div><div><br></div><div>varis <- do.call('cbind',vars)</div><div><br></div><div>da <- expand.grid(FatA=rep(LETTERS[1:6],4),</div><div>                  FatB=letters[1:2])</div><div><br></div><div>dad <- data.frame(da,varis)</div><div><br></div><div>tabs <- apply(dad[,-c(1:2)],</div><div>      2,</div><div>      function(x) aov(x ~ FatA*FatB,dad))</div><div>lapply(tabs,summary) </div></div></div>