Fatorial com reptições perdidas
Boa Páscoa a todos, estou analizando um experimento fatorial como variáveis respostas altura, diâmetro do peito, número de galhos, mdida da copa da árvore, os fatores são tipo de plantio (puro e mista) e aração (arado e aarado). Cada medida tem 3 repetições, o problema é que em alguns casos, a árvore morreu, estou com um impasse, não sei se considero zero, fato que infuenciaria na média ou resposta perdida "missing" (NA). Fiz uma tentativa usando o NA, rodei o modelo com aov, verifiquei os pressupostos, fiz o desdobramento das interações, quando tentei rodar o tukey, o R me deu a mensagem dizendo que x e y não tem a mesma dimensão. Alguem poderia me sugerir o que fazer com a arvore morta? Caso tenha que dropar os NA's como proceder? Edson Lira Estatístico Manaus-Amazonas
Edson, Você tá usando a agricolae::HSD.test()? Bem, internamente ela não tem proteção contra os NA. O que eu faço é já remover os NA da base para iniciar a análise. Faça dados <- dados[complete.cases(dados),] Lembre-se a HSD.test() vai usar a média harmônica do número de repetições que só é recomendado para desbalanceamentos pouco severos. Outra coisa implícita na analises é a de que os NA surgem ao acaso, ou seja, as árvores que morreram não podem ter sido por função dos tratamentos. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Walmes, estou usando o TukeyHSD, não estou usando o agricolae. Usei o complete.cases e zerou o banco. Edson Lira Estatístico Manaus-Amazonas --- Em dom, 24/4/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Fatorial com reptições perdidas Para: r-br@listas.c3sl.ufpr.br Data: Domingo, 24 de Abril de 2011, 17:40 Edson, Você tá usando a agricolae::HSD.test()? Bem, internamente ela não tem proteção contra os NA. O que eu faço é já remover os NA da base para iniciar a análise. Faça dados <- dados[complete.cases(dados),] Lembre-se a HSD.test() vai usar a média harmônica do número de repetições que só é recomendado para desbalanceamentos pouco severos. Outra coisa implícita na analises é a de que os NA surgem ao acaso, ou seja, as árvores que morreram não podem ter sido por função dos tratamentos. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Edson, A TukeyHSD() não desdobra interação. Ela compara os efeitos principais e comparar todos os efeitos de interação entre si. Isso não é o desdobramento usual que fazemos em estatística experimental. Na situação de desbalanceamento e presença de interação, os contrastes da TukeyHSD() são viesados. Desdobrar essa interação na presença de desbalanceamento requer muita algebra e código R. Até onde eu conheço não tem função que retorne isso pronto. Coloquei um procedimento exemplo de como eu tenho feito isso nas Rídiculas do LEG, http://www.leg.ufpr.br/doku.php/ridiculas. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Valeu Walmes. Edson Lira Estatístico Manaus-Amazonas --- Em dom, 24/4/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Fatorial com reptições perdidas Para: r-br@listas.c3sl.ufpr.br Data: Domingo, 24 de Abril de 2011, 19:39 Edson, A TukeyHSD() não desdobra interação. Ela compara os efeitos principais e comparar todos os efeitos de interação entre si. Isso não é o desdobramento usual que fazemos em estatística experimental. Na situação de desbalanceamento e presença de interação, os contrastes da TukeyHSD() são viesados. Desdobrar essa interação na presença de desbalanceamento requer muita algebra e código R. Até onde eu conheço não tem função que retorne isso pronto. Coloquei um procedimento exemplo de como eu tenho feito isso nas Rídiculas do LEG, http://www.leg.ufpr.br/doku.php/ridiculas. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (2)
-
Edson Lira -
Walmes Zeviani