ANOVA por medidas repetidas

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

Evandro, Algumas duvidas: 1) Os analisadores seriam uma restrição a casulização em meu ponto de vista o popular BLOCO, caso seja isso as interações não fazem sentido com esse fator 2) quando vc fala em medidas repetidas, quer dizer no tempo? 3) magnificações, foi um fator de interesse ou mais uma restrição? Att Em 3 de junho de 2011 10:59, Evandro Malanski <evanmal@hotmail.com>escreveu:
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@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Vou explicar melhor o caso: O objetivo deste trabalho foi comparar duas técnicas de medição. Isso pois em um outro trabalho eu uso a medição por duas técnicas diferentes, então preciso saber se as técnicas produzem o mesmo resultado. Contudo, aproveitei para verificar se o analisador é um fator que pode causar erros. Também, organismos menores podem estar mais sujeitos a erros de medição que organismos maiores (usando o microscópio), então verifiquei se o fator magnificação também poderia causar erros. Foi assim que "desenhei" o experimento. Então, as medidas repetidas não são ao longo do tempo, mas sim entre analisadores. A restrição da magnificação é somente para separar organismos maiores de menores, tanto é que tenho 30 + 30 + 30 = 90. Neste caso seria óbvio que haveriam diferenças significativas nas medições, como mostrou naquele sumário. Ficou melhor explicado? ____________________ 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 Date: Fri, 3 Jun 2011 11:12:52 -0300 From: eder@leg.ufpr.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] ANOVA por medidas repetidas Evandro,Algumas duvidas:1) Os analisadores seriam uma restrição a casulização em meu ponto de vista o popular BLOCO, caso seja isso as interações não fazem sentido com esse fator2) quando vc fala em medidas repetidas, quer dizer no tempo? 3) magnificações, foi um fator de interesse ou mais uma restrição?Att Em 3 de junho de 2011 10:59, Evandro Malanski <evanmal@hotmail.com> escreveu: 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@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Então, esta melhor explicado. Parece que esta correto então sua análise. 90 unidades experimentais vc teria se quisesse comparar as magnificações desconsiderando os outros fatores. Justifique aos editores da revista que vc tem realmente 540 unidades experimentais. Se quiser pode retirar as interações das magnificações, mas os resultados deste modelo não mudará suas conclusões. Vc também pode avaliar uma correlação de Pearson entre os valores medidos das tecnicas e dos avaliadores para cada magnificação. --- Em sex, 3/6/11, Evandro Malanski <evanmal@hotmail.com> escreveu: De: Evandro Malanski <evanmal@hotmail.com> Assunto: Re: [R-br] ANOVA por medidas repetidas Para: "R-BR Fórum" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 3 de Junho de 2011, 14:27 Vou explicar melhor o caso: O objetivo deste trabalho foi comparar duas técnicas de medição. Isso pois em um outro trabalho eu uso a medição por duas técnicas diferentes, então preciso saber se as técnicas produzem o mesmo resultado. Contudo, aproveitei para verificar se o analisador é um fator que pode causar erros. Também, organismos menores podem estar mais sujeitos a erros de medição que organismos maiores (usando o microscópio), então verifiquei se o fator magnificação também poderia causar erros. Foi assim que "desenhei" o experimento. Então, as medidas repetidas não são ao longo do tempo, mas sim entre analisadores. A restrição da magnificação é somente para separar organismos maiores de menores, tanto é que tenho 30 + 30 + 30 = 90. Neste caso seria óbvio que haveriam diferenças significativas nas medições, como mostrou naquele sumário. Ficou melhor explicado? ____________________ Evandro Malanski, Oc. Curso de Mestrado em Oceanografia Biológica Laboratório de Ecologia do Ictioplâncton - LEI Instituto de Oceanografia - IOUniversidade 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 Date: Fri, 3 Jun 2011 11:12:52 -0300 From: eder@leg.ufpr.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] ANOVA por medidas repetidas Evandro, Algumas duvidas: 1) Os analisadores seriam uma restrição a casulização em meu ponto de vista o popular BLOCO, caso seja isso as interações não fazem sentido com esse fator 2) quando vc fala em medidas repetidas, quer dizer no tempo? 3) magnificações, foi um fator de interesse ou mais uma restrição? Att Em 3 de junho de 2011 10:59, Evandro Malanski <evanmal@hotmail.com> escreveu: 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 - IOUniversidade 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@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Evandro, boa tarde Se tirar as interações (quais quer sejam) os resultados mudam sim... saem efeitos do modelo, aumentam os graus de liberdade e a variância do erro que é usada para testar os efeitos que ficam (teste F)! Não entendi "bem" como e experimento foi concebido, mas pelo meu pobre entendimento, alguns dos efeitos podem ser aninhados aos outros (não por interação) e outra que pelas esperanças pode ser que alguns testes não devam ser feitos pelo erro e sim pelas interações, isso é especificado pelo argumento Error() que o Benilton citou! att, FH 2011/6/3 Emmanuel Arnhold <emmanuelarnhold@yahoo.com.br>
Então, esta melhor explicado. Parece que esta correto então sua análise. 90 unidades experimentais vc teria se quisesse comparar as magnificações desconsiderando os outros fatores. Justifique aos editores da revista que vc tem realmente 540 unidades experimentais. Se quiser pode retirar as interações das magnificações, mas os resultados deste modelo não mudará suas conclusões. Vc também pode avaliar uma correlação de Pearson entre os valores medidos das tecnicas e dos avaliadores para cada magnificação.
--- Em *sex, 3/6/11, Evandro Malanski <evanmal@hotmail.com>* escreveu:
De: Evandro Malanski <evanmal@hotmail.com> Assunto: Re: [R-br] ANOVA por medidas repetidas Para: "R-BR Fórum" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 3 de Junho de 2011, 14:27
Vou explicar melhor o caso: O objetivo deste trabalho foi comparar duas técnicas de medição. Isso pois em um outro trabalho eu uso a medição por duas técnicas diferentes, então preciso saber se as técnicas produzem o mesmo resultado. Contudo, aproveitei para verificar se o analisador é um fator que pode causar erros. Também, organismos menores podem estar mais sujeitos a erros de medição que organismos maiores (usando o microscópio), então verifiquei se o fator magnificação também poderia causar erros. Foi assim que "desenhei" o experimento. Então, as medidas repetidas não são ao longo do tempo, mas sim entre analisadores. A restrição da magnificação é somente para separar organismos maiores de menores, tanto é que tenho 30 + 30 + 30 = 90. Neste caso seria óbvio que haveriam diferenças significativas nas medições, como mostrou naquele sumário. Ficou melhor explicado? ____________________ *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*
------------------------------ Date: Fri, 3 Jun 2011 11:12:52 -0300 From: eder@leg.ufpr.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] ANOVA por medidas repetidas
Evandro, Algumas duvidas: 1) Os analisadores seriam uma restrição a casulização em meu ponto de vista o popular BLOCO, caso seja isso as interações não fazem sentido com esse fator 2) quando vc fala em medidas repetidas, quer dizer no tempo? 3) magnificações, foi um fator de interesse ou mais uma restrição? Att
Em 3 de junho de 2011 10:59, Evandro Malanski <evanmal@hotmail.com<http://br.mc431.mail.yahoo.com/mc/compose?to=evanmal@hotmail.com>
escreveu:
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@listas.c3sl.ufpr.br<http://br.mc431.mail.yahoo.com/mc/compose?to=R-br@listas.c3sl.ufpr.br> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-----Anexo incorporado-----
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br<http://br.mc431.mail.yahoo.com/mc/compose?to=R-br@listas.c3sl.ufpr.br> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

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@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@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)

