本发明涉及数据处理技术，尤其涉及一种基于医保大数据测算疾病患病率(prevalence)的方法，具体针对汇总数据格式下高效计数患病率对应分子、分母进行统计运算。
背景技术：
：医保数据(claimsdata)，是医保业务过程中获得的数据，由于无需抽样调查，天然覆盖较大规模的人群，并较为完整真实的记录该人群在一定时间范围内的就医信息，被越来越多的用于疾病的患病率研究。常规的流行病学研究中，计算患病率需要明确分母的人口规模，计数某一时间内对应的病例数量。但对于现在积累的海量医保数据，无论数据存储的期别、格式和体量，还是数据指标的跨度、缺失和个体脱保等，都会导致直接实施上述思路遇到一些问题，从而制约医保大数据的实际应用。美国的medicaid、medicare联合办公中心总结开发了适用于两者数据库的常见慢性病症的患者界定规则，但同时明确指出目标患者的界定规则因数据库而异、因病而异。美国jensen等对影响分母的界定参数进行了探讨，提出就诊频率不高的病症，最短观察时间设定较短会导致患病率低估，设定较长又会导致高估，具体需要因病而异。上述探索均集中在国外，而我们国家医保报销模式决定了医保数据的实际存储以月为单位，患者参保信息与就诊信息之间以就诊号关联。国外研究都基于个体原始数据，数据量与运算资源性能相对匹配，而我国医保人口超过13亿，多人次多年份的海量数据情境下，传统的构建个体面板数据的运算策略将会导致数据稀疏、信息冗余、存储和运算困难等问题。技术实现要素：本发明的目的是提供一种基于医保大数据测算疾病患病率的新方法，基于优化数据中间存储格式，通过汇总每月医保数据的多个关键参数，包括：每月参保个体总数、每月新增参保个体数、每月就诊记录总数、每月就诊记录诊断缺失总数，再推导患病率计算所需的分子与分母信息，进而计算得到患病率。本发明可测算的疾病要求不能彻底治愈，一旦诊断，终身罹患。本发明的原理是：基于人月的概念，计数每月参保个体总数、每月新增参保个体数、每月就诊记录总数、每月就诊记录诊断缺失总数，通过目标疾病定义抽提目标患者，并推导诊断随机缺失情境下的“隐形患者”，再根据患病率公式进行推导，计算得到患病率。本发明方法可测算的疾病包括多发性骨髓瘤、浆细胞白血病、浆细胞病、男性乳腺癌、黑色素瘤等。本发明提供的技术方案是：一种基于医保大数据测算疾病患病率的方法，通过汇总每月医保数据的多个关键参数(包括：每月参保个体总数、每月新增参保个体数、每月就诊记录总数、每月就诊记录诊断缺失总数)，推导患病率计算所需的分子与分母信息，进而计算得到患病率；患病率计算的分子指特定时间内总人口中的所有患有某病的新旧病例，分母即特定时间内的总人口；包括如下步骤：a1.确定数据库范围(如时间跨度、地域分布、门诊/住院)；a2.数据库的基本清洗和目标疾病的定义；a3.分母信息的汇总具体分为三组：参保但从未报销的个体、参保且有报销记录但未出现目标诊断的个体和参保有报销记录且出现目标诊断的个体。根据每位观察对象在每个月的参保状态、纳入参保人次，剔除未参保人次。a4.分子信息的汇总；明确目标疾病的定义后，进行对应的分子信息抽提，具体分为两组：直接观测到的患者和需填补测算的患者。后者测算基于就诊信息的诊断缺失与是否罹患目标疾病的关联不存在统计学意义。对之前已经判断为目标疾病的个体，则直接填补该个体为患病状态，即为出现目标诊断的个体。a5.分子患者的基本特征核查与统一，如年龄、性别、民族、户籍等a6.患病率的计算，将上述汇总的分子信息与分母信息求商，计算患病率。本发明的有益效果是：本发明提供一种基于医保大数据测算疾病患病率的新方法，针对汇总格式数据，高效地计数得到患病率计算对应的分子和分母，进而计算得到患病率。通过本发明所提供的以医保大数据为基础的患病率测算方法，简便快速有效，可以替代耗费人财物的现有的传统流行病学人群调查，同时转变以面板数据为基础的患病率计算策略，优化数据存储、运算的执行效率，降低隐私泄露风险和数据共享门槛，促进医保大数据的转化应用，为切实填补我国罕见病的患病率资料空白提供一种简便快速的方法。附图说明图1是本发明提供的计算患病率的方法的流程框图。图2是本发明具体实施例提供的患病率分子的基本特征核查与统一的流程示意图。具体实施方式下面结合附图，通过实施例进一步描述本发明，但不以任何方式限制本发明的范围。本发明提供一种基于医保大数据测算疾病患病率的新方法，基于优化数据中间存储格式，通过汇总每月医保数据的多个关键参数，推导患病率计算所需的分子与分母信息，进而计算得到患病率。患病率计算的分子指特定时间内总人口中的所有患有某病的新旧病例，分母即特定时间内的总人口。本发明的具体实施方式如下：a.患病率对应分母的计算a1.确定数据库范围(如时间跨度、地域分布、门诊/住院)；a2.数据库的基本清洗和目标疾病的定义；数据库的基本清洗包括以下基本步骤：(1)数据库中变量的完整性和逻辑性核查；(2)数据库中文本内容的编码标准化和自然语言处理；(3)数据库中国际疾病分类(internationalclassificationofdiseases，icd)的版本确定和统一。本发明中，目标疾病的定义以医保数据库中出现对应疾病的名称或icd编码为准，具体需要充分考虑文本和icd编码的多种表达形式，构建尽量全面的字典库。a3.分母信息的汇总；患病率的分母具体分为三组第一组：参保但从未报销的个体该部分患者从未因病就诊，只有参保记录，无报销记录，在患病率/发病率计算时仅用作分母。具体需要统计观察时间内，每个观察对象在每个月的参保状态(1＝参保，0＝未参保)，然后把未参保的人月剔除(如附图1)，把参保的人月合计放入分母。以每个月为例，第一组分母的人月总和对应计算公式如式1：其中，t表示第t个月份；insurancet,n为第n个该组个体在第t个月份的参保状态。n代表第一组分母的人月总和。第二组：参保且有报销记录但未出现目标诊断的个体该部分患者曾因病就诊，但未出现目标诊断，同时拥有参保记录、报销记录，同样在患病率/发病率计算时仅用作分母。具体需要统计观察时间内，每个观察对象在每个月的参保状态(1＝参保，0＝未参保)，然后同样把未参保的人月剔除(如附图1)，但参保的人月不能直接放入分母，而是根据诊断状态分为三种情况：第一种情况：未因病就诊无报销记录的人月直接纳入分母(如附图1)，对每月而言，即未因病就诊无报销记录的当月人数m1,1；第二种情况：因病就诊且诊断完整的人月应纳入分母计算(如附图1)，对每月而言即因病就诊且诊断完整的当月人数m1,2；第三种情况：有就诊但诊断缺失的人月应考虑后续填补(如附图1)，提取因病就诊但诊断缺失的当月人数m1,3。以每个月为例，第二组分母的人月总和对应计算公式如式2：其中，t表示第t个月份；insurancet,m为第m个该组个体在第t个月份的参保状态。m代表第二组分母的人月总和。第三组：参保有报销记录且出现目标诊断的个体该部分患者有因病就诊，且出现过目标诊断，同时拥有参保记录、报销记录，在患病率/发病率计算时用作分子和分母。就分母而言，具体需要统计观察时间内，每个观察对象在每个月的参保状态(1＝参保，0＝未参保)，然后同样把未参保的人月剔除(如附图1)，但参保的人月仍不能直接放入分母，而是根据诊断状态分为三种情况：第一种情况：未因病就诊无报销记录的人月直接纳入分母(如附图1)，对每月而言即未因病就诊无报销记录的当月人数k1,1；第二种情况：因病就诊且诊断完整的人月应纳入分母计算(如附图1)，对每月而言即因病就诊且诊断完整的当月人数k1,2；第三种情况：有就诊但诊断缺失的人月应考虑后续填补(如附图1)，对每月而言即因病就诊但诊断缺失的当月人数k1,3。以每个月为例，第三组分母的人月总和对应计算公式如式3：其中，t表示第t个月份；insurancet,k为第k个该组个体在第t个月份的参保状态；k代表示第三组分母的人月总和。a4.分子信息的汇总；根据目标疾病的定义后，进行对应的分子信息抽提，具体分为两组：分子第一组：直接观测到的患者所有观测到患者均计入，需要计算每个月诊断为目标疾病的人数之和，即∑tcase_real，其中，t表示第t个月份；case_real表示每个月诊断为目标疾病的人数。分子第二组：需填补测算的患者部分就诊记录存在诊断缺失，但该部分记录不宜直接剔除，如表1所示，表1患病率计算时分子填补示意图理想状态下的患病率prevalence的计算公式为式4：其中，a代表参保且有就诊且就诊记录诊断不缺失的个体中可以抓取到的目标疾病患者人数，b代表参保且有就诊且就诊记录诊断不缺失的个体中可以抓取到的非目标疾病患者人数，c代表参保且有就诊且就诊记录诊断缺失的个体中理论上可抓取到的目标疾病患者人数，d代表参保且有就诊且就诊记录诊断缺失的个体中可以抓取到的非目标疾病患者人数，e代表参保但从未就诊的个体人数。而直接剔除后的患病率计算公式表示为式5：理想状态下的患病率与直接剔除后的患病率两者显然不等。因此，需要对诊断缺失部分，基于一定假设，进行适当估算，以获取c和d的数值。本发明采用的初始假设是就诊信息的诊断缺失与是否罹患某种罕见病的关联不存在统计学意义，即若满足此假设，则c表示为式6：其中，c+d为诊断缺失的总记录数，可直接计数得到。注意，需要填补测算的诊断缺失部分包括两部分：参保且有报销记录但未出现目标诊断的个体中当月有就诊但诊断缺失的人数m1,3和参保有报销记录且既往已经出现目标诊断的个体中当月有就诊但诊断缺失的人数k1,3，其中，m1,3按照上述假设进行填补，而k1,3则按照是否既往已经判断为目标疾病分别处理。若之前未判断为目标疾病，则同样按照上述假设进行填补，若之前已经判断为目标疾病，则直接填补为患病状态。填补后的目标疾病人数可分别为两部分：根据参保且有报销记录但未出现目标诊断的个体中当月有就诊但诊断缺失的人数所估算的目标患者总数根据参保有报销记录且既往已经出现目标诊断的个体中当月有就诊但诊断缺失的人数所估算的目标患者总数其中，t表示第t个月份，caseimpute_m表示每个月内，根据参保且有报销记录但未出现目标诊断的个体中当月有就诊但诊断缺失的人数所估算的目标患者数，caseimpute_k表示每个月内，根据参保且有报销记录且既往已经出现目标诊断的个体中，当月有就诊但诊断缺失的人数所估算的目标患者数。a5.分子患者的基本特征核查与统一；如附图2，医保数据分为“参保人员信息表”、“普通门(急)诊费用及结算信息表”和“门诊大病、门诊统筹、住院、家庭病床费用及结算信息表”3张表单，各表单之间的关联变量为aac001，对年龄、性别、民族、户籍等需要进行多部核查和统一，以达到每个关联变量对应唯一的身份识别id(如身份证)，同时每个唯一的身份识别id对应的年龄、性别、民族和户籍等信息内部一致。a6.患病率的计算当各个月份的患病率波动较小(无显著的季节或其他时间趋势、随访人群分母规模较大时)，某一整年的患病率计算公式如式7(时间跨度为一年，以人月为单位)：其中，每个月的患者总数∑tcase＝case_real+caseimpute_m+caseimpute_k，t表示第t个月份，case_real表示当月内实际抓取到的患者总数，caseimpute_m表示当月内根据参保且有报销记录但未出现目标诊断的个体中，根据有就诊但诊断缺失的人数所估算的目标患者总数，caseimpute_k表示根据参保有报销记录且既往已经出现目标诊断的个体中，根据当月有就诊但诊断缺失的人数所估算的目标患者数。每个月的分母总数计算∑tpersonmonth＝∑tpersonmonth1+∑tpersonmonth2+∑tpersonmonth3。其中，t表示第t个月份，∑tpersonmonth1对应参保但从未报销的个体所贡献的人月，∑tpersonmonth2代表参保且有报销记录但未出现目标诊断的个体所贡献的人月，∑tpersonmonth3代表参保有报销记录且出现目标诊断的个体。当各个月的患病率波动较大时，如部分月份的参保人群基数较小(如少于10万人)，则可采用月患病率的中位数，或者以年为单位计数唯一患者数，公式如式8：其中，t表示第t个月份，每个月的新增患者总数∑tnewcase＝newcasereal+newcaseimpute_m。其中，newcasereal为当月实际观察到新增患者数，newcaseimpute_m为在利用既往月份诊断信息对当前月份诊断信息进行填补后，剩余的参保有报销记录但诊断缺失的记录按照诊断未缺失部分的新增患者比例，估算出的新增患者数。下面通过实例来对本发明做进一步说明。假设基于某省2012年的医保数据库计算多发性骨髓瘤的患病率，则数据库范围为2012年的城镇职工参保人员(143,670,014人)、城镇居民参保人员(22,890,524人)。完成基本的数据清洗后(如报销日期、就诊日期变量缺失、异常等)，多发性骨髓瘤的临床诊断方式结合文本、icd编码表述如表2：表2多发性骨髓瘤的诊断描述和icd编码列举icd-9icd-10icd-o-3多发性骨髓瘤骨髓瘤病203.0c90.051m97320/3浆细胞性骨髓瘤203.0c90.002m97320/3多发性骨髓瘤203.0c90.001m97320/3骨髓瘤肾病203.0c90.003+m97320/3数据库中包含诊断信息的字段名称共有6个，分别是主要诊断名称、主要诊断编码、次要诊断名称1、次要诊断编码1、次要诊断名称2和次要诊断编码2，则实际定义按照数据库内字段结构定义如下：上述各字段必须包含字段取值(字段之间为“或者”关系)：骨髓瘤，卡勒，骨髓癌/骨髓ca，骨髓病，c90，m9732，203.0；上述各字段必须排除字段取值(字段之间为“或者”关系)：浆细胞，孤立性，c90.1，c90.2具体分子抓取时，针对主要诊断、次要诊断1、次要诊断2、主要诊断编码、次要诊断编码1和次要诊断编码2共六个字段展开，要求全部六个字段中，必须具备至少一个“必须包含字段”，但必须不含有“必须排除字段”。然后进行分母信息的汇总，其中参保但从未报销的个体、参保且有报销记录但未出现目标诊断的个体和参保有报销记录且出现目标诊断的个体在每个月的分布，如表3所示：表3分母信息的汇总表之后进行分子的提取和信息核查与统一，其中直接观察到的患者信息如表4：表4患者信息汇总表每个月的患病率计算结果如表5：表5每个月的患病率计算结果分子(人数)分母(人数或人月)患病率(/10万)120126621660.95180122869091.46199132531031.50225134052491.68292155608651.88508153317033.31469158337212.96491157132053.12404154140552.62463157320582.9445697105854.70210116569191.8040171665605382.41各个月份的患病率相对稳定，因此最终计算为:即患病率为2.41/10万人月，若直接采用月患病率的中位数为2.25/10万人月。需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。当前第1页12