
André, Veja se este exemplo pode te ajudar(pelo menos de início para suas análises). Como não deu maiores detalhes suponho que, quando usou as palavras negativo e positivo na sua tabela, estava querendo indicar ausência e presença dos parasitas. Primeiro pode verificar se existe diferença entre as proporções de negativo e positivo, assim pode verificar se a presença ou não de cada parasita tem diferença significativa em seus dados: ##Diferença entre Proporções e IC95%(d) dados<-matrix(c(250,15,34,14),ncol=2,byrow=T) dados ## Calculando as proporções entre Cryptosporidium negativo e Cryptosporidium positivo p11<-(dados[1,1]/(sum(dados[1,]))) p11 p21<-(dados[2,1]/(sum(dados[2,]))) p21 d<-p11-p21 d vd<-((p11*(1-p11))/(sum(dados[1,])-1)) + ((p21*(1-p21))/(sum(dados[2,])-1)) dvd<-sqrt(vd) z<-qnorm(0.975) li<- d - (z*dvd) ls<- d + (z*dvd) cbind(d,li,ls) # Aqui você verificou que a diferença entre Cryptosporidium negativo e positivo é significativa #Abaixo o famoso Qui-Quadrado. Note que o argumento sim=200 faz com que o p-valor seja calculado por simulações para que o pressupostos da validade do teste Qui-quadrado não sejam violados. chisq.test(dados,sim=500) ##Há evidências de se rejeitar H0,logo Cryptosporidium e Giardia estão associados #Pode também ter verificar qual a chance de ausência ou não dos parasitas. ##Razão de Chances ou Odds Ratio (OR) e IC95%(OR) OR<-(dados[1,1]*dados[2,2])/(dados[1,2]*dados[2,1]) vf<-(1/dados[1,1])+(1/dados[1,2])+(1/dados[2,1]+(1/dados[2,2])) dpf<-sqrt(vf) z<-qnorm(0.975) li<-exp(log(OR)-z*dpf) ls<-exp(log(OR)+z*dpf) cbind(OR,li,ls) ## A chance de ausência Cryptosporidium e Giardia é 6,8 vezes maior que a presença podendo variar entre 3 e 15,4 vezes ao nível de confiança de 95%. -- Atenciosamente Felipe E. Barletta Mendes Estatístico - Conre3 9766-A +55 (41)-92077191 +55 (41)-33287216