sas系统数据管理(编辑修改稿)

sas系统数据管理(编辑修改稿)内容摘要：

proc copy in=lib1 out=lib2 MTYPE=(data catalog)。 run。 例 将逻辑库 LIB1的所有 SAS文件移动到 LIB2中。 proc copy in=lib1 out=lib2 move。 run。 例 选择以 DAT开头的所有文件。 data (alter=zsw)。 var=39。 move39。 run。 proc copy in=ResDat out=work move。 select dat: (alter=zsw)。 run。 第八章 对数极大似然估计 极大似然估计法 (maximum likelihood, ML)，是不同于最小二乘法的另一种参数估计方法，是从极大似然原理发展起来的其他估计方法的基础。 虽然其应用没有最小二乘法普遍，但在计量经济学理论上占据很重要的地位，因为极大似然原理比最小二乘原理更本质地揭示了通过样本估计母体参数的内在机理，计量经济学理论的发展更多的是以极大似然估计原理为基础的，对于一些特殊的计量经济学模型，只有极大似然方法才是很成功的估计方法。 EViews包含了一些常用方法，如最小二乘法、非线性最小二乘法、加权最小二乘法、 TSLS、 GMM、 ARIMA、ARCH、 GARCH等方法，这些方法可以解决可能遇到的大多数估计问题。 但是，我们在研究中也可能会碰到一些不在上述之列的特殊的模型，这些模型可能是现存方法的一个扩展，也可能是一类全新的问题。 为了能解决这些特殊的问题 ， EViews提供了 对数极大似然估计对象 这一工具来估计各种不同类型的模型。 对数极大似然估计对象提供了一个一般的 ， 开放的工具 ， 可以通过这个工具极大化相关参数的似然函数对一大类模型进行估计。 使用对数极大似然估计对象估计时 ， 我们用 EViews的序列生成器 ， 将样本中各个观测值的对数似然贡献描述为一个未知参数的函数。 可以给出似然函数中一个或多个参数的解析微分 ， 也可以让 EViews自动计算数值微分。 EViews将寻找使得指定的似然函数最大化的参数值 ， 并给出这些参数估计的估计标准差。 在本章 ， 我们将详细论述对数极大似然估计对象 ，说明其一般特征。 并给出了一些可以使用该方法的具体的例子。 167。 对数极大似然估计的基本原理 167。 极大似然估计的基本原理 设总体的概率密度函数为 P， 其类型是已知的 ， 但含有未知参数 （ 向量 ） 。 我们的目的就是依据从该总体抽得的随机样本 y1, y2, … , yT ， 寻求对  的估计。 观测值 y1, y2, … , yT 的联合密度函数被给定为 （ ） 其中： y = ( y1, y2, … , yT )。 将这一联合密度函数视为参数  的函数 ， 称为样本的似然函数 （ likelihood function）。 TttyPL1)()。 ( ψy 极大似然原理就是寻求参数的估计值 ， 使得所给样本值的概率密度 （ 即似然函数 ） 的值在这个参数值之下 ， 达到最大。 在当前的情形下 ， 就是寻求  的估计值 ， 使得似然函数 L(y。 ) 相对于给定的观测值 y1, y2, … , yT 而言达到最大值 ， 就被称为极大似然估计量。 ψˆψˆ 在 L(y。 ) 关于 i（ i =1, 2, … , n， n是未知参数的个数 ）的偏导数存在时 ， 要使 L(y。 ) 取最大值 ，  必须满足 , i =1, 2, …, n （ ） 由上式可解得 n1 向量  的极大似然估计值 ，而式 ()也被称为似然函数。 0)。 (   ψyLiψˆ 因为 L(y。 ) 与 ln[L(y。 ))] 在同一点处取极值 ， 所以也可以由 , i =1, 2, …, n （ ） 求得，因为对数可将乘积变成求和，所以，式 ()往往比直接使用式 ()来得方便。 式 ()也被称为对数似然函数。 0)。 (ln   ψyLi 考虑多元线性回归模型的一般形式 , t =1, 2 , … , T () 其中 k 是解释变量个数 ， T 是观测值个数 ， 随机扰动项 ～ ， 那么 yt 服从如下的正态分布： ～ 其中 () tktkttt uxxxy   22110tu ),0( 2Nty ),(2 tNktkttt xxx   22110 y 的随机抽取的 T 个样本观测值的联合概率函数为 （ ） 这就是变量 y的似然函数 ， 未知参数向量 ={1, 2,… k, 2}。 对似然函数求极大值和对数似然函数求极大值是等价 的 ， 式 ()的 对数似然函数形式 为： （ ） TtttyTTTttTeyPyyyPL122)(2121212π)2(1)(),(),( βTtttTtttyyTL122212222)(21)π2l n (21)(21)π2l n (2),(lnβ 注意 ， 可以将对数似然函数写成 t 时刻所有观测值的对数似然贡献和的形式 ， （ ） 这里对数似然的单个贡献 （ 用小写字母表示 ） 由下面的式子给出： （ ） ),(),(ln 212  ββ TttlL2222 )(21)π2l n (21),(ttt yl  βTt ,2,1  式 （ ） 也可用标准正态分布的密度函数  表示 （ ） 式中 标准正态分布的对数似然函数  为 （ ） 这里对数似然函数每个观测值的贡献式 ()又可以由下面的式子给出： （ ） TtttTtyTL122122 )(21)l n (21π)2l n (2),(ln β    Tttty12 )l n (21)(ln Tttt zTz1221π)2l n (2)(ln   tttyz )l n (21ln),( 2    tttyl βTt ,2,1  167。 EViews极大似然对象概述 用对数极大似然估计来估计一个模型，主要的工作是建立用来求解似然函数的说明文本。 用 EViews指定对数极大似然函数的说明是很容易的，因为似然函数的说明只是一系列对序列的赋值语句，这些赋值语句在极大化的过程中被反复的计算。 我们所要做的只是写下一组语句，在计算时，这些语句将描述一个包含每个观测值对似然函数贡献的序列。 注意到 ， 我们能将对数似然函数写成所有观测值 t 的对数似然贡献和 的形式 ， 这里 单个贡献 由下面的式子给出：。