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

[bac3917]

#Code below throws error: subscript out of bound

library(RCurl) H112.raw<-getURL("https://raw.githubusercontent.com/bac3917/Cauldron/master/Huitema11_1.csv") H112<-read.csv(text=H112.raw) head(H112)

J-N approach helps to define regions of significance and non-significance

use probemod package

library(probemod) adj2<-lm(y~x+D,data=H112) ls(adj2) JMresult1 <- jn(adj2, dv='fitted.values', iv='xlevels', mod='qr',alpha = .01)