<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:12pt"><div><span>Evandro, pelo que entendi seu experimento quer comparar a concordância das leituras dos analisadores,</span></div><div><span>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.<br></span></div><div><span><br></span></div><div><span>Quanto as técnicas: existe alguma padrão? Caso haja,tem um pacote no R chamado </span><span>DiagnosisMed</span><span> que calcula ROC´s, sensibilidade e especificidade da técnica. Instale o pacote </span><span>DiagnosisMed e veja ?ROC.<br></span></div><div><br><span></span></div><div>Como em todos os pacotes do R, existem exemplo de rotinas. <br><span></span></div><div><span><br></span></div><div>Espero ter te ajudado.<br></div><div>Edson
Lira<br>Estatístico<br>Manaus-Amazonas<br></div><div style="font-family: Courier New, courier, monaco, monospace, sans-serif; font-size: 12pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><font size="2" face="Arial"><hr size="1"><b><span style="font-weight:bold;">De:</span></b> Evandro Malanski <evanmal@hotmail.com><br><b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br<br><b><span style="font-weight: bold;">Enviadas:</span></b> Sexta-feira, 3 de Junho de 2011 9:59<br><b><span style="font-weight: bold;">Assunto:</span></b> [R-br] ANOVA por medidas repetidas<br></font><br><div id="yiv2022262164">
<style><!--
#yiv2022262164 .yiv2022262164hmmessage P
{
margin:0px;padding:0px;}
#yiv2022262164 .yiv2022262164hmmessage
{
font-size:10pt;font-family:Tahoma;}
--></style>
Pessoal, tenho algumas dúvidas quanto a ANOVA. Se puderem me ajudar, ficarei muito grato.<br><br>Montei o seguinte experimento com medição de organismos:<br>3 analisadores diferentes<br>3 magnificações no microscópio possíveis<br>2 técnicas de medição diferentes<br>90 organismos<br><br>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.<br>Ou seja, para cada organismo, tenho 6 medidas.<br>Eu montei a planilha da seguinte forma, mostrando somente as 10 primeiras medições:<br><br>> compara<br> Larva Técnica Analisador Magnificação Medida<br>1 1 1 1 20x
3.8303<br>2 2 1 1 20x 4.2159<br>3 3 1 1 20x 4.4470<br>4 4 1 1 20x 4.9549<br>5 5 1 1 20x
4.0265<br>6 6 1 1 20x 3.9235<br>7 7 1 1 20x 3.7697<br>8 8 1 1 20x 3.9166<br>9 9 1 1 20x 3.5777<br>10
10 1 1 20x 3.7515<br>...<br><br>Como eu tenho medidas repetidas, usei a seguinte análise do R:<br><br>compara.aov<- aov(Medida~Analisador+Magnificação+Técnica+<br> Analisador:Magnificação+Analisador:Técnica+Magnificação:Técnica+<br> Analisador:Magnificação:Técnica)<br><br>Ou seja, até mesmo a interação entre os fatores analiso.<br>Contudo, agora que enviei para correção o meu trabalho, questionaram que a análise está incorreta.<br>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:<br><br>>
summary(compara.aov)<br> Df Sum Sq Mean Sq F value Pr(>F) <br>Analisador 2 0.43 0.21 0.2253 0.7984 <br>Magnificação 2 1278.94 639.47 672.4928 <2e-16
***<br>Técnica 1 0.03 0.03 0.0288 0.8654 <br>Analisador:Magnificação 4 0.13 0.03 0.0330 0.9979 <br>Analisador:Técnica 2 0.44 0.22 0.2336 0.7918 <br>Magnificação:Técnica 2 0.86 0.43 0.4521 0.6366
<br>Analisador:Magnificação:Técnica 4 0.10 0.03 0.0269 0.9986 <br>Residuals <font style="" color="#FF0000"><b>522 </b></font>496.37 0.95 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 <br><br>Esta análise eu fiz para observar se haveriam diferenças significativas entre técnicas.<br>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.<br>Alguém saberia me explicar melhor isso?<br>E, como eu poderia
solucionar esse erro? Modificando a organização da planilha dos dados? Modificando a análise?<br>Agradeço a ajuda.<br><br><div><font style="" color="#0099ff" face="Arial"><font color="#0033ff">____________________</font></font></div>
<div><font style="" color="#0099ff" face="Arial"><font color="#0033ff"><b>Evandro Malanski, Oc. </b><br><i>Curso de Mestrado em Oceanografia Biológica</i></font></font></div>
<div><font style="" color="#0099ff" face="Arial"><font color="#0033ff"><i>Laboratório de Ecologia do Ictioplâncton - LEI</i></font></font></div>
<div><font style="" color="#0033ff" face="Arial"><i>Instituto de Oceanografia - IO</i></font></div>
<font style="" face="Arial"><i><font color="#0099ff"><font color="#0033ff">Universidade Federal de Rio Grande - FURG<br></font></font><font color="#0099ff"><font color="#0033ff">+55(53)3233-6529 / +55(53)8439-6755<br><br>Master Course in Biological Oceanography<br>Ichthyoplankton Ecology Laboratory - LEI<br>Oceanography Institute - IO<br>Federal University of Rio Grande - FURG<br></font></font></i></font><font style="" face="Arial"><i><font color="#0099ff"><font color="#0033ff">+55(53)3233-6529 / +55(53)8439-6755</font></font></i></font><br><br>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br><br><br></div></div></div></body></html>