我有以下概率模型，其中i和a之间具有交互作用：

probit5 <- glm(y ~ a + b + c + d + e + f + g + h + i + j + k + l + i*a,family = binomial(link = "probit"),data = data)

但是当我使用intEff函数时：

intEff(probit5,c("i","a"),data = data)

我收到以下错误：

Error in dimnames(x) <- dn : 
  length of 'dimnames' [2] not equal to array extent
In addition: Warning message:
In cbind(deriv1,deriv2,deriv3,nn,deriv0) :
  number of rows of result is not a multiple of vector length (arg 4)

数据框“ data”看起来像这样，而类说这是一个数据框。

可复制的示例：

a <- c(0,0)
b <- c(18,18,18)
c <- c(3.667,3.667,3.667)
d <- c(4.8,4.8,4.8)
e <- c(1,1,1)
f <- c(4.54,4.54,4.54)
g <- c(3.6364,3.6364,3.6364)
h <- c(2.8,2.8,2.8)
i <- c(1,0)
j <- c(1,1)
k <- c(0,0)
l <- c(0,0)
y <- c(0,0)

data <- data.frame(y,a,b,c,d,e,f,g,h,i,j,k,l)

probit <- glm(y ~ a + b + c + d + e + f + g  + h + i + j + k + l + i*a,data = data)

intEff(probit,data = data)

但这显示了以下错误：

Error in apply(X[,vars],2,table) : dim(X) must have a positive length

也许重要的是要注意data.csv是与read.csv2一起导入的

实际示例的系数：