R语言如何在生存分析与Cox回归中计算IDI,NRI指标

Posted 拓端数据部落

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了R语言如何在生存分析与Cox回归中计算IDI,NRI指标相关的知识,希望对你有一定的参考价值。

原文链接: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*5x<-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 0

M1表示IDI

M2表示NRI

M3表示中位数差异


图形演示


点击标题查阅往期内容



更多内容,请点击左下角“阅读原文”查看

R语言如何在生存分析与Cox回归中计算IDI,NRI指标

R语言如何在生存分析与Cox回归中计算IDI,NRI指标

R语言如何在生存分析与Cox回归中计算IDI,NRI指标


案例精选、技术干货 第一时间与您分享

长按二维码加关注

更多内容,请点击左下角“阅读原文”查看

以上是关于R语言如何在生存分析与Cox回归中计算IDI,NRI指标的主要内容,如果未能解决你的问题,请参考以下文章

R语言生存分析生成仿真数据构建COX回归分析模型并计算C-index

R语言生存分析之COX比例风险模型构建及C-index计算示例

R语言 | 生存分析之R包survival的单变量和多变量Cox回归

R语言生存分析COX回归分析实战:放疗是否会延长胰脏癌症患者的生存时间

R语言生存分析COX回归分析实战:以乳腺癌数据为例

R语言生存分析COX回归分析实战:两种治疗方法发生肾功能损害的情况