Let\'s consider data following with probit regression performed.
data(\'UnionWage\', package = \'pglm\') anb <- pglm(union ~ wage + exper + rural, UnionWag