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.