关于时间序列分析Word文档格式.docx

关于时间序列分析Word文档格式.docx

《关于时间序列分析Word文档格式.docx》由会员分享，可在线阅读，更多相关《关于时间序列分析Word文档格式.docx（17页珍藏版）》请在冰豆网上搜索。

（用因变量的散点图和直方图及其包含的正态分布检验随机性，大多数服从正态分布。

）

（2）平稳性：

样本序列的自相关函数在某一固定水平线附近摆动，即方差和数学期望稳定为常数。

样本序列的自相关函数只是时间间隔的函数，与时间起点无关。

其具有对称性，能反映平稳序列的周期性变化。

特征识别利用自相关函数ACF：

ρk=γk/γ0

其中γk是yt的k阶自协方差，且ρ0=1、-1<

ρk<

1。

平稳过程的自相关系数和偏自相关系数都会以某种方式衰减趋近于0，前者测度当前序列与先前序列之间简单和常规的相关程度，后者是在控制其它先前序列的影响后，测度当前序列与某一先前序列之间的相关程度。

实际上，预测模型大都难以满足这些条件，现实的经济、金融、商业等序列都是非稳定的，但通过数据处理可以变换为平稳的。

4.预测类型

（1）点预测：

确定唯一的最好预测数值，其给出了时间序列未来发展趋势的一个简单、直接的结果。

但常产生一个非零的预测误差，其不确定程度为点预测值的置信区间。

（2）区间预测：

未来预测值的一个区间，即期望序列的实际值以某一概率落入该区间范围内。

区间的长度传递了预测不确定性的程度，区间的中点为点预测值。

（3）密度预测：

序列未来预测值的一个完整的概率分布。

根据密度预测，可建立任意置信水平的区间预测，但需要额外的假设和涉及复杂的计算方法。

5.基本步骤

（1）分析数据序列的变化特征。

（2）选择模型形式和参数检验。

（3）利用模型进行趋势预测。

（4）评估预测结果并修正模型。

3.3.2随机时间序列

系统中某一因素变量的时间序列数据没有确定的变化形式，也不能用时间的确定函数描述，但可以用概率统计方法寻求比较合适的随机模型近似反映其变化规律。

（自变量不直接含有时间变量,但隐含时间因素）

1．自回归AR（p）模型

（R：

模型的名称P：

模型的参数）（自己影响自己，但可能存在误差，误差即没有考虑到的因素）

（1）模型形式（εt越小越好，但不能为0：

ε为0表示只受以前Y的历史的影响不受其他因素影响）

yt=φ1yt-1+φ2yt-2+……+φpyt-p+εt

式中假设：

yt的变化主要与时间序列的历史数据有关，与其它因素无关；

εt不同时刻互不相关，εt与yt历史序列不相关。

式中符号：

p模型的阶次，滞后的时间周期，通过实验和参数确定；

yt当前预测值，与自身过去观测值yt-1、…、yt-p是同一序列不同时刻的随机变量，相互间有线性关系，也反映时间滞后关系；

yt-1、yt-2、……、yt-p同一平稳序列过去p个时期的观测值；

φ1、φ2、……、φp自回归系数，通过计算得出的权数，表达yt依赖于过去的程度，且这种依赖关系恒定不变；

εt随机干扰误差项，是0均值、常方差σ2、独立的白噪声序列，通过估计指定的模型获得。

（2）识别条件

当k>

p时，有φk=0或φk服从渐近正态分布N（0,1/n）且（|φk|>

2/n1/2）的个数≤4.5%，即平稳时间序列的偏相关系数φk为p步截尾，自相关系数rk逐步衰减而不截尾，则序列是AR（p）模型。

实际中，一般AR过程的ACF函数呈单边递减或阻尼振荡，所以用PACF函数判别（从p阶开始的所有偏自相关系数均为0）。

（3）平稳条件

一阶：

|φ1|<

二阶：

φ1+φ2<

1、φ1-φ2<

1、|φ2|<

φ越大，自回归过程的波动影响越持久。

（4）模型意义

仅通过时间序列变量的自身历史观测值来反映有关因素对预测目标的影响和作用，不受模型变量相互独立的假设条件约束，所构成的模型可以消除普通回归预测方法中由于自变量选择、多重共线性等造成的困难。

2．移动平均MA（q）模型

（1）模型形式

yt=εt-θ1εt-1-θ2εt-2-……-θpεt-p

（2）模型含义

用过去各个时期的随机干扰或预测误差的线性组合来表达当前预测值。

AR（p）的假设条件不满足时可以考虑用此形式。

总满足平稳条件，因其中参数θ取值对时间序列的影响没有AR模型中参数p的影响强烈，即这里较大的随机变化不会改变时间序列的方向。

（3）识别条件

q时，有自相关系数rk=0或自相关系数rk服从N（0,1/n（1+2∑r2i）1/2）且（|rk|>

2/n1/2（1+2∑r2i）1/2）的个数≤4.5%，即平稳时间序列的自相关系数rk为q步截尾，偏相关系数φk逐步衰减而不截尾，则序列是MA（q）模型。

实际中，一般MA过程的PACF函数呈单边递减或阻尼振荡，所以用ACF函数判别（从q阶开始的所有自相关系数均为0）。

（4）可逆条件

|θ1|<

|θ2|<

1、θ1+θ2<

当满足可逆条件时，MA（q）模型可以转换为AR（p）模型

3．自回归移动平均ARMA（p,q）模型

（1）模型形式

yt=φ1yt-1+φ2yt-2+……+φpyt-p+εt-θ1εt-1-θ2εt-2-……-θpεt-p

p和q是模型的自回归阶数和移动平均阶数；

φ和θ是不为零的待定系数；

εt独立的误差项；

yt是平稳、正态、零均值的时间序列。

（2）模型含义

使用两个多项式的比率近似一个较长的AR多项式，即其中p+q个数比AR（p）模型中阶数p小。

前二种模型分别是该种模型的特例。

一个ARMA过程可能是AR与MA过程、几个AR过程、AR与ARMA过程的迭加，也可能是测度误差较大的AR过程。

（3）识别条件

平稳时间序列的偏相关系数φk和自相关系数rk均不截尾，但较快收敛到0，则该时间序列可能是ARMA（p,q）模型。

实际问题中，多数要用此模型。

因此建模解模的主要工作是求解p、q和φ、θ的值，检验εt和yt的值。

（4）模型阶数

AIC准则：

最小信息准则，同时给出ARMA模型阶数和参数的最佳估计，适用于样本数据较少的问题。

目的是判断预测目标的发展过程与哪一随机过程最为接近。

因为只有当样本量足够大时，样本的自相关函数才非常接近母体的自相关函数。

具体运用时，在规定范围内使模型阶数从低到高，分别计算AIC值，最后确定使其值最小的阶数是模型的合适阶数。

模型参数最大似然估计时AIC=（n-d）logσ2+2（p+q+2）

模型参数最小二乘估计时AIC=nlogσ2+（p+q+1）logn

式中：

n为样本数，σ2为拟合残差平方和，d、p、q为参数。

其中：

p、q范围上线是n较小时取n的比例，n较大时取logn的倍数。

实际应用中p、q一般不超过2。

4．自回归综合移动平均ARIMA（p,d,q）模型

（1）模型识别

平稳时间序列的偏相关系数φk和自相关系数rk均不截尾，且缓慢衰减收敛，则该时间序列可能是ARIMA（p,d,q）模型。

模型形式类似ARMA（p,q）模型，但数据必须经过特殊处理。

特别当线性时间序列非平稳时，不能直接利用ARMA（p,q）模型，但可以利用有限阶差分使非平稳时间序列平稳化，实际应用中d一般不超过2。

若时间序列存在周期性波动，则可按时间周期进行差分，目的是将随机误差有长久影响的时间序列变成仅有暂时影响的时间序列。

即差分处理后新序列符合ARMA（p,q）模型，原序列符合ARIMA（p,d,q）模型。

3.3.3建模解模过程

1.数据检验

检验时间序列样本的平稳性、正态性、周期性、零均值，进行必要的数据处理变换。

（1）作直方图：

检验正态性、零均值。

按图形Graphs—直方图Histogram的顺序打开如图3.15所示的对话框。

图3.15

将样本数据送入变量Variable框，选中显示正态曲线Displaynormalcurve项，点击OK运行，输出带正态曲线的直方图，如图3.16所示。

图3.16

从图中看出：

标准差不为1、均值近似为0，可能需要进行数据变换。

（2）作相关图：

检验平稳性、周期性。

按图形Graphs—时间序列TimeSeries—自相关Autocorrelations的顺序打开如图3.17所示的对话框。

图3.17

将样本数据送入变量Variable框，选中自相关Autocorrelations和偏自相关PartialAutocorrelations项，暂不选数据转换Transform项，点击设置项Options，出现如图3.18所示对话框。

图3.18

因为一般要求时间序列样本数据n>

50，滞后周期k<

n/4，所以此处控制最大滞后数值MaximumNumberofLags设定为12。

点击继续Continue返回自相关主对话框后，点击OK运行系统，输出自相关图如图3.19所示。

图3.19

从图中看出；

样本序列数据的自相关系数在某一固定水平线附近摆动，且按周期性逐渐衰减，所以该时间序列基本是平稳的。

（3）数据变换：

若时间序列的正态性或平稳性不够好，则需进行数据变换。

常用有差分变换（利用transform—CreateTimeSeries）和对数变换（利用Transform—Compute）进行。

一般需反复变换、比较，直到数据序列的正态性、平稳性等达到相对最佳。

2.模型识别

分析时间序列样本，判别模型的形式类型，确定p、d、q的阶数。

（1）判别模型形式和阶数

①相关图法：

运行自相关图后，出现自相关图（图3.19）和偏自相关图（图3.20）。

图3.20

自相关系数和偏相关系数具有相似的衰减特点：

衰减快，相邻二个值的相关系数约为0.42，滞后二个周期的值的相关系数接近0.1，滞后三个周期的值的相关系数接近0.03。

所以，基本可以确定该时间序列为ARMA（p,q）模型形式，但还不能确定是ARMA（1,1）或是ARMA（2,2）模型。

但若前四个自相关系数分别为0.40、0.16、0.064、0.0256，则可以考虑用AR

（1）模型。

另外，值得说明的是：

只是ARMA模型需要检验时间序列的平稳性，若该序列的偏自相关函数具有显著性，则可以直接选择使用AR模型。

实际上，具体应用自相关图进行模型选择时，在观察ACF与PACF函数中，应注意的关键问题是：

函数值衰减的是否快；

是否所有ACF之和为-0.5，即进行了过度差分；

是否ACF与PACF的某些滞后项显著和容易解释的峰值等。

但是，仅依赖ACF图形进行时间序列的模型识别是比较困难的。

②参数估计：

从（m,m-1）开始试验，一般到m=p+q=1/n。

实际应用中，往往从（1,1）、……、（2,2），逐个计算比较它们的AIC值（或SBC值），取其值最小的确定为模型。

（2）建立时间序列新变量

无论是哪种模型形式，时间序列总是受自身历史数据序列变化的影响，因此需将历史数据序列作为一个新的时间序列变量。

按数据转换transform—建立时间序列CreateTimeSeries的顺序展开对话框，图3.21。

图3.21

①在功能Function下拉框中选择变量转换的函数，其中：

非季节差分Differences:

计算时间序列连续值之间的非季节性差异。

季节性差分SeasonalDifferences:

计算时间序列跨距间隔恒定值之间的季节性差异，跨距根据定义的周期确定。

领先移动平均Priormovingaverage:

计算先前的时间序列数值的平均值。

中心移动平均Centeredmovingaverage:

计算围绕和包括当前值的时间序列数值的平均值。

中位数Runningmedians:

计算围绕和包括当前值的时间序列的中位数。

累积和Cumulativesum:

计算直到包括当前值的时间序列数值的累计总数。

滞后顺序Lag:

根据指定的滞后顺序，计算在前观测量的值。

领先顺序Lead:

根据指定的领先顺序，计算连续观测量的值。

平滑Smoothing:

以混合数据平滑为基础，计算连续观测量的值。

以上各项主要用在生成差分变量、滞后变量、平移变量，并且还要关注差分、滞后、平移的次数，以便在建立模型、进行参数估计时，使方程达到一致。

②在顺序Order框中填入在前或在后的时间序列数值间隔的数目。

在新变量NewVariable框中接受左边框移来的源变量。

在名称Name框中定义新变量的名称，但必单击改变Change方能成立。

③单击OK运行系统，在原数据库中出现新变量列。

另外，若需产生周期性时间序列的日期型变量，则按数据Data—定义日期DefineDates的顺序展开如图3.22所示对话框。

图3.22

在样本CasesAre栏中选择定义日期变量的时间间隔，在起始日期FirstCaseIs栏中设定日期变量第一个观测量的值，单击OK完成定义。

3.参数估计

采用最大似然估计或最小二乘估计等方法估计φ、θ参数值，并进行显著性检验。

按分析Analyze—时间序列Timeseries—ARIMA模型的顺序展开如图3.23对话框。

图3.23

在图3.23中：

选择原时间序列变量进入因变量框；

根据模型识别结果和建立的新时间变量，选择一个或多个变量进入自变量框；

暂时不进行因变量的数据转换；

与自变量的选择对应，根据模型识别结果或实验的思路设定p、（d）、q的值；

选择模型中包含常数项；

分别单击保存和设置按钮，展开如图3.24和3.25对话框。

图3.24

图3.24中：

在建立变量CreateVariable栏选择新建变量结果暂存原数据文件Addtofile项，也可选择用新建变量代替原数据文件中计算结果Replaceexisting项；

在设定置信区间百分比%ConfidenceIntervals下拉框选择95；

在预测样本PredictCases栏选择根据时期给出预测结果的方法。

图3.25

图3.25中：

在收敛标准ConvergenceCriteria栏选择迭代次数Maximumiterations、参数变化精度Parameterchange、平方和变化精度Sumofsquareschange，当运算达到其中一个参数的设定，则迭代终止；

在估计初始值InitialValuesforEstimation栏选择由过程自动选择Automatic或由先前模型提供Applyfrompreviousmodel，一般默认前者；

在预测方法ForecastingMethod栏选择无条件Unconditional或有条件最小二乘法Conditionalleastsquares；

在输出控制Display栏选择最初和最终参数的迭代摘要Initialandfinalparameterswithiterationsummary或详细资料details、或只显示最终参数Finalparametersonly。

单击OK，系统立即执行，输出信息如下：

MODEL:

MOD_1

Splitgroupnumber:

1Serieslength:

48

Nomissingdata.

Melard'

salgorithmwillbeusedforestimation.

Conclusionofestimationphase.

Estimationterminatedatiterationnumber7because:

Sumofsquaresdecreasedbylessthan.001percent.

FINALPARAMETERS:

Numberofresiduals48

Standarderror1.1996949

Loglikelihood-75.463915

AIC156.92783

SBC162.54143

AnalysisofVariance:

DFAdj.SumofSquaresResidualVariance

Residuals4565.0999231.4392678

VariablesintheModel:

BSEBT-RATIOAPPROX.PROB.

AR1.02318739.31945836.0725835.94245925

MA1-.44871554.28829314-1.5564558.12660552

CONSTANT-.02421308.25505018-.0949346.92478827

Thefollowingnewvariablesarebeingcreated:

NameLabel

FIT_1Fitfor样本数据fromARIMA,MOD_1CON

ERR_1Errorfor样本数据fromARIMA,MOD_1CON

LCL_195%LCLfor样本数据fromARIMA,MOD_1CON

UCL_195%UCLfor样本数据fromARIMA,MOD_1CON

SEP_1SEoffitfor样本数据fromARIMA,MOD_1CON

各个输出统计量的意义：

常数项：

认为是取值恒为1的常数变量，其系数就是自变量为0时因变量的最优预测值，也称为预测基准值。

系数：

反映自变量对因变量影响的权重。

标准误：

表明样本数据的可靠性。

在（残差）参数近似服从正态分布条件下，系数加减两倍的标准误差近似等于总体参数95%的置信区间。

其值越小，置信区间越窄；

并且其对于系数的相对值越小，估计结果越精确。

t统计量：

估计系数与标准误差的比值，检验变量的不相关性。

一般给定5%显著水平，则拒绝原假设的0值位于95%的置信区间外，其绝对值必大于2。

t概率值：

其值越小，则拒绝原假设不相关性的证据越充分。

其值接近0.05与t统计量接近2相对应。

均值：

度量变量的集中度，传递随机变量的位置信息。

标准差：

度量变量的离散度，传递随机变量的规模信息。

平方和：

残差平方和是许多统计量的组成部分，孤立考察无太大价值。

准则：

信息准则AIC和SBC用于模型的选择，越小越好，但受自由度约束较为严重。

R2校正：

是模型中自变量对因变量变动的解释比例，度量方程预测因变量的成功程度，其是回归标准误差与因变量标准差比较的结果。

另一个比较方法是回归标准误差不超过因变量均值的10%则为好的模型。

DW统计：

用于检验随机误差项是否存在序列相关。

LN似然：

用于模型比较和假设检验，越大越好。

残差图：

4.模型检验

检验新建模型的合理性。

若检验不通过，则调整（p,q）值，重新估计参数和检验，反复进行直到接受为止。

但模型识别、参数估计、检验修正三个过程之间相互作用、相互影响，有时需要交叉进行、反复实验，才能最终确定模型形式。

（1）相关图检验残差白噪声：

因为白噪声过程是序列无关的，所以白噪声过程的自相关函数和偏自相关函数在自相关图中均为等于0的水平直线。

（2）散点图检验残差独立性：

以误差值为纵坐标、以预测值为横坐标，观察散点分布的均匀性、随机性。

理想预测模型的预测误差一定是不可预测的、无规律的、序列无关的。

相应的DW统计量仅适用检验一阶序列。

（3）直方图检验残差零均值：

零均值仅检验残差序列无关，若正态分布则检验独立性。

（4）概率图检验残差自相关：

以显著性水平0.05计算χ2（）概率值，。

（5）均方差检验预测的效果：

以预测误差的均方差最小为标准，注意预测误差仅与预测周期有关，而与起始时刻无关。

5.模型预测

预测系统研究对象的未来某时刻状态。

列出预测模型，计算预测值。