Salida del modelo logístico en R

Question

Estoy tratando de interpretar los siguientes tipos de modelo logístico:

mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial)

¿Es la salida del predict(mdl) la probabilidad esperada de éxito para cada punto de datos? ¿Hay una manera simple de tabular las probabilidades para cada nivel de factor del modelo, en lugar de todos los puntos de datos?

Answer 1

Las páginas de ayuda de

predict.glm

estado: "de modo que para un valor predeterminado modelo binomial el valor predeterminado predicciones son de log-odds (probabilidades en la escala logit) y" tipo = "respuesta"' da la predicción de las probabilidades". Así, predict(mdl) devuelve el registro de(probabilidades), y el uso de "tipo = "respuesta" devuelve la predicción de probabilidades. Usted puede encontrar este juguete ejemplo instructivo:

> y <- c(0,0,0,1,1,1,1,1,1,1)
> prop.table(table(y))
y
  0   1 
0.3 0.7 
> glm.y <- glm(y~1, family = "binomial")
> ## predicted log(odds)
> predict(glm.y)
        1         2         3         4         5         6         7         8 
0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 
        9        10 
0.8472979 0.8472979 
> ## predicted probabilities (p = odds/(1+odds))
> exp(predict(glm.y))/(1+exp(predict(glm.y)))
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 
> predict(glm.y, type = "response")
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 

Respecto a su segunda pregunta, usted puede ser que desee comprobar hacia fuera los efectos de un paquete de http://socserv.socsci.mcmaster.ca/jfox/Misc/effects/index.html por John Fox; véase también su JSS artículo "Efecto Muestra en R para Modelos Lineales Generalizados" (páginas 8 a 10).