微信公众号:医学统计与R语言
如果你觉得对你有帮助,欢迎转发
输入1:
"rio")
结果1:
Gender Race Age Income Cigarettes Smoking_status
输入2:
m3 <- lrm(Smoking_status~., x =T,y = T,data = bilog)newdata <- bilog[101:200,]pred.lg<- predict(m3,newdata )newdata$prob <- 1/(1+exp(-pred.lg))val.prob(newdata$prob, newdata$Smoking_status ,m=10, cex=1)
结果2: