9.1 Ajuste: función glm

Para el ajuste (estimación de los parámetros) de un modelo lineal generalizado a un conjunto de datos (por máxima verosimilitud) se emplea la función glm:

ajuste <- glm(formula, family = gaussian, datos, ...)

El parámetro family indica la distribución y el link. Por ejemplo:

• gaussian(link = "identity"), gaussian(link = "log")

• binomial(link = "logit"), binomial(link = "probit")

• poisson(link = "log")

• Gamma(link = "inverse")

Para cada distribución se toma por defecto una función link (mostrada en primer lugar; ver help(family) para más detalles).

Muchas de las herramientas y funciones genéricas disponibles para los modelos lineales son válidas también para este tipo de modelos: summary, coef, confint, predict, anova, ….

Veremos con más detalle el caso particular de la regresión logística.