运用stata进行时间序列分析

1、STATA 从入门到精通第十一章第十一章 时间序列分析时间序列分析Page 2STATA从入门到精通从入门到精通11.1 基本时间序列模型的估计基本时间序列模型的估计n 在许多情况下，人们用时间序列的观测时期代表的时间作为模型的解在许多情况下，人们用时间序列的观测时期代表的时间作为模型的解释变量，用来表示被解释变量随时间的自发变化趋势。这种变量称为释变量，用来表示被解释变量随时间的自发变化趋势。这种变量称为时间变量，也叫做趋势变量。时间变量，也叫做趋势变量。n 时间变量通常用时间变量通常用t表示，其在用时间序列构建的计量经济模型中得到表示，其在用时间序列构建的计量经济模型中得到广泛的应用，它可

2、以单独作为一元线性回归模型中的解释变量，也可广泛的应用，它可以单独作为一元线性回归模型中的解释变量，也可以作多元线性回归模型中的一个解释变量，其对应的回归系数表示被以作多元线性回归模型中的一个解释变量，其对应的回归系数表示被解释变量随时间变化的变化趋势，时间变量也经常用在预测模型中。解释变量随时间变化的变化趋势，时间变量也经常用在预测模型中。Page 3STATA从入门到精通从入门到精通11.1.1 定义时间序列定义时间序列在在stata中的实现中的实现n 在进行时间序列的分析之前，首先要定义变量为时间序列数据。只有定义之在进行时间序列的分析之前，首先要定义变量为时间序列数据。只有定义之后，才

3、能对变量使用时间序列运算符号，也才能使用时间序列分析的相关命后，才能对变量使用时间序列运算符号，也才能使用时间序列分析的相关命令。定义时间序列用令。定义时间序列用tsset命令，其基本命令格式为：命令，其基本命令格式为：n tsset timevar , options n 其中，其中， timevar为时间变量。为时间变量。Options分为两类，或者定义时间单位，或者分为两类，或者定义时间单位，或者定义时间周期（即定义时间周期（即timevar两个观测值之间的周期数）。两个观测值之间的周期数）。Options的相关描述的相关描述如表如表11-1所示。所示。Page 4STATA从入门到精通

4、从入门到精通n注：（注：（1）units表示时间单位，对于表示时间单位，对于%tc，允许的时间单位包括：，允许的时间单位包括：second、seconds、secs、secs、minutes、minute、mine、min、hours、hour、days、weeks、week。对于其他。对于其他%t的格式，的格式，Stata自动自动获得其时间单位，获得其时间单位，delta选项经常与选项经常与%tc格式一起使用。格式一起使用。时间单位时间单位格式说明格式说明Clocktimetimevar的格式为%tc,0=1jan1960 00:00:00.000，1=1jan1960 00:00:00.0

5、01即0代表1960年1月1日的第一秒，1为1960年1月1日的第二秒，依次后推。dailytimevar的格式为%td，0=1jan1960，1=2jan1960；即0为1960年第一天，1为1960年第二天，依次后推。weeklytimevar的格式为%tw，0=1960w1，1=1960w2；即0为1960年第一周，1为1960年第二周，依次后推。monthlytimevar的格式为%tm，0=1，1=；即0为1960年第一月，1为1960年第二月，依次后推。quarterlytimevar的格式为%tq，0=1960q1，1=1960q2；即0为1960年第一季，1为1960年第二季，

