Quedarnos con un elemento específico de un objeto en R. Viene de una duda planteada por un lector que surgía a raíz de un post dedicado al $p$-valor y al tamaño muestral. José Ignacio desea almacenar en un objeto el $p$-valor asociado a un test de Wald para la independencia entre los niveles de un factor. Una sugerencia para poder trabajar con este elemento:
y <- c(rep(1, 200), rep(0, 100))
x <- c(rep(1, 32), rep(0, 168), rep(1, 15), rep(0, 85))
datos <- data.frame(y, x)
table(datos)
modelo.1 <- glm(y ~ x, data = datos, family = binomial)
summary(modelo.1)
Para la realización del test de Wald con R vamos a emplear la librería lmtest:
library(lmtest)
guarda <- waldtest(modelo.1)
str(guarda)
Y ahora la “parte interesante” de esta entrada: cómo nos quedamos con el $p$-valor asociado al test de independencia de Wald que está en el objeto guarda:
pvalor_test <- guarda[2, 4]
Recordad siempre que str es la clave cuando realizamos este tipo de tareas. Saludos.