spss时间序列分析步骤?spss如何进行时间序列预测

牵着乌龟去散步 生活 19 0

大家好,今天小编来为大家解答spss时间序列分析步骤这个问题,spss如何进行时间序列预测很多人还不知道,现在让我们一起来看看吧!

本文目录

  1. 如何使用SPSS做时间序列分析
  2. 如何用SPSS进行时间序列分析
  3. SPSS的时间序列分析怎么做

一、如何使用SPSS做时间序列分析

1、我们在使用SPSS做数据分析的时候,有时需要利用SPSS做时间序列分析,那么时间序列分析应该注意什么和具体该如何去操作呢?

2、首先,我们在SPSS里面导入Excel里面的一组测试数据用来做时间序列分析。在如图所示的对话框中“打开现有数据源”下面选择图示的excel文件。

3、然后在弹出的“打开Excel数据源”框内,“工作表”下面选择你输入数据的Excel sheet表格,单击“确定”。

4、接着,我们需要查看我们导入的数据,比如是否有缺失数据,数据的分布是怎么样的。 *** 一:点击左下角“数据视图”,查看原数据(使用数据不多的情况); *** 二:依次点击“分析-描述统计-描述“查看数据情况(数据多的情况下推荐)。

5、在完成上面的步骤后,做时间序列分析前需要对数据进行一个预处理,即为数据定义日期。

6、首先,我们在如图所示的菜单上依次点击“数据--定义日期”。

7、接着,我们在弹出的“定义日期”对话框内,设置日期的格式。在图示的案例中,我们现在“年份,月”作为日期格式。

8、确定日期格式后,我们在SPSS数据表格里面的“数据视图”可以看到新插入的日期“Year”“Month”“Date”(新变量默认名称)。

9、首先,我们利用指数平滑法时间序列分析。指数平滑法的使用特点是将较大的权数放在最近的资料。

10、我们依次点击之一排菜单栏里面的“分析-预测-创建模型”,弹出“时间序列建模器”。

11、现在进行一些设置:在“变量”选项下,将需要进行时间序列预测的变量拖入图示的“因变量”框内;在 *** 中,选择“指数平滑法”。

12、其他的一些设置包括【统计量】,我们勾选“平稳的R方”“拟合优度”“显示预测值”;在【图表】中选“观察值”“预测值”“拟合值”;在【保存】中勾选“预测值”;在【选项】下填写我们需要预测到的指定日期。

13、全部设置完成后,点击【确定】,即可在输出文档里面看到时间序列建模程序显示的结果(右击结果图表可以复制,导出到Excel等操作)。

14、我们还可以使用ARIMA模型来进行时间序列的预测,使用的特点是将非平稳时间序列转化为平稳时间序列,然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA模型中ARIMA(p,d,q)称为差分自回归移动平均模型,AR是自回归, p为自回归项; MA为移动平均,q为移动平均项数,d为时间序列成为平稳时所做的差分次数。

15、使用的 *** 和上面的类似,依次点击之一排菜单栏里面的“分析-预测-创建模型”,弹出“时间序列建模器”。

16、由于在指数平滑法中我们做了设置,这里就不需要再次设置。这里需要的设置是把【变量】下面的自变量拖入【因变量】框内和【自变量】框内,然后在 *** 中选择“ARIMA”即可。

17、全部设置完成后,单击确定,即可在刚才的输出文档里面看到使用ARIMA模型的预测结果(同样,这些结果右击可以进行复制导出等操作)。

二、如何用SPSS进行时间序列分析

1、建立工作文件,创建并编辑数据。结果如下图所示。

2、在命令行输入lsycx,然后回车。

3、弹出equation窗口,如图所示。观察t统计量、可决系数等,可知模型通过经济意义检验,查表与X的t统计量比较发现,t检验值显著。模型对Y的解释程度高达99.3%。

4、将样本期范围从1978-2003年扩展为1978-2004年:在workfile窗口中依次点击proc->Structure。

5、弹出WorkfileStructure窗口,将2003改为2004,然后点击ok,如图所示。

6、在Group窗口中输入2004年X的值,如图所示。

7、在equation窗口中点击Forecast。

9、在workfile窗口中会生成一个yf,双击打开它,如图所示,即可看到对2004年的预测值。

三、SPSS的时间序列分析怎么做

(1)一般概念:系统中某一变量的观测值按时间顺序(时间间隔相同)排列成一个数值序列,展示研究对象在一定时期内的变动过程,从中寻找和分析事物的变化特征、发展趋势和规律。它是系统中某一变量受其它各种因素影响的总结果。

