基于贝叶斯模型平均的基准剂量估计及其在白银市人群镉暴露评估中的应用

[背景]最优模型法估计基准剂量(BMD)时未考虑模型选择的不确定性,国内暂无使用贝叶斯模型平均(BMA)估计BMD的研究.[目的]将BMA应用于我国镉污染的暴露评估研究,讨论BMA在基于剂量-反应模型的BMD估计中的作用,为有害物质健康风险评估提供方法学支持.[方法]模拟研究基于甘肃省白银市某镉污染区数据获得的不同剂量-反应模型(Gamma、Log-logistic、Log-probit、Two-stage和Weibull模型)的参数和中国5个镉污染区的尿镉范围,模拟不同正确模型、不同剂量组数(5、8)、不同样本量(50、100、200)的模拟数据,比较BMA和最优模型法的优劣.实例分析使用甘肃省白银市某镉污染区的镉暴露数据.模拟研究和实例分析均以尿镉为镉暴露指标,β2-微球蛋白异常率为效应指标,基准反应设定为10％,分别使用正确模型(模拟数据时使用的模型)、最优模型法[赤池信息准则(AIC)最小的模型]和BMA估计BMD及其95％置信区间下限(BMDL);比较不同方法估计的BMD、BMDL及相对偏差.[结果]对于模拟研究,随着样本量或剂量组数增加,BMD的第5百分位数和第90百分位数区间大致呈现变窄趋势;当正确模型为单一模型时,BMA估计的BMD相对偏差大于最优模型法;当正确模型为等权重混合模型时,BMA估计的BMD相对偏差小于最优模型法.基于实例数据,最优模型为Log-probit模型(AIC=1814.46),其次为Log-logistic模型(AIC=1814.57);Log-probit模型、Log-logistic模型和BMA估计的BMD(BMDL)分别为3.46(2.68)、3.16(2.33)和2.92(2.07)μg·g?1.[结论]在正确模型已知时,最优模型法依然是值得推荐的方法;但在有害物质剂量-反应关系高度不确定或来源人群、暴露分组等不同时,BMA与最优模型法相比,理论上更充分地考虑了可能存在的多个备选模型,能提供更为稳定的BMD和BMDL估计.