<div dir="ltr">Eu recorri à função que o Leonardo sugeriu... e rodei..<div><br></div><div># Sugestão do Leonardo:</div><div><br></div><div><div style="font-size:12.8000001907349px"><div style="font-size:12.8000001907349px">A1= as.integer(c("1",NA,"1","0","1","0","1","1","0","1","1","0","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1",NA,"1",NA,"0","1","0","0","1","1","1","1","1","1","0","0","1","1","1","1","0","1","0","1","0","0","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","0","1","1","1","1","1","1","0"))</div><div style="font-size:12.8000001907349px">A2= as.integer(c("1","1","1","0","1","0","1","1","0","1","1","0","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","1","0","1","1","1","1","1","1","1","1",NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,"1","0","1","1","1","1","0","0","1","1","1","1","0","1","0","1",NA,NA,"1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","0","1","1","1","1","1","1","1"))</div><div style="font-size:12.8000001907349px"><br></div></div><div><div style="font-size:12.8000001907349px"> sucA1=sum(A1, na.rm = TRUE)</div><div style="font-size:12.8000001907349px"> sucA2=sum(A2, na.rm = TRUE)</div><div><div style="font-size:12.8000001907349px"> tenA1=length(na.omit(A1))</div><div style="font-size:12.8000001907349px"> tenA2=length(na.omit(A2))</div><div><span style="font-size:12.8000001907349px">prop.test(c(sucA1,sucA2),c(tenA1,tenA2)).</span><br></div></div></div><div><span style="font-size:12.8000001907349px"><br></span></div><div># De acordo com a função:</div><div><br></div><div><div>smokers  <- c( 83, 90, 129, 70 )</div><div>patients <- c( 86, 93, 136, 82 )</div><div>prop.test(smokers, patients) </div></div><div><br></div><div># Então, semelhantemente, rodei:</div><div><br></div><div>prop.test(A1, A2) <br></div><div><br></div><div># Output:</div><div><br></div><div>Erro em prop.test(A1, A2) : elements of 'n' must be positive.<br></div><div><br></div><div>Comentário: acredito que a função não admite dados perdidos 'NA'.</div><div><br></div><div>É isto mesmo?</div><div><br></div><div><br></div><div><br></div><div><br></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div>Luiz Roberto Martins Pinto<br>Prof. Pleno/DCET/UESC</div><div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div><div>Ilhéus-Bahia-Brasil<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div></div></div></div>
<br><div class="gmail_quote">Em 19 de agosto de 2015 18:57, Luiz Roberto Martins Pinto <span dir="ltr"><<a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Semelhança no perfil.<div><br></div><div>O perfil de cada acesso é descrito por 101 marcadores moleculares.</div><div><br></div><div>Pensando, agora, na sugestão do Leonardo, e pensando no seu questionamento, talvez avaliar a proporção de marcadores iguais fosse uma boa análise deste perfil.  Não tinha pensado nisto...</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Luiz Roberto.</div><div><br></div></font></span></div><div class="gmail_extra"><span class=""><br clear="all"><div><div><div dir="ltr"><div><div>Luiz Roberto Martins Pinto<br>Prof. Pleno/DCET/UESC</div><div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div><div>Ilhéus-Bahia-Brasil<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div></div></div></div>
<br></span><div><div class="h5"><div class="gmail_quote">Em 19 de agosto de 2015 18:36, Cesar Rabak <span dir="ltr"><<a href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Luiz Roberto,<div><br></div><div>Frente a sua resposta para o Leonardo, é imperioso propor-lhe a questão:</div><div><br></div><div>Qual seria a definição de semelhança entre os dois acessos?</div><div><br></div><div>Com essa informação podemos tentar encontrar um teste estatístico com mais potência para a comparação...</div><div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">2015-08-19 16:10 GMT-03:00 Luiz Roberto Martins Pinto <span dir="ltr"><<a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Leonardo,<div><br></div><div>Agradeço a sugestão. Mas é possível que as proporções sejam as mesmas entre dois acessos diferentes, se os marcadores estiverem em posições diferentes.</div><span><font color="#888888"><div><br></div><div>Luiz Roberto.</div></font></span></div><div class="gmail_extra"><span><br clear="all"><div><div><div dir="ltr"><div><div>Luiz Roberto Martins Pinto<br>Prof. Pleno/DCET/UESC</div><div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div><div>Ilhéus-Bahia-Brasil<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div></div></div></div>
<br></span><div><div><div class="gmail_quote">Em 19 de agosto de 2015 15:07, Leonardo Ferreira Fontenelle <span dir="ltr"><<a href="mailto:leonardof@leonardof.med.br" target="_blank">leonardof@leonardof.med.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>




<div><div>Luiz Roberto, como estou sem R ao alcance vou fornecer uma solução aproximada.<br></div>
<div> </div>
<div>Em primeiro lugar, você vai querer transformar esses vetores de texto em vetores de números inteiros, utilizando a função <a href="https://stat.ethz.ch/R-manual/R-devel/library/base/html/integer.html" target="_blank">as.integer()</a>. Quanto à comparação dos vetores propriamente dita, a solução mais simples é com a função <a href="https://stat.ethz.ch/R-manual/R-devel/library/stats/html/prop.test.html" target="_blank">prop.test()</a>. Utilize <a href="https://stat.ethz.ch/R-manual/R-devel/library/base/html/sum.html" target="_blank">sum(x, na.rm = TRUE)</a> para calcular o número de sucessos em cada vetor, e <a href="https://stat.ethz.ch/R-manual/R-devel/library/base/html/length.html" target="_blank">length</a>(na.omit(x)) ou length(x[!<a href="http://is.na" target="_blank">is.na</a>(x)]) para calcular o número de tentativas em cada vetor. Então utilize a função prop.test, informando como primeiro argumento c(sucessosA1, sucessosA2) e segundo argumento c(tentativasA1, tentativasA2).<br></div>
<div> </div>
<div><div><a href="http://lattes.cnpq.br/9234772336296638" title="Currículo Lattes" target="_blank">Leonardo Ferreira Fontenelle</a><br></div>
</div><div><div>
<div> </div>
<div> </div>
<div>Em Qua 19 ago. 2015, às 14:46, Luiz Roberto Martins Pinto escreveu:<br></div>
</div></div><blockquote type="cite"><div><div><div dir="ltr"><div>Caros companheiros da Lista R,<br></div>
<div> </div>
<div><div>#  Os dados abaixo referem-se a 2 acessos de mandioca (A1 e A2), e 101 marcadores RAPD (1=presença da marca; 0=ausência da marca; NA=dado perdido); portanto os dados são binários.<br></div>
<div> </div>
<div># Preciso de uma função que me retorne a probabilidade (p-value) do acesso A1 ser semelhante ao acesso A2.<br></div>
<div> </div>
<div>A1=c("1",NA,"1","0","1","0","1","1","0","1","1","0","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1",NA,"1",NA,"0","1","0","0","1","1","1","1","1","1","0","0","1","1","1","1","0","1","0","1","0","0","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","0","1","1","1","1","1","1","0")<br></div>
<div>A2=c("1","1","1","0","1","0","1","1","0","1","1","0","1","1","1","1","1","0","1","1","1","1","1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","1","0","1","1","1","1","1","1","1","1",NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,"1","0","1","1","1","1","0","0","1","1","1","1","0","1","0","1",NA,NA,"1","1","1","1","1","1","1","1","0","0","1","1","1","1","1","1","0","1","1","1","1","1","1","1")<br></div>
</div>
<div> </div>
<div>Agradeço a boa-vontade de quem puder ajudar.<br></div>
<div> </div>
<div>Abraços,<br></div>
<div> </div>
<div>Luiz Roberto.<br></div>
<div> </div>
<div><div><div dir="ltr"><div><div><div>Luiz Roberto Martins Pinto<br></div>
<div>Prof. Pleno/DCET/UESC<br></div>
</div>
<div>Laboratório de Estatística Computacional<br></div>
<div>Universidade Estadual de Santa Cruz<br></div>
<div><div>Ilhéus-Bahia-Brasil<br></div>
<div> </div>
<div><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br></div>
<div>skype: lrmpinto<br></div>
</div>
<div><div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br></div>
<div> </div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div><div><u>_______________________________________________</u><br></div>
<div>R-br mailing list<br></div>
<div><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br></div>
<div><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></div>
<div>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div>
</blockquote><div> </div>
</div>

<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div></div></div>
</blockquote></div><br></div>