<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Mesmo com o aviso use o comando summary(m0) e o quadro será formado.<BR> <BR>Att.<BR> <BR>Tiago.<br><br><BR><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div><div> </div><div>################################################################# </div><br> <BR><div>> Date: Wed, 4 Dec 2013 18:20:19 -0200<br>> From: joseclaudio.faria@gmail.com<br>> To: r-br@listas.c3sl.ufpr.br<br>> Subject: [R-br] Parcela subdividida: problema na declaração do erro associado a parcela<br>> <br>> Pessoal,<br>> <br>> O CMR abaixo é da apostila do Walmes.<br>> <br>> #------------------------------------------------------------------------------------------<br>> # dados<br>> ps <- expand.grid(BL=c("I", "II", "III", "IV"),<br>> ES=c("e1", "e2"),<br>> AD=c("a1", "a2", "a3"))<br>> 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)<br>> str(ps)<br>> #<br>> #------------------------------------------------------------------------------------------<br>> # análise de variância (erro A = BL x AD)<br>> m0 <- aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data=ps) #<br>> termo Error para declarar erro A<br>> m0 <- aov(alt ~ BL + AD * ES + Error(BL:AD), data=ps) # o mesmo com<br>> fórmula comprimida<br>> <br>> summary(m0)<br>> #------------------------------------------------------------------------------------------<br>> <br>> Encontrei muitos outros exemplos (na web e em livros) com essa<br>> notação: Error(BL:AD)<br>> <br>> Contudo o R (3.0.2patched) não está aceitando essa notação!<br>> Fornece sempre a mensagem de erro:<br>> In aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data = ps) :<br>> modelo Error() é singular<br>> <br>> Contudo o R está aceitando esse modelo:<br>> m0 <- aov(alt ~ BL + AD + Error(BL/AD) + ES + AD:ES, data=ps)<br>> <br>> Mas nesse caso o resíduo A (associado a bloco e parcela - AD) não fica correto.<br>> <br>> Alguém que anda trabalhando com esses modelos poderia ajudar?<br>> <br>> Ab,<br>> -- <br>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>> Jose Claudio Faria<br>> Estatistica<br>> UESC/DCET/Brasil<br>> joseclaudio.faria at gmail.com<br>> Telefones:<br>> 55(73)3680.5545 - UESC<br>> 55(73)9100.7351 - TIM<br>> 55(73)8817.6159 - OI<br>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>> _______________________________________________<br>> R-br mailing list<br>> R-br@listas.c3sl.ufpr.br<br>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<br>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.<br></div> </div></body>
</html>