[R-br] Chi-quadrado demonstra associação, mas e daí?

[Felipe]

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