Source: https://www.math.univ-toulouse.fr/~besse/Wikistat/pdf/tp_ozone1_ancova_logit.pdf