[问答] 【求助】各位牛人帮忙看一下我写的状态空间模型哪里出了问题？ [推广有奖]

0关注
0粉丝

小学生

57%

还不是VIP/贵宾

威望: 0 级
论坛币: 31 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 125 点
帖子: 6
精华: 0
在线时间: 0 小时
注册时间: 2008-3-4
最后登录: 2014-5-6

楼主

brusselee 发表于 2008-5-8 01:45:00 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我用卡尔曼滤波估计三因子广义高斯仿射模型（参考文献：范龙振，上交所利率期限结构的三因子广义高斯仿射模型），在eviews中定义状态空间模型如下：

n1=C(1)-(C(51)^2/(C(41)^2*1))*(1-2*(1-EXP(-C(41)*1))/C(41)+(1-EXP(-2*C(41)*1))/(2*C(41)))-(C(52)^2/(C(42)^2*1))*(1- 2*(1-EXP(-C(42)*1))/C(42)+(1-EXP(-2*C(42)*1))/(2*C(42)))-(C(53)^2/(C(43)^2*1))*(1-2*(1-EXP(-C(43)*1))/C(43)+(1- EXP(-2*C(43)*1))/(2*C(43)))+((1-EXP(-C(41)*1))/(C(41)*1))*SV1+((1-EXP(-C(42)*1))/(C(42)*1))*SV2+((1-EXP(- C (43)*1))/(C(43)*1))*SV3+[VAR=C(61)]

n2=C(1)-(C(51)^2/(C(41)^2*2))*(2-2*(1-EXP(-C(41)*2))/C(41)+(1-EXP(-2*C(41)*2))/(2*C(41)))-(C(52)^2/(C(42)^2*2))*(2- 2*(1-EXP(-C(42)*2))/C(42)+(1-EXP(-2*C(42)*2))/(2*C(42)))-(C(53)^2/(C(43)^2*2))*(2-2*(1-EXP(-C(43)*2))/C(43)+(1- EXP(-2*C(43)*2))/(2*C(43)))+((1-EXP(-C(41)*2))/(C(41)*2))*SV1+((1-EXP(-C(42)*2))/(C(42)*2))*SV2+((1-EXP(- C(43)*2))/(C(43)*2))*SV3+[VAR=C(62)]

n3=C(1)-(C(51)^2/(C(41)^2*3))*(3-2*(1-EXP(-C(41)*3))/C(41)+(1-EXP(-2*C(41)*3))/(2*C(41)))-(C(52)^2/(C(42)^2*3))*(3- 2*(1-EXP(-C(42)*3))/C(42)+(1-EXP(-2*C(42)*3))/(2*C(42)))-(C(53)^2/(C(43)^2*3))*(3-2*(1-EXP(-C(43)*3))/C(43)+(1- EXP(-2*C(43)*3))/(2*C(43)))+((1-EXP(-C(41)*3))/(C(41)*3))*SV1+((1-EXP(-C(42)*3))/(C(42)*3))*SV2+((1-EXP(- C(43)*3))/(C(43)*3))*SV3+[VAR=C(63)]

n4=C(1)-(C(51)^2/(C(41)^2*4))*(4-2*(1-EXP(-C(41)*4))/C(41)+(1-EXP(-2*C(41)*4))/(2*C(41)))-(C(52)^2/(C(42)^2*4))*(4- 2*(1-EXP(-C(42)*4))/C(42)+(1-EXP(-2*C(42)*4))/(2*C(42)))-(C(53)^2/(C(43)^2*4))*(4-2*(1-EXP(-C(43)*4))/C(43)+(1- EXP(-2*C(43)*4))/(2*C(43)))+((1-EXP(-C(41)*4))/(C(41)*4))*SV1+((1-EXP(-C(42)*4))/(C(42)*4))*SV2+((1-EXP(- C(43)*4))/(C(43)*4))*SV3+[VAR=C(64)]

n5=C(1)-(C(51)^2/(C(41)^2*5))*(5-2*(1-EXP(-C(41)*5))/C(41)+(1-EXP(-2*C(41)*5))/(2*C(41)))-(C(52)^2/(C(42)^2*5))*(5- 2*(1-EXP(-C(42)*5))/C(42)+(1-EXP(-2*C(42)*5))/(2*C(42)))-(C(53)^2/(C(43)^2*5))*(5-2*(1-EXP(-C(43)*5))/C(43)+(1- EXP(-2*C(43)*5))/(2*C(43)))+((1-EXP(-C(41)*5))/(C(41)*5))*SV1+((1-EXP(-C(42)*5))/(C(42)*5))*SV2+((1-EXP(- C(43)*5))/(C(43)*5))*SV3+[VAR=C(65)]

@STATE SV1=EXP(-(C(41)-C(11)))*SV1(-1)+C(10)/(C(41)-C(11))-C(10)*EXP(-(C(41)-C(11)))/(C(41)-C(11))+[VAR=C(51)^2*(1-EXP(-2*(C(41)-C(11))))/(2*(C(41)-C(11)))]
@STATE SV2=EXP(-(C(42)-C(22)))*SV2(-1)+C(20)/(C(42)-C(22))-C(20)*EXP(-(C(42)-C(22)))/(C(42)-C(22))+[VAR=C(52)^2*(1-EXP(-2*(C(42)-C(22))))/(2*(C(42)-C(22)))]
@STATE SV3=EXP(-(C(43)-C(33)))*SV3(-1)+C(30)/(C(43)-C(33))-C(30)*EXP(-(C(43)-C(33)))/(C(43)-C(33))+[VAR=C(53)^2*(1-EXP(-2*(C(43)-C(33))))/(2*(C(43)-C(33)))]

其中，n1,n2,n3,n4,n5为可观测变量（分别有849个样本量），C(41)，C(42)，C(43)，C(51)，C(52)，C(53)，C(10)，C(20)，C(30)，C(11)，C(22)，C(33)，C(61)，C(62)，C(63)，C(64)，C(65) 都为待估计参数。

模型我已经检查好几遍，并没有发现问题，但就是在估计时总是会出现“Missing value found in state innovation vector”提示，不知道这是什么意思？到底是什么地方出了问题？还望各位牛人不吝赐教！

多谢！！！