(2)研究实质:通过处理预测目标本身的时间序列数据,获得事物随时间过程的演变特性与规律,进而预测事物的未来发展。它不研究事物之间相互依存的因果关系。

(3)假设基础:惯性原则。即在一定条件下,被预测事物的过去变化趋势会延续到未来。暗示着历史数据存在着某些信息,利用它们可以解释与预测时间序列的现在和未来。

近大远小原理(时间越近的数据影响力越大)和无季节性、无趋势性、线性、常数方差等。

(4)研究意义:许多经济、金融、商业等方面的数据都是时间序列数据。

时间序列的预测和评估技术相对完善,其预测情景相对明确。

尤其关注预测目标可用数据的数量和质量,即时间序列的长度和预测的频率。

(1)趋势性:某个变量随着时间进展或自变量变化,呈现一种比较缓慢而长期的持续上升、下降、停留的同性质变动趋向,但变动幅度可能不等。

(2)周期性:某因素由于外部影响随着自然季节的交替出现高峰与低谷的规律。

(3)随机性:个别为随机变动,整体呈统计规律。

(4)综合性:实际变化情况一般是几种变动的叠加或组合。预测时一般设法过滤除去不规则变动,突出反映趋势性和周期性变动。

认识时间序列所具有的变动特征,以便在系统预测时选择采用不同的 *** 。

(1)随机性:均匀分布、无规则分布,可能符合某统计分布。(用因变量的散点图和直方图及其包含的正态分布检验随机性,大多数服从正态分布。)

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

样本序列的自相关函数只是时间间隔的函数,与时间起点无关。其具有对称性,能反映平稳序列的周期性变化。

特征识别利用自相关函数ACF:ρk=γk/γ0

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

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

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

(1)点预测:确定唯一的更好预测数值,其给出了时间序列未来发展趋势的一个简单、直接的结果。但常产生一个非零的预测误差,其不确定程度为点预测值的置信区间。

(2)区间预测:未来预测值的一个区间,即期望序列的实际值以某一概率落入该区间范围内。区间的长度传递了预测不确定性的程度,区间的中点为点预测值。

(3)密度预测:序列未来预测值的一个完整的概率分布。根据密度预测,可建立任意置信水平的区间预测,但需要额外的假设和涉及复杂的计算 *** 。

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

(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、独立的白噪声序列,通过估计指定的模型获得。

当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)。

一阶:|φ1|<1。二阶:φ1+φ2<1、φ1-φ2<1、|φ2|<1。φ越大,自回归过程的波动影响越持久。

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

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

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

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

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

当k>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)。

一阶:|θ1|<1。二阶:|θ2|<1、θ1+θ2<1。

当满足可逆条件时,MA(q)模型可以转换为AR(p)模型

3.自回归移动平均ARMA(p,q)模型

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

式中符号: p和q是模型的自回归阶数和移动平均阶数;

φ和θ是不为零的待定系数;εt独立的误差项;

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

使用两个多项式的比率近似一个较长的AR多项式,即其中p+q个数比AR(p)模型中阶数p小。前二种模型分别是该种模型的特例。

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

平稳时间序列的偏相关系数φk和自相关系数rk均不截尾,但较快收敛到0,则该时间序列可能是ARMA(p,q)模型。实际问题中,多数要用此模型。因此建模解模的主要工作是求解p、q和φ、θ的值,检验εt和yt的值。

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的倍数。

4.自回归综合移动平均ARIMA(p,d,q)模型

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

模型形式类似ARMA(p,q)模型,但数据必须经过特殊处理。特别当线性时间序列非平稳时,不能直接利用ARMA(p,q)模型,但可以利用有限阶差分使非平稳时间序列平稳化,实际应用中d一般不超过2。

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

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

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

(1)作直方图:检验正态性、零均值。

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

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

从图中看出:标准差不为1、均值近似为0,可能需要进行数据变换。

(2)作相关图:检验平稳性、周期性。

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

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

因为一般要求时间序列样本数据n>50,滞后周期k<n/4,所以此处控制更大滞后数值Maximum Number of Lags设定为12。点击继续Continue返回自相关主对话框后,点击OK运行系统,输出自相关图如图3.19所示。

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

若时间序列的正态性或平稳性不够好,则需进行数据变换。常用有差分变换(利用transform—Create Time Series)和对数变换(利用Transform—Compute)进行。一般需反复变换、比较,直到数据序列的正态性、平稳性等达到相对更佳。

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

运行自相关图后,出现自相关图(图3.19)和偏自相关图(图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值),取其值最小的确定为模型。

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

按数据转换transform—建立时间序列Create Time Series的顺序展开对话框,图3.21。

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

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

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

