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. Jose 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(cbind(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.