6、依次后推。harfyearlytimevar的格式为%th，0=1960h1，1=1960h2；即0为从1960起的第一个半年，1为从1960年起第二个半年，依次后推。yearlytimevar的格式为%ty，1960=1960，1961=1960generictimevar的格式为%tgformat(%fmt)用户定义的其他时间周期时间周期 例子例子delta(#)例如delta(1)或delta(2)delta(exp)例如delta(7*24)delta(#units)例如delta(7 days)或delta(15 minutes)或delta(7 days 15 minutes)。见

7、注（1）delta(exp)units)例如delta(2+3) weeks)Page 5STATA从入门到精通从入门到精通n 可以通过以下三种方式来定义时间序列。例如，想要生成格式为可以通过以下三种方式来定义时间序列。例如，想要生成格式为%td的时间序列，并定义该时间序列为的时间序列，并定义该时间序列为t，则可以用以下三种方法：，则可以用以下三种方法： 方法方法1 方法方法2 方法方法3format t %td tsset ttsset t,dailytsset t, format(%td)Page 6STATA从入门到精通从入门到精通n 【例【例11.1】使用】使用文件文件“cpi.dta

8、”的数据来对的数据来对tsset命令的应用进行说命令的应用进行说明。该例子是我国明。该例子是我国1983年年1月年至月年至2007年年8月的居民消费价格指数月的居民消费价格指数CPI。部分数据如表。部分数据如表11-2所示：所示：n 表表11-2 我国居民消费价格指数我国居民消费价格指数CPIYear Year monthmonthcpicpi19831100.619832100.919833100.919834100.419835101.219836101.919837100.9Page 7STATA从入门到精通从入门到精通11.1.2 对时间序列进行修匀对时间序列进行修匀n 时间序列的形成

9、是各种不同的因素对事物的发展变化共同起作用的结时间序列的形成是各种不同的因素对事物的发展变化共同起作用的结果。这些因素概括起来可以归纳为四类：长期趋势因素、季节变动因果。这些因素概括起来可以归纳为四类：长期趋势因素、季节变动因素、循环变动因素和不规则变动因素。素、循环变动因素和不规则变动因素。n 时间序列构成分析就是要观察现象在一个相当长的时期内，由于各个时间序列构成分析就是要观察现象在一个相当长的时期内，由于各个影响因素的影响，使事物发展变化中出现的长期趋势、季节变动、循影响因素的影响，使事物发展变化中出现的长期趋势、季节变动、循环变动和不规则变动。环变动和不规则变动。n 通过测定和分析过去

10、一段时间之内现象的发展趋势，可以认识和掌握通过测定和分析过去一段时间之内现象的发展趋势，可以认识和掌握现象发展变化的规律性，为统计预测提供必要的条件，同时也可以消现象发展变化的规律性，为统计预测提供必要的条件，同时也可以消除原有时间序列中长期趋势的影响，更好地研究季节变动和循环变动除原有时间序列中长期趋势的影响，更好地研究季节变动和循环变动等问题。测定和分析长期趋势的主要方法是对时间序列进行修匀等问题。测定和分析长期趋势的主要方法是对时间序列进行修匀。Page 8STATA从入门到精通从入门到精通n 数据数据=修匀部分修匀部分+粗糙部分，运用粗糙部分，运用Stata进行修匀使用进行修匀使用ts

11、smooth命令，其基本命令，其基本命令格式如下所示：命令格式如下所示：n tssmooth smoothertype newvar = exp if in , .n 其中其中smoothertype有一系列目录，如下表有一系列目录，如下表11-4所示：所示：平滑的种类平滑的种类smoothertype移动平均不加权ma加权ma递归单指数过滤器exponential双指数过滤器dexponential非季节性Holt-Winters修匀hwinters季节性Holt-Winters修匀shwinters非线性过滤器nlPage 9STATA从入门到精通从入门到精通n 【例【例11.2】继续使用

12、】继续使用上例上例的数据来对的数据来对tssmooth命令的应用进行说明。命令的应用进行说明。在本例中对该组数据进行修匀，以便消除不规则变动的影响，得到时在本例中对该组数据进行修匀，以便消除不规则变动的影响，得到时间序列长期趋势，本例修匀的方法是利用之前的间序列长期趋势，本例修匀的方法是利用之前的1个月和之后的个月和之后的2个月个月及本月进行平均。及本月进行平均。Page 10STATA从入门到精通从入门到精通11.2 ARIMA模型的估计、单位根与协整模型的估计、单位根与协整n 时间序列模型一般分为四类，分别是自回归过程、移动平均过程、自时间序列模型一般分为四类，分别是自回归过程、移动平均过

