The RATS Software Forum

Code: Select all

calendar 2000 1 1
allocate 2008:5
open data Sspace.xls
data(format=xls, org=obs) / *
table
compute kstart=2001:1,kend=2008:5

equation kalmeq returns
# constant vol_g dif_inf dif_tx var_re libor abcp swap stan_poo_vol gdpre_rate ratio_cap_gdp
declare vector b(11) bmean(11)
ewise bmean(i)=(i==2)
compute bmean=b
associate(perm) kalmeq b

linreg returns
# constant vol_g dif_inf dif_tx var_re libor abcp swap stan_poo_vol gdpre_rate ratio_cap_gdp
compute rseesq=.9*%seesq

system kalmeq
kfset(constant,noscale,likelihood=likely) xxx
#rseesq
tvarying tvx
end(system)

compute pi5 = .20**2
compute pi2 = .5**2
compute pi3 = 10000.
compute pi7 = 0.00000001
compute pi8 = 1.00

dim xxx(11,11) tvx(11,11)
*
compute likely=%const(0.0)
compute b=bmean
compute xxx=%const(0.0)
compute xxx(1,1)=pi5*pi3*rseesq
do i=1,4
  compute decayfac=i**(2*.5)
  compute xxx(i+1,i+1)=pi5*decayfac
  compute xxx(i+5,i+5)=pi5*pi2*decayfac*rseesq
end do i
compute tvx=pi7*xxx
*
do time=kstart,kend
  if time==kstart
     kalman(start=time)
  else {
     compute b=pi8*b+(1-pi8)*bmean
     compute xxx=(pi8**2)*xxx
     kalman
     }
end do time

Code: Select all

ewise bmean(i)=(i==2)

Code: Select all

dlm(y=returns,c=||1.0,stan_poo_vol||,sv=1.0,var=concentrate,exact,free=2,type=filter) / xstates vstates

Code: Select all

DLM Setup Options 

A=RECTANGULAR or FRML[RECT] [identity matrix]
This gives the A matrices. The expression should evaluate to an NxN matrix.

C=RECTANGULAR or FRML[RECT] [zero matrix]
This gives the C matrices—it should evaluate to an NxM matrix.

Y=VECTOR or FRML[VECT]  [no measurement equation]
This gives the Y vectors—it should evaluate to an M vector

Y=VECTOR or FRML[VECT]  [no measurement equation]
This gives the Y vectors—it should evaluate to an M vector

SV=SYMMETRIC or FRML[SYMMETRIC] [zero]
This gives the covariance matrices of the v’s—it should evaluate to an MxM SYMMETRIC.

SW=SYMMETRIC or FRML[SYMMETRIC] [zero]
This gives the covariance matrices of the w’s—it should evaluate to an NxN SYMMETRIC array.

SX0=SYMMETRIC or FRML[SYMMETRIC] [zero]
This gives the covariance of X(0)—it should evaluate to an NxN SYMMETRIC array.

X0=VECTOR or FRML[VECTOR] [zero vector]
This gives the initial X vector—it should evaluate to an N vector.

Z=VECTOR or FRML[VECTOR] [zero]
This allows the inclusion of exogenous shifts in the state equation—it should evaluate to an N VECTOR.

B=RECTANGULAR or FRML[RECT] [NxN identity matrix]
This gives the loadings from the w's to the X's. It should evaluate to an NxL matrix. 

F=RECTANGULAR or FRML[RECT] [zero matrix]
This gives the B matrices. The expression should evaluate to an NxP matrix.

Q=SYMMETRIC or FRML[SYMMETRIC] [zero]
This gives the Q matrices—it should evaluate to an NxN SYMMETRIC.

R=SYMMETRIC or FRML[SYMMETRIC] [zero]
This gives the R matrices—it should evaluate to a PxP SYMMETRIC.

Code: Select all

## MAT13. Store into Out-of-Range Matrix or Series Element
The Error Occurred At Location 0113 of loop/block
Line 4 of loop/block