事件发生率:比较频次前累计风险暴露时间
两份数据可以有相同企业数和相同事件数,却代表完全不同的观察时长。跟踪三个月与跟踪一年贡献的证据不同。比较频次之前,应先逐对象建立暴露账,说明它何时进入风险集、何时停止贡献时间。
交付是事件次数除以累计风险暴露时间,并附清楚单位。这个发生率不能自动解释为一年内发生概率。首次事件与重复事件还需要不同计数规则,若只换分子、不改相应暴露分母,就可能制造错误比较。
事件与计时规则同时定义
假设一个研究登记库跟踪企业入组后的首次经核实重大报告更正。合格条件要求研究规定历史内尚无符合标准的更正。事件日期采用预先确定的生效日期,不使用研究者碰巧下载通知的时间。
单位为企业日,每家从入组起最多跟踪三百六十五天;首次事件或更早失去可观察覆盖时停止相应贡献。每行保存真正观察开始与终止点,不能因为大家采用同一最长规则,就把实际时间都填为一年。不同企业的日历入组日期可以不同。
CDC归档的风险测量课程在发生率部分解释按指定终点累计对象时间。本文企业日账是对这一分母逻辑的自拟非临床应用,不代表CDC研究或认证这些公司数据,也不把企业日与人数分母混用。
先逐行建立暴露账
假设四家合格企业的记录如下。时间均为入组后经过天数,事件在该段终点计入。
| 企业 | 开始 | 终止 | 终点原因 | 首次事件 | 风险日 |
|---|---|---|---|---|---|
| A | 0 | 365 | 研究结束 | 0 | 365 |
| B | 0 | 365 | 首次事件 | 1 | 365 |
| C | 0 | 180 | 首次事件 | 1 | 180 |
| D | 0 | 90 | 失去覆盖 | 0 | 90 |
累计为365+365+180+90,等于1000企业日,事件两次。观察发生率为2/1000,即每天每企业0.002次,也可写每一千企业日两次。
C在第一百八十天首次事件后不再贡献风险时间;D在第九十天失去覆盖后不再贡献,因为后来事件状态未知,而不是已经发生事件。补上它们剩余的一百八十五天及二百七十五天,会加入不合格或未观察暴露。2/4是观察事件比例,D在三百六十五天的状态未解决,不能直接确定完整队列的统一期间概率。
本例A确认至三百六十五天无事件,B则恰在该日发生合格事件。如果所有无事件对象都完整观察至同一期间,事件对象虽然较早停表,事件比例仍可描述该统一期间;不能笼统说事件停表就使固定期比例无效。这里问题来自D在九十天右删失。应保留终点原因及同日事件优先规则,不能只留下一个停止日期。
每组用自己的暴露分母
另设可比登记库E,同一事件定义下,在五千企业日观察到五次事件,发生率为每一千企业日一次。原登记库相对E的率比为2,率差为每一千企业日一次。
如果只比两次与五次,会得到与发生频率相反的印象,暴露分母解释了差别。率比在两边采用同一时间单位时可以抵消单位,展示的原发生率和率差仍必须带单位,不能把每年与每天直接相除。
Cochrane手册第六点七节区分次数、率与二元结果,并要求事件总数和风险时间对应。虽然分母不是企业数,仍应保留企业数量用于描述研究:一千天集中在极少企业,与分散在很多企业,依赖结构可能完全不同。
从率换概率必须另加模型
0.002乘365得到每企业年0.73次,只是改了时间单位,不是年概率百分之七十三。率可以超过每单位一次,概率却不能超过一。
若另作严格示意,假设首次事件的日hazard恒定为0.002,没有竞争性退出,则三百六十五天概率为1−exp(−0.002乘365),约0.5181,与0.73不同。指数换算加入了关于时间结构及事件过程的假设,不能说四企业两事件的账本已经验证了这种恒定机制。
若报告季或企业年龄影响风险,单一总体率会丢失时序信息,需要适当生存分析或可论证的时变模型。收购导致原定义事件不再可发生,也不能与数据商临时失联等同;应识别竞争终态,而非把所有停止都当同一种无害删失。
重复事件必须规定重新入集
换一个独立目标:可在恢复后再次发生的临时运营故障。某企业E在第九十和第一百八十天发生故障,分别到第一百天和第一百九十天恢复,研究在第三百六十五天结束。按假设规则,恢复前没有资格发生下一次新故障。
风险区间为(0,90]、(100,180]、(190,365],长度90+80+175=345天。重复事件两次,率为2/345,约每天每企业0.005797次。两段十天恢复期按既定资格规则排除,不是为了调整结论临时减掉。
同一历史若采用首次事件目标,只用一次事件和九十天,得到该单企业记录的1/90。两种分子不能互相粘贴到对方分母。若实际过程允许故障后立即再次入集,则暴露规则又不同,需要由研究定义决定。
同企业重复事件通常存在依赖,拥有时间分母不等于独立泊松模型已经成立。计算区间前,应检查复发、集中程度及聚类,更多事件行没有创造更多独立企业。永久吸收性的倒闭也不能照搬“恢复后重新入集”规则。
处理日期不完整和信息性观察
年度复核若只知道某一年内发生更正,确切日期仍未知。按年中分配是一项假设,不是恢复出的时间戳。应保留已知区间,核对不同合理日期对暴露合计的影响;精确事件时序重要时,应使用适合区间信息的方法。
失去观察还可能与即将出问题有关。如果困难公司总在更正前消失,正确累计已知天数,也不能保证总体率无偏。需要追查覆盖中断机制、保留原因,并比较组间跟踪完整性;准确分母不能自行补回未见事件。
| 步骤 | 输入与计算 | 缺失冲突的处理 |
|---|---|---|
| 固定资格 | 事件、首次或重复目标 | 不合并不相容事件 |
| 建立时钟 | 原点、进入、终止、单位 | 未知日期仍列未知 |
| 累计暴露 | 每企业不重叠合格区间 | 清理重复及重叠时间 |
| 对齐事件 | 合格日期与企业标识 | 合并指向同一事实的重复通知 |
| 比较组别 | 相同定义和率单位 | 覆盖变化与频次分开 |
| 评估前提 | 复发、删失、竞争退出 | 限定解释或修订模型 |
两个来源日期不同时,先区分法律生效日、公告日与研究发现日,再选择研究时钟。应保存采用规则与备选日期,静默替换可能同时改变事件顺序及暴露时间。
新增漏报事件被核实、终点日期纠正、企业被发现不合格或复发规则修订,都会触发重算。事件数量桥和时间桥应分别能复现,才不会把延长覆盖或更换风险集造成的变化,误归因于企业真实行为。
时间暴露也要去重
两个提供方可能同时覆盖一家企业。另设没有资格中断的观察账:一方提供零至一百八十天,另一方提供九十至二百七十天。直接相加一百八十加一百八十,会报三百六十企业日;区间并集实际只有二百七十天,重叠九十天。这是独立暴露例,不是在主例首次事件后继续补时间。
应有明确覆盖合并规则,并保留各段来源。同一事件也可能在两个数据流重复出现,需要唯一事实标识。只去事件重、不去时间重,会让率偏小;只去时间重、两条通知却都计事件,会让率偏大。数量桥和暴露桥必须一起核对。
可将每家企业的合格区间排序,检查重叠、间隙及对应事件,最后再汇总。无法确认间隙期间是否仍可观察时,不能用最早开始到最晚结束一减了之,把中间缺口自动填满。覆盖是否连续属于证据问题,不是日期格式问题。
同时固定时间换算口径。研究按实际天数累计,与把每个日历年简单换成三百六十五天,并非在所有日期下都相同;可保留原始日期和计算日数,展示单位另作转换。不能只因两个表都写企业年,就默认它们累计相同暴露。
分组率比较还应保存每组暴露构成。如果新增覆盖主要来自原本高风险时期,率变化可能同时反映观察组合变化。至少列出主要时期、进入批次及对象构成,让读者知道比较是同一目标的频率变化,还是采集范围已改变。
常见问题
每一千企业日两次是否就是年度概率?
不是。它是含时间分母的率,固定期间概率还需要事件过程信息或额外模型。
重复故障能否直接加入首次事件分子?
不能。重复事件需要相应合格暴露及依赖处理,不能只增加次数而沿用旧分母。