提高诊断试验效率:多变量综合诊断模型的应用——SPSS软件实现

发布于 2022年6月6日 星期一 10:22:05 浏览:2549
原创不易,转载请注明来源,感谢!
附件下载:
1654482163-多变量综合诊断模型的应用.zip 请勿重复点击,如无响应请耐心等待或稍后再试。

在多变量综合诊断文章中以logistic回归模型为例简要介绍了多变量综合诊断分析基于SPSS软件(提高诊断试验效率:多变量综合诊断——SPSS软件实现)及Medcalc软件的实现过程(提高诊断试验效率:多变量综合诊断——Medcalc软件实现),在本篇推文中将仍以该以logistic回归模型为例,进一步介绍如何更好地利用该模型做出临床诊断。

关键词:诊断试验; 多变量综合诊断; logistic回归; L评分

在多变量综合诊断分析的知识小贴士中提到,logistic回归模型需要通过将就诊者各参数信息代入公式后计算得到Logit(P),再通过与模型截断值比较,做出是否患病的判断。本文便介绍该方法的具体运用。

一、案例简述

某课题组通过就诊者的基本特征、实验室检查等多个指标构建了评价该就诊者是否患病的logistic回归模型为:Logit(P) = -4.068 + 1.617×年龄(60岁以下=0,60岁及以上=1) + 1.935×是否吸烟(不吸烟=0,吸烟=1)+ 0.420×肿瘤标志物A。现欲探索如何将该模型便捷地用于临床诊断。构建模型时共纳入424例就诊者,包括就诊者的年龄(60岁以下=0,60岁及以上=1)、是否吸烟、是否饮酒、肿瘤标志物A;肺癌诊断的金标准为穿刺活检病理检查(非肺癌=0,肺癌=1)。部分数据见图1,本文案例可从“附件下载”处下载。

表1

二、问题分析

本案例的研究目的是将综合诊断模型便捷地应用于临床诊断,可以将模型简化得到L评分计算公式,利用该公式计算原始数据中每位就诊者的L评分,然后进行ROC分析找到模型最优时对应的L评分,即截断值,后面只要通过计算就诊者的L评分就可以进行疾病诊断。

三、软件操作及结果解读

(一) 计算L评分

1.软件操作

从(提高诊断试验效率:多变量综合诊断——SPSS软件实现,链接)中得到诊断“是否患肺癌”的危险得分Logit(P) = -4.068 + 1.617×年龄(60岁以下=0,60岁及以上=1) + 1.935×是否吸烟(不吸烟=0,吸烟=1)+ 0.420×肿瘤标志物A。为了便于应用,可以把模型系数同时除以一个可以让模型变简单的数值(一般选择除以最小的系数),Logit(P)变为L评分,本例中除以0.42,得到L评分= -9.7 + 3.9×年龄(60岁以下=0,60岁及以上=1) + 4.6×是否吸烟(不吸烟=0,吸烟=1) +肿瘤标志物A。如果原危险得分公式已经很简单易算,此步骤可以省去,直接将Logit(P)计算公式作为L评分的计算公式。

然后将所有参与模型构建的就诊者年龄(60岁以下=0,60岁及以上=1)、是否吸烟情况和肿瘤标志物A数据代入L评分计算公式,得到新变量“L评分”,具体操作如下:

  • 点击“转换”—“计算变量”(图2)。
图2
  • 在“计算变量”对话框进行L评分的计算:“目标变量”下录入生成的新变量名称“L评分”,右侧数字表达式中录入L评分的计算公式:“-9.7 + 3.9 * age + 4.6 * smoking+ marker_A”,然后点击“确定”(图3)。
图3

2.结果解读

在“数据视图”中可以看到新生成变量“L评分”(图4)。

图4

(二) 确定截断值

通过ROC曲线可以选择最佳的诊断截断值,ROC曲线越靠近左上角,试验的准确性就越高。

1.软件操作

  • 选择“分析”—“分类”—“ROC曲线”(图5)。
图5
  • 在“ROC曲线”对话框中,“检验变量”中选择“L评分”,“状态变量”中选择“金标准诊断结果”,并将其下方的“状态变量值”设为“1”。在下方“显示”部分全部勾选,最后点击“确定”(图6)。
图6

2.结果解读

图7是通过“L评分”生成的ROC曲线。

图7

图8显示ROC曲线下面积为0.860 (95%CI:0.823~0.896)。

图8

图9是不同L评分下的“敏感度”(一般写为“灵敏度”)和“1-特异性”(一般写为“特异度”),即ROC曲线上的坐标点。

图9

要想获得诊断截断值,首先要计算约登指数(又称为正确指数,Youden index),约登指数越高,该诊断方法发现真正病人与非病人的能力越强。约登指数=灵敏度+特异度-1=灵敏度-(1-特异度),可以将图9中数值复制到Excel表中进行计算,然后对算出的约登指数做降序排列,找到约登指数最大为0.586,此时对应的点就是图7曲线中最接近左上角(1.0,1.0)的点,此时灵敏度为0.730,特异度为1-0.144=0.856,L评分为2.715,即为通过L评分诊断就诊者是否患病的截断值(图10)。

图10

(三) 模型应用

由以上分析结果可知,当L评分<2.715时,可以认为该就诊者未患病,当L评分>2.715时,可以认为该就诊者患病。假设临床中遇到1位就诊者,其年龄为45岁,有吸烟史,肿瘤标志物A的数值为7。则该就诊者的L评分= -9.7 + 3.9×0 + 4.6×1 + 7 = 1.9 < 2.715,因此可以诊断该就诊者未患病,即可快速实现模型的应用。

四、结论

通过对综合诊断分析得到的logistic回归模型进行简化,得到L评分= -9.7 + 3.9×年龄(60岁以下=0,60岁及以上=1) + 4.6×是否吸烟(不吸烟=0,吸烟=1) + 肿瘤标志物A。将案例中所有就诊者的年龄(60岁以下=0,60岁及以上=1)、是否吸烟和肿瘤标志物A的数据代入该公式可以计算得到L评分。通过L评分的ROC曲线分析,得到曲线下面积为0.860 (95%CI:0.823~0.896),也得到每个L评分对应的曲线上坐标(“敏感度”和“1-特异性”),通过坐标值计算得到约登指数最大为0.586,对应L评分为2.715。

综上,只要根据就诊者的年龄、是否吸烟及肿瘤标志物A算出的L评分>2.715,就诊断为该就诊者患病;若L评分<2.715,可认为该就诊者未患病。这种诊断方法的灵敏度为0.730,特异度为0.856。

End
文章目录 沉浸式阅读