今天我们来讲一下如何用SPSS和Medcalc来绘制ROC曲线,并找到截断值。
假设某实验室人员选用一种新方法对同一批患者进行了诊断,患者的真实状态已知。actual代表真实患病状态,即以金标准(如病理诊断)判断是否患病的结果,其中1代表患病,0代表不患病;assay1代表一种新检验方法测量的结果。
一、用SPSS进行ROC曲线的绘制
1. 选择Analyze→ROC Curve…
2. 主对话框设置:将已知的疾病情况actual送入State Variable框中,assay1送入Test Variable中,并在Value of State Variable框中填1(即代表患病状态),在Display中选择所有,其中Coordinate points of the ROC Curve会给出曲线上的坐标点→OK。
3. 结果
(1) ROC曲线
(2) ROC曲线下面积:从Area Under the Curve的结果可知,assay1的ROC曲线下面积为0.856(95%置信区间:0.825-0.886,P<0.001)。
(3) ROC曲线上的坐标点:如下图所示,我们可以根据Coordinates of the Curve的结果可以得到一系列灵敏度和1-特异度的值。
要想获得截断值,就是最接近左上角(0,1.0)的点所对应的坐标点,我们可以将这两列值复制到Excel表中,根据正确指数最大选出最佳临界点。
正确指数又称约登指数(Youden’s index),表示检验方法发现真正病人与非病人的总能力,是灵敏度与特异度之和减去1,即约登指数=灵敏度+特异度-1,在Excel中,用灵敏度-(1-特异度)得到的就是约登指数,对相减的结果进行排序,可以得到正确指数的最大值,即最佳临界点。
确认删除