求助，做SFA时用超越对数生产函数，为何STATA和frontier4.1估计结果差异非常大？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请教，用超越对数生产函数做SFA时候，用的是Battese and Coelli(1992)模型，为啥stata和Frontier4.1的估计结果差异非常大，请高手指点，送3000论坛币，谢谢
啊！！！
（1）stata命令如下：. xtfrontier  logy logk logl logklogk logllogl logklogl tlogk tlogl t tt,tvd
（2）Frontier4.1的ins.txt文件内容是
1       1=ERROR COMPONENTS MODEL, 2=TE EFFECTS MODEL
bc92.dta.txt   DATA FILE NAME
bc94.out.txt    OUTPUT FILE NAME
1       1=PRODUCTION FUNCTION, 2=COST FUNCTION
y        LOGGED DEPENDENT VARIABLE (Y/N)
29         NUMBER OF CROSS-SECTIONS
18       NUMBER OF TIME PERIODS
522        NUMBER OF OBSERVATIONS IN TOTAL
9       NUMBER OF REGRESSOR VARIABLES (Xs)
y       MU (Y/N) [OR DELTA0 (Y/N) IF USING TE EFFECTS MODEL]
y     ETA (Y/N) [OR NUMBER OF TE EFFECTS REGRESSORS (Zs)]
n      STARTING VALUES (Y/N)
        IF YES THEN
            BETA0

        BETA1 TO
   
BETAK



SIGMA SQUARED

GAMMA

MU
[OR DELTA0

ETA
DELTA1 TO

DELTAP]



NOTE: IF YOU ARE SUPPLYING STARTING VALUES

AND YOU HAVE RESTRICTED MU [OR DELTA0] TO BE

ZERO THEN YOU SHOULD NOT SUPPLY A STARTING

VALUE FOR THIS PARAMETER.

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏45 回帖

关键词：frontier 超越对数生产函数 frontie Stata front 结果 函数 Stata 对数 SFA

没用过frontier4.1
你看看算法是否相同

很奇怪，如果没有加入t，则两种软件结果完全一样；如果加入t，就完全不一样。问题可能就出在t上，不知是否这样？
另外，想问一下，用STATA来做随机前沿，采用BC(1992)的模型，并基于超越对数生产函数，t是否从1，2，3，...18.依次类推呢？感谢牛人回答啊

其实，STATA是输出不来无效函数的！
STATA用的是BC（1992）程序，最好还是用frontier软件比较合适。
我也和你遇到了同样的问题！

好多钱呀，我想要，可惜我不会

最好给出数据，别人才可能帮忙验证的

xhw88184004 发表于 2009-9-21 16:57
请教，用超越对数生产函数做SFA时候，用的是Battese and Coelli(1992)模型，为啥stata和Frontier4.1的估计结果差异非常大，请高手指点，送3000论坛币，谢谢
啊！！！
（1）stata命令如下：. xtfrontier  logy logk logl logklogk logllogl logklogl tlogk tlogl t tt,tvd
（2）Frontier4.1的ins.txt文件内容是
1       1=ERROR COMPONENTS MODEL, 2=TE EFFECTS MODEL
bc92.dta.txt   DATA FILE NAME
bc94.out.txt    OUTPUT FILE NAME
1       1=PRODUCTION FUNCTION, 2=COST FUNCTION
y        LOGGED DEPENDENT VARIABLE (Y/N)
29         NUMBER OF CROSS-SECTIONS
18       NUMBER OF TIME PERIODS
522        NUMBER OF OBSERVATIONS IN TOTAL
9       NUMBER OF REGRESSOR VARIABLES (Xs)
y       MU (Y/N) [OR DELTA0 (Y/N) IF USING TE EFFECTS MODEL]
y     ETA (Y/N) [OR NUMBER OF TE EFFECTS REGRESSORS (Zs)]
n      STARTING VALUES (Y/N)
        IF YES THEN
            BETA0

        BETA1 TO
   
BETAK



SIGMA SQUARED

GAMMA

MU
[OR DELTA0

ETA
DELTA1 TO

DELTAP]



NOTE: IF YOU ARE SUPPLYING STARTING VALUES

AND YOU HAVE RESTRICTED MU [OR DELTA0] TO BE

ZERO THEN YOU SHOULD NOT SUPPLY A STARTING

VALUE FOR THIS PARAMETER.

至今这个问题依然悬而未决，是不愿意作答，还是没有发现呢？

光这样说不好判断，发数据来看一下。


已回复，见12楼。

蓝色 发表于 2009-10-7 09:27
最好给出数据，别人才可能帮忙验证的

相关数据附上了，其中：logY,logK,logL分别是总产出，资本存量，劳动投入量（均取对数）；h,insti,gov,fdi,open分别表示影响技术无效方程的变量。
  现在恳请各位高人帮忙验证一下，如果用Battese and Coelli(1995)模型，且采用超越对数生产函数，如下：
    logY=a0+a1*logK+a2*logL+a3*(logK*logK)+a4*(logL*logL)+a5*(logK*logL)+a6*(logK*t)+a7*(logL*t)+a8*t+a9*(t*t)
感谢！！！

根据楼主提供的数据在stata上跑了一下，发现根本就不收敛呵。
我也发现STATA和frontier4.1估计结果有差异，估计应该从算法上找原因。楼主可查一下xtfrontier的命令详解和frontier4.1的使用手册。
另有一问题顺便请教：在技术无效方程中，因变量是否就是无效项 u_it ？能否对效率值te回归？还有技术无效方程的影响变量是取原值呢还是取对数值？
呵呵，没帮助解决问题，倒是提了一对问题。先谢了！