适用于个体配对资料的加权分位数合回归WQS?
WQS是用来分析混合物效应的一种方法,但是已有的gWQS包不能用于分析个体配对资料(时间分层病例交叉设计),有文献建立了一种方法可以解决这个问题,我下载了文献和作者写的代码,但是不知如何带入自己的数据和变量。方法比较新,目前百度上没有教程。
图中是部分代码的截图
我想知道应该如何代入自己的数据和变量,代码什么位置需要更改。
可以接受500以下的付费咨询
参考文献:Lee M, Rahbar MH, Samms-Vaughan M, et al. A generalized weighted quantile sum approach for analyzing correlated data in the presence of interactions. Biometrical Journal. 2019;61:934–954. https://doi.org/10.1002/bimj.201800259
这段代码看起来是在进行加权分位数回归分析。代码需要更改的地方是数据和变量的部分。可能需要修改的地方有:
data.Nt[ind,],这里应该改成你自己的数据,ind是索引
data.Nt.bs$Case,这里应该改成你自己的病例变量名
data.Nt.bs$q.meta11, data.Nt.bs$q.meta12等,这里应该改成你自己的变量名
data.Nt.bs$GENOTYPE.2, data.Nt.bs$GENOTYPE.3, 这里应该改成你自己的基因型变量名
在更改变量名之后,还需要确保这些变量与数据的结构是匹配的。