Evandro, pelo que entendi seu experimento quer comparar a concordância das leituras dos analisadores, nesse caso tem um pacote no R chamado psych, que tem uma rotina chamada wkappa que dá a medida de concordância entre os analisadores. Instale o pacote psych e veja ?wkappa. Quanto as técnicas: existe alguma padrão? Caso haja,tem um pacote no R chamado DiagnosisMedque calcula ROC´s, sensibilidade e especificidade da técnica. Instale o pacote DiagnosisMed e veja ?ROC. Como em todos os pacotes do R, existem exemplo de rotinas. Espero ter te ajudado. Edson Lira Estatístico Manaus-Amazonas ________________________________ De: Evandro Malanski <evanmal@hotmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Sexta-feira, 3 de Junho de 2011 9:59 Assunto: [R-br] ANOVA por medidas repetidas 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 - IOUniversidade 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@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Ta um pouco difícil de entender. Mas me parece que vc tem realmente 540 unidades experimentais. 90 x 2 tecnicas x 3 avaliadores. As técnicas e avaliações fora feitas em "momentos" ou "tempos" diferentes e isso pode ser um problema dependendo da maneira como vc casualizou seu experimento. Se vc casualizou a distribuição de seus 18 tratamentos as suas 540 unidades experimentais não deve ter problema. Se vc fez técnicas e avaliadores em tempos diferentes (sistematizou em parte e não casualisou) aí seu erro experimental pode ter problema. --- Em sex, 3/6/11, Evandro Malanski <evanmal@hotmail.com> escreveu: De: Evandro Malanski <evanmal@hotmail.com> Assunto: [R-br] ANOVA por medidas repetidas Para: r-br@listas.c3sl.ufpr.br Data: Sexta-feira, 3 de Junho de 2011, 13:59 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 - IOUniversidade 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 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Acabei de pegar um livro aqui na biblioteca da universidade, e tem um capítulo que tem um exemplo semelhante ao meu e fala sobre os erros, como os que relatei. O livro é: "Statistics: an introduction using R" de Michael Crawley. Vou dar uma estudada aqui. Quando solucionar o problema, escreverei para avisar. Muito obrigado aos que estão ajudando e questionando! ____________________ 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 Date: Fri, 3 Jun 2011 07:53:44 -0700 From: emmanuelarnhold@yahoo.com.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] ANOVA por medidas repetidas Ta um pouco difícil de entender. Mas me parece que vc tem realmente 540 unidades experimentais. 90 x 2 tecnicas x 3 avaliadores. As técnicas e avaliações fora feitas em "momentos" ou "tempos" diferentes e isso pode ser um problema dependendo da maneira como vc casualizou seu experimento. Se vc casualizou a distribuição de seus 18 tratamentos as suas 540 unidades experimentais não deve ter problema. Se vc fez técnicas e avaliadores em tempos diferentes (sistematizou em parte e não casualisou) aí seu erro experimental pode ter problema. --- Em sex, 3/6/11, Evandro Malanski <evanmal@hotmail.com> escreveu: De: Evandro Malanski <evanmal@hotmail.com> Assunto: [R-br] ANOVA por medidas repetidas Para: r-br@listas.c3sl.ufpr.br Data: Sexta-feira, 3 de Junho de 2011, 13:59 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 - IOUniversidade 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 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Evandro, A unidade experimental onde incide o efeito de magficação é o organismo. Dentro de organismo foram obtidas as medidas indexadas por outros fatores, analizador(3)*técnica(2). Então dadas as informações fornecidas, o modelo estatístico correspondente ao experimento (não planejado, porque se você você saberia o modelo) em fórmula R é y~magn*aval*tecn+Error(org) Error: org Df Sum Sq Mean Sq F value Pr(>F) magn 1 0.165 0.16500 0.2386 0.6265 Residuals 88 60.863 0.69163 Error: Within Df Sum Sq Mean Sq F value Pr(>F) aval 2 0.66 0.33009 0.3000 0.7409 tecn 1 0.03 0.03458 0.0314 0.8594 magn:aval 2 0.22 0.10797 0.0981 0.9065 magn:tecn 1 0.08 0.07777 0.0707 0.7905 aval:tecn 2 1.20 0.60220 0.5474 0.5789 magn:aval:tecn 2 2.12 1.06026 0.9637 0.3823 Residuals 440 484.06 1.10014 À 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 (7)
-
Benilton Carvalho
-
Eder David Borges da Silva
-
Edson Lira
-
Emmanuel Arnhold
-
Evandro Malanski
-
Fernando Henrique Toledo
-
Walmes Zeviani