一种孤立性肺结节恶性概率预测模型的建立方法

摘要

本发明公开了一种孤立性肺结节恶性概率预测模型的建立方法，所述方法具体如下：采集患者基本信息和术前1‑7天血清肿瘤标记物水平，再结合患者的CT影像学报告，将患者病例按毛玻璃样病变GGO比例分为GGO比例大于等于50％和小于50％两组，每组病例按3：1比例均设置实验组和验证组，将实验组病例相关资料进行单因素分析，初步筛选出独立危险因素，带入多因素分析，获得判断SPN良、恶性的独立危险因素，应用Logistics回归获得SPN恶性概率预测模型，再应用验证组病例资料带入模型进行验证。本发明提供的模型简单易用，并且所用的指标均是常规检查就能得出来的指标，易于采用，根据这种模型可以给医生进一步诊断治疗提供有效的中间参考信息。

说明书

技术领域

本发明涉及医学技术领域，具体涉及一种孤立性肺结节恶性概率预测模型的建立方法。

背景技术

目前学术界公认将直径小于或等于3cm，周围为肺组织包绕的圆形影像学结节影定义为孤立性肺结节(SPN)。文献统计，大约每500张X线胸片中可能有1张(0.2％)显示SPN，90％的SPN在健康人群体检中被无意发现，在美国每年可经健康体检发现超过15万个SPN。常见的SPN病理诊断结果包括良性的非特异性肉芽肿、炎性肉芽肿、错构瘤、硬化型血管瘤等，以及恶性的原发性肺癌及肺转移瘤等。国际上已有大样本、前瞻性临床研究证实对偶发的孤立性肺结节在必要时进行外科干预，明确良恶性并行规范治疗，可以降低与肺癌有关的病死率。

目前临床上选择对SPN病例进一步检查甚至外科干预的时机主要基于医师的临床经验，受主观因素影响巨大，不同专业方向的医师往往意见相左(如同样的孤立性肺结节，外科医生通常倾向手术活检，内科医生通常倾向药物治疗并随诊观察，放射科医生通常建议CT引导下穿刺活检)，并无客观、可量化指标供参考。

发明内容

本发明的目的在于提供一种孤立性肺结节恶性概率预测模型的建立方法，该方法准确性高，简单易行，能够给进一步的诊断治疗提供有效的中间信息。

为实现上述目的，本发明方法采用如下技术方案：

一种孤立性肺结节恶性概率预测模型的建立方法，所述方法具体如下：采集患者基本信息和术前1-7天血清肿瘤标记物水平，再结合患者的CT影像学报告，将患者病例按GGO比例分为GGO比例大于等于50％和小于50％两组，每组病例按3：1比例均设置实验组和验证组，将实验组病例相关资料进行单因素分析，初步筛选出独立危险因素，带入多因素分析，获得判断SPN良、恶性的独立危险因素，应用Logistics回归获得SPN恶性概率预测模型，再应用验证组病例资料带入模型进行验证。

具体地，所述患者基本信息包括年龄、性别、吸烟情况(是否吸烟、吸烟量、戒烟时间等)、家族肿瘤病史、既往肿瘤病史等。

所述血清肿瘤标记物水平包括血清癌胚抗原CEA和血清骨胶素Cyfra21-1水平。

所述CT影像学报告中包括结节位置、CT最大径、毛玻璃样病变GGO所占比例、毛刺征、分叶征、钙化、空洞、边界清楚及胸膜牵拉。

所述独立危险因素包括年龄、CT最大径、有无毛刺、是否钙化、有无胸膜牵拉、有无肿瘤史、边界是否清楚、Cyfra21-1和CEA。

