我在R中创建了一个glm模型,其中包含3个独立变量和一个与二进制相关的变量:
model <- glm(y ~ x1 + x2 + x3, data = train, family = binomial)
在给定自变量特定值的情况下,使用简单的公式来预测结果非常有用:
newdata <- data.frame(x1 = 10, x2 = 8, x3 = 5)
predict(model, newdata, type = "response").
但是考虑到我们想要找出x1呢?
y = 1
x2 = 10
x3 = 2
我听说for循环是获取具有可能值组合的列表的不错选择。但是不知道如何编码。 有人知道吗? 还是解决此问题的任何内置函数?
答案 0 :(得分:0)
认为我已经掌握了
x <- seq(1,100,by=1)
r <- vector(mode="numeric", length=100)
for (x1 in v) {
r[x1] <- 1/(1+exp(-(CoefX1*X1+CoeffX2*10+CoeffX3*2)))
}
cbind(x,r)