<!DOCTYPE html>
<html>
<head>
<title></title>
</head>
<body><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">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">prop.test()</a>. Utilize <a href="https://stat.ethz.ch/R-manual/R-devel/library/base/html/sum.html">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">length</a>(na.omit(x)) ou length(x[!is.na(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 id="sig4487139"><div class="signature" id="signature"><a href="http://lattes.cnpq.br/9234772336296638" title="Currículo Lattes">Leonardo Ferreira Fontenelle</a><br></div>
</div>
<div> </div>
<div> </div>
<div>Em Qua 19 ago. 2015, às 14:46, Luiz Roberto Martins Pinto escreveu:<br></div>
<blockquote type="cite"><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">luizroberto.uesc@gmail.com</a><br></div>
<div>skype: lrmpinto<br></div>
</div>
<div><div><a href="http://lattes.cnpq.br/2732314327604831">http://lattes.cnpq.br/2732314327604831</a> <br></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">R-br@listas.c3sl.ufpr.br</a><br></div>
<div><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">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">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></div>
</blockquote><div> </div>
</body>
</html>