有坛友反映有些双边随机模型代码包没办法安装，而且命令包不对，请大家在文件夹中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
        }
      
        tabstat  Neg -NI, stat(mean sd p25 p50 p75) format(%6.3f) c(s)  //算出U,W的双边效应估计拟合值
4.根据分组要求重新把bug修订为 **时间特征**和**地区特征**，
logout, save(hlh_z0) excel replace:  /// //输出 Excel 表格
      tabstat  Pos 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 表格
         tabstat  NI,  ////
             by(year) stat(mean sd p25 p50 p75) format(%6.3f) c(s)        //净效应时间特征