Luiz Roberto, como estou sem R ao alcance vou fornecer uma solução aproximada. Em primeiro lugar, você vai querer transformar esses vetores de texto em vetores de números inteiros, utilizando a função as.integer()[1]. Quanto à comparação dos vetores propriamente dita, a solução mais simples é com a função prop.test()[2]. Utilize sum(x, na.rm = TRUE)[3] para calcular o número de sucessos em cada vetor, e length[4](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). Leonardo Ferreira Fontenelle[5] Em Qua 19 ago. 2015, às 14:46, Luiz Roberto Martins Pinto escreveu:
Caros companheiros da Lista R,
# 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.
# Preciso de uma função que me retorne a probabilidade (p-value) do # acesso A1 ser semelhante ao acesso A2.
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") 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")
Agradeço a boa-vontade de quem puder ajudar.
Abraços,
Luiz Roberto.
Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil
luizroberto.uesc@gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831
_________________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Links: 1. https://stat.ethz.ch/R-manual/R-devel/library/base/html/integer.html 2. https://stat.ethz.ch/R-manual/R-devel/library/stats/html/prop.test.html 3. https://stat.ethz.ch/R-manual/R-devel/library/base/html/sum.html 4. https://stat.ethz.ch/R-manual/R-devel/library/base/html/length.html 5. http://lattes.cnpq.br/9234772336296638