R中快速绘制PR曲线并计算AUPR-modEvA包

前言:

绘制ROC曲线的R包有很多,例如pROC、ROCR等,并可以获得AUC,有的顺便可以绘制PR曲线,但不能获得AUPR。最近需要获取AUPR,但查询网上几乎没有关于AUPR的,这里查找到一个简单的R包-modEvA,写下使用的简单案例。

术语:

PR曲线:precision-recall curve,以recall作为横坐标轴,precision作为纵坐标轴

AUPR:Area under the precision-recall curve,PR曲线下的面积

参考:

Barbosa A.M., Real R., Muñoz A.R. & Brown J.A. (2013) New measures for assessing model equilibrium and prediction mismatch in species distribution models. Diversity and Distributions 19: 1333-1338 (DOI: 10.1111/ddi.12100)

https://cran.r-project.org/web/packages/modEvA/modEvA.pdf</

  • 4
    点赞
  • 21
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
在Python,你可以使用 lifelines 来进行生存数据的Cox分析,建立Cox生存模型,并计算AUC、MCC、AUPR、DCA、KM曲线和校准曲线。 首先,需要安装 lifelines ,可以使用以下命令进行安装: ``` pip install lifelines ``` 接下来,我们将详细解释并提供示例代码来完成你的需求。 1. 导入所需的库和模块: ```python import pandas as pd import numpy as np from lifelines import CoxPHFitter from lifelines.utils import concordance_index from lifelines.utils import k_fold_cross_validation from sklearn.metrics import matthews_corrcoef, roc_auc_score, average_precision_score, precision_recall_curve import matplotlib.pyplot as plt ``` 2. 加载数据并进行预处理: ```python # 假设你的生存数据保存在名为 'data.csv' 的CSV文件,其含特征和生存时间 data = pd.read_csv('data.csv') # 对数据进行预处理,例如删除缺失值等 data = data.dropna() ``` 3. 创建Cox生存模型并拟合数据: ```python # 创建CoxPHFitter对象 cph = CoxPHFitter() # 拟合Cox生存模型 cph.fit(data, 'time', event_col='event') ``` 4. 计算Cox模型的性能指标: ```python # 预测生存概率 data['predicted_proba'] = cph.predict_survival_function(data).T.iloc[:, 0] # 计算AUC auc = roc_auc_score(data['event'], data['predicted_proba']) # 计算MCC mcc = matthews_corrcoef(data['event'], data['predicted_proba'] > 0.5) # 计算AUPR aupr = average_precision_score(data['event'], data['predicted_proba']) # 计算DCA dca = k_fold_cross_validation(cph, data, 'time', event_col='event', k=5).mean() # 绘制KM曲线 cph.plot() # 绘制校准曲线 cph.plot_calibration_curve(data, 'time', event_col='event') ``` 以上代码提供了一个简单的示例,你可以根据你的数据和需求进行适当的修改。请确保在使用之前,先根据你的实际情况加载和预处理数据。 希望这个示例能帮助到你!
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值