R语言如何在生存分析与Cox回归中计算IDI,NRI指标
原文链接:http://tecdat.cn/?p=6095 原文出处:拓端数据部落公众号 读取样本数据
D=D\[!is.na(apply(D,1,mean)),\] ;
dim(D)
## \[1\] 4167
查询部分数据(结果和预测因子)
head(D)
##time statusage albumin edema protime bili
## 14001 58.765232.601.012.2 14.5
## 2 45000 56.446274.140.010.61.1
## 3 10121 70.072553.480.512.01.4
## 4 19251 54.740592.540.510.31.8
## 5 15040 38.105413.530.010.93.4
## 6 25031 66.258733.980.011.00.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
## 14001
## 2 45000
## 3 10121
## 4 19251
## 5 15040
## 6 25031
``````
head(covs0)
##age albumin edema protime
## 1 58.765232.601.012.2
## 2 56.446274.140.010.6
## 3 70.072553.480.512.0
## 4 54.740592.540.510.3
## 5 38.105413.530.010.9
## 6 66.258733.980.011.0
``````
head(covs1)
##age albumin edema protime bili
## 1 58.765232.601.012.2 14.5
## 2 56.446274.140.010.61.1
## 3 70.072553.480.512.01.4
## 4 54.740592.540.510.31.8
## 5 38.105413.530.010.93.4
## 6 66.258733.980.011.00.8
推理
t0=365*5
x<-IDI (outcome, covs0, covs1, t0, npert=200) ;
输出
##Est. Lower Upper p-value
## M1 0.090 0.052 0.1190
## M2 0.457 0.340 0.5660
## M3 0.041 0.025 0.0620
M1表示IDI
M2表示NRI
M3表示中位数差异
图形演示
文章图片
文章图片
参考文献
1.R语言绘制生存曲线估计|生存分析|如何R作生存曲线图
【R语言如何在生存分析与Cox回归中计算IDI,NRI指标】2.R语言生存分析可视化分析
3.R语言如何在生存分析与Cox回归中计算IDI,NRI指标
4.r语言中使用Bioconductor 分析芯片数据
5.R语言生存分析数据分析可视化案例
6.r语言ggplot2误差棒图快速指南
7.R 语言绘制功能富集泡泡图
8.R语言如何找到患者数据中具有差异的指标?(PLS—DA分析)
9.R语言中的生存分析Survival analysis晚期肺癌患者4例
推荐阅读
- 任时光绽放成六月繁花
- 我从来不做坏事
- 考研英语阅读终极解决方案——阅读理解如何巧拿高分
- 樱花雨
- 如何寻找情感问答App的分析切入点
- 拍照一年啦,如果你想了解我,那就请先看看这篇文章
- mybatisplus如何在xml的连表查询中使用queryWrapper
- MybatisPlus使用queryWrapper如何实现复杂查询
- 【生信技能树】R语言练习题|【生信技能树】R语言练习题 - 中级
- 人如果没梦想,和咸鱼有什么区别(自媒体时代把握住就能咸鱼翻身)