双边随机前沿模型，SFA-two-tier程序+案例数据+论文

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

   包括双边随机
前沿模型，SFA,Two-tier Stochastic Frontier Mo
del：案例数据+命令+程序命令源代码+案例文章****面板改进版***
*--
---------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.dro
p if year==.  //去掉缺失值，要多打一个=号，去除空缺值；

2. 
 
    rename u_hat       Neg  //  u代表了抑制
效应（或称为技术效率损失）
        rename w_hat      
 Pos  //W或V则代表了促进效应（或称为技术效率增益）
    renam
e uw_diff_exp   NI   //净效应
        repla
ce NI=Pos-Neg     //可能算的有差异可以备用
*-转换成百分数
，便于统计        
        foreach v of varli
st Neg-NI{
          replace `v' = `v'*1
00
        }
      
        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 p
50 p75) format(%
6.3f) c(s)  //分区域统计
**时间
特征**
    *-w_hat u_hat uw_diff (Panel 1-
2)
logout, save(hlh_z1) excel replace:  
/// //输出 Excel 表格   
        tabstat w_h
at u_hat uw_diff,  
             by(
year) stat(mean sd p25 p50 p75) format(%

6.3f) c(s)  //促进效应、抑制效应时间特征
            
   
    *-NI (Panel 3-4)
logout, save(hl
h_z2) excel replace:  /// //输出 Excel 表格

         tabstat  NI,  
            
 by(year) stat(mean sd p25 p50 p75) form
at(%
6.3f) c(s)        //净效应时间特征
   

下载链接：https://download.csdn.net/download/weixin_45892228/89139350

点击下载：双边随机前沿模型，SFA-two-tier程序+案例数据+论文