：详解的区分度和校准度。基于风险预测模型的预后研究一直以来都是研究者关注的热点，各种各样的预测模型质量参差不齐，常常让人眼花缭乱，那么如何去评价一个模型的好坏，或者说当你构建出一个疾病风险预测模型后，它到底靠不靠谱，值不值得去推广和使用呢?这是一个我们需要去好好考量的问题。
一个好的疾病风险预测模型，它不只是简单的因变量和自变量的数学组合，它背后的实际临床意义才是我们所要把握的重点，这就要求预测模型不仅要有很好的区分度(discrimination)，同时还要具备良好的校准度(calibration)。
discrimination和calibration是我们在评价预测模型时最常用到的一对指标，但是2015年circ cardiovasc qual outcomes杂志(：4.5)上发表的一项关注心血管疾病预测模型的系统发现，63%的研究报告了预测模型的discrimination信息，但仅36%的研究报告了calibration信息，使得预测模型的质量成为研究泛滥的重灾区。
本文详细介绍一下的区分度和校准度这两个重要指标，尤其是常常被人忽略的calibration。
区分度(discrimination).
介绍calibration之前，我们先简单介绍一下discrimination。顾名思义，一个好的疾病风险预测模型，它能够把未来发病风险高、低不同的人群正确地区分开来，预测模型通过设置一定的风险界值，高于界值判断为发病，低于界值则判断为不发病，从而正确区分个体是否会发生结局事件，这就是预测模型的区分度(discrimination)。
评价预测模型区分能力的指标，最常用的就是大家非常熟悉的下面积(auc)，也叫c统计量(c-statistics)。auc越大，说明预测模型的判别区分能力越好。一般auc<0.6认为区分度较差，0.6-0.75认为模型有一定的区分能力，>0.75认为区分能力较好。
校准度(calibration).
预测模型的校准度(calibration)，是评价一个疾病风险模型预测未来某个个体发生结局事件概率准确性的重要指标，它反映了模型预测风险与实际发生风险的一致程度，所以也可以称作为一致性。校准度好，提示预测模型的准确性高，校准度差，则模型有可能高估或低估疾病的发生风险。
在实际的应用中，通常用hosmer-lemeshow good of fit test(拟合优度检验)来评价预测模型的校准度。hosmer-lemeshow检验的基本思路如下：
1. 首先根据预测模型来计算每个个体未来发生结局事件的预测概率;
2. 根据预测概率从小到大进行排序，并按照十分位等分成10组;
3. 分别计算各组的实际观测数和模型预测数，其中模型预测数，即每个人的预测概率*人数，再求总和，这里人数即为1，最后总和就相当于每个个体预测概率的直接加和;
4. 根据每组实际观测数和模型预测数计算卡方值(自由度=8)，再根据卡方分布得到对应的p值。
若所得的统计量卡方值越小，对应的p值越大，则提示预测模型的校准度越好。若检验结果显示有统计学显著性(p<0.05)，则表明模型预测值和实际观测值之间存在一定的差异，模型校准度差。
区分度和校准度的操作:
一、建立数据库.
某研究人员拟建立一个关于冠心病患者支架介入术后再次发生mace事件(major adverse cardiovascular events，主要心血管不良事件)的风险预测模型，并对该风险模型的预测能力进行评价。数据库格式如下图所示。
其中因变量(结局事件)为event，自变量(影响因素)为性别(gender)、年龄(age)、收缩压(sbp)、吸烟(smoking)、低密度脂蛋白胆固醇(ldl)及冠脉病变syntax评分(syntax)。
二、构建预测模型.
本研究利用构建预测模型(若研究为含有时间变量的生存数据，则可采用cox回归模型)。的操作步骤对大家来说应该早就是小case了，操作方法如下：
1. analyze → regression → binary logistic regression
2. 将因变量event选入dependent框中，将各个自变量选入covariates框中
3. 点击save，在predicted values下勾选probabilities，目的是为了在数据库中新生成一个概率值，用于绘制和校准曲线图。
4. 点击options，勾选hosmer-lemeshow goodness-of-fit，用于输出hosmer-lemeshow拟合优度检验的结果。
三、结果.
variable in the equation中输出了每个影响因素的回归系数(β)、or值、95% ci以及p值等信息。回归方程如下：
logit (p)= -8.713 - 0.899*gender 0.05*age 0.021*sbp 0.912*smoking 0.438*ldl 0.07*syntax
