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

ちょっとめんどくさいがマニュアルでブートストラップができる。 解説はhttp://www.ats.ucla.edu/stat/stata/faq/ownboot.htm 要は、rclassのプログラムを組んで、simulateコマンドでシミュレートすると、bstatコマンドで組み込みのbootstrapと同じように結…