昨日对文章Wu C, Yan Y, Chen X, Gong J, Guo Y, Zhao Y, Yang N, Dai J, Zhang F, Xiang H. Short-term exposure to ambient air pollution and type 2 diabetes mortality: A population-based time series study. Environ Pollut. 2021 Nov 15;289:117886. doi: 10.1016/j.envpol.2021.117886. Epub 2021 Jul 31. PMID: 34371265
进行解读,然后绘制效应曲线后,文章的作者联系了我
我借机问作者要了原始的绘图代码,作者也很大方,得到作者的代码后,我比较了一下,和我的方法一基本相同,我用作者代码按我的数据跑了一下
这是作者代码出的图:
下面这个是我的代码出的图:
感觉我的代码曲线更加顺滑一点,因为他是用R自带的函数取的系数,我是用mgcv包自带函数提取系数。
所以大家就放心的用吧!
文章《R语言广义可加模型在环境污染方面的应用 (2)》之后续----文章作者联系我啦
![](https://img-home.csdnimg.cn/images/20240711042549.png)