[R-br] ANOVA por medidas repetidas
Benilton Carvalho
beniltoncarvalho em gmail.com
Sexta Junho 3 11:16:45 BRT 2011
Voce precisa levar em consideracao as repeticoes e sua fonte de erro.
Veja o uso da funcao Error() dentro do aov() - leia com carinho a ajuda.
O item acima te levara muito provavelmente para modelos mistos, uma
estrategia muito comum para o ajuste de modelos assim. Para tanto,
veja o pacote lme4 e sua extensiva documentacao.
b
2011/6/3 Evandro Malanski <evanmal em hotmail.com>:
> Pessoal, tenho algumas dúvidas quanto a ANOVA. Se puderem me ajudar, ficarei
> muito grato.
>
> Montei o seguinte experimento com medição de organismos:
> 3 analisadores diferentes
> 3 magnificações no microscópio possíveis
> 2 técnicas de medição diferentes
> 90 organismos
>
> Dos 90 organismos, 30 foram medidos pela magnificação 1, outros 30 na
> magnificação 2, e os últimos 30 na magnificação 3, por cada um dos 3
> analisadores, e usando as 2 técnicas.
> Ou seja, para cada organismo, tenho 6 medidas.
> Eu montei a planilha da seguinte forma, mostrando somente as 10 primeiras
> medições:
>
>> compara
> Larva Técnica Analisador Magnificação Medida
> 1 1 1 1 20x 3.8303
> 2 2 1 1 20x 4.2159
> 3 3 1 1 20x 4.4470
> 4 4 1 1 20x 4.9549
> 5 5 1 1 20x 4.0265
> 6 6 1 1 20x 3.9235
> 7 7 1 1 20x 3.7697
> 8 8 1 1 20x 3.9166
> 9 9 1 1 20x 3.5777
> 10 10 1 1 20x 3.7515
> ...
>
> Como eu tenho medidas repetidas, usei a seguinte análise do R:
>
> compara.aov<- aov(Medida~Analisador+Magnificação+Técnica+
> Analisador:Magnificação+Analisador:Técnica+Magnificação:Técnica+
> Analisador:Magnificação:Técnica)
>
> Ou seja, até mesmo a interação entre os fatores analiso.
> Contudo, agora que enviei para correção o meu trabalho, questionaram que a
> análise está incorreta.
> Segue o sumário da análise que fiz, e destaco o grau de liberdade do
> resíduo, que foi a partir daí que me questionaram:
>
>> summary(compara.aov)
> Df Sum Sq Mean Sq F
> value Pr(>F)
> Analisador 2 0.43 0.21 0.2253
> 0.7984
> Magnificação 2 1278.94 639.47 672.4928
> <2e-16 ***
> Técnica 1 0.03 0.03 0.0288
> 0.8654
> Analisador:Magnificação 4 0.13 0.03 0.0330 0.9979
> Analisador:Técnica 2 0.44 0.22 0.2336
> 0.7918
> Magnificação:Técnica 2 0.86 0.43 0.4521 0.6366
> Analisador:Magnificação:Técnica 4 0.10 0.03 0.0269 0.9986
> Residuals 522 496.37
> 0.95
> ---
> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> Esta análise eu fiz para observar se haveriam diferenças significativas
> entre técnicas.
> Falaram que minha comparação não está baseada nas 90 larvas medidas, mas sim
> em todas as 540 medidas que eu tenho, como verificado nos graus de liberdade
> dos resíduos.
> Alguém saberia me explicar melhor isso?
> E, como eu poderia solucionar esse erro? Modificando a organização da
> planilha dos dados? Modificando a análise?
> Agradeço a ajuda.
>
> ____________________
> Evandro Malanski, Oc.
> Curso de Mestrado em Oceanografia Biológica
> Laboratório de Ecologia do Ictioplâncton - LEI
> Instituto de Oceanografia - IO
> Universidade Federal de Rio Grande - FURG
> +55(53)3233-6529 / +55(53)8439-6755
>
> Master Course in Biological Oceanography
> Ichthyoplankton Ecology Laboratory - LEI
> Oceanography Institute - IO
> Federal University of Rio Grande - FURG
> +55(53)3233-6529 / +55(53)8439-6755
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
--
Successful people ask better questions, and as a result, they get
better answers. (Tony Robbins)
Mais detalhes sobre a lista de discussão R-br