预测糖尿病患者心血管并发症动态发病风险的方法及系统与流程

本发明涉及预测发病风险的技术领域，尤其是指一种预测糖尿病患者心血管并发症动态发病风险的方法及系统。

背景技术：

从20世纪70年代开始，多个国家和地区开展了大规模的心血管疾病风险因素的研究，并先后推出了多个心血管疾病风险评估方法，如英国前瞻性糖尿病预测模型，该模型用cox比例风险模型建立，通过糖尿病患者的性别、糖尿病诊断年龄、吸烟情况，以及收缩压、胆固醇、低密度脂蛋白、糖化血红蛋白的值，预测2型糖尿病患者的冠心病发病风险，健康经济学家、临床医生和糖尿病患者等可免费使用模型的软件包。

现有的糖尿病并发症预测模型使用患者的横断面信息进行预测，没有考虑患者预测变量的动态变化(如多次测量的血糖值)和累计效应，对于变化的生理指标，预测风险不能动态的更新。另外，国外这些预测模型的建立基于已有的研究队列人群，由于这些人群和中国人群有较大的种族差异和并发症流行率的差异，导致人群的基线发病风险不同，使这些模型不适用于中国糖尿病人群的并发症预测。

技术实现要素：

为此，本发明所要解决的技术问题在于克服现有技术中无法对个体进行心血管疾病风险动态预测的问题，从而提供一种对个体进行心血管疾病风险动态预测，以达到预防并阻止病情发展目的的预测糖尿病患者心血管并发症动态发病风险的方法及系统。

为解决上述技术问题，本发明的一种预测糖尿病患者心血管并发症动态发病风险的方法，包括：从多个卫生系统中收集、整合糖尿病人群的公共卫生数据；构建贝叶斯多变量联合模型并对所述模型进行评估；利用评估后的模型对个体进行预测。

在本发明的一个实施例中，所述的公共卫生数据包括不同时间多次重复测量的血糖、血压、血脂、尿酸等纵向指标，冠心病、脑卒中等心血管疾病是否发病及发病时间的结局变量。

在本发明的一个实施例中，所述数据将通过身份证号利用软件sas进行匹配、整合。

在本发明的一个实施例中，所述构建贝叶斯多变量联合模型的方法为：假设在多个时间点收集的糖尿病患者的纵向变量，与发生心血管并发症有关，使用贝叶斯多变量联合模型拟合纵向数据和结局数据，所述联合模型包括纵向子模型和生存子模型，其中所述纵向子模型采用线性混合效应模型，所述生存子模型采用cox比例风险回归模型。

在本发明的一个实施例中，所述纵向子模型公式：yki(t)＝xki(t)βk+zki(t)bki，其中yki(t)表示在时刻t第i个糖尿病患者的第k个变量的变量值，xki(t)和zki(t)分别是固定效应βk和随机效应bki的系数。

在本发明的一个实施例中，所述生存子模型公式：hi(t)＝h0(t)exp[γ·wi(t)+α·fi^*(ti)]，其中hi(t)表示第i个糖尿病患者在时刻t患心血管疾病的风险函数，h0(t)表示基线风险函数，wi(t)是协变量，相应的系数是γ，fi^*(ti)是与纵向变量相关得函数，关联参数α将两个子模型联系起来，假设生存结局与纵向变量在时点的期望值、在时点的变化趋势、在全部时间的累计效应有关。

在本发明的一个实施例中，所述模型的拟合用软件的程序包进行。

在本发明的一个实施例中，对所述模型进行评估的方法为：模型拟合优度评价采用对数似然函数值、aic信息准则和bic信息准则，当所述aic信息准则和bic信息准则数值越小，而所述对数似然函数值越大，表示模型拟合越好。

在本发明的一个实施例中，所述对个体进行预测时，对糖尿病患者心血管并发症发病风险有影响的是全部时间空腹血糖的累计效应、全部时间收缩压的累计效应、最近一次检查的尿酸值，以及糖尿病患者的性别、糖尿病诊断时间、身体质量指数。

