基于g-computation联合混合效应模型控制未测混杂因素的因果推断方法模拟研究及实例验证

目的通过模拟实验和实例验证探讨基于g-computation的联合混合效应模型(joint mixed-effects model,JMM)控制纵向研究未测混杂因素进行因果推断时的效果及性能特点。方法通过计算机模拟产生包含基线及两次随访时点的纵向数据,模拟条件包括样本含量、有无未测混杂因素及未测混杂效应大小,分别利用基于g-computation的JMM、线性混合效应模型、固定效应模型和纵向目标极大似然估计方法估计因果效应,通过平均绝对偏差(mean absolute deviation,MAD)、标准误、均方根误差(root mean square error,RMSE)、95%置信区间覆盖率(95%confidence interval coverage,95%CI coverage)评价比较各方法因果推断的效果。利用绝经期女性队列体检数据,应用四类模型分别估计绝经期女性血清卵泡刺激素(follicle-stimulating hormone,FSH)水平与腰椎骨密度间因果关系,对各模型在真实纵向数据中的因果推断效果进行验证。结果JMM控制未测混杂因素的因果推断准确性最佳,但稳定性略差。当研究中存在较强未测混杂效应时,仅JMM可准确估计因果效应,且其在大样本量时估计的精确性和真实性较好。结论基于g-computation的JMM可有效控制纵向研究中未测混杂因素进行近似无偏因果推断。

前列腺癌患者经放射治疗后发生第二原发膀胱癌的风险分析

目的探究经放射治疗后前列腺癌患者发生第二原发膀胱癌(SPBC)的风险并评估结果对未测得混杂的稳健性。方法选取美国“监测、流行病学和结果(SEER)”数据库中2004年至2009年病理诊断明确的前列腺癌患者资料,根据是否接受放射治疗分为放疗组和未放疗组,应用倾向性评分匹配(PSM)对两组患者进行匹配;使用竞争风险模型评估放疗后前列腺癌患者SPBC的发生风险,并采用未测得混杂敏感性分析,评估结果和结论的稳健性。结果本研究共纳入83064名前列腺癌根治术患者,其中接受放疗的为5141例(88例发生SPBC),未接受放疗的为77923例(753例发生SPBC)。匹配前后竞争风险模型显示,放疗增加前列腺癌患者发生SPBC的风险(匹配前SHR为1.438,95%置信区间1.133~1.824,P=0.003;匹配后SHR为1.586,95%置信区间1.076~2.215,P=0.007)。敏感性分析结果显示,未测得混杂与放疗和SPBC发生的关联为OR=1.65,SHR=1.84可以使观察到的放疗-SPBC发生风险关联无统计学意义。结论对放疗导致前列腺癌患者SPBC发生风险增高应谨慎判断,某个与结局或放疗关联较低的混杂因素就可能将结果推至无效。