所述预测模型具体为：P＝e^X/(1+e^X)；当GGO成分小于50％时：X＝-10.6+(0.064×年龄)+(1.358×最大径)+(1.698×毛刺)-(2.492×钙化)+(1.446×胸膜牵拉)+(2.101×肿瘤史)+(0.639×Cyfra21-1)+(0.633×CEA)。当GGO成分大于等于50％时：X＝-1.180+(0.057×年龄)-(2.30×边界清楚)+(1.032×CEA)；其中e是自然对数；“年龄”的单位是“岁”；“最大径”指CT肺窗上结节最大直径，单位为cm；“毛刺”、“钙化”、“胸膜牵拉”、“肿瘤史”、“边界清楚”为定性指标，“是”计为1，“否”记为0；CEA和Cyfra21-1单位均为ng/mL，所述检测应用医疗单位统一试剂盒。

本发明方法具有如下优点：

以上指标均是通过对病人的询问或者常规检查得出来的，不涉及有创检查或操作，为患者最大限度减少病痛。

利用Logistics回归获得的模型，这种方法基于临床“大数据”方法的应用，可靠性高。

本发明提供的模型简单易用，并且所用的指标均是常规检查就能得出来的指标，易于采用，根据这种模型可以给医生进一步诊断治疗提供有效的中间参考信息，具有较高的参考价值。

附图说明

图1是GGO成分小于50％病例的ROC曲线(左为实验组，右为验证组)；

图2是GGO成分大于等于50％病例的ROC曲线(左为实验组，右为验证组)。

具体实施方式

下面将参照附图更详细地描述本发明的具体实施例。虽然附图中显示了本发明的具体实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。

如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语，故应解释成“包含但不限定于”。说明书后续描述为实施本发明的较佳实施方式，然所述描述乃以说明书的一般原则为目的，并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。

实施例1一种孤立性肺结节恶性概率预测模型的建立方法

申请人采用1234例孤立性肺结节病例，全部病例均获得明确的病理诊断。采集患者一般状况，包括年龄、性别、吸烟情况(是否吸烟、吸烟量、戒烟时间等)、家族肿瘤病史、既往肿瘤病史。采集术前1-7天血清肿瘤标记物水平，包括血清癌胚抗原(carcinoembryonicantigen,CEA)和血清骨胶素(Cyfra21-1)水平。由胸外科、呼吸科、放射科医师共同阅片评估CT影像学特点，包括结节位置、CT最大径(肺窗)、毛玻璃样病变(ground-glass opacity,GGO)所占比例、毛刺征、分叶征、钙化、空洞、边界清楚及胸膜牵拉征。

考虑到结节GGO比例不同会严重影响良恶性判断，将全部病例按GGO比例分为两组(GGO比例大于等于50％或小于50％)。每组病例按3：1比例均设置实验组和验证组。将实验组病例相关资料进行单因素分析，初步筛选出独立危险因素，带入多因素分析，获得判断SPN良恶性的独立危险因素，应用Logistics回归获得SPN恶性概率预测模型如下：

P＝e^X/(1+e^X)；当GGO成分小于50％时：X＝-10.6+(0.064×年龄)+(1.358×最大径)+(1.698×毛刺)-(2.492×钙化)+(1.446×胸膜牵拉)+(2.101×肿瘤史)+(0.639×Cyfra21-1)+(0.633×CEA)。当GGO成分大于等于50％时：X＝-1.180+(0.057×年龄)-(2.30×边界清楚)+(1.032×CEA)；其中e是自然对数；“年龄”的单位是“岁”；“最大径”指CT肺窗上结节最大直径，单位为cm；“毛刺”、“钙化”、“胸膜牵拉”、“肿瘤史”、“边界清楚”为定性指标，“是”计为1，“否”记为0；CEA和Cyfra21-1单位均为ng/mL，所述检测应用医疗单位统一试剂盒。

为了进一步证实本发明所述预测模型的预测效果，申请人对其进行了如下对比验证，具体如下：

我们对该模型预测效力进行了初步检测，绘制受试者工作曲线(receiveroperating characteristic curve,ROC)并对曲线下面积(area under the ROC curve,AUC)进行比较。

对GGO成分小于50％的病例，AUC值在实验组为0.960(95％CI,0.946-0.973)，验证组为0.948(95％CI,0.913-0.983)。对GGO成分大于等于50％的病例，AUC值在实验组为0.847(95％CI,0.784-0.909)，验证组为0.833(95％CI,0.673-0.993)。如图1、图2所示。

根据统计学理论，ROC曲线下的面积值在1.0和0.5之间。在AUC>0.5的情况下，AUC越接近于1，说明试验效果越好。AUC在0.5～0.7时有较低准确性，AUC在0.7～0.9时有一定准确性，AUC在0.9以上时有较高准确性，能够协助选择最佳的识别界限值。ROC曲线越靠近左上角，试验的准确性就越高。最靠近左上角的ROC曲线的点是错误最少的最好阈值，其假阳性和假阴性的总数最少。应用MedCalc软件直接得出，或通过比较各点的约登指数(Youden index)＝敏感性-(1-特异性)，得出对GGO成分小于50％的病例，当选择0.521作为数据截点(Cutoff)的时候，试验效果最佳，敏感度为92.7％，特异性为86.0％。对GGO成分大于等于50％的病例，当选择0.867作为数据截点的时候，敏感度为76.2％，特异性为72.7％。因此，该模型兼备较高的灵敏度和特异性。

实施例2

病例1：男，55岁，体检发现左下叶结节，GGO成分小于50％，CT最大径1.8cm，无毛刺、钙化、胸膜牵拉，既往无肿瘤史，血清肿瘤标记物Cyfra21-1值3.34ng/mL，CEA值2.82ng/mL。

带入预测模型：P＝e^X/(1+e^X)，其中X＝-10.6+(0.064×年龄)+(1.358×最大径)+(1.698×毛刺)-(2.492×钙化)+(1.446×胸膜牵拉)+(2.101×肿瘤史)+(0.639×Cyfra21-1)+(0.633×CEA)。计算得到P＝0.328，小于界值0.521。

经过多方会诊，最后对病人进行手术，病人术后进行病理检测，诊断为肉芽肿性炎，考虑为陈旧性结核可能，系良性结节。

实施例3

病例2，男，79岁，体检发现右下叶结节，GGO成分大于50％，边界不清，血清肿瘤标记物CEA值5.33ng/mL。

带入预测模型：P＝e^X/(1+e^X)，其中X＝-1.180+(0.057×年龄)-(2.30×边界清楚)+(1.032×CEA)。计算得到P＝0.998，大于界值0.867。

在此基础上，继续CT引导下穿刺病理明确诊断为肺腺癌，系恶性结节，接受了后续手术治疗。

实施例4

病例3，女，64岁，体检发现右下叶结节，GGO成分小于50％，CT最大径1.2cm，无毛刺、钙化、胸膜牵拉，既往宫颈癌病史，血清肿瘤标记物Cyfra21-1值2.82ng/mL，CEA值2.14ng/mL。

带入预测模型：P＝e^X/(1+e^X)，其中X＝-10.6+(0.064×年龄)+(1.358×最大径)+(1.698×毛刺)-(2.492×钙化)+(1.446×胸膜牵拉)+(2.101×肿瘤史)+(0.639×Cyfra21-1)+(0.633×CEA)。计算得到P＝0.5946，大于界值0.521。

经过多方会诊，参考上述信息，患者接受手术治疗，术后病理诊断为浸润性腺癌，系恶性结节。

实施例5

病例4，女，71岁，体检发现右上叶结节，GGO成分大于50％，边界清，血清肿瘤标记物CEA值2.78g/mL。

带入预测模型：P＝e^X/(1+e^X)，其中X＝-1.180+(0.057×年龄)-(2.30×边界清楚)+(1.032×CEA)。计算得到P＝0.969，大于界值0.867。

参考上述信息，并且经过多方会诊，患者接受手术治疗，术后病理诊断为微浸润腺癌，系恶性结节。

虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。