[R-br] Fatorial tiplo com dados adicionais

[Tiago Souza Marçal]

Te indico o pacote ExpDes se você quer uma rotina simples.

Att.

Tiago.

Alisson Lucrécio <alisson.lucrecio em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150509/7a348621/attachment.html>
-------------- Próxima Parte ----------
_______________________________________________
R-br mailing list
R-br em 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.