《对数极大似然估计》PPT课件

《《对数极大似然估计》PPT课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第十八章对数极大似然估计EViews包含了一些常用方法，如最小二乘法、非线性最小二乘法、加权最小二乘法、TSLS、GMM、ARIMA、ARCH、GARCH等方法，这些方法可以解决可能遇到的大多数估计问题。但是，我们在研究中也可能会碰到一些不在上述之列的特殊的模型，这些模型可能是现存方法的一个扩展，也可能是一类全新的问题。为了能解决这些特殊的问题，EViews提供了对数极大似然估计这一工具来估计各种不同类型的模型。对数极大似然估计提供了一个一般的，开放的工具，可以通过这个工具极大化相关参数的似然函数对一大类模型进行估计。使用对数极大似然估计时，我们用EViews的序列生

2、成器，将样本中各个观测值的对数似然贡献描述为一个未知参数的函数。可以给出似然函数中一个或多个参数的解析微分，也可以让EViews自动计算数值微分。EViews将寻找使得指定的似然函数最大化的参数值，并给出这些参数估计的估计标准差。在本章，我们将详细论述对数极大似然估计，说明其一般特征。并给出了一些可以使用该方法的具体的例子。§18.1概论用对数极大似然估计来估计一个模型，主要的工作是建立用来求解似然函数的说明文本。用EViews指定对数极大似然函数的说明是很容易的，因为似然函数的说明只是一系列对序列的赋值语句，这些赋值语句在极大化的过程中被反复的计算。我们所要做的只是

3、写下一组语句，在计算时，这些语句将描述一个包含每个观测值对似然函数贡献的序列。首先，我们简单地回顾一下线性回归模型的对数极大似然估计方法。考虑多元线性回归模型的一般形式t=1,2,…,T(1)其中k是解释变量个数，T是观测值个数，随机扰动项～，那么服从如下的正态分布：～(2)其中Y的随机抽取的T个样本观测值的联合概率为(3)这就是变量Y的似然函数。对似然函数求极大值和对对数似然函数求极大值是等价的，对数似然函数为(4)注意到，我们能将对数似然函数写成所有观测值t的对数似然贡献和的形式，(5)这里单个贡献由下面的式子给出：(6)以只含一个解释变量的方程为例。假定知道模型

4、参数的真实值，并且想用EViews产生一个包含每个观测值的贡献的序列。可以将已知的参数赋值给系数向量的c(1)到c(3)元素，然后把下面的赋值语句作为EViews的命令或程序来执行。Seriesres=y-c(1)-c(2)*xSeriesvar=c(3)SerieslogL1=-log(2*3.14159*var)/2-(res^2/var)/2前面两行语句描述了用来存储计算时的中间结果的序列。第一个语句创建了残差序列：res，而第二个语句创建了方差序列：var。而序列logL1包含了每个观测值的对数似然贡献的集合。下面考虑稍复杂的例子，假设数据是由条件异方差回归模

5、型生成的：(18.1)这里，x,y,z是观测序列，而1,2,3,,是模型的参数。有T个观测值的样本的对数似然函数可以写成：(18.2)这里，是标准正态分布的密度函数。Eviews中的标准正态分布的对数似然函数为，将对数似然函数写成所有观测值t的对数似然贡献的和的形式：(18.3)这里每个观测值的贡献由下面的式子给出：(18.4)将这一例子的对数极大似然函数过程写成下面的赋值语句：Seriesres=y-c(1)-c(2)*x-c(3)*zSeriesvar=c(4)*z^c(5)SerieslogL1=log(@dnorm(res/@sqrt(var)))

6、-log(var)/2前面两行语句创建了残差序列res和方差序列var，参数c(1),c(2),c(3)代表了回归系数，c(4)代表了，c(5)代表了，序列logL1包含了每个观测值的对数似然贡献的集合。现在假定不知道模型参数的真实值，而想使用数据来估计它。参数的极大似然估计被定义为：使得样本中所有随机抽取的一组观测值的联合概率密度，即似然函数取最大值的那组参数值。而对数极大似然方法使得寻找这些极大似然估计变得容易了。只需创建一个新的对数似然对象，把上面的赋值语句输入到logL的说明窗口，然后让EViews来估计这个模型。在输入赋值语句时，只需对上面的文本做两处微小的

7、改动就可以了。首先，把每行开头的关键字series删掉（因为似然说明暗含了假定序列是当前的）。第二，必须在说明中加入额外的一行（它为包含似然贡献的序列命名）。这样，要在logL说明窗口输入下面的内容：@logLlogL1res=y-c(1)-c(2)*x-c(3)*zvar=c(4)*z^c(5)logL1=log(@dnorm(res/@sqrt(var)))-log(var)/2对数似然函数的第一行，@logLlogL1，告诉EViews用logL1序列来存储似然贡献。余下的行定义了中间结果的计算和实际的似然贡献的计算。当用EViews估计模型参