我第一次接触GBD全球疾病负担研究数据是在做某个人群健康评估的时候。那时候团队需要回答一个特别现实的问题资源有限的情况下到底该优先干预哪些健康问题翻来覆去找数据最后发现没有任何一份现成的统计年鉴能直接回答直到有人甩给我一份GBD的Excel导出文件说“这是全球疾病负担研究的结果你自己按需拉数吧”。当时我对着那堆指标缩写其实是一头雾水的DALY、YLD、YLL、UI区间、SDI分层每个词都见过但组合到一起就不知道该怎么解读。后来啃了不少方法学资料也亲手跑了不知道多少轮数据提取才慢慢摸清楚这套体系的门道。这篇就当作一个踩过坑的从业者给你们做的GBD入门拆解不讲空话尽量把“这东西是什么、怎么算、怎么用、有什么坑”讲透。1. GBD到底在解决什么问题1.1 GBD不是一本统计年鉴而是一套数据治理框架很多第一次接触的人会把GBD理解成“一个很大的健康数据库”用的时候感觉像查字典点开网页选好病种导出数据完事。这么理解不算全错但它会带来一个很大隐患你并不清楚手里的数字是怎么来的也就很难判断这些数字到底能不能用在自己的研究场景里。GBD本质上是一套全球统一的健康数据测量和估算体系由健康指标与评估研究机构IHME牵头联合全球上千名研究人员持续运作。它要解决的核心矛盾非常简单每个地区都缺数据但决策者不能等数据齐了才做判断。于是GBD采用了一整套多层次的数据整合、交叉校验和模型估算方法把全球每个地区、每个年龄段、两种性别、几百种疾病和伤害的健康损失统一换算成可比的指标。这个思路让我想到做天气预报地面观测站再多也有盲区但你不会因为某个地方缺气象站就拒绝发布天气预报总得靠观测加模型把整片天空拼出来。GBD干的就是这件事。所以项目标题里的“介绍”两个字实际意味着要理解两件事第一这套数据怎么读第二这套数据是怎么被生产出来的。如果只学会点网页下载那你拿到的只是一堆数字而不是研究结果。真正会用GBD的人通常看得懂估算链路知道哪些数字是实测扎实的哪些数字是模型硬推出来的。1.2 谁在用、主要用在哪GBD的数据使用场景比想象中广。做流行病学研究的会拿它来分析某一类疾病在不同人群中的发病和死亡趋势做卫生经济学的会用它给干预措施算成本效益做临床研究的会用疾病负担数据来论证某个研究方向的重要性方便在基金申请里讲清楚立项价值做得更细的还有在做康复、精神卫生、伤害预防这类“在传统死亡统计里看不见”的领域因为单看死亡率会严重低估这些健康问题的影响只有把失能维度也拉进来才说得清楚。我自己感受最深的一点是GBD特别适合用来回答“横向比较”和“纵向趋势”这两类问题。横向比较指的是不同地区之间、或者同一个地区里不同疾病之间的负担排序纵向趋势则是指看十年、二十年的变化轨迹。这两类分析在医院内部报告、科研论文、卫生规划报告里出现频率特别高。你要是能把GBD数据从原始下载到最终成图全流程跑通基本就掌握了一项通用的技能换到多少种疾病、多少层级的地区套路都是同一套。2. 核心逻辑一个负担数字是怎么算出来的2.1 DALY把“早死”和“失能”放到同一杆秤上GBD最核心的产出是DALY伤残调整生命年它的计算是这套体系的中轴。DALY这个概念是1990年代初由世界银行和世界卫生组织的专家们推动成型的后来在GBD项目里得到系统应用和迭代。DALY的设计初衷很简单健康损失不只有死亡这一种表现有些人因病长期卧床有些人带着严重症状活了几十年这些人并没有“死亡”但他们同样损失了大量健康生命年如果指标只衡量死亡就会把这些损失完全漏掉。DALY把两类损失加起来一类是过早死亡损失的生命年YLL另一类是带病生存期间损失的健康年YLD。用一句话可以概括成一个人因某种疾病死亡了算他损失了多少本该拥有的健康寿命一个人没有死但因为疾病处于某种失能状态就需要把他带病生活的年头折算成相当于损失了多少完整健康的年头。两者相加得到的就是这种疾病导致的总负担单位统一为“年”。这个思路相当于把所有健康损失换算成同一种货币然后才能拿不同疾病之间进行直接比较。举个例子可能更容易理解。某种癌症的负担可能主要来自YLL因为确诊后很多人寿命显著缩短而某种慢性腰背痛负担则几乎全部来自YLD因为很少有人直接因为这个疾病死亡但受到影响的人非常多每个人损失的健康年虽不致命乘以庞大的人群基数以后总量非常惊人。2.2 YLL与YLD的估算逻辑YLL的计算需要两个核心输入死亡人数和期望寿命。GBD在做估算的时候针对不同年龄段的死亡都有一套对应的标准期望寿命作为参照一个人越年轻死亡损失的YLL就越大。这其实很公平一个十岁孩子和一个七十岁老人因同种疾病死亡对应的YLL完全不同因为各自剩余的本应健康生存的年数差异很大反映在DALY上就是完全不同的负担量级。YLD部分的计算则比YLL复杂得多因为它不只是看有多少人患病还要看这些患者具体处于什么状态、这种状态带来了多大程度的健康损失。GBD为每一种疾病阶段定义了一套健康状态的失重权重范围在0到1之间0代表完全健康1代表等同于死亡。比如某个疾病导致的轻度症状权重可能只有0.05说明患者损失了5%的健康而终末期疾病权重可能接近0.8意味着那个状态非常接近“完全没有健康生命”。然后用患病率乘以权重再乘以持续时间得到YLD。这个环节是整个GBD估算体系里最容易产生误差的地方因为失能权重本身来自全球多国的调查数据不同人群对“什么是好健康状态”的感受其实存在差异模型已经做了大量标准化处理但不确定性依然存在这也是为什么GBD每个指标都会给出95%不确定性区间。2.3 年龄标准化和两个关键思路除了DALYGBD还提供发病率、患病率、死亡率、期望寿命、健康期望寿命等多种指标不同研究场景侧重点不一样。做病因分析可能更关注发病率和患病率做最终健康结果比较就看DALY和死亡率想看人群整体健康水平则用期望寿命和健康期望寿命HALE。这几个指标各有各的视角单独用哪一个都无法覆盖全部信息组合起来才能讲完整故事。在跨地区和跨时间比较时年龄标准化是绕不开的一个操作。因为不同地区的人口年龄结构可能差异很大一个老龄化严重的地区癌症死亡人数天然会高于年轻人口为主的地区这不是因为它的医疗或环境更差纯粹是年龄结构造成的。GBD官方下载工具里会同时提供粗率和标化率如果你要做区域比较请务必选年龄标准化率否则得出的结论很可能是错的。另外需要提一下的是GBD 2010版本之后官方做了一个重要调整取消了以前采用的年龄权重和时间贴现。早年的DALY计算里中间年龄段健康年的权重会被认为更高未来年份的损失也会被打折折现。后来经过大量争论学界普遍认为这种做法不够透明并且带有明显的价值判断所以现在GBD报告的都是未经过年龄权重和时间贴现的“纯”健康年损失。这个调整让不同年代的结果之间有可比性你在看老文献时如果发现数值和现在下载的不一致很可能就是算法版本导致的不一定是数据错误。3. 数据从哪来GBD喂给模型的“食材清单”3.1 上游数据源普查、登记、监测和文献GBD之所以能覆盖全球几乎所有地区前提是它搜集了大量来源的数据。简单分一下类的话上游数据大概有四种人口普查和调查数据、人口动态登记系统数据、疾病监测系统数据以及已发表的流行病学研究文献。第四种可能最容易被忽略但在一部分地区正式的死亡登记系统覆盖并不完整大量信息来自科学家发表的研究论文GBD团队会做系统检索把那些报告了患病率、发病率、死亡率的研究提取进数据库再参与后续的估算。这带来的结果就是数据的“底层质量”参差不齐。有的地区死亡登记系统非常完善每个死亡个案都有明确死因有的地区可能连人口普查都是多年以前的数据上面的数字还未必可靠。GBD的能力不在于让这些数据变完美而在于用统计手段把不完整、不一致的数据整合到一起并准确告诉使用者“我们的估计有多大的不确定性”。所以我在看GBD每次更新的版本说明时都会重点观察某个地区或疾病的数据来源数量有没有增加数据来源多了估计就会更扎实UI区间也会明显变窄。3.2 中游标准化死因推断与ICD编码映射原始数据进来以后首先要做标准化处理这里面有两个特别关键的步骤死因推断和编码映射。死因推断针对的是没有医学诊断死因的地区。很多死亡发生在医疗体系覆盖不到的地方家人只能大概描述死者生前的症状这时候需要通过死因推断方法把描述性文本对应到可能的死因类别。GBD设计了专门的算法来做这件事很大程度上提高了死因数据的覆盖范围但推断本身带有概率性所以对应的死亡负担数据通常不确定性更大。编码映射则是把不同来源、不同年份、不同编码系统下的诊断和死因名称统一到GBD自身的一套疾病层级目录里。国际疾病分类ICD从第9版、第10版到第11版不同地区的使用习惯和扩展码都不一样某些旧数据用的还是非常粗略的分类名称GBD团队必须一个个做映射。这个环节听着琐碎实际上是最影响数据质量的环节之一如果映射错了后续模型吃得再准也白搭。3.3 下游模型DisMod-MR 2.1和CODEm标准化之后的数据不会直接变成输出指标中间还要经过估算模型的处理其中两个模型一定要知道DisMod-MR 2.1和CODEm。DisMod-MR 2.1是个贝叶斯元回归工具主要用来估计非致命健康结果也就是发病率、患病率、持续时间、缓解率这些参数。它厉害的地方在于可以做多层次的跨地区信息借用数据充足的地区更多依赖本地数据数据缺乏的地区则更多借用周边或全球层面的信息然后综合得出估计值。这套机制的假设是“同一类疾病的流行病学特征在全球范围内存在一定的相似性”虽然不完全对但在缺数据场景下已经是目前最可行的办法。CODEm则是处理死因数据的模型专门用来估计各类疾病和伤害的死亡人数。它会同时考虑多种与死亡相关的因素在不同候选模型里做交叉验证选出预测表现最好的组合再综合多个模型的结果。用大白话说就是在“根据已知死因数据往回推未知死因”的过程中让多套模型互相投票而不是赌某一个模型的运气。这也是GBD数据中死亡估计相对比较稳定的原因之一。4. 实操从GBD官网下载数据到完成第一张趋势分析表4.1 获取数据的主要途径搞懂指标体系之后真正的动手环节就来了。现在获取GBD数据主要有几种方式我用下来觉得各有优缺点。最常规的是GBD Results Tool也就是IHME官网的在线下载工具适合做按需取数。选好地区范围、年龄组、性别、年份、指标和疾病原因它会生成一张交叉表可以导出CSV格式的数据。这个工具适合不需要太复杂变量的常规分析我自己大多数时候都用它。如果下载需求比较批量或者需要做自动化分析更推荐直接从GHDx全球健康数据交换平台下载完整结果数据集里面的文件通常很大但字段更全适合用Python或R进行灵活拆分。还有一个工具是GBD Compare它本质是个可视化交互界面适合快速看图找方向不适合作为研究中提取最终数据的来源。做严肃分析我还是建议用Results Tool或GHDx因为导出的数据才是结构化、带完整字段的。4.2 下载界面的核心选项怎么选以Results Tool为例进去以后要设置的那套条件看起来很多但只要想清楚分析目标就不难。Years的选项直接决定你能做多长时间跨度的趋势分析官方提供的数据覆盖从1990年到最近发布年份通常会选择完整区间至少得有十年才能看出趋势变化。Location选择分析的地理单元这个体系是多层级的全球、区域层级、国家、国家以下行政区逐级向下。你一定要先想清楚自己的研究问题聚焦在哪一层比如想做全球层面的整体比较选Global或者区域层级就够了想做某个特定国家或地区内部的差异就得选更细的行政区划。每下拉一级背后可用的数据和模型稳定性都会不同这很正常。Cause选择疾病原因这是GBD的疾病层级目录分为三级。Level 1是最粗的分类比如传染病、孕产妇、新生儿和营养性疾病非传染性疾病伤害三大类往下细分到Level 2比如肿瘤、心血管疾病、慢性呼吸系统疾病再到Level 3就是具体的单个病种比如肺癌、缺血性心脏病。建议选Level 2或Level 3作为分析粒度这样既有足够的区分度又不会碎片化到难以解读。Metrics选项里有Number例子数、Percent构成比、Rate每10万人率三种维度实际研究中Number适合看总量Rate则适合比较不同人口规模的地区。我一般两个都会下载一个看结构一个看强度。4.3 下载后的数据结构与预处理要点导出的CSV打开以后你会看到有很多行列反复重复的组合字段年份、地区、性别、年龄组、疾病、指标呈长表格式排列。首先要做的就是了解字段含义然后根据需要做数据透视或者宽表转换。在做数据清洗的时候有几个小点容易被忽视。第一是年龄组编码“All Ages”和各个具体年龄组是分开的行如果你自己在做年龄标准化只按“每10万人率”里的年龄标准化率这一列取数就好不要自己用粗率再去套什么权重。第二是性别字段“Both”这种合计行和Male、Female分开存在做合计分析时记得选Both别把男女行相加再除以二那完全是多此一举。第三就是UI区间上下界字段官方会给出低值和高值做成图时可以考虑画误差线这样审稿人会有更多安全感。实际分析过程中我最常用的操作是从长表里筛选出自己关心的疾病子集和年份范围然后计算DALY构成比也就是每种疾病占该地区总负担的比例再排一个序。这个表做出来以后你通常能一眼看出来该地区最主要的几个健康问题是什么再叠加一个分年份的趋势折线就能看到哪些疾病负担在上升、哪些在下降整个逻辑就闭环了。4.4 要让结果可信必须学会看UI区间这可能是新手最容易忽略的地方。GBD是一个模型估算系统所以每个输出值都带有一个95%不确定性区间UI这个区间的宽度直接反映了一个估计值的可信程度。如果两个疾病的点估计排第一第二但它们的UI区间完全重叠那你就不能很自信地说谁负担更重正确的说法应该描述为“两者负担相近实际上尚无证据表明存在显著排序差异”。我还记得自己第一次给某区域做疾病负担排序分析的时候看到DALY最高的两三个病种确实很清晰但排第四和第五的差距只有零点几个百分点我当时直接按点估计写了顺序被审稿人质疑了。后来重新拉出UI区间一看果然后两个病种的区间几乎完全重叠。做排序分析时尽量把UI区间体现出来这样既显得你懂数据也避免得出过度确定的结论。5. 新手必须知道的坑和排查心得5.1 指标口径最容易搞混我见过太多人在分析里把“发病率”和“患病率”混用两者是完全不同的东西。发病率指的是某时期内新发病例数反映的是风险患病率指的是时点或期间内所有现存病例数反映的是存量。对某些发病率低但生存期长的慢性病来说患病率会远远高于发病率而对病死率高的疾病来说则可能反过来。GBD结果表里这两列是同时存在的如果你的研究问题关注“新发风险”用患病率就错了要描述“现患人群规模”用发病率就不合适。动笔之前先把这个问题想清楚比学会多少统计方法都重要。年龄组选择也是一个常见坑。GBD的结果下载默认会提供非常多的年龄组选项包括什么早新生儿期、晚新生儿期、婴幼儿期、5到9岁、10到14岁等等。有些疾病只在特定年龄段有意义比如某些儿童传染病你只需要看0到14岁数据而有些老年相关疾病集中在中老年组。如果你统一跑一个“All Ages”的合计虽然方便但会掩盖年龄段内部的巨大差异。更好的做法是先用全体年龄的粗率做一个总体判断然后再按年龄段分层看分布这样既能看见宏观轮廓又不会丢失细节。5.2 版本更替和断点问题GBD每隔一段时间会发布新版本每次发布都会把历史年份的数据重新估算。这意味着你在某一年下载的1990年数据和明年下载的1990年数据未必完全一致因为模型、输入数据和死因归类都在持续优化。做纵向研究的时候一定要一次性用完同一个版本的数据并且始终明确标注“基于GBD 2021版数据”。如果把不同版本导出的数表拼在一起可能会在版本边界上看到莫名其妙的断点那不是真实趋势的突变而是估计方法换代导致的假象。关于版本还有一个细节是GBD的疾病列表每版本都会微调某些病种可能会有拆分或合并。如果你的研究跨版本比较数据先确认两次下载使用了完全相同的疾病定义和代码体系否则很容易在分析里引入系统性误差。5.3 模型估计不等于实测数据最后想强调一个认知层面上的坑GBD所有输出本质上都是估算值不是某个地区官方报告的直接统计结果。哪怕某个地区有非常完备的健康数据登记系统GBD也不会直接照搬当地官方死亡统计它需要统一处理去重、错分和口径差异最终输出的是经过模型校正的估计。这意味着当某个区域的GBD数据和那个区域自己的官方统计数据不完全一致时这未必代表GBD错了只能说明两种数据服务于不同的比较目的。做本地化政策分析时比较稳妥的做法是先看GBD的估计值再看当地官方报告两者差异如果在UI区间范围内一般问题不大如果差异很大就要去看看GBD用的数据来源是否包含了该地最新的公开数据有时候只是数据更新时滞造成的。这种情况下写清楚结论的适用范围比强行让两套数字对齐更有价值。5.4 分析过程中的小习惯养成记录版本信息、下载日期和筛选条件的习惯这会让你在写方法学部分时省很多力气。GBD的分析结果必须在方法部分写明使用了哪个版本、哪个指标、年龄标准化与否、UI区间是否提取这些信息缺一项审稿人或阅读者就很难判断结论的可靠性。我自己现在会在每个分析项目目录下放一个README文件专门记录下载时间、版本号、所有筛选条件和字段映射关系看起来很笨但半年后回来重新看项目时你会庆幸自己当时记过这些。6. 技能落地与延伸方向数据下载下来、能做趋势分析只是入门。这套技能的延伸路径其实很长而且每条路径的实用价值都非常高。最常见的一个延伸方向是疾病负担分解分析也就是把DALY变化拆解成人口增长、人口老龄化、流行病学变化三个部分。这个方法可以回答一类特别有价值的问题某个疾病负担上升到底是因为人变多了、人变老了还是疾病本身风险增加了三种原因对应完全不同的干预策略。学术界现在有专门的R包和操作指南在做这个事但核心输入数据就是GBD的结果表学会以后在慢性病、伤害、心理健康研究里都能派上用场。还有就是结合GBD里的SDI社会人口指数分层数据做关联分析。SDI是GBD官方的一个人群分层变量综合了收入水平、受教育程度和生育情况等信息。沿着这个维度看疾病负担能看到很多有意思的结构性差异这在全球健康趋势分析里非常常见。分析这类问题时记得区分“关联”和“因果”SDI分层数据更多反映的是相关关系不是因果关系。往数据工程方向走还可以做自动化的数据更新流程。因为GBD数据每年或每两年更新一次如果项目里需要持续跟踪趋势建议把下载脚本、清洗脚本和绘图脚本封装成流水线新版本一发布改个版本号就能全部重跑一遍。GBD的导出接口逐渐在开放API化但现阶段最常见的方式还是半自动下载再本地处理一个能稳定跑通的脚本能节省大量重复劳动。说到底GBD是一把屠龙刀但它本身并不会告诉你该砍哪里。它的价值取决于你能不能提出好问题、能不能理解估算背后的逻辑、能不能正确处理不确定性。工具会持续更新但“理解数据从哪来、敢于质疑数字、永远把不确定性放在心里”这套方法论是比任何数据版本都更持久的东西。
全球疾病负担研究GBD入门：从DALY到数据应用全解析
我第一次接触GBD全球疾病负担研究数据是在做某个人群健康评估的时候。那时候团队需要回答一个特别现实的问题资源有限的情况下到底该优先干预哪些健康问题翻来覆去找数据最后发现没有任何一份现成的统计年鉴能直接回答直到有人甩给我一份GBD的Excel导出文件说“这是全球疾病负担研究的结果你自己按需拉数吧”。当时我对着那堆指标缩写其实是一头雾水的DALY、YLD、YLL、UI区间、SDI分层每个词都见过但组合到一起就不知道该怎么解读。后来啃了不少方法学资料也亲手跑了不知道多少轮数据提取才慢慢摸清楚这套体系的门道。这篇就当作一个踩过坑的从业者给你们做的GBD入门拆解不讲空话尽量把“这东西是什么、怎么算、怎么用、有什么坑”讲透。1. GBD到底在解决什么问题1.1 GBD不是一本统计年鉴而是一套数据治理框架很多第一次接触的人会把GBD理解成“一个很大的健康数据库”用的时候感觉像查字典点开网页选好病种导出数据完事。这么理解不算全错但它会带来一个很大隐患你并不清楚手里的数字是怎么来的也就很难判断这些数字到底能不能用在自己的研究场景里。GBD本质上是一套全球统一的健康数据测量和估算体系由健康指标与评估研究机构IHME牵头联合全球上千名研究人员持续运作。它要解决的核心矛盾非常简单每个地区都缺数据但决策者不能等数据齐了才做判断。于是GBD采用了一整套多层次的数据整合、交叉校验和模型估算方法把全球每个地区、每个年龄段、两种性别、几百种疾病和伤害的健康损失统一换算成可比的指标。这个思路让我想到做天气预报地面观测站再多也有盲区但你不会因为某个地方缺气象站就拒绝发布天气预报总得靠观测加模型把整片天空拼出来。GBD干的就是这件事。所以项目标题里的“介绍”两个字实际意味着要理解两件事第一这套数据怎么读第二这套数据是怎么被生产出来的。如果只学会点网页下载那你拿到的只是一堆数字而不是研究结果。真正会用GBD的人通常看得懂估算链路知道哪些数字是实测扎实的哪些数字是模型硬推出来的。1.2 谁在用、主要用在哪GBD的数据使用场景比想象中广。做流行病学研究的会拿它来分析某一类疾病在不同人群中的发病和死亡趋势做卫生经济学的会用它给干预措施算成本效益做临床研究的会用疾病负担数据来论证某个研究方向的重要性方便在基金申请里讲清楚立项价值做得更细的还有在做康复、精神卫生、伤害预防这类“在传统死亡统计里看不见”的领域因为单看死亡率会严重低估这些健康问题的影响只有把失能维度也拉进来才说得清楚。我自己感受最深的一点是GBD特别适合用来回答“横向比较”和“纵向趋势”这两类问题。横向比较指的是不同地区之间、或者同一个地区里不同疾病之间的负担排序纵向趋势则是指看十年、二十年的变化轨迹。这两类分析在医院内部报告、科研论文、卫生规划报告里出现频率特别高。你要是能把GBD数据从原始下载到最终成图全流程跑通基本就掌握了一项通用的技能换到多少种疾病、多少层级的地区套路都是同一套。2. 核心逻辑一个负担数字是怎么算出来的2.1 DALY把“早死”和“失能”放到同一杆秤上GBD最核心的产出是DALY伤残调整生命年它的计算是这套体系的中轴。DALY这个概念是1990年代初由世界银行和世界卫生组织的专家们推动成型的后来在GBD项目里得到系统应用和迭代。DALY的设计初衷很简单健康损失不只有死亡这一种表现有些人因病长期卧床有些人带着严重症状活了几十年这些人并没有“死亡”但他们同样损失了大量健康生命年如果指标只衡量死亡就会把这些损失完全漏掉。DALY把两类损失加起来一类是过早死亡损失的生命年YLL另一类是带病生存期间损失的健康年YLD。用一句话可以概括成一个人因某种疾病死亡了算他损失了多少本该拥有的健康寿命一个人没有死但因为疾病处于某种失能状态就需要把他带病生活的年头折算成相当于损失了多少完整健康的年头。两者相加得到的就是这种疾病导致的总负担单位统一为“年”。这个思路相当于把所有健康损失换算成同一种货币然后才能拿不同疾病之间进行直接比较。举个例子可能更容易理解。某种癌症的负担可能主要来自YLL因为确诊后很多人寿命显著缩短而某种慢性腰背痛负担则几乎全部来自YLD因为很少有人直接因为这个疾病死亡但受到影响的人非常多每个人损失的健康年虽不致命乘以庞大的人群基数以后总量非常惊人。2.2 YLL与YLD的估算逻辑YLL的计算需要两个核心输入死亡人数和期望寿命。GBD在做估算的时候针对不同年龄段的死亡都有一套对应的标准期望寿命作为参照一个人越年轻死亡损失的YLL就越大。这其实很公平一个十岁孩子和一个七十岁老人因同种疾病死亡对应的YLL完全不同因为各自剩余的本应健康生存的年数差异很大反映在DALY上就是完全不同的负担量级。YLD部分的计算则比YLL复杂得多因为它不只是看有多少人患病还要看这些患者具体处于什么状态、这种状态带来了多大程度的健康损失。GBD为每一种疾病阶段定义了一套健康状态的失重权重范围在0到1之间0代表完全健康1代表等同于死亡。比如某个疾病导致的轻度症状权重可能只有0.05说明患者损失了5%的健康而终末期疾病权重可能接近0.8意味着那个状态非常接近“完全没有健康生命”。然后用患病率乘以权重再乘以持续时间得到YLD。这个环节是整个GBD估算体系里最容易产生误差的地方因为失能权重本身来自全球多国的调查数据不同人群对“什么是好健康状态”的感受其实存在差异模型已经做了大量标准化处理但不确定性依然存在这也是为什么GBD每个指标都会给出95%不确定性区间。2.3 年龄标准化和两个关键思路除了DALYGBD还提供发病率、患病率、死亡率、期望寿命、健康期望寿命等多种指标不同研究场景侧重点不一样。做病因分析可能更关注发病率和患病率做最终健康结果比较就看DALY和死亡率想看人群整体健康水平则用期望寿命和健康期望寿命HALE。这几个指标各有各的视角单独用哪一个都无法覆盖全部信息组合起来才能讲完整故事。在跨地区和跨时间比较时年龄标准化是绕不开的一个操作。因为不同地区的人口年龄结构可能差异很大一个老龄化严重的地区癌症死亡人数天然会高于年轻人口为主的地区这不是因为它的医疗或环境更差纯粹是年龄结构造成的。GBD官方下载工具里会同时提供粗率和标化率如果你要做区域比较请务必选年龄标准化率否则得出的结论很可能是错的。另外需要提一下的是GBD 2010版本之后官方做了一个重要调整取消了以前采用的年龄权重和时间贴现。早年的DALY计算里中间年龄段健康年的权重会被认为更高未来年份的损失也会被打折折现。后来经过大量争论学界普遍认为这种做法不够透明并且带有明显的价值判断所以现在GBD报告的都是未经过年龄权重和时间贴现的“纯”健康年损失。这个调整让不同年代的结果之间有可比性你在看老文献时如果发现数值和现在下载的不一致很可能就是算法版本导致的不一定是数据错误。3. 数据从哪来GBD喂给模型的“食材清单”3.1 上游数据源普查、登记、监测和文献GBD之所以能覆盖全球几乎所有地区前提是它搜集了大量来源的数据。简单分一下类的话上游数据大概有四种人口普查和调查数据、人口动态登记系统数据、疾病监测系统数据以及已发表的流行病学研究文献。第四种可能最容易被忽略但在一部分地区正式的死亡登记系统覆盖并不完整大量信息来自科学家发表的研究论文GBD团队会做系统检索把那些报告了患病率、发病率、死亡率的研究提取进数据库再参与后续的估算。这带来的结果就是数据的“底层质量”参差不齐。有的地区死亡登记系统非常完善每个死亡个案都有明确死因有的地区可能连人口普查都是多年以前的数据上面的数字还未必可靠。GBD的能力不在于让这些数据变完美而在于用统计手段把不完整、不一致的数据整合到一起并准确告诉使用者“我们的估计有多大的不确定性”。所以我在看GBD每次更新的版本说明时都会重点观察某个地区或疾病的数据来源数量有没有增加数据来源多了估计就会更扎实UI区间也会明显变窄。3.2 中游标准化死因推断与ICD编码映射原始数据进来以后首先要做标准化处理这里面有两个特别关键的步骤死因推断和编码映射。死因推断针对的是没有医学诊断死因的地区。很多死亡发生在医疗体系覆盖不到的地方家人只能大概描述死者生前的症状这时候需要通过死因推断方法把描述性文本对应到可能的死因类别。GBD设计了专门的算法来做这件事很大程度上提高了死因数据的覆盖范围但推断本身带有概率性所以对应的死亡负担数据通常不确定性更大。编码映射则是把不同来源、不同年份、不同编码系统下的诊断和死因名称统一到GBD自身的一套疾病层级目录里。国际疾病分类ICD从第9版、第10版到第11版不同地区的使用习惯和扩展码都不一样某些旧数据用的还是非常粗略的分类名称GBD团队必须一个个做映射。这个环节听着琐碎实际上是最影响数据质量的环节之一如果映射错了后续模型吃得再准也白搭。3.3 下游模型DisMod-MR 2.1和CODEm标准化之后的数据不会直接变成输出指标中间还要经过估算模型的处理其中两个模型一定要知道DisMod-MR 2.1和CODEm。DisMod-MR 2.1是个贝叶斯元回归工具主要用来估计非致命健康结果也就是发病率、患病率、持续时间、缓解率这些参数。它厉害的地方在于可以做多层次的跨地区信息借用数据充足的地区更多依赖本地数据数据缺乏的地区则更多借用周边或全球层面的信息然后综合得出估计值。这套机制的假设是“同一类疾病的流行病学特征在全球范围内存在一定的相似性”虽然不完全对但在缺数据场景下已经是目前最可行的办法。CODEm则是处理死因数据的模型专门用来估计各类疾病和伤害的死亡人数。它会同时考虑多种与死亡相关的因素在不同候选模型里做交叉验证选出预测表现最好的组合再综合多个模型的结果。用大白话说就是在“根据已知死因数据往回推未知死因”的过程中让多套模型互相投票而不是赌某一个模型的运气。这也是GBD数据中死亡估计相对比较稳定的原因之一。4. 实操从GBD官网下载数据到完成第一张趋势分析表4.1 获取数据的主要途径搞懂指标体系之后真正的动手环节就来了。现在获取GBD数据主要有几种方式我用下来觉得各有优缺点。最常规的是GBD Results Tool也就是IHME官网的在线下载工具适合做按需取数。选好地区范围、年龄组、性别、年份、指标和疾病原因它会生成一张交叉表可以导出CSV格式的数据。这个工具适合不需要太复杂变量的常规分析我自己大多数时候都用它。如果下载需求比较批量或者需要做自动化分析更推荐直接从GHDx全球健康数据交换平台下载完整结果数据集里面的文件通常很大但字段更全适合用Python或R进行灵活拆分。还有一个工具是GBD Compare它本质是个可视化交互界面适合快速看图找方向不适合作为研究中提取最终数据的来源。做严肃分析我还是建议用Results Tool或GHDx因为导出的数据才是结构化、带完整字段的。4.2 下载界面的核心选项怎么选以Results Tool为例进去以后要设置的那套条件看起来很多但只要想清楚分析目标就不难。Years的选项直接决定你能做多长时间跨度的趋势分析官方提供的数据覆盖从1990年到最近发布年份通常会选择完整区间至少得有十年才能看出趋势变化。Location选择分析的地理单元这个体系是多层级的全球、区域层级、国家、国家以下行政区逐级向下。你一定要先想清楚自己的研究问题聚焦在哪一层比如想做全球层面的整体比较选Global或者区域层级就够了想做某个特定国家或地区内部的差异就得选更细的行政区划。每下拉一级背后可用的数据和模型稳定性都会不同这很正常。Cause选择疾病原因这是GBD的疾病层级目录分为三级。Level 1是最粗的分类比如传染病、孕产妇、新生儿和营养性疾病非传染性疾病伤害三大类往下细分到Level 2比如肿瘤、心血管疾病、慢性呼吸系统疾病再到Level 3就是具体的单个病种比如肺癌、缺血性心脏病。建议选Level 2或Level 3作为分析粒度这样既有足够的区分度又不会碎片化到难以解读。Metrics选项里有Number例子数、Percent构成比、Rate每10万人率三种维度实际研究中Number适合看总量Rate则适合比较不同人口规模的地区。我一般两个都会下载一个看结构一个看强度。4.3 下载后的数据结构与预处理要点导出的CSV打开以后你会看到有很多行列反复重复的组合字段年份、地区、性别、年龄组、疾病、指标呈长表格式排列。首先要做的就是了解字段含义然后根据需要做数据透视或者宽表转换。在做数据清洗的时候有几个小点容易被忽视。第一是年龄组编码“All Ages”和各个具体年龄组是分开的行如果你自己在做年龄标准化只按“每10万人率”里的年龄标准化率这一列取数就好不要自己用粗率再去套什么权重。第二是性别字段“Both”这种合计行和Male、Female分开存在做合计分析时记得选Both别把男女行相加再除以二那完全是多此一举。第三就是UI区间上下界字段官方会给出低值和高值做成图时可以考虑画误差线这样审稿人会有更多安全感。实际分析过程中我最常用的操作是从长表里筛选出自己关心的疾病子集和年份范围然后计算DALY构成比也就是每种疾病占该地区总负担的比例再排一个序。这个表做出来以后你通常能一眼看出来该地区最主要的几个健康问题是什么再叠加一个分年份的趋势折线就能看到哪些疾病负担在上升、哪些在下降整个逻辑就闭环了。4.4 要让结果可信必须学会看UI区间这可能是新手最容易忽略的地方。GBD是一个模型估算系统所以每个输出值都带有一个95%不确定性区间UI这个区间的宽度直接反映了一个估计值的可信程度。如果两个疾病的点估计排第一第二但它们的UI区间完全重叠那你就不能很自信地说谁负担更重正确的说法应该描述为“两者负担相近实际上尚无证据表明存在显著排序差异”。我还记得自己第一次给某区域做疾病负担排序分析的时候看到DALY最高的两三个病种确实很清晰但排第四和第五的差距只有零点几个百分点我当时直接按点估计写了顺序被审稿人质疑了。后来重新拉出UI区间一看果然后两个病种的区间几乎完全重叠。做排序分析时尽量把UI区间体现出来这样既显得你懂数据也避免得出过度确定的结论。5. 新手必须知道的坑和排查心得5.1 指标口径最容易搞混我见过太多人在分析里把“发病率”和“患病率”混用两者是完全不同的东西。发病率指的是某时期内新发病例数反映的是风险患病率指的是时点或期间内所有现存病例数反映的是存量。对某些发病率低但生存期长的慢性病来说患病率会远远高于发病率而对病死率高的疾病来说则可能反过来。GBD结果表里这两列是同时存在的如果你的研究问题关注“新发风险”用患病率就错了要描述“现患人群规模”用发病率就不合适。动笔之前先把这个问题想清楚比学会多少统计方法都重要。年龄组选择也是一个常见坑。GBD的结果下载默认会提供非常多的年龄组选项包括什么早新生儿期、晚新生儿期、婴幼儿期、5到9岁、10到14岁等等。有些疾病只在特定年龄段有意义比如某些儿童传染病你只需要看0到14岁数据而有些老年相关疾病集中在中老年组。如果你统一跑一个“All Ages”的合计虽然方便但会掩盖年龄段内部的巨大差异。更好的做法是先用全体年龄的粗率做一个总体判断然后再按年龄段分层看分布这样既能看见宏观轮廓又不会丢失细节。5.2 版本更替和断点问题GBD每隔一段时间会发布新版本每次发布都会把历史年份的数据重新估算。这意味着你在某一年下载的1990年数据和明年下载的1990年数据未必完全一致因为模型、输入数据和死因归类都在持续优化。做纵向研究的时候一定要一次性用完同一个版本的数据并且始终明确标注“基于GBD 2021版数据”。如果把不同版本导出的数表拼在一起可能会在版本边界上看到莫名其妙的断点那不是真实趋势的突变而是估计方法换代导致的假象。关于版本还有一个细节是GBD的疾病列表每版本都会微调某些病种可能会有拆分或合并。如果你的研究跨版本比较数据先确认两次下载使用了完全相同的疾病定义和代码体系否则很容易在分析里引入系统性误差。5.3 模型估计不等于实测数据最后想强调一个认知层面上的坑GBD所有输出本质上都是估算值不是某个地区官方报告的直接统计结果。哪怕某个地区有非常完备的健康数据登记系统GBD也不会直接照搬当地官方死亡统计它需要统一处理去重、错分和口径差异最终输出的是经过模型校正的估计。这意味着当某个区域的GBD数据和那个区域自己的官方统计数据不完全一致时这未必代表GBD错了只能说明两种数据服务于不同的比较目的。做本地化政策分析时比较稳妥的做法是先看GBD的估计值再看当地官方报告两者差异如果在UI区间范围内一般问题不大如果差异很大就要去看看GBD用的数据来源是否包含了该地最新的公开数据有时候只是数据更新时滞造成的。这种情况下写清楚结论的适用范围比强行让两套数字对齐更有价值。5.4 分析过程中的小习惯养成记录版本信息、下载日期和筛选条件的习惯这会让你在写方法学部分时省很多力气。GBD的分析结果必须在方法部分写明使用了哪个版本、哪个指标、年龄标准化与否、UI区间是否提取这些信息缺一项审稿人或阅读者就很难判断结论的可靠性。我自己现在会在每个分析项目目录下放一个README文件专门记录下载时间、版本号、所有筛选条件和字段映射关系看起来很笨但半年后回来重新看项目时你会庆幸自己当时记过这些。6. 技能落地与延伸方向数据下载下来、能做趋势分析只是入门。这套技能的延伸路径其实很长而且每条路径的实用价值都非常高。最常见的一个延伸方向是疾病负担分解分析也就是把DALY变化拆解成人口增长、人口老龄化、流行病学变化三个部分。这个方法可以回答一类特别有价值的问题某个疾病负担上升到底是因为人变多了、人变老了还是疾病本身风险增加了三种原因对应完全不同的干预策略。学术界现在有专门的R包和操作指南在做这个事但核心输入数据就是GBD的结果表学会以后在慢性病、伤害、心理健康研究里都能派上用场。还有就是结合GBD里的SDI社会人口指数分层数据做关联分析。SDI是GBD官方的一个人群分层变量综合了收入水平、受教育程度和生育情况等信息。沿着这个维度看疾病负担能看到很多有意思的结构性差异这在全球健康趋势分析里非常常见。分析这类问题时记得区分“关联”和“因果”SDI分层数据更多反映的是相关关系不是因果关系。往数据工程方向走还可以做自动化的数据更新流程。因为GBD数据每年或每两年更新一次如果项目里需要持续跟踪趋势建议把下载脚本、清洗脚本和绘图脚本封装成流水线新版本一发布改个版本号就能全部重跑一遍。GBD的导出接口逐渐在开放API化但现阶段最常见的方式还是半自动下载再本地处理一个能稳定跑通的脚本能节省大量重复劳动。说到底GBD是一把屠龙刀但它本身并不会告诉你该砍哪里。它的价值取决于你能不能提出好问题、能不能理解估算背后的逻辑、能不能正确处理不确定性。工具会持续更新但“理解数据从哪来、敢于质疑数字、永远把不确定性放在心里”这套方法论是比任何数据版本都更持久的东西。