大家新年快乐!!!!
我难得元旦有个假期,还来给各位录视频,各位不觉得感动吗。。。。。。
知乎视频www.zhihu.com这期视频主要是我之前帮助的那个同门非要我把上次给他做的回归的C-index弄出来,我跟他解释了说对现在做的这个二元logistic回归模型,我给你画的ROC曲线的曲线下面积AUC值的大小就可以评价了。
然鹅他不听
然鹅他之前请我吃了顿大餐
所以我还是给他把C-index弄出来了
顺带还送了一张calibration curve
以下是分步的代码
第一步:加载包,读取数据,数据修饰
library(openxlsx)
library(survival)
library(lattice)
library(Formula)
library(ggplot2)
library(Hmisc)
library(rms)
lxrdata<-read.xlsx("cknormo.xlsx", sheet =1)
dd=datadist(lxrdata)
options(datadist="dd")
第二步:做logistic回归,绘制nomogram图
ckf <- lrm(type ~ age + gender +area+ sym ,data = lxrdata,x=T,y=T)
nom <-