简体   繁体   English

针对费率运行GEE模型,并使用也是费率的协变量进行调整(GENMOD SAS)

[英]Running a GEE model for a rate, and adjusting with a covariate that is also a rate (GENMOD SAS)

I want to run a GEE for clustered data - I am trying to get incidence rate ratios (IRR) for antibiotic reactions between two drugs. 我想为聚类数据运行GEE-我试图获取两种药物之间抗生素反应的发生率比(IRR)。 I have searched for information on constructing GEE models (GENMOD in SAS, xtgee in Stata) but I can't find criteria on what type of variables can be included as covariates. 我已经搜索了有关构造GEE模型的信息(SAS中的GENMOD,Stata中的xtgee),但是我找不到关于可以将哪些类型的变量包括为协变量的标准。 My model is this: 我的模型是这样的:

 proc genmod data = mydata;
  class Pt fev1_cat;
  model rate_pip = cumulative_dose_before fev1_cat Average_Dose_Admis mero_rate  / 
  type3 dist=poisson link=log;
  repeated subject=Pt;
 run;

rate_pip is the rate of adverse events (AE) for antibiotic in question, mero_rate is the rate of AE for a different antibiotic. rate_pip是所讨论抗生素的不良事件(AE)发生率,mero_rate是其他抗生素的AE发生率。 The other variables are either categorical or continuous. 其他变量是分类变量或连续变量。

if I adjust the GEE with a covariate that is a rate, is it 1) a correct use of the GEE model, and 2) would the interpretation of the exp(coef) be the IRR between the two rates of AE, or is it interpreted as: for each unit increase in rate of mero_rate, the IRR of rate_pip is x times higher/lower? 如果我用比率作为协变量来调整GEE,是1)正确使用GEE模型,2)exp(coef)的解释是两个AE比率之间的IRR,还是解释为:对于mero_rate速率的每个单位增加,rate_pip的IRR是高/低x倍?

Thanks - I'd be happy to clarify/provide more info if needed. 谢谢-如果需要,我很乐意澄清/提供更多信息。

I can't say whether this is a correct use of a GEE model without knowing a little more about the data structure, but I don't know of anything that's special about GEE models that would preclude the use of rate variables as predictors (as compared to say a ordinary least squares regression model). 我不能在不了解更多数据结构的情况下说这是否正确使用了GEE模型,但是我不知道GEE模型有什么特别之处,它会排除使用速率变量作为预测变量(如比起普通的最小二乘回归模型而言)。 If the model is okay without the mero_rate predictor, it would probably be okay with it too. 如果没有mero_rate预测变量,该模型还可以,那么它也可能也可以。 Maybe the caveat is that it can't be too correlated with the other predictors. 可能需要注意的是,它不能与其他预测变量过于相关。

As far as interpretation goes, I think you've pretty much got it. 就解释而言,我认为您已经理解了。 The log of the incidence rate increases by beta units for a mero_rate value of x+1 events per unit time, compared to a mero_rate value of x events per unit time, all other things equal. 在所有其他条件相同的情况下,与每单位时间x个事件的mero_rate值相比,每单位时间x + 1事件的mero_rate值的发生率的对数增加了beta单位。

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM