probemod индекс вне границ

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. Даже при восстановлении этой строки появляется другая ошибка.

Я бы написал сопроводителю пакета.