领先移动平均Prior moving average:计算先前的时间序列数值的平均值。

中心移动平均Centered moving average:计算围绕和包括当前值的时间序列数值的平均值。

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

累积和Cumulative sum:计算直到包括当前值的时间序列数值的累计总数。

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

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

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

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

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

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

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

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

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

在样本Cases Are栏中选择定义日期变量的时间间隔,在起始日期First Case Is栏中设定日期变量之一个观测量的值,单击OK完成定义。

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

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

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

根据模型识别结果和建立的新时间变量,选择一个或多个变量进入自变量框;暂时不进行因变量的数据转换;

与自变量的选择对应,根据模型识别结果或实验的思路设定p、(d)、q的值;选择模型中包含常数项;

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

在建立变量Create Variable栏选择新建变量结果暂存原数据文件Add to file项,也可选择用新建变量代替原数据文件中计算结果Replace existing项;

在设定置信区间百分比%Confidence Intervals下拉框选择95;

在预测样本Predict Cases栏选择根据时期给出预测结果的 *** 。

在收敛标准Convergence Criteria栏选择迭代次数Maximum iterations、参数变化精度Parameter change、平方和变化精度Sum of squares change,当运算达到其中一个参数的设定,则迭代终止;

在估计初始值Initial Values for Estimation栏选择由过程自动选择Automatic或由先前模型提供Apply from previous model,一般默认前者;

在预测 *** Forecasting Method栏选择无条件Unconditional或有条件最小二乘法Conditional least squares;

在输出控制Display栏选择最初和最终参数的迭代摘要Initial and final parameters with iteration summary或详细资料details、或只显示最终参数Final parameters only。

单击OK,系统立即执行,输出信息如下:

Split group number: 1 Series length: 48

Melard's algorithm will be used for estimation.

Conclusion of estimation phase.

Estimation terminated at iteration number 7 because:Sum of squares decreased by less than.001 percent.

DF Adj. Sum of Squares Residual Variance

Residuals 45 65.099923 1.4392678

AR1.02318739.31945836.0725835.94245925

MA1-.44871554.28829314-1.5564558.12660552

CONSTANT-.02421308.25505018-.0949346.92478827

The following new variables are being created:

FIT_1 Fit for样本数据 from ARIMA, MOD_1 CON

ERR_1 Error for样本数据 from ARIMA, MOD_1 CON

LCL_1 95% LCL for样本数据 from ARIMA, MOD_1 CON

UCL_1 95% UCL for样本数据 from ARIMA, MOD_1 CON

SEP_1 SE of fit for样本数据 from ARIMA, MOD_1 CON

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

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

标准误:表明样本数据的可靠性。在(残差)参数近似服从正态分布条件下,系数加减两倍的标准误差近似等于总体参数95%的置信区间。其值越小,置信区间越窄;并且其对于系数的相对值越小,估计结果越精确。

t统计量:估计系数与标准误差的比值,检验变量的不相关性。一般给定5%显著水平,则拒绝原假设的0值位于95%的置信区间外,其绝对值必大于2。

t概率值:其值越小,则拒绝原假设不相关性的证据越充分。其值接近0.05与t统计量接近2相对应。

均值:度量变量的集中度,传递随机变量的位置信息。

标准差:度量变量的离散度,传递随机变量的规模信息。

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

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

R2校正:是模型中自变量对因变量变动的解释比例,度量方程预测因变量的成功程度,其是回归标准误差与因变量标准差比较的结果。另一个比较 *** 是回归标准误差不超过因变量均值的10%则为好的模型。

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

LN似然:用于模型比较和假设检验,越大越好。

spss时间序列分析步骤?spss如何进行时间序列预测-第1张图片-

检验新建模型的合理性。若检验不通过,则调整(p,q)值,重新估计参数和检验,反复进行直到接受为止。但模型识别、参数估计、检验修正三个过程之间相互作用、相互影响,有时需要交叉进行、反复实验,才能最终确定模型形式。

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

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

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

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

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

(4)概率图检验残差自相关:以显著性水平0.05计算χ2()概率值,。

(5)均方差检验预测的效果:以预测误差的均方差最小为标准,注意预测误差仅与预测周期有关,而与起始时刻无关。

预测系统研究对象的未来某时刻状态。列出预测模型,计算预测值。

好了,文章到这里就结束啦,如果本次分享的spss时间序列分析步骤和spss如何进行时间序列预测问题对您有所帮助,还望关注下本站哦!

标签: 时间序列 spss 步骤 预测 进行

抱歉,评论功能暂时关闭!