Comparações Múltiplas em modelo linear com co-variável.

Prezados, Estou realizando uma análise estatística através de um modelo de análise de covariância: y ~ sítio + bloco + tratamento + densidade + sitio:bloco + sitio:tratamento onde: sitio, bloco e tratametno são qualitativos densidade é quantitativo, foi inserido no modelo para reduzir a variação dos dados Os resultados da análise apontam efeitos significativos de tratamento e densidade. O desafio, agora, é fazer comparações múltiplas entre os tratamentos, no entanto, fixando-se densidades médias que variam entre os tratamentos, ou seja, descontando os efeitos de sítios e blocos que não são significativos. Tenho densidades médias para os quatro tratamentos: (400, 1200, 1150 e 900 para os tratamentos A, B, C e D), desta forma, utilizando os parâmetros estimados para cada tratamento + efeitos da densidade, gostaria de saber como eu poderia comparar as estimativas e verificar se os tratamentos diferem entre si. Eu poderia usar o tukey? Lembrando que a variável resposta é contínua... Caso não tenha sido claro, trago alguns resultados para exemplificar o problema. Desde já agradeço o empenho de todos e Feliz Natal! Att. -- Danilo Scorzoni Ré Engenheiro Florestal Mestrando em Ciência Florestal FCA / UNESP - Botucatu (14) 8180-2494

Você pode usar a dobradinha contrast::contrast() com multcomp::glht(). O correto não seria comparar os tratamentos sob a mesma densidade? Veja este script da análise de covariância em experimento zootecnico, pode ser um começo. http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap14ancova-iso.R À 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 ==========================================================================
participantes (2)
-
Danilo Scorzoni Ré
-
Walmes Zeviani