基准剂量估计R软件包bmd的改进与应用

摘要

目的 本研究旨在构建一个更加便捷的R扩展包简化基准剂量估计(benchmark dose, BMD)的计算和操作过程。方法 使用R 4.1.3对原有R软件包bmd进行修改和整合，构建一个R扩展包。采用研究中封装完成的R扩展包将微囊藻毒素LR(microcystin-LR,MCLR)的数据与18种剂量-反应模型拟合，根据Akaike信息标准(Akaike information criterion, AIC)最小原则，得到最佳拟合模型下MCLR的BMD和基准剂量下限(benchmark dose lower confidence limit, BMDL)。结果 封装后的R扩展包将原有bmd包的运算过程进行简化。MCLR的基准剂量估计中2个参数的Weibull type 2模型AIC值最小，在该模型下，MCLR的BMD和BMDL分别为1.39 mg/kg和0.49 mg/kg。结论 封装后的R扩展包简化了BMD的计算和操作过程，可以同时得到18种剂量-反应模型的AIC、BMD和BMDL,直观选择AIC最小的模型相对应的BMD及BMDL。

关键词

基准剂量估计 / 剂量-反应 / 微囊藻毒素