0

我想计算一个协方差矩阵,其中矩阵分量是通过指数移动平均值(即S(t)=a*Y(t)+(1-a)*S(t-1))获得的。

在我看来mewma应该这样做,但是当我为 60,400 矩阵调用函数时,我得到一个“订阅越界错误”。

我相信它来自这条线:

ucl <- h4[m1[l], m2[a[2] - 1]]

可以mewma为 dim 大于 9 的矩阵调用该函数吗?

否则有替代包吗?

我已关注: https ://stats.stackexchange.com/questions/19328/how-to-compute-an-exponentially-weighted-covariance-matrix-function-in-r

4

0 回答 0