SAS学习系列37. 时间序列分析Ⅰ—平稳性及纯随机性检验

《SAS学习系列37. 时间序列分析Ⅰ—平稳性及纯随机性检验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、37.时间序列分析Ⅰ—平稳性及纯随机性检验（一）基本概念一、什么是时间序列？为了研究某一事件的规律，依据时间发生的顺序将事件在多个时刻的数值记录下来，就构成了一个时间序列。对时间序列进行观察、研究，找寻它变化发展的规律，预测它将来的发展趋势就是时间序列分析。例如，国家或地区的年度财政收入，股票市场的每日波动，气象变化，工厂按小时观测的产量等等。注：随温度、高度等变化而变化的离散序列，也可以看作时间序列。二、时间序列的特点（1）顺序性；（2）随机性；（3）前后时刻（不一定相邻）的依存性；（4）整体呈趋势性和周期性。三、时间序列

2、的分类按研究对象的数目：一元时间序列、多元时间序列；按序列统计特性：平稳时间序列、非平稳时间序列；按分布规律：高斯时间序列、非高斯时间序列。四、研究方法1.平稳时间序列分析；2.非平稳时间序列分析（确定性分析、随机性分析）。五、其它任何时间序列经过合理的函数变换后都可以被认为是由下列三部分叠加而成：（1）趋势项部分；（2）周期项部分；（3）随机项部分（随机信号、随机噪声）图1.四种趋势：线性、二次、指数增长、S型例如，手机销售的月记录按年增长（趋势项）；按季节周期波动（周期项）；随机信号和随机噪声。时间序列分析的主要任务就是

3、：上面三部分分解出来，是研究平稳随机过程的变化规律，建立特定的ARIMA模型（要求大体平稳、可能含有周期但不能有规则性的线性指数等类型趋势项）。六、方法性工具1.差分运算（1）k步差分间隔k期的观察值之差：Δk=xt-xt-k（2）p阶差分Δxt=xt-xt-1称为一阶差分；称为p阶差分；SAS函数实现：diffn(x)2.延迟算子延迟算子作用于时间序列，时间刻度减小1个单位（序列左移一位）：Bxt=xt-1,……,Bpxt=xt-p.SAS函数实现：lagn(x)用延迟算子表示k步差分和p阶差分为：Δk=xt-xt-k=(

4、1-Bk)xt（二）平稳时间序列一、概念平稳时间序列按限制条件的严格程度，分为严平稳时间序列：序列所有的统计性质都不会随着时间的推移而发生变化；宽平稳时间序列：序列的主要性质近似稳定，即统计性质只要保证序列的二阶矩平稳，即对任意的时间t，s，k，序列Xt满足：二、平稳时间序列的统计性质（1）均值为常数；（2）自协方差只依赖于时间跨度；若定义自协方差函数为γ(t,s)=E(Xt-μt)(Xs-μs)则可由二元函数简化为一元函数γ(t-s)，得延迟k自协方差函数：γ(k)=γ(t,t+k)由此易知平稳时间序列必具有常数方差：D(

5、Xt)=E(Xt-μt)2=γ(t,t)=γ(0)时间序列自相关函数：延迟k自相关函数：基本性质：（1）ρ(0)=1;（2）ρ(-k)=ρ(k);（3）自相关阵为对称负定阵；（4）非唯一性。注意：协方差函数和相关函数——度量两个不同事件（Xt，Yt）彼此之间的相互影响的程度。自协方差函数和自相关函数——度量用一事件（Xt）在两个不同时期之间的相互影响的程度。三、样本估计值总体均值的估计值：延迟k自协方差函数的估计值：总体方差的估计值：延迟k自相关函数的估计值：四、平稳性检验（1）时序图检验若无明显的趋势性和周期性，则平稳；（

6、2）自相关图检验零均值平稳序列的自相关函数要么截尾要么拖尾；若时间序列零均值化后出现缓慢衰减或周期性衰减，则说明存在趋势性和周期性（非平稳）；（3）单位根检验就是通过检验时间序列自回归特征方程的特征根是在单位圆内（平稳）还是在单位圆及单位圆外（非平稳）。通常用ADF检验法。Dickey和Fuller(1979)利用如下的广义自回归模型其中，Δxj,t表示x的一阶差分；xj,t-1表示延迟一期；Δxj,t-k表示延迟k期再一阶差分；εk,t表示扰动项。上述回归模型生成的xj,t-1的t值正好对应ADF统计量，做假设检验：H0:

7、非平稳；H1：平稳。t值在1%,5%,10%置信水平的临界值分别为：-3.524233,-2.902358,-2.588587.以此判断序列是否平稳。注：若Xt不平稳，可以依次对Xt做一阶、二阶…差分，直到序列平稳。例1.平稳性检验——ADF检验的SAS实现。代码：datasimulation;doi=1to100;x=rannor(1234);output;end;run;datatimeseries;setsimulation;x_1st_lag=lag1(x);x_1st_diff=dif1(x);x_1st_diff

8、_1st_lag=dif1(lag1(x));x_1st_diff_2nd_lag=dif1(lag2(x));x_1st_diff_3rd_lag=dif1(lag3(x));x_1st_diff_4th_lag=dif1(lag4(x));x_1st_diff_5th_lag=dif1