把你的code贴出来。

这个code比较长，我又研究了一下，大概问题出在下面这句话上：

rndns(nobs,rows(sigma_res),seed1)*chol(sigma_res);

这里的chol函数我如果单独运行chol(sigma_res)，它会说矩阵是非正定的。可是总体运行的时候它就是正常运行，不显示这个问题。因为程序是我借用的，所以我不知道这一块是不是有问题，是不是应该改？

谁来帮帮忙？

单看这一句程序，确实看不出问题........

From: Pasquale Della Corte

Hi,

may be you have problems with the inverse matrix.

In place of inv( ), use invswp( ) that allows you to get the inverse
of a singular matrix.

Best

Pasquale