R语言如何在生存分析与Cox回归中计算IDI，NRI指标
R语言如何在生存分析与Cox回归中计算IDI，NRI指标
拓端
发布于 2020-07-16 15:10:49
发布于 2020-07-16 15:10:49
代码可运行
运行总次数：0
代码可运行
原文链接：http://tecdat.cn/?p=6095
读取样本数据
D=D[!is.na(apply(D,1,mean)),] ; dim(D)
## [1] 416 7查询部分数据（结果和预测因子）
head(D)
## time status age albumin edema protime bili
## 1 400 1 58.76523 2.60 1.0 12.2 14.5
## 2 4500 0 56.44627 4.14 0.0 10.6 1.1
## 3 1012 1 70.07255 3.48 0.5 12.0 1.4
## 4 1925 1 54.74059 2.54 0.5 10.3 1.8
## 5 1504 0 38.10541 3.53 0.0 10.9 3.4
## 6 2503 1 66.25873 3.98 0.0 11.0 0.8模型0和模型1的结果数据和预测变量集
outcome=D[,c(1,2)]
covs1<-as.matrix(D[,c(-1,-2)])
covs0<-as.matrix(D[,c(-1,-2, -7)])
head(outcome)
## time status
## 1 400 1
## 2 4500 0
## 3 1012 1
## 4 1925 1
## 5 1504 0
## 6 2503 1head(covs0)
## age albumin edema protime
## 1 58.76523 2.60 1.0 12.2
## 2 56.44627 4.14 0.0 10.6
## 3 70.07255 3.48 0.5 12.0
## 4 54.74059 2.54 0.5 10.3
## 5 38.10541 3.53 0.0 10.9
## 6 66.25873 3.98 0.0 11.0head(covs1)
## age albumin edema protime bili
## 1 58.76523 2.60 1.0 12.2 14.5
## 2 56.44627 4.14 0.0 10.6 1.1
## 3 70.07255 3.48 0.5 12.0 1.4
## 4 54.74059 2.54 0.5 10.3 1.8
## 5 38.10541 3.53 0.0 10.9 3.4
## 6 66.25873 3.98 0.0 11.0 0.8推理
t0=365*5
x<-IDI (outcome, covs0, covs1, t0,npert=200) ;输出
## Est. Lower Upper p-value
## M1 0.090 0.052 0.119 0
## M2 0.457 0.340 0.566 0
## M3 0.041 0.025 0.062 0M1表示IDI
M2表示NRI
M3表示中位数差异
图形演示
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2020-07-09，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
暂无评论
推荐阅读
编辑精选文章
换一批
推荐阅读
相关推荐
如何在生存分析与Cox回归中计算IDI，NRI指标|附代码数据
更多 >
💥开发者 MCP广场重磅上线！
精选全网热门MCP server，让你的AI更好用 🚀
腾讯云开发者
