# NOT RUN {
data("ProgramEffectiveness")
## Greene (2003), Table 21.1, col. "Probit"
fm_probit <- glm(grade ~ average + testscore + participation,
data = ProgramEffectiveness, family = binomial(link = "probit"))
summary(fm_probit)
# }
Run the code above in your browser using DataLab