*stata的命令(z为0-1变量): treatreg y x1-xn, tr(z=w1-wm) two
*等价于以下命令: prob z w1-wm predict gw, xb g lambda=normalden(gw)/normal(gw) if z==1 replace lambda=-normalden(gw)/normal(-gw) if z==0 reg y z x1-xn lambda
*stata的命令(z为0-1变量): heckman y x1-xn, two sel(z=w1-wm)
*等价于以下命令: prob z w1-wm predict gw, xb g lambda=normalden(gw)/normal(gw) if z==1 reg y x1-xn lambda if z==1