[R-br] Parcela subdividida: problema na declaração do erro associado a parcela

[walmes .]

Essa mensagem de erro sempre aparece. Eu suponho que seja algo relacionado
a matriz criada com cada especificação para o termo de erro. No CMR abaixo
você vê o mesmo modelo declarado de duas diferentes formas. No final tem a
matriz para os fatores usados no termo de R e como elas são matrizes
diferentes (em conteúdo e dimensões) eu suponho que a origem do erro esteja
atrelado à isso. Aí tem que abrir a fonte para ver o que de fato tá pegando.

ps <- expand.grid(BL=c("I", "II", "III", "IV"),
                  ES=c("e1", "e2"),
                  AD=c("a1", "a2", "a3"))
ps$alt <- c(58,77,38,52,44,59,30,34,85,90,73,77,
            59,68,45,55,66,93,67,64,54,75,53,48)
str(ps)

## modelo para experimento de parcelas subsdividas
m0 <- aov(alt~BL+AD*ES+Error(BL:AD), data=ps)

## criando o fator parcela fora da aov
ps$parcela <- with(ps, interaction(BL, AD))

## mesmo modelo
m1 <- aov(alt~BL+AD*ES+Error(parcela), data=ps)

## a mesma anova
summary(m0)
summary(m1)

a <- model.matrix(~BL:AD, ps)
dim(a)
b <- model.matrix(~parcela, ps)
dim(b)

À disposição.
Walmes.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131204/dc1861f2/attachment.html>