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

[Tiago Souza Marçal]

Mesmo com o aviso use o comando summary(m0) e o quadro será formado.
 
Att.
 
Tiago.


################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ################################################################# 
 
> Date: Wed, 4 Dec 2013 18:20:19 -0200
> From: joseclaudio.faria em gmail.com
> To: r-br em listas.c3sl.ufpr.br
> Subject: [R-br] Parcela subdividida: problema na declaração do erro associado a parcela
> 
> 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
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
> _______________________________________________
> 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.
 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131204/98b1c903/attachment.html>