Te indico o pacote ExpDes se você quer uma rotina simples. Att. Tiago. Alisson Lucrécio <alisson.lucrecio@ifgoiano.edu.br> escreveu: Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 níveis de solo e 1 nível de solução (Água). Segue um exemplo abaixo. fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, 0.50))) Rep <- apply(data.frame(Solo = rep(1:5, each = 18), Sol = rep(rep(c(0.5, 1), each = 9), 5), Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), 1, sum) fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), Sol = rep("Água", 15)) Rep <- rep(1:5, each = 3) adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) Muito obrigado. -- Alisson Lucrecio da Costa _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.