随机前沿模型(SFA)-原理解读

本教学视频由合工大经济学院张王飞制作，请尊重作者的知识产权，购买后使用，勿用于商业用途，或上传至网络！谢谢随机前沿模型（SFA）原理和软件实现一、SFA原理在经济学中，常常需要估计生产函数或者成本函数。生产函数f(x)的定义为：在给定投入x情况下的最大产出。但现实中的产商可能达不到最大产出的前沿，为了，假设产商i的产量为：iiiyf(x,)（1）其中，为待估参数；i为产商i的水平，满足i01。如果i=1，则产商i正好处于效率前沿。同时，考虑生产函数还会受到随机冲击，故将方程（1）改写成：iviiiyf(x,)e（2）其中，ive0为随机冲击。方程（2）意味着生产函数的前沿ivif(x,)e是随机的，故此类模型称为“随机前沿模型”（stochasticfrontiermodel）。随机前沿模型最早由Aigner,LovellandSchmidt(1977)提出，并在实证领域运用广泛，KumbhakarandLovell(2000)为该领域的研究写了一本著作，有兴趣的同学可以去参考。假设ok1i1ikif(x,)exx（柯布道格拉斯生产函数，共有K个投入品），则对方程（2）取对数可得：Ki0kkiiik1lny=+lnxln（3）由于i01，故iln0。定义iiu=-ln0，则方程3可以写成：Ki0kkiiik1lny=+lnx-u其中，iu0为“无效率”项，反映产商i距离效率前沿面的距离。混合扰动项iii分布不对称，使用OLS估计不能估计无效率项iu。为了估计无效率项iu，必须对ii、的分布作出假设，并进行更有效率的MLE（最大似然估计）估计。一般，无效率项的分布假设有如下几种：(1)半正态分布(2)截断正态分布(3)指数分布在一般的论文中，使用的最多的是半正态分布随机前沿模型可以很容易地用于估计成本函数，经过与生产函数的随机前沿模型类似的推导可得：Ki0yikkiiik1lnc=+lnylnP+u本教学视频由合工大经济学院张王飞制作，请尊重作者的知识产权，购买后使用，勿用于商业用途，或上传至网络！谢谢其中，ic为产商i的成本，iy为产出，kiP为要素K的价格，iu为无效率项，i为成本函数的随机冲击。注意混合误差项的形式（符号）。对于成本函数，iu=0意味着产商达到最低成本的效率前沿；反之，如果iu0，则产商需付出更高的成本。iu是否存在的检验使用随机前沿模型的前提是无效率项iu存在，此假定可以通过检验“220u1uH:=0vsH:0”来判断是否成立。使用单边的广义似然比检验。二、软件实现Frontier4.1软件是由TimCoelli开发的一款专门用于完成随机前沿分析的软件，它可以用最大似然估计随机前沿成本模型和随机前沿生产模型，下面简单介绍一下该软件的使用方法，更加详细的说明可以参考英文指导《AGuidetoFRONTIERVersion4.1:AComputerProgramforStochasticFrontierProductionandCostFunctionEstimation》Eg1.DTA用于输入数据，是一个纯文本文件，数据文件的格式必须是3+K[+p]列。第一列是评价体系的序号；第二列是时期t；第三列是因变量；第四列之后是K个自变量；[+p]仅当选择TEEFFECTSMODEL模型输入。EG1.INS设置命令11=ERRORCOMPONENTSMODEL,2=TEEFFECTSMODEL选择模型eg1.dtaDATAFILENAME数据文件eg1.outOUTPUTFILENAME结果存储文件21=PRODUCTIONFUNCTION,2=COSTFUNCTION选择生产模型（1）还是成本模型（2）nLOGGEDDEPENDENTVARIABLE(Y/N)变量是不是已经进行了对数运算25NUMBEROFCROSS-SECTIONS评价体系数目1NUMBEROFTIMEPERIODS本教学视频由合工大经济学院张王飞制作，请尊重作者的知识产权，购买后使用，勿用于商业用途，或上传至网络！谢谢时期数目25NUMBEROFOBSERVATIONSINTOTAL总记录数目2NUMBEROFREGRESSORVARIABLES(Xs)自变量个数YMU(Y/N)[ORDELTA0(Y/N)IFUSINGTEEFFECTSMODEL]假设U的分布。Y表示截断分布，N表示半正态分布nETA(Y/N)[ORNUMBEROFTEEFFECTSREGRESSORS(Zs)]y表示时变模型，n表示非时变模型。nSTARTINGVALUES(Y/N)选择n其他设置保持不变