Contingency Tables: Fisher & Chi Squared

Tenemos tres poblaciónes (N1, N2, N3) con parámetros p (Proporciones poblacionales) para una enfermedad X

Obtenemos una muestra de 100 individuos para cada población.

Y la cantidad (proporción) de positivos.

Aquí vemos la estimación de la proporción de positivos para cada población, con los límites de confianza de Wilson (Agresti, 2007)

Efecto del grupo:

H_O : P_N1 = P_N2 = P_N3

Pr.Positive <- table(resultados.n$group,resultados.n$resultado)
mosaicplot(Pr.Positive, color = c("darkred", "gold"), xlab ="Group", ylab = "Pr.", main= "Proportion of Positive",
           legend= TRUE)

Aproximación con χ²

###
chi2 <- chisq.test(cbind(grouped$positivos,grouped$negativos))
chi2

    Pearson's Chi-squared test

data:  cbind(grouped$positivos, grouped$negativos)
X-squared = 111.07, df = 2, p-value < 2.2e-16

grouped$expected <- chi2$expected[c(1,2,3)]
grouped$std.res <- chi2$stdres[c(1,2,3)]
grouped