13、程、自回归移动平均过程、单整自回归移动平均过程。回归移动平均过程、单整自回归移动平均过程。n 1、 自回归过程自回归过程n 如果一个剔出均值和确定性成分的线性过程可表达为如果一个剔出均值和确定性成分的线性过程可表达为n xt = 1xt-1 + 2 xt-2 + + p xt-p + utn 其中其中 i, i = 1, p 是自回归参数，是自回归参数，ut 是白噪声过程，则称是白噪声过程，则称xt为为p阶自回归过程，阶自回归过程，用用AR(p)表示。表示。xt是由它的是由它的p个滞后变量的加权和以及个滞后变量的加权和以及ut相加而成。相加而成。n 2、移动平均过程、移动平均过程n 如果一个剔

14、出均值和确定性成分的线性随机过程可用下式表达如果一个剔出均值和确定性成分的线性随机过程可用下式表达n xt = ut + 1 ut 1 + 2 ut -2 + + q ut q n 其中其中 1, 2, , q是回归参数，是回归参数，ut为白噪声过程，则上式称为为白噪声过程，则上式称为q阶移动平均过阶移动平均过程，记为程，记为MA(q) 。Page 11STATA从入门到精通从入门到精通n 3、自回归移动平均过程、自回归移动平均过程n 由自回归和移动平均两部分共同构成的随机过程称为自回归移动平均过程，由自回归和移动平均两部分共同构成的随机过程称为自回归移动平均过程，记为记为ARMA(p, q)

15、, 其中其中p, q分别表示自回归和移动平均部分的最大阶数。分别表示自回归和移动平均部分的最大阶数。ARMA(p, q) 的一般表达式是的一般表达式是n xt = 1xt-1 + 2xt-2 + p xt-p + ut + 1ut-1 + 2 ut-2 + .+ q ut-qn 4、单整自回归移动平均过程、单整自回归移动平均过程n 对于对于ARMA过程（包括过程（包括AR过程），如果特征方程过程），如果特征方程 (L) = 0 的全部根取值在单的全部根取值在单位圆之外，则该过程是平稳的；如果若干个或全部根取值在单位圆之内，则位圆之外，则该过程是平稳的；如果若干个或全部根取值在单位圆之内，则该过

16、程是强非平稳的。除此之外还有第三种情形，即特征方程的若干根取值该过程是强非平稳的。除此之外还有第三种情形，即特征方程的若干根取值恰好在单位圆上。这种根称为单位根，这种过程也是非平稳的。恰好在单位圆上。这种根称为单位根，这种过程也是非平稳的。n 若若随机过程随机过程yt 经过经过d 次差分之后可变换为一个以次差分之后可变换为一个以 (L)为为p阶自回归算子，阶自回归算子， (L)为为q阶移动平均算子的平稳、可逆的随机过程，则称阶移动平均算子的平稳、可逆的随机过程，则称yt 为（为（p, d, q）阶单）阶单整整(单积单积)自回归移动平均过程，记为自回归移动平均过程，记为ARIMA (p, d,

17、q)。Page 12STATA从入门到精通从入门到精通11.2.1 时间序列时间序列相关性相关性检验的检验的stata实现实现n 在进行在进行arima分析前，对序列的特征应该有相应的了解。包括自相关分析前，对序列的特征应该有相应的了解。包括自相关图，偏自相关图和图，偏自相关图和Q统计量。统计量。n 自相关刻画它序列自相关刻画它序列 的邻近数据之间存在多大程度的相关性。的邻近数据之间存在多大程度的相关性。n 偏自相关度量的是偏自相关度量的是k期间距的相关而不考虑期间距的相关而不考虑k -1期的相关。期的相关。n p阶滞后的阶滞后的Q-统计量的原假设是：序列不存在统计量的原假设是：序列不存在p阶

