APP下载

p维TVPAR模型中参数的极大似然估计

2013-12-03施三支王泽升

吉林大学学报(理学版) 2013年6期
关键词:模拟计算游动估计值

施三支,闫 丽,王泽升

(1.长春理工大学 理学院,长春 130022;2.吉林大学 数学学院,长春 130012)

时间序列模型,特别是随时间变化参数的向量自回归(TVP-VAR)模型,在经济、金融等领域应用广泛[1-4].对随时间变化参数的自回归模型(TVPAR模型),目前采用的普遍方法是将所有的参数(包括随机误差)都视为随机游动过程.由于随机游动过程是非平稳的,因此,当随机误差被视为随机游动过程在分析金融数据的长期变化时是不可靠的.文献[5-6]用Markov链Monte Carlo(MCMC)方法讨论了TVP-VAR模型中随机误差不变的情形.本文主要考虑p维TVPAR模型,当变系数仍为系数已知的时间序列,而方差为同方差,且与p维时间序列模型的误差方差均为已知时,采用极大似然估计方法,对p维TVPAR模型中的未知参数进行估计,给出了估计的显式表达式,并进行了模拟计算.

p维TVPAR模型具有如下形式:

(1)

这里{xt,t=1,2,…,n}为可观察数据.本文要求模型满足下列条件:

2)βi,ui是系数参数,且-1<βi<1(i=1,2,…,p)未知,-1

3)βt,i(i=1,2,…,p;t=1,2,…,n)为不可观测数据;βt,1,βt,2,…,βt,p不相关;

4)β0,i=0(i=1,2,…,p);x0=0.

1 模型的极大似然估计

由模型(1)知:

记β=(β1,β2,…,βp)T,u=(u1,u2,…,up)T.由假设条件1)知

即有

(2)

根据式(2)知,似然函数为

(3)

用类似的方法可以讨论模型(1)的平稳性:记Xt=(xt,xt-1,…,xt-p+1)T,Xt是p×1维随机向量.记Bt=(βt,1,βt,2,…,βt,p),则Bt为1×p维AR(1)过程,Bt各元素之间相互独立.记U=diag(u1,u2,…,up),ηt=(ηt,1,ηt,2,…,ηt,p),将模型(1)改写为如下形式:

Xt=(M+Dt)Xt-1+εt,Dt=Dt-1U+μt,

(4)

2 模拟计算

表1 λ=0.1时不同p值下参数β的真值与估计值的比较Table 1 True values and estimated values of β under different p when λ=0.1

表2 λ=0.01时不同p值下参数β的真值与估计值的比较Table 2 True values and estimated values of β under diffrent p when λ=0.01

[1] Primiceri G E.Time Varying Structural Vector Autoregressions and Monetary Policy [J].Review of Economic Studies,2005,72(3):821-852.

[2] Benati L,Mumtaz H U S.Evolving Macroeconomic Dynamics:A Structural Investigation [R].Working Paper 746.Frankfurt:European Central Bank,2007.

[3] Baumeister C,Durinck E,Peersman G.Liquidity,Inflation and Asset Prices in a Time-Varying Framework for the Euro Area [J].National Bank of Belgium in Its Series Research Series,2008:10-17.

[4] Nakajima J,Kasuya M,Watanabe T.Bayesian Analysis of Time-Varying Parameter Vector Autoregressive Model for the Japanese Economy and Monetary Policy [J].Journal of the Japanese and International Economics,2011,25(3):225-245.

[5] Canova F,Gambetti L,Pappa E.The Structural Dynamics of Output Growth and Inflation:Some International Evidence [J].The Economic Journal,2007,117(3):167-191.

[6] Gambetti L,Pappa E,Canova F.The Structural Dynamics of US Output and Inflation:What Explains the Changes [J].Journal of Money,Credit,and Banking,2006,40(2/3):369-388.

[7] Weiss A A.The Stability of the AR(1) Process with an AR(1) Coefficient [J].Journal of Time Series Analysis,1985,6(3):181-186.

[8] SHI San-zhi,WANG De-hui,SONG Li-xin,et al.Bayesian Estimation of Time Vary Parameter Autoregression Model [J].Journal of Jilin University:Science Edition,2011,49(5):857-860.(施三支,王德辉,宋立新,等.随时间变化系数参数AR模型的Bayes估计 [J].吉林大学学报:理学版,2011,49(5):857-860.)

(责任编辑:赵立芹)

研究简报

猜你喜欢

模拟计算游动估计值
R1234ze PVTx热物性模拟计算
球轴承用浪型保持架径向游动量的测量
把手放进袋子里
一道样本的数字特征与频率分布直方图的交汇问题
2018年4月世界粗钢产量表(续)万吨
挤出发泡片材褶皱分析及模拟计算
父亲
2014年2月世界粗钢产量表
实际发射工况下底排药柱结构完整性的模拟计算
2014年5月世界粗钢产量表万吨