Ihavenothing
您好,我尝试按照上面的代码画了图,因为发现纵坐标pred都是负值(范围大约是-4~0),觉得有些奇怪,不太确定我是否正确理解了您分享的方法,非常希望能和您确认一下方法。
在已有x(连续)和y变量(已编码为0,1)的情况下是否直接打下面这部分代码就可以了?
library(splines)
model = glm(y ~ bs(x), family = binomial())
pred = predict(model, bs(x))
plot(x, pred)
另外,平滑分析的方法是不是无法处理协变量,比如全模型是
fit<-lm(y~age+gender+edu+marr1+marr2+SES+x,data=d)的情况下,要看x与“y为1的概率”的关系,也是用上面的代码吗?
衷心感谢您的帮助!