[R-SIG-Finance] DLM and matrices with 0 eigenvalues

Rebecca Sela rsela at stern.nyu.edu
Thu Feb 19 14:39:17 CET 2009


I am using DLM to fit a state space model.  The covariance matrix of states (W) is given by:
a 0 a 0
0 0 0 0
a 0 a 0
0 0 0 0
where a is a parameter to be estimated.  Even though the matrix is positive semidefinite, sometimes DLM gives me an error that W is not a valid variance matrix.  As far as I can tell, the reason is that one of R's computed eigenvalues is very slightly negative (something like -5E-17).  Is there a way to work around this?

Thanks!

Rebecca



More information about the R-SIG-Finance mailing list