有坛友反映有些双边随机模型代码包没办法安装，而且命令包不对，请大家在文件夹中S里面的程序包复制到plus文件下的S文件夹应该可以的。

包括双边随机前沿模型，SFA,Two-tier Stochastic Frontier Model：案例数据+命令+程序命令源代码+案例文章****面板改进版***
*-----------Panel SFA 一般化设定-----------------
*
*y[it] = a + x[it]B + v[it] - u[it] (1)
*ε[it]=v[it]− ui //这是我们算出vit的关键
*v[it] ~ N(0, vsigma^2[it]) (2)
*u[it] ~ N(mu[it], usigma^2[it]) (3)
*mu[it] = q[it]B (4)
*vsigma^2[it] = w[it]B (5)
*usigma^2[it] = z[it]B (6)
*
*-----------------------------------------------



补充内容 (2023-8-11 16:02):
有几处代码小bug解释一下，1.drop if year==.//去掉缺失值，要多打一个=号，去除空缺值；
2.
rename u_hat Neg//u代表了抑制效应（或称为技术效率损失）
rename w_hat Pos//W或V则代表了促进效应（或称为技术效率增益）
rename uw_diff_exp NI //净效应
replace NI=Pos-Neg //可能算的有差异可以备用
3.添加完善
*-转换成百分数，便于统计
foreach v of varlist Neg-NI{
replace `v' = `v'*100
}

tabstatNeg -NI, stat(mean sd p25 p50 p75) format(%6.3f) c(s)//算出U,W的双边效应估计拟合值
4.根据分组要求重新把bug修订为 **时间特征**和**地区特征**，
logout, save(hlh_z0) excel replace:/// //输出 Excel 表格
tabstatPos Neg uw_diff,////
by(zone) stat(mean sd p25 p50 p75) format(%6.3f) c(s)//分区域统计
**时间特征**
*-w_hat u_hat uw_diff (Panel 1-2)
logout, save(hlh_z1) excel replace:/// //输出 Excel 表格
tabstat w_hat u_hat uw_diff,////
by(year) stat(mean sd p25 p50 p75) format(%6.3f) c(s)//促进效应、抑制效应时间特征

*-NI (Panel 3-4)
logout, save(hlh_z2) excel replace:/// //输出 Excel 表格
tabstatNI,////
by(year) stat(mean sd p25 p50 p75) format(%6.3f) c(s) //净效应时间特征