data("SwissLabor")
### Gerfin (1996), Table I.
fm_probit <- glm(participation ~ . + I(age^2), data = SwissLabor,
family = binomial(link = "probit"))
summary(fm_probit)
### alternatively
fm_logit <- glm(participation ~ . + I(age^2), data = SwissLabor,
family = binomial)
summary(fm_logit)
Run the code above in your browser using DataLab