本发明还提供了一种预测糖尿病患者心血管并发症动态发病风险的系统，包括：数据收集模块，用于从多个卫生系统中收集、整合糖尿病人群的公共卫生数据；构建模型模块，用于构建贝叶斯多变量联合模型并对所述模型进行评估；构建模型模块，用于构建贝叶斯多变量联合模型并对所述模型进行评估。

本发明的上述技术方案相比现有技术具有以下优点：

本发明所述的预测糖尿病患者心血管并发症动态发病风险的方法及系统，通过糖尿病不良结局心血管疾病的动态发病风险预测，社区糖尿病人群能够及时获知自己的心血管疾病发生风险，使其意识到疾病的严重性和易感性，促使糖尿病患者主动寻求自我管理及提高遵医嘱依从性，进而预防并阻止病情的发展，是降低其发生、促使其向良好的结局转归、减少疾病负担的关键。

附图说明

为了使本发明的内容更容易被清楚的理解，下面根据本发明的具体实施例并结合附图，对本发明作进一步详细的说明，其中

图1是本发明预测糖尿病患者心血管并发症动态发病风险的方法流程图；

图2是本发明某糖尿病患者的心血管疾病的发生风险的示意图。

具体实施方式

实施例一

如图1所示，本实施例提供一种预测糖尿病患者心血管并发症动态发病风险的方法，包括如下步骤：步骤s1：从多个卫生系统中收集、整合糖尿病人群的公共卫生数据；步骤s2:构建贝叶斯多变量联合模型并对所述模型进行评估；步骤s3：利用评估后的模型对个体进行预测。

本实施例所述预测糖尿病患者心血管并发症动态发病风险的方法，所述步骤s1中，从多个卫生系统中收集、整合糖尿病人群的公共卫生数据，有利于适用于中国糖尿病人群的并发症预测；所述步骤s2中，构建贝叶斯多变量联合模型并对所述模型进行评估，从而有利于对个体进行心血管疾病风险动态预测，以达到预防并阻止病情发展的目的；所述步骤s3中，利用评估后的模型对个体进行预测，从而可以预测糖尿病患者心血管疾病的发生风险，由于通过基于建立的联合模型对个体进行心血管疾病风险动态预测，因此可以达到预防并阻止病情发展的目的。

为了实现对中国糖尿病人群的并发症预测，所述的公共卫生数据包括不同时间多次重复测量的血糖、血压、血脂、尿酸等纵向指标，冠心病、脑卒中等心血管疾病是否发病及发病时间的结局变量，由于使用预测人群重复测量的体检数据，考虑了预测变量的动态变化和累积效应，可以动态的更新预测风险。所述数据将通过身份证号利用软件sas进行匹配、整合。

下面详细论述如何从多个卫生系统中收集、整合糖尿病人群的公共卫生数据：

对某市某区65周岁以上参加体检的糖尿病人群进行心血管并发症的动态风险预测：用到的卫生系统有体检数据库、慢性病管理数据库、糖尿病健康管理档案、死亡报告数据库。从体检数据库中提取不同时间多次重复测量的血糖、血压、血脂、尿酸等指标；糖尿病健康管理档案中提取血糖、用药等数据。需要收集的数据还包括慢性病发病数据库中的冠心病、脑卒中等是否发病及发病时间，以及死亡报告数据库中死因、发病时间、死亡时间等。所述数据将通过身份证号利用软件sas进行匹配、整合。以上收集的数据分为纵向变量和结局变量进行分析，纵向变量包括动态变量(比如血糖、收缩压、尿酸等)和静态变量(比如性别、吸烟、饮酒)，结局变量是脑卒中、冠心病等心血管疾病是否发病和发病时间。

