完成期末作业的同时来更一下博客

背景介绍

对中国储蓄存款总额(Y,亿元)与GDP(亿元)两个变量进行一元线性回归,检验并消除自相关。

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言


自相关分为一阶自相关和高阶自相关,通常使用DW检验来查看是否存在自相关,然后再用LM检验确定自相关的阶数。

自相关检验

DW检验

自相关检验dwr语言 自相关检验dw查表_统计学_02


在拟合方程的时候就会显示出DW统计量,本案例中是0.14

自相关检验dwr语言 自相关检验dw查表_统计学_03


自相关检验dwr语言 自相关检验dw查表_最小二乘_04


需要查找DW统计量分布表来确定接受和拒绝的上下限。

本案例拒绝H0,认为存在一阶自相关

LM检验

DW统计量的局限在于只能检验一阶自相关,因此再进行LM检验

LM统计量既可以检验一阶自相关也可以检验高阶自相关

同样,H0为不存在自相关,H1为存在自相关

Eviews操作如下:先建立一个方程,然后在view中选择LM检验

自相关检验dwr语言 自相关检验dw查表_最小二乘_05

选择2阶检验

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_06


查看统计量的值以及对应的p值

p值小于0.05说明存在二阶自相关

但是到底存在多少阶自相关要反复检验

如何判定自相关阶数:

使用AIC或SIC准则判断

选择使得AIC最小的阶数

自相关检验dwr语言 自相关检验dw查表_最小二乘_07

FGLS消除自相关

自相关的解决:可行广义最小二乘(FGLS)

一阶

经济计量模型中误差项自相关的最常见形式是一阶线性自回归形式。

自相关检验dwr语言 自相关检验dw查表_最小二乘_08


其中自相关检验dwr语言 自相关检验dw查表_差分_09为残差,自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_10是非自相关的,并满足一定假设条件。

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_11


由于

自相关检验dwr语言 自相关检验dw查表_统计学_12



自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_13

根据DW统计量的值计算出一阶自回归系数 自相关检验dwr语言 自相关检验dw查表_差分_14为0.93

然后构造广义差分变量gdx(自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_15)和gdy(自相关检验dwr语言 自相关检验dw查表_统计学_16

自相关检验dwr语言 自相关检验dw查表_差分_17


自相关检验dwr语言 自相关检验dw查表_统计学_18


自相关检验dwr语言 自相关检验dw查表_差分_19

输入gdy和gdx的表达式,其中(-1)表示t-1期

自相关检验dwr语言 自相关检验dw查表_最小二乘_20


自相关检验dwr语言 自相关检验dw查表_统计学_21

估计新的模型

自相关检验dwr语言 自相关检验dw查表_最小二乘_22

再查看DW统计量检验是否消除了自相关

自相关检验dwr语言 自相关检验dw查表_差分_23

高阶

对于高阶自相关的处理以二阶为例

自相关检验dwr语言 自相关检验dw查表_最小二乘_24


对应的广义差分变换为

自相关检验dwr语言 自相关检验dw查表_最小二乘_25

因此需要先对残差序列进行回归,得到广义差分变换的系数。

系统会自动保存最后一次回归的残差,保存在resid里面。

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_26

可以看出第一项是缺失的,因为我之前做了二阶LM检验,残差是LM检验的残差。但我们要得到原模型的残差,所以要对原模型再回归一次。

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_27

这时的残差才是原模型的残差。
这一点需要特别小心!

自相关检验dwr语言 自相关检验dw查表_统计学_28

新建一个序列res把这个残差保存下来

自相关检验dwr语言 自相关检验dw查表_差分_29

对残差进行回归

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_30

得到一阶和二阶自回归系数

自相关检验dwr语言 自相关检验dw查表_回归_31

然后进行二阶广义差分
定义变量gdx2和gdy2

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_32


自相关检验dwr语言 自相关检验dw查表_最小二乘_33


自相关检验dwr语言 自相关检验dw查表_最小二乘_34

然后进行新模型的回归

自相关检验dwr语言 自相关检验dw查表_自相关检验dwr语言_35

通过DW统计量可以看出此时已经消除了自相关

自相关检验dwr语言 自相关检验dw查表_最小二乘_36

或者进行LM检验也可看出消除了自相关

自相关检验dwr语言 自相关检验dw查表_最小二乘_37

最后将广义最小二乘模型的系数估计量作为原模型的系数估计量即可

自相关检验dwr语言 自相关检验dw查表_差分_38