18、自相关；备选假设为：序列阶自相关；备选假设为：序列存在存在p阶自相关。阶自相关。n 在在Stata中实现相关性检验的基本命令格式如下所示：中实现相关性检验的基本命令格式如下所示：n 命令格式命令格式1（做出自相关和偏自相关图）：（做出自相关和偏自相关图）：n corrgram varname if in , corrgram_optionsn 命令格式命令格式2（做出自相关图）：（做出自相关图）：n ac varname if in , ac_optionsn 命令格式命令格式3（做偏自相关图）：（做偏自相关图）：n pac varname if in , pac_optionsPage 13

19、STATA从入门到精通从入门到精通n 以上三个命令格式的选项的相关描述分别如表以上三个命令格式的选项的相关描述分别如表11-5、11-6、11-7所示：所示：n表表11-5 corrgram_options的相关描述的相关描述 表表11-6 ac_options的相关描述的相关描述表表11-7 ac_options的相关描述的相关描述主要选项主要选项描述描述lags(#)*滞后阶数noplot不进行作图yw通过Yule-Walker方程组，计算偏自相关PAC主要选项主要选项描述描述lags(#)*滞后阶数generate(newvar)生成新变量，默认不做图level(#)置信度，默认95%f

20、ft通过傅里叶转化计算AC主要选项主要选项描述描述lags(#)*滞后阶数generate(newvar)level(#)生成新变量，默认不做图置信度，默认95%yw通过Yule-Walker方程组，计算偏自相关PACPage 14STATA从入门到精通从入门到精通n 【例【例11.3】使用表】使用表11-8的数据来对的数据来对Stata中自相关与偏自相关的应用进行说明。中自相关与偏自相关的应用进行说明。该数据给出了中国该数据给出了中国1953-1984年的国民生产总值年的国民生产总值GNP、私人国内总投资、私人国内总投资I、GNP的隐性价格折算因子的隐性价格折算因子P（以（以1972为基期）

21、、半年期商业票据利率为基期）、半年期商业票据利率R。在。在本例中我们对本例中我们对GNP时间序列进行分析，观察期相关图和自相关图，从而得到时间序列进行分析，观察期相关图和自相关图，从而得到GNP时间序列的类型。时间序列的类型。部分部分数据说明数据说明下下表所示。表所示。年份年份中国中国GNP私人国内总私人国内总投资投资GNPGNP的隐性的隐性价格折算因价格折算因子（子（1972=1）半年期商业半年期商业票据利率票据利率1953623.685.30.5882.521954616.183.10.5961.591955657.5103.80.6082.191956671.6102.60.6283.3

22、11957683.8970.6493.821958680.987.50.662.471959721.71080.6763.96Page 15STATA从入门到精通从入门到精通11.2.2 时间序列稳定时间序列稳定性性检验的检验的stata实现实现n 检验序列的平稳性，可以用检验序列的平稳性，可以用phillips-perron检验，检验，dickey-fuller检验，检验，以及应用以及应用GLS扩展的扩展的dickey-fuller检验。其基本命令格式如下：检验。其基本命令格式如下：n 命令格式命令格式1（dickey-fuller检验）：检验）：n dfuller varname if i

23、n ,optionn 命令格式命令格式2（GLS扩展的扩展的dickey-fuller检验）：检验）：n dfgls varname if in , optionsn 命令格式命令格式3（phillips-perron检验）：检验）：n pperron varname if in , optionsn 以上三个命令格式的选项的相关描述分别如表以上三个命令格式的选项的相关描述分别如表11-10、11-11、11-12所示所示：Page 16STATA从入门到精通从入门到精通n 表表11-10 dickey-fuller检验检验options的相关描述的相关描述n 表表11-11 GLS扩展的扩展

24、的dickey-fuller检验检验options的相关描述的相关描述n 表表11-12 phillips-perron检验检验检验检验options的相关描述的相关描述主要选项主要选项描述描述noconstant没有截据项trend包括时间趋势drift包括漂移项regress 显示回归结果lags(#) 滞后阶数主要选项主要选项描述描述maxlag(#)最大滞后阶数notrend没有时间趋势ers利用插值法计算临界值主要选项主要选项描述描述noconstant没有截据项trendregress 有趋势项显示回归结果lags(#)最大滞后阶数Page 17STATA从入门到精通从入门到精通n

25、 【例【例11.4】继续使用】继续使用上例的上例的数据来对数据来对Stata中平稳性检验的相关应用中平稳性检验的相关应用进行说明。这里要求使用进行说明。这里要求使用dickey-fuller检验、检验、GLS扩展的扩展的dickey-fuller检验和检验和phillips-perron检验三种方法，对检验三种方法，对GNP的一阶差分进行的一阶差分进行平稳性检验。平稳性检验。Page 18STATA从入门到精通从入门到精通11.2.3 ARIMA模型模型的的stata实现实现n 时间序列的自回归移动平均法可是通过使用时间序列的自回归移动平均法可是通过使用arima命令来实现。其基本命命令来实现

26、。其基本命令格式如下：令格式如下：n arima depvar indepvars if in weight , optionsn 在使用在使用arima模型前，需要先检验数据的平稳性和相关性，然后经过判断才能使用。模型前，需要先检验数据的平稳性和相关性，然后经过判断才能使用。主要选项主要选项描述描述noconstant没有截据项Arima(#p,#d,#q)Arima(p,d,q)模型Ar(numlist)Ar的滞后阶数Ma(numlist) Ma的滞后阶数Constraints(constraints)线性约束collinear保留多重共线性变量Sarima(#p,#d,#q,#s)季节a

27、rima模型Mar(numlist,#s)季节ar的滞后阶数Mma( numlist,#s)季节ma的滞后阶数Page 19STATA从入门到精通从入门到精通n 【例【例11.5】使用表】使用表11-14的数据来对的数据来对Stata中中ARIMA模型的相关应用进模型的相关应用进行说明。该表给出了某地区每年的年度总人口数。行说明。该表给出了某地区每年的年度总人口数。部分部分数据如数据如下下：年份年份年底总人口数年底总人口数(万人万人)194954167195055196195156300195257482195358796195460266195561465195662828195764653

28、195865994195967207Page 20STATA从入门到精通从入门到精通11.3 VAR与与VEC模型模型的估计及解释的估计及解释n 1、VAR模型的阶数选择模型的阶数选择n 在在Stata中中VAR模型阶数选择的实现，是通过如下基本命令来实现的：模型阶数选择的实现，是通过如下基本命令来实现的：n depvarlist if in , preestimation_options主要选项主要选项描述描述maxlag(#)最高滞后阶数; 默认是滞后4期exog(varlist)外生变量constraints(constraints)对外生变量的线性约束noconstant 没有常数项l

29、evel(#) 置信度，默认95%separator(#)分割线Page 21STATA从入门到精通从入门到精通n 2、构建、构建VAR模型模型n 在在Stata中构建中构建VAR模型的实现，是通过如下基本命令来实现的：模型的实现，是通过如下基本命令来实现的：n var depvarlist if in , options主要选项主要选项描述描述模型1noconstant 没有常数项lags(numlist)VAR滞后阶数 exog(varlist) 外生变量模型2 constraints(numlist)线性约束 nolog 不显示迭代过程 noisure一步迭代dfk自由度调节small小

30、样本t,f统计量报告结果 level(#)置信度Page 22STATA从入门到精通从入门到精通n 3、平稳性条件考察、平稳性条件考察n 在在Stata中中VAR模型平稳性条件考察的实现，是通过如下基本命令来实现的：模型平稳性条件考察的实现，是通过如下基本命令来实现的：n varstable , options主要选项主要选项描述描述estimates(estname)考察VAR（estname）的平稳性graph对伴随矩阵的特征值作图dlabel将特征值标记为到单位圆的距离Page 23STATA从入门到精通从入门到精通n 4、残差的正态性和自相关检验、残差的正态性和自相关检验n 在在Sta

31、ta中中VAR模型残差的正态性和自相关检验的实现，是通过如下基本命令模型残差的正态性和自相关检验的实现，是通过如下基本命令来实现的：来实现的：n varnorm , options主要选项主要选项描述描述jbera statistics Jarque-Bera 统计量skewness偏度kurtosis峰度estimates(estname)cholesky 已估计的var名称使用Cholesky 分解separator(#)分割线Page 24STATA从入门到精通从入门到精通n 5、格兰杰因果检验、格兰杰因果检验n 在在Stata中中VAR模型格兰杰因果检验的实现，是通过如下基本命令来模型

32、格兰杰因果检验的实现，是通过如下基本命令来实现的：实现的：n vargranger , estimates(estname) separator(#)Page 25STATA从入门到精通从入门到精通n 6、脉冲分析、脉冲分析n（1）irf文件的创建、显示、激活和清除文件的创建、显示、激活和清除nVAR模型脉冲分析的实现，首先是要创建模型脉冲分析的实现，首先是要创建irf文件。在文件。在Stata中是通过如下基本命令来实现的：中是通过如下基本命令来实现的：n命令格式命令格式1（VAR模型的模型的irf文件创建）：文件创建）：irf create irfname , var_optionsn命令格

33、式命令格式2（SVAR模型的模型的irf文件创建）：文件创建）：irf create irfname , svar_optionsn命令格式命令格式3（VEC模型的模型的irf文件创建）：文件创建）：irf create irfname , vec_optionsn创建创建irf文件之后，显示处于当下活动状态的文件之后，显示处于当下活动状态的irf，输入以下命令：，输入以下命令：irf setn激活激活irf文件，可以输入以下命令：文件，可以输入以下命令：irf set ifr_namen清除活动的清除活动的irf文件，可以输入以下命令：文件，可以输入以下命令：irf set, clear主要

34、选项主要选项描述描述set(filename, replace)创建文件replace如果文件已存在，则替换文件order(varlist)Cholesky排序estimates(estname) 以估计的VAR名称Page 26STATA从入门到精通从入门到精通n （2）Irf作图作图n Irf文件作图，可以输入以下命令：文件作图，可以输入以下命令：irf graph stat , optionsn stat的相关描述的相关描述 options的相关描述的相关描述主要选项主要选项描述描述irf irfoirf正交irfdm动态乘子cirf 累计irfcoirf 累计正交irfcdm累计同台乘

35、子fevdCholesky 方差分解sirf结构IRFsfevd结构 Cholesky 方差分解主要选项主要选项描述描述set(filename) 使文件激活irf(irfnames)IRF 结果名称impulse(impulsevar)脉冲变量response(endogvars)响应变量Page 27STATA从入门到精通从入门到精通n 6 johansen检验检验n 当变量之间同阶单整时，可以运用当变量之间同阶单整时，可以运用johansen检验查看变量之间是否协整。检验查看变量之间是否协整。Stata中中VAR模型模型johansen检验的实现，是通过如下基本命令来实现的：检验的实现，

36、是通过如下基本命令来实现的：n vecrank depvar if in , options 主要选项主要选项描述描述lags(#) VAR模型的最高滞后阶数trend(constant)VAR模型有常数项，协整方程有常数项trend(rconstant)VAR模型有常数项，协整方程无常数项trend(trend)VAR模型有趋势项，协整方程有趋势项trend(rtrend)VAR模型有趋势项，协整方程无趋势项trend(none)VAR模型无常数项，协整方程无常数项Page 28STATA从入门到精通从入门到精通n 【例【例11.6】表】表11-10给出了我国给出了我国CPI、利率、利率R、

37、狭义货币供应量、狭义货币供应量M1经过修匀后的数据。经过修匀后的数据。其中狭义货币供应量增长率经过其中狭义货币供应量增长率经过SAR修匀后记为修匀后记为M1sar，贷款利率记为，贷款利率记为r,cpi经过经过sa修匀后记为修匀后记为cpisa。数据区间是从。数据区间是从1994年年1月月2007年年12月。本例中将要建立一个关月。本例中将要建立一个关于变量于变量m1sar 、变量、变量cpisa和变量和变量r的的VAR模型，部分数据如表模型，部分数据如表11-23所示：所示：monthyearm1sarcpisar119940.190123392 20.93511929 12.24219940

38、.166035575 23.36645208 12.24319940.154944021 22.50509823 12.24419940.145357052 21.894884512.24519940.155072278 21.58277167 9619940.181137027 22.78261176 9719940.233693845 24.00026113 9819940.284786294 25.76049934 9919940.291513079 27.16382803 9Page 29STATA从入门到精通从入门到精通11.4 ARCH与与GARCH的估计及解释的估计及解释n 1、

39、ARCH模型模型n 若一个平稳随机变量若一个平稳随机变量xt可以表示为可以表示为AR(p) 形式，其随机误差项的方差可用误形式，其随机误差项的方差可用误差项平方的差项平方的q阶分布滞后模型描述，阶分布滞后模型描述，n xt = 0 + 1 xt -1 + 2 xt -2 + + p xt - p + utn t2 = E(ut2) = 0 + 1 ut -1 2 + 2 ut -22 + + q ut - q2n 则称则称ut 服从服从q阶的阶的ARCH过程，记作过程，记作ut ARCH (q)。其中。其中第一第一式称作均值方式称作均值方程，程，第二第二式称作式称作ARCH方程。方程。n 2、

40、GRACH模型模型n ARCH模型中的第二式模型中的第二式是关于是关于 t2的分布滞后模型。为避免的分布滞后模型。为避免ut2的滞后项过多，的滞后项过多，可采用加入可采用加入 t2的滞后项的方法（回忆可逆性概念）。对于的滞后项的方法（回忆可逆性概念）。对于第二第二式，可给出如式，可给出如下形式，下形式，n t2 = 0 + 1 ut 1 2 + 1 t -12n 此模型称为广义自回归条件异方差模型，用此模型称为广义自回归条件异方差模型，用GARCH (1, 1) 表示。其中表示。其中ut 1称称为为ARCH项，项， t -1称为称为GARCH项。项。Page 30STATA从入门到精通从入门到

41、精通n 在在Stata中中ARCH模型的实现，是通过如下基本命令来实现的：模型的实现，是通过如下基本命令来实现的：n arch depvar indepvars if in weight , optionsModel noconstant没有常数项 arch(numlist)ARCH 滞后阶数 garch(numlist)GARCH 滞后阶数 saarch(numlist)简单非对称 ARCH 模型 tarch(numlist)门限ARCH 模型 aarch(numlist)非对称 ARCH模型 narch(numlist)非线性ARCH模型 narchk(numlist)带有位移的非线性AR

42、CH模型 abarch(numlist)绝对值ARCH 模型 atarch(numlist)绝对门限ARCH模型 sdgarch(numlist)garch项的滞后项 earch(numlist)Nelsons EGARCH 模型的信息项egarch(numlist)log（garch）的滞后项 parch(numlist)幂ARCH 模型 tparch(numlist)门限幂 ARCH模型 aparch(numlist)非对称幂ARCH 模型 nparch(numlist)非线性幂ARCH 模型 nparchk(numlist)带有位移的非线性幂 ARCH 模型 pgarch(numlist)幂 GARCH 模型 constraints(constraints)线性约束Model 2archm均值方程加入方差项archmlags(numlist)均值方程加入滞后阶数 archmexp(exp)将exp转换为ARCH-IN-MEAN的形式arima(#p, #d, #q)ARIMA(p,d,q) 模型ar(numlist)ar模型 ma(numlist)ma模型Model 3 het(varlist)条件方差估计中带有外生变