Question

我正在研究R项目。数据存储在3个动物园对象中，每个对象有9列。我想对每个列进行格兰杰测试并将结果（Pr（＆gt; F）值）存储在3个矩阵9x9中。

我尝试使用for（）循环并遇到以下问题

for(i in 1:ncol(dfvola.pre)) {
   grangertest(order=4, 
       dfvola.pre[,2], 
       dfvola.pre[,i])$'Pr(>F)'[2]
 }

产生错误Error in solve(vc[ovar, ovar]) : subscript out of bounds

同时

grangertest(order=4, dfvola.pre[,2], 
            dfvola.pre[,3])$'Pr(>F)'[2]

工作正常，即我可以将Pr（> F）存储在变量或矩阵单元格中。但要制作3个矩阵9x9，我必须编写81 * 3行代码，我想知道R中是否有更优雅的解决方案。我是R的新手并且有一种感觉，也许应该有一个，甚至更优雅的嵌套循环：）

以防万一，这是我的数据（动物园对象）：

head(dfvola.pre)
                    CH           HK          SI          RU         EN
2002-02-08 -0.04105356 -0.021797675 -0.03411227 -0.06821175 -0.1728407
2002-02-11 -0.07383230 -0.034095749 -0.05364423 -0.12425684 -0.1656535
2002-02-12 -0.09345455  0.011086639 -0.06846780 -0.14499508 -0.0169436
2002-02-13 -0.09492381  0.009393133 -0.08113240 -0.13592678  2.4487733
2002-02-14 -0.08875820  0.026173698 -0.08273877 -0.07034832  0.4948738
2002-02-15 -0.07895699 -0.041345747 -0.07501171 -0.01333353  0.2248262
                     US          GR          FR          JP
2002-02-08  0.015261541 -0.02886118 -0.01236761 -0.03934752
2002-02-11 -0.102323758 -0.11697662 -0.12982389 -0.08438561
2002-02-12  0.001706797 -0.11818056 -0.14477277 -0.09312417
2002-02-13 -0.085467795 -0.10278221 -0.07810301 -0.11467290
2002-02-14 -0.082738612 -0.13175354 -0.12782784 -0.12296171
2002-02-15 -0.098760674 -0.08439060 -0.08903754 -0.11871091

Answer 1

i变为等于2时出现问题。然后，您正在测试两个相同数据集之间的差异，从而导致错误：

x = runif(100)
grangertest(order=4, x, x)
Error in solve(vc[ovar, ovar]) : subscript out of bounds

我不熟悉grangertest，但可能是协方差矩阵变得单一。解决方案是跳过i = 2：

for(i in 1:ncol(dfvola.pre)) {
 if(!i == 2) {
   grangertest(order=4, 
       dfvola.pre[,2], 
       dfvola.pre[,i])$'Pr(>F)'[2]
  }
 }

循环中的格兰杰测试

1 个答案: