割と簡単。
ただのブートストラップは、
bootstrap, reps(100) seed(1): reg y x
ここでreps()の引数はサンプリング回数
ちょい複雑なやつマージナルエフェクト
bootstrap margf = (expression), reps(100) seed(1): reg y x
このうちexpressionの中には、マージナルエフェクトの公式を推定値を_b[x]として書いてやればいい。
２つ以上マージナルエフェクトが必要なら
bootstrap margf1 = (expression)
margf2 + (expression), reps(100) seed(1): reg y x
と続ければいい。
割と簡単だし、結構計算も早い。
