作者tncarol (^^)
看板Statistics
标题[程式] RATS 正交化(orthogonal)问题
时间Tue Apr 20 23:29:28 2010
----------------------------------------------------------------------------
[软体程式类别]:
RATS
[程式问题]:
VEC模型中两外生变数之正交化处理
[软体熟悉度]:
低(1~3个月)
[问题叙述]:
在做VEC模型时,必须将模型中之两个外生变数(rstoxx,rsp)做正交化处理,官网上提供以
下程式,
declare rect f(n,n)
ewise f(i,j)=sigma(i,j)/sqrt(sigma(j,j))
impulse(decomp=f)
但跑出下面的错误:
## MAT12. Syntax: EWISE array(I,J)=expression or array(I)=expression
想请教各位哪里写错了?
非常非常非常感谢回答!!!
[程式范例]:
calendar(irregular)
allocate 2919
open data "d:\data.txt"
data(format=free,org=cols) / lnftse lndax lnftsemib lnibex lncac rftse rdax
rftsemib ribex rcac rstoxx rsp
dec vect[series] y(5) u(5)
dec vect[frm1] resid(5)
declare rect f(2,2)
ewise f(rstoxx,rsp)=sigma(rstoxx,rsp)/sqrt(sigma(rsp,rsp))
impulse(decomp=f)
set Z = lnftse - 0.698855 * lndax -0.543291 * lnftsemib + 0.130312 * lnibex +
0.502418 * lncac - 2.368344
nonlin(parmset=meanparms) c1 a1 b111 b112 b113 b114 b115 b116 b117 $
b121 b122 b123 b124 b125 b126 b127 $
b131 b132 b133 b134 b135 b136 b137 $
b141 b142 b143 b144 b145 b146 b147 $
b151 b152 b153 b154 b155 b156 b157 $
b161 $
b171 $
c2 a2 b211 b212 b213 b214 b215 b216 b217 $
b221 b222 b223 b224 b225 b226 b227 $
.
.
.
----------------------------------------------------------------------------
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 58.114.240.210
1F：推 virginfish:sigma有declare吗？ 04/21 12:32
2F：→ tncarol:谢谢提醒! 05/10 21:26