stata时间序列分析_第七章 时间序列分析（3）

案例应用

在本小节中，我们通过一个案例来将本章所涉及的时间序列分析过程一一串连并予以展现，旨在令读者能够更深刻地体会到时间序列数据的分析特点及要点。

以1984年第1季度到2009年第3季度的美国联邦基金回报率(f)与三年期债券回报率(b)为例进行分析，具体流程如下：

使用stata16打开在目录“D:\stata16\shuju\chap07”中的“0701.dta”数据文件，命令如下：

use“D:\stata16\shuju\chap07\0701.dta”,clear

然后创造一个时间变量，并使用tsset命令使得数据成为时间序列数据。在Command界面输入以下命令：

gen date=tq(1984q1)+_n-1

format% tq date /*将日期格式转化为字符串格式*/

tsset date /*时间序列数据进行回归前，需要使用tsset命令进行设定*/

其中，tq(1984q1)是一种拟函数，表示输入所有的数据都会转化为整数当量。以上述为例，Stata是以1960年基础时间原点，所以1984q1表示第一季度距离1960年间隔96个季度，对应生成的整数当量为96，_n-1类似于一个步长值为1的循环语句，Stata中以增量为1进行叠加到date中。

(1)时间序列平稳性检验

由于对于时间序列的检验，一般检验其三种形式，第一种形式对包含截距项形式的回归模型检验。我们先对f执行ADF检验，在检验之前，我们需要对其滞后阶数进行选择，由于Stata中没有直接的方式进行判断滞后阶数，而是需要自己进行编程，选择滞后的阶数。这里选用Schwert(1989)建议的最大滞后阶数P_max=[12*(T/100)^1/4]的法则确定最大阶数，然后利用AIC准则，找出合适的阶数。在Command界面输入如下命令：

di 12*(104/100)^(1/4) /*104为本数据集的样本个数*/