Question

我正尝试使用R包rms尝试间隔验证。在做这个comnand：

fit.glm<-glm(y.0 ~ x1 + rcs(x2, 3) + rcs(x4, 3) + log(x5) + x7 + x9 + rcs(x2,3):log(x5) + rcs(x2, 3):x7 + rcs(x4, 3)*log(x5), x = TRUE, y = TRUE)
summary(fit.glm)
summary(fitted(fit.glm))

后面的命令给了我

> summary(fitted(fit.glm))
    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
-0.27450  0.02072  0.04557  0.07351  0.08772  1.00600

为什么它是负值...，它可能是R错误，还是我遗漏了某些东西......？

Answer 1

如果在R中执行逻辑回归，fitted.values的范围应为0到1.但是，在您提供的示例中，您只执行了普通的线性回归。要执行逻辑回归，您需要在glm函数中指定错误分布，在您的情况下，family=binomial。例如：

glm.fit <- glm(y~var1+var2,data=data,family=binomial)

可以拟合。来自R的值给出负值吗？

1 个答案: