本篇是GBD相关的第13篇推文
接着上次APC模型的介绍,我们介绍了APC网页版模型的使用,但网页版其实的核心是用R语言来编写的,于是我尝试联系APC网页版模型的通讯作者,很高兴通讯作者回复了,并将原代码分享给了我。我根据GBD研究的实际情况进行了微调,这里分享给大家。
首先,我们需要载入分析需要的R包
后面大部分代码和前面的第12篇推文类似,获取和网页版差不多的格式
获取人口学数据
获取发病数据
两者年龄层取交集后整理成我们需要的数据格式
接着,我们用prepare_rates包将incidence_population转变成list格式,用于APC模型
接着我们运行APC模型,代码很简洁,只要一行就可以
后面作图也只要1行代码即可搞定
最后一张图上可出现几乎所有的结果图
结果解读基本上在上一个推文中已经讲过了,这里就不赘述了。