r roc曲线 语言_调整ROC曲线的R实现

  前言  

近日有朋友在《统计咨询》公众号咨询：如何绘制调整混杂变量C后，X预测Y的ROC曲线？这个问题马上引起了小编的注意，因为在ROC曲线方面，常见的是联合多个指标对结局进行预测，而调整混杂变量的情况比较少见，但具有重要的应用价值。比如，我们希望在校正性别的情况下，评估年龄对肾移植术后生存与否的诊断效能，就需要计算调整ROC曲线下的面积 (AUC)。

  调整ROC的原理  

感兴趣的朋友可以查看发表在 Statistical Methods inMedical Research 原文： Standardized and weighted time-dependent receiver operating characteristic curves to evaluate the intrinsic prognostic capacities of a marker by taking into account confounding factors 。文章介绍的是采用逆概率权重的原理对标准化后的预测因子进行诊断效能评估的，而预测因子的标准化是根据对照组中混杂变量的分布情况来对全部样本进行标准化的。

逆概率权重在倾向评分加权中的应用更为常见，感兴趣的朋友可以查看往期文章《三组的逆概率权重IPTW》。

  R语言实例  

首先介绍一下本文用到的数据 (dataDIVAT) ，这是一个来自法国的肾移植队列研究，样本量为4267，我们用到的变量为①结局death：随访结束时的生存状态，