Resultados Prueba chi cuadrado frente a logit binomial

Question

mis datos tienen este aspecto:

Tengo dos grupos de personas: estudiantes y ancianos.

59 de 266 estudiantes dijeron "SÍ" 23 de 127 ancianos dijeron "SÍ"

El resto dijo "no"

Intento averiguar si esta diferencia en sus respuestas es estadísticamente significativa.

Mi primera forma fue la prueba de chi cuadrado que dice que son diferentes con una p inferior a 0,05. Sin embargo, mi logit binomial con una sola variable independiente trae este resultado resultado de logit binomial, lo que significa que la diferencia no es significativa.

¿Qué me estoy perdiendo?

Gracias.

Answer 1

Acabo de hacer esto yo mismo y obtuve resultados totalmente diferentes a los tuyos.

x_student <- rep("student", 266)
y_student <- c(rep(0, 266 - 59), rep(1, 59))
x_elderly <- rep("elderly", 127)
y_elderly <- c(rep(0, 127 - 23), rep(1, 23))
x <- c(x_student, x_elderly)
y <- c(y_student, y_elderly)
L <- glm(y ~ x, family = binomial)
summary(L)

m <- matrix(c(
    23, 137 - 23,
    59, 266 - 59
), 2, 2)
chisq.test(m)

Mi prueba chi-cuadrado devuelve un valor p de $0.253$ mientras que mi coeficiente de regresión logística tiene un valor p de $0.354$ ninguno de los cuales coincide con su trabajo. (Podría creer que cometí un error en mi código, pero creo que es correcto). Aunque estos p-valores son diferentes, están lo suficientemente cerca, y ninguno señala datos que sean inconsistentes con la hipótesis nula implícita de que las personas mayores y los estudiantes responden de forma diferente a la encuesta. Dado que ambos grupos tienen aproximadamente un $80/20$ dividido, esto no es sorprendente, sobre todo porque sólo se dispone de unos cientos de observaciones y, por tanto, no de mucha potencia para detectar una pequeña diferencia.