所述构建贝叶斯多变量联合模型的方法为：假设在多个时间点收集的糖尿病患者的纵向变量，与发生心血管并发症有关，使用贝叶斯多变量联合模型拟合纵向数据和结局数据，所述联合模型包括纵向子模型和生存子模型，其中所述纵向子模型采用线性混合效应模型，所述生存子模型采用cox比例风险回归模型。所述纵向子模型公式：yki(t)＝xki(t)βk+zki(t)bki，其中yki(t)表示在时刻t第i个糖尿病患者的第k个变量的变量值，xki(t)和zki(t)分别是固定效应βk和随机效应bki的系数。所述生存子模型公式：hi(t)＝h0(t)exp[γ·wi(t)+α·fi^*(ti)]，其中hi(t)表示第i个糖尿病患者在时刻t患心血管疾病的风险函数，h0(t)表示基线风险函数，wi(t)是协变量，相应的系数是γ，fi^*(ti)是与纵向变量相关得函数，关联参数α将两个子模型联系起来，假设生存结局与纵向变量在时点的期望值、在时点的变化趋势、在全部时间的累计效应有关。模型的拟合和个体预测用软件r的程序包“jmbayes”进行。

下面详细说明如何构建贝叶斯多变量联合模型：

假设在多个时间点收集的糖尿病患者的纵向变量，与发生糖尿病并发症心血管疾病有关，使用贝叶斯多变量联合模型拟合纵向数据和结局数据。模型的拟合用软件r的程序包“jmbayes”进行。先对空腹血糖、收缩压、舒张压、胆固醇、甘油三酯、高密度脂蛋白、低密度脂蛋白、尿酸、心率等变量拟合单变量联合模型，考虑变量的时点期望值、时点变化趋势和在全部时间的累计效应，结合单变量联合模型的结果拟合贝叶斯多变量联合模型。拟合好多变量联合模型，对糖尿病患者心血管并发症发病风险有影响的是全部时间空腹血糖的累计效应、全部时间收缩压的累计效应、最近一次检查的尿酸值，以及糖尿病患者的性别、糖尿病诊断时间、身体质量指数，公式为：h＝h0exp(γ1·qzage+γ2·bmi+γ3·sex+α1·y血糖+α2·y收缩压+α3·y尿酸)，h表示糖尿病患者患心血管疾病的风险函数，h0表示基线风险函数，qzage表示糖尿病诊断时间，bmi表示身体质量指数，sex表示性别，γ1、γ2、γ3是相应的系数，y血糖表示空腹血糖值，y收缩压表示收缩压的值，y尿酸表示尿酸的值，α1、α2、α3分别为对应的关联参数。

对所述模型进行评估的方法为：模型拟合优度评价采用对数似然函数值、aic信息准则和bic信息准则，当所述aic信息准则和bic信息准则数值越小，而所述对数似然函数值越大，表示模型拟合越好。用roc曲线下面积、brier分数、平均绝对误差评估模型的预测能力。

所述对个体进行预测时，对糖尿病患者心血管并发症发病风险有影响的是全部时间空腹血糖的累计效应、全部时间收缩压的累计效应、最近一次检查的尿酸值，以及糖尿病患者的性别、糖尿病诊断时间、身体质量指数。具体地，使用累计三个纵向变量(空腹血糖、收缩压和尿酸水平)数据和糖尿病患者心血管并发症的发生情况，建立了糖尿病患者心血管并发症动态发生风险的预测模型，利用所述模型预测1名糖尿病患者心血管疾病的发生风险，从图2中可以看到该名糖尿病患者的心血管疾病的发生风险。

实施例二

基于同一发明构思，本实施例提供了一种预测糖尿病患者心血管并发症动态发病风险的系统，其解决问题的原理与所述预测糖尿病患者心血管并发症动态发病风险的方法类似，重复之处不再赘述。

本实施例所述的预测糖尿病患者心血管并发症动态发病风险的系统包括：

数据收集模块，用于从多个卫生系统中收集、整合糖尿病人群的公共卫生数据；

构建模型模块，用于构建贝叶斯多变量联合模型并对所述模型进行评估；

预测模块，利用评估后的模型对个体进行预测。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

显然，上述实施例仅仅是为清楚地说明所作的举例，并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。