Cláudio... isso não seria um warnings! É apenas para avisar que o a solução do sistema está vindo por uma inversa generalizada pois sua matriz de incidência não tem posto completo!
2013/12/4 Jose Claudio Faria <[email protected]> > 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. >
_______________________________________________ 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.
