R 中的 fitted.values 可以给出负值吗?
Can fitted.values from R give negative values?
我正在尝试使用 R
包 rms
进行间隔验证。在执行此命令时:
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
错误,还是我遗漏了什么...?
如果您在 R 中执行逻辑回归,fitted.values
的范围应为 0 到 1。但是,在您提供的示例中,您只是执行了普通的线性回归。要执行逻辑回归,您需要在 glm
函数中指定误差分布,在您的例子中为 family=binomial
。例如:
glm.fit <- glm(y~var1+var2,data=data,family=binomial)
我正在尝试使用 R
包 rms
进行间隔验证。在执行此命令时:
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
错误,还是我遗漏了什么...?
如果您在 R 中执行逻辑回归,fitted.values
的范围应为 0 到 1。但是,在您提供的示例中,您只是执行了普通的线性回归。要执行逻辑回归,您需要在 glm
函数中指定误差分布,在您的例子中为 family=binomial
。例如:
glm.fit <- glm(y~var1+var2,data=data,family=binomial)