r - 如何更改 R 代码以将 cox 比例风险模型与 BDLIM 相结合？

Question

我读过论文“贝叶斯分布式滞后交互模型，以确定儿童健康中的围产期脆弱性窗口”。我还尝试了 Github 网站上的 R 代码。我还阅读了另一篇论文“Prenatal fineparticle exposure and early child孢子气对母体压力和胎儿性别的影响”，作者还应用 BDLIMs 来调用 OR 的效果修改。我一直在寻找如何将 BDLIM 中的 y 轴转换为 OR 的方法，但一直没有找到。此外，我想知道如何将 y 轴的标签更改为 OR 并将 x 轴的标签更改为孕周。但是，我想使用 cox 比例风险模型来计算我的数据集中的 HR。同时，我们还想应用 BDLIMs 来检查性别效应修饰。所以，我可以问你如何更改 R 代码以将 cox 比例风险模型与 BDLIM 结合起来吗？我的数据：y 轴是一种疾病（有和没有自闭症 (1,0)） x 轴是每周大气汞暴露。我想通过带有 BDLIM 的 cox 比例风险模型来研究暴露与疾病之间的关联。

R code:
fit <- bdlim(Y=data$autism,X=lag, Z=data$ID_SEX, G=data$ID_SEX, 
             inter.model="all", niter=100, nthin=5,
             basis.opts = list(df=7, type="bs"),seed=1234)

我希望y轴是HR，但实际输出是估计效果

标签： r