第七章自回归模型
- 格式:ppt
- 大小:719.00 KB
- 文档页数:41
自回归(AR)模型自回归模型(Autoregressive model)的形式为:1122n n n p n p n X X X X ϕϕϕε---=++++ (5.1.1)式中p ϕϕ,,1 为模型参数;n X 为因变量,12,,...,n n n p X X X ---为“自”变量。
这里“自”变量是同一(因此称为“自”)变量,但属于以前各个时期的数值,所谓自回归即是此含义。
最后,},1,0,{ ±=n n ε是白噪声序列,即ko k n n n E E δσεεεε2)(,0)(==+,也就是说随机序列}{n ε的均值为零,方差为2εσ,且互不相关,它代表不能用模型说明随机因素。
假定()0,()t n E X t n ε=<,即随机影响与数据值无关。
p 为模型的阶数。
用)(p AR 来简记此模型。
引入向后推移算子B : k n n k B X X -=,C C B k =, ,1,0=k (C 为常数)并记)1()(221p p p B B B B ϕϕϕ----=Φ则式(5.1.1)可重写为()p n n B X εΦ= (5.1.2)称多项式方程()0p λΦ=为)(p AR 模型的特征方程,它的p 个根p λλλ,,,21 称为模型的特征根。
特征根可能是实数,也可能是复数。
如果这p 个特征根都在单位圆外,即1,1,2,...,i i p λ>=则称)(p AR 模型是稳定的或平稳的。
称上式为平稳性条件。
这里应引起读者注意的是,平稳时间序列{}n X 是指n X 的均值为常数(我们设其为零)且自相关函数为齐次的随机时间序列;而平稳的()AR p )则指它满足平稳性条件:()0p λΦ=的根均在单位圆外。
这两种“平稳”是两个不同的概念。
如,对于)1(AR 模型,其特征方程为011=-x ϕ特征根111-=ϕλ,从而)1(AR 的平稳性条件是11<ϕ。
在条件11<ϕ下,有111110N k N n n n n k n N k X X X εϕϕεϕ----==+==+∑ ∑∞=-=01k k n k εϕ (5.1.3)由于k ε表示第k 期的预测误差,因此上式表示对平稳的)1(AR 模型,n X 可由过去各期的误差线性表示。
自回归模型一、 预测方法综述预测方法大体上分为定性预测法、时间序列预测法和因果模型预测法。
定性预测法是在数据资料掌握不多的情况下,依靠人的经验和分析能力,用系统的、逻辑的思维方法,把有关资料加以综合、进行预测的方法。
定性预测法包括特尔斐法、主观概率预测法、判断预测法等方法。
时间序列预测法是依据预测对象过去的统计数据,找到其随时间变化的规律,建立时序模型,以判断未来数值的预测方法。
其基本思想是:过去的变化规律会持续到未来,即未来是过去的延伸。
时间序列预测法包括时间序列平滑法、趋势外推法、季节变动预测法等确定型时间序列的预测方法和马尔可夫法、随机型时间序列的预测方法。
因果模型预测法是把所要预测的对象同其他有关因素联系起来进行分析,制定出揭示因果关系的模型,然后根据模型进行预测。
因果模型预测法包括回归分析预测法、经济计量模型法、投入产出预测法等。
由于时间序列预测法和因果模型预测法都是以统计资料为依据,应用统计方法进行预测的,所以有时两者统称为统计预测。
到目前为止,已有近二百种预测方法。
1987年,Ledes和Farbor首次将神经网络引入到预测领域中,无论是从思想上、还是技术上都是一种拓宽和突破。
常用的分析和预测方法有下面几种:(1) 投资分析方法。
这是市场分析家常用的方法。
(2) 时间序列分析法。
这种方法主要是通过建立综合指数之间的时间序列相关辩识模型,如自回归移动平均模型(ARMA)、齐次非平稳模型(ARIMA)等来预测未来变化。
(3) 神经网络预测法。
神经网络是一种最新的时间序列分析方法。
(4) 其他预测方法。
如专家评估法和市场调查法等定性方法、季节变动法、马尔柯夫法和判别分析法等定量预测方法。
传统的预测方法大都采用线性模型来近似地表达预测对象的发展规律。
如最常用的AR模型预测,就是在时间序列平稳的假设基础之上,对其建立线性模型,然后采用模型外推的方法预测其未来值。
然而这些方法只适用于平稳时间序列的预测。
自回归模型(Autoregressive Model)是一种经典的时间序列预测模型,在许多领域中都有着广泛的应用。
它的核心思想是利用过去时间点的观测值来预测未来的观测值。
在本文中,我将介绍自回归模型的概念,并使用Python实现一个简单的自回归模型。
1.自回归模型概述自回归模型是建立在时间序列数据上的统计模型。
它假设当前时刻的观测值是过去时刻的观测值的线性组合,其中线性关系由模型的参数确定。
自回归模型可以被表示为如下形式:X_t = c + Σ(φ_i *X_(t-i)) + ε_t 其中,X_t是当前时刻的观测值,c是常数项,φ_i是参数,ε_t是误差项。
根据历史观测值和参数的不同,自回归模型可以分为不同阶数的自回归模型,如一阶自回归模型(AR(1))、二阶自回归模型(AR(2))等。
2.自回归模型的Python实现为了实现自回归模型,我们需要借助Python中的统计分析库statsmodels。
我们需要安装statsmodels库,可以使用以下命令进行安装: pip install statsmodels接下来,我们使用一个示例数据集来演示自回归模型的实现。
假设我们有一个包含100个观测值的时间序列数据,可以使用以下代码生成一个随机的时间序列数据:import numpy as np生成随机时间序列数据np.random.seed(0) data = np.random.randn(100)我们可以使用statsmodels库中的AR模型来建立自回归模型,并进行参数估计和预测。
以下是一个简单的自回归模型的实现代码示例: fromstatsmodels.tsa.ar_model import AutoReg构建AR模型model = AutoReg(data, lags=1)拟合模型model_fit = model.fit()打印模型系数print(model_fit.params)进行单步预测predictions = model_fit.predict(start=len(data), end=len(data))print(predictions)在上述代码中,我们首先使用AutoReg类构建了一个自回归模型,其中lags参数指定了模型的阶数,这里我们选择了一阶自回归模型(lags=1)。
练习题7.1参考解答(1)先用第一个模型回归,结果如下:22216.4269 1.008106 t=(-6.619723) (67.0592)R 0.996455 R 0.996233 DW=1.366654 F=4496.936PCE PDI =-+==利用第二个模型进行回归,结果如下:122233.27360.9823820.037158 t=(-5.120436) (6.970817) (0.257997)R 0.996542 R 0.996048 DW=1.570195 F=2017.064t t t PCE PDI PCE -=-++==(2)从模型一得到MPC=1.008106;从模型二得到,短期MPC=0.982382,长期MPC=0.982382+(0.037158)=1.01954练习题7.2参考答案(1)在局部调整假定下,先估计如下形式的一阶自回归模型:*1*1*0*t t t t u Y X Y +++=-ββα估计结果如下:122ˆ15.104030.6292730.271676 se=(4.72945) (0.097819) (0.114858)t= (-3.193613) (6.433031) (2.365315)R =0.987125 R =0.985695 F=690.0561 DW=1.518595t t t Y X Y -=-++根据局部调整模型的参数关系,有****1 1 t tu u αδαβδββδδ===-=将上述估计结果代入得到:*1110.2716760.728324δβ=-=-=*20.738064ααδ==-*0.864001ββδ==故局部调整模型估计结果为:*ˆ20.7380640.864001t tY X =-+经济意义解释:该地区销售额每增加1亿元,未来预期最佳新增固定资产投资为0.864001亿元。
运用德宾h检验一阶自相关:(121(1 1.34022d h =-=-⨯=在显著性水平05.0=α上,查标准正态分布表得临界值,由于,则接收21.96h α=21.3402 1.96h h α=<=原假设0=ρ,说明自回归模型不存在一阶自相关。