22.3 Predicción

Las predicciones se obtienen también con la función predict:

p.est <- predict(modelo, type = "response")

El parámetro type = "response" permite calcular las probabilidades estimadas de la segunda categoría.

Podríamos obtener una tabla de clasificación:

cat.est <- as.numeric(p.est > 0.5)
tabla <- table(datos$nsatisfa, cat.est)
tabla
##       cat.est
##         0  1
##   bajo 44 10
##   alto  7 38
print(100*prop.table(tabla), digits = 2)
##       cat.est
##           0    1
##   bajo 44.4 10.1
##   alto  7.1 38.4

Por defecto predict obtiene las predicciones correspondientes a las observaciones (modelo$fitted.values). Para otros casos hay que emplear el argumento newdata.