如何用stata下计算三阶段DEA效率？

frontier a b c d
predict u,u  \\可直接求出管理无效率
predict xb,xb  \\截距项
        gen v=slack11-xb-u   //求出技术无效率
        egen vmax=max(v)
        egen xbmax=max(xb)
        gen x1= inv1+xbmax-xb+vmax-v  //算出inv1调整后的值X1，其它投入变量同理
        replace x1=inv1 if x1==.    //缺失值用原始值替代
        drop xb u v vmax xbmax

已经给出要审核