逆ミルズ

probitを推定した後、
predict y_hat, xb
でlatent variableの予測値を計算し、
gen IMR = normalden(y_hat)/normal(y_hat)
で計算。