作者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