Question

我正在尝试探讨审核关系，但我不断收到以下错误：

Error in cov[interactionterm, interactionterm] : subscript out of bounds

有人能告诉我我做错了什么吗？可重复的例子：

library(probemod)
set.seed(42)
DV<-sample(1:10,20,replace=T)
IV<-sample(1:10,20,replace=T)
MOD<-sample(1:10,20,replace=T)
model<-lm(DV~IV+MOD)
jn(model, dv='DV', iv='IV', mod='MOD')

Answer 1

对我而言，它看起来像包中的一个错误。如果查看jn函数，则索引协方差矩阵时会出错：

   interactionterm <- paste(jnret$iv, jnret$mod, sep = ":")
    jna <- tcrit^2 * cov[interactionterm, interactionterm] - 
      beta.hat[interactionterm]^2

interactionterm将在您的案例中分配字符串"IV:MOD"。这显然不是正确的索引，并且在jna帮助中运行示例时会出现相同的错误。即使修复此行，也会出现另一个错误。

我会写信给包的维护者。

Answer 2

用model<-lm(DV~IV+MOD)替换model<-lm(DV~IV*MOD)，可能会解决此问题。

probemod下标超出范围

2 个答案: