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

[Jose Claudio Faria]

Pessoal,

O CMR abaixo é da apostila do Walmes.

#------------------------------------------------------------------------------------------
# dados
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)
#
#------------------------------------------------------------------------------------------
# análise de variância (erro A = BL x AD)
m0 <- aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data=ps)   #
termo Error para declarar erro A
m0 <- aov(alt ~ BL + AD * ES + Error(BL:AD), data=ps)   # o mesmo com
fórmula comprimida

summary(m0)
#------------------------------------------------------------------------------------------

Encontrei muitos outros exemplos (na web e em livros) com essa
notação: Error(BL:AD)

Contudo o R (3.0.2patched) não está aceitando essa notação!
Fornece sempre a mensagem de erro:
In aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data = ps) :
  modelo Error() é singular

Contudo o R está aceitando esse modelo:
m0 <- aov(alt ~ BL + AD + Error(BL/AD) + ES + AD:ES, data=ps)

Mas nesse caso o resíduo A (associado a bloco e parcela - AD) não fica correto.

Alguém que anda trabalhando com esses modelos poderia ajudar?

Ab,
-- 
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
Jose Claudio Faria
Estatistica
UESC/DCET/Brasil
joseclaudio.faria at gmail.com
Telefones:
55(73)3680.5545 - UESC
55(73)9100.7351 - TIM
55(73)8817.6159 - OI
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\