У меня есть следующая пробитная модель с эффектом взаимодействия между 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" выглядит следующим образом, и класс говорит, что это data.frame
Воспроизводимый пример:
a <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
b <- c(18, 18, 18, 18, 18, 18, 18, 18, 18, 18)
c <- c(3.667, 3.667, 3.667, 3.667, 3.667, 3.667, 3.667, 3.667, 3.667, 3.667)
d <- c(4.8, 4.8, 4.8, 4.8, 4.8, 4.8, 4.8, 4.8, 4.8, 4.8)
e <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
f <- c(4.54, 4.54, 4.54, 4.54, 4.54, 4.54, 4.54, 4.54, 4.54, 4.54)
g <- c(3.6364, 3.6364, 3.6364, 3.6364, 3.6364, 3.6364, 3.6364, 3.6364, 3.6364, 3.6364)
h <- c(2.8, 2.8, 2.8, 2.8, 2.8, 2.8, 2.8, 2.8, 2.8, 2.8)
i <- c(1, 0, 1, 0, 1, 0, 1, 0, 1, 0)
j <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
k <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
l <- c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0)
y <- c(0, 0, 1, 0, 0, 0, 0, 0, 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,
family = binomial(link = "probit"),
data = data)
intEff(probit, c("i", "a"), data = data)
Но это показывает следующую ошибку:
Error in apply(X[, vars], 2, table) : dim(X) must have a positive length
Возможно, важно отметить, что data.csv был импортирован с read.csv2
Коэффициенты реального примера:
