背景

Mixture analysis是近几年环境流行病学最热门的方法,即“混合物分析”,其中BKMR是非常非常重要的方法，无法忽略它强大的解释能力。小编上一期推文混合物分析新思路：BKMR套路解析中详细介绍过。陆陆续续有小伙伴咨询，大样本(>1000人)跑的时候太慢了，有的人甚至说要跑5天。那有没有更快速跑BKMR的解决方案呢？也许本章节就有答案。

bkmrhat包解析

bkmr包有4个缺点，1)没有并行同时运行多个MC链2）没有跨多个MC链的推断3）参数的有限的后验总结4）有限的诊断。其实，bkmrhat就是用future包的多线程写入原来的包，达到经典bkmr包速度的4-24倍(根据CPU核心线程不同不同)，甚至更快。Bkmrhat包地址 https://github.com/kdevick/bkmr-cma。主要函数如下，相当于bkmr增加了paralled功能。

R实战

案例与BKMR案例相同，dataset1数据集(文末获取附件)。结果解释参考BKMR章节，一样。

Tips

目前未见比较好的文章有详细运用BKMRhat包各项功能的案例，有见到的小伙伴还请文章底部留言告知一下。谢谢。

小结

BKMRhat包提供了相对BKMR包高达20倍以上速度提升，可计算总效应趋势(多污染物模式)、相对重要性(PIP)；也可分MC链展示分链多污染物模式、分链单污染物模式。

数据和代码索要请公众号后台发送信息 bkmrhat 

本公众号建立了学术交流群(群)，仅供SCI学术交流，人数有限需要实名制。入群请加笔者微信 popnie，请备注说明：姓名-学校(单位)-专业。

参考文献：1https://github.com/alexpkeil1/bkmrhat 2 https://cran.r-project.org/web/packages/bkmrhat/vignettes/bkmrhat-vignette.html

本文使用 文章同步助手 同步