STLS分两种情况,其实只要弄清楚原理了就很好写了
第一种Truncation,并且假设从OLS开始第一步,那么就是
Reg y x
Predict yhat, xb
drop if y > 2*yhat
然后不断循环,最终能得到趋于相同的结果。
第二种情况Censoring,也假设从OLS开始第一步,那么就是
Reg y x
Predict yhat, xb
replace y = yhat if y > 2*yhat
drop if yhat < 0
然后不断循环重复上述步骤,最终也能得到相同的结果