Verdade! Como deu erro eu não tentei rodar o summary.
Muito obrigado Tiago! 2013/12/4 Tiago Souza Marçal <[email protected]>: > 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: [email protected] >> To: [email protected] >> 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 >> [email protected] >> 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. > > _______________________________________________ > R-br mailing list > [email protected] > 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. -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ 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 [email protected] 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.
