本征（同伴）错误：“ x”中的值无穷或缺失

Question

VAR估计中的错误可能在vars pachange的根函数中

亲爱的参与者， 我有一个包含16个参数和2515个观测值的数据集。我的问题是相关矩阵。 我已经检查了空单元格，NA和零 （Y [！is.na（Y）]）

运行以下代码时出现错误：

  Y <-csv中的mydata
t <-增长（Y）
k <-ncol（Y）
空间<-200
CV <-array（NA，c（k，k，（t-space）））
colnames（CV）=行名（CV）= colnames（Y）
对于（i in 1：dim（CV）[3]）{
   var1 <-VAR（Y [i：（space + i-1），]，p = nlag，type =“ const”）
   if（any（roots（var1）> 1））{
      CV [,, i] <-CV [,,（i-1）]
   }其他{
     CV [,, i] <-gfevd（var1，n.ahead = nfore）$ fevd
   }
   如果（i %% 500 == 0）{print（i）}
}

 

我得到的错误如下：

  

“本征（同伴）错误：'x'中的值无穷或缺失”。

它应该通过k矩阵为k中的每个i和j元素进行计算。