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

时间:2016-08-12 15:38:52

标签: sas rates irr gee

我想为集群数据运行GEE - 我试图获得两种药物之间抗生素反应的发生率比率(IRR)。我已经搜索了有关构建GEE模型的信息(SAS中的GENMOD,Stata中的xtgee),但我无法找到可以作为协变量包含哪些类型的变量的标准。我的模型是这样的:

 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是所讨论的抗生素的不良事件(AE)率,mero_rate是不同抗生素的AE率。其他变量要么是分类的,要么是连续的。

如果我用一个速率的协变量调整GEE,是1)正确使用GEE模型,2)exp(coef)的解释是两个AE速率之间的IRR,或者它被解释为:对于每个单位增加的mero_rate率,rate_pip的IRR是高出/低于x倍?

谢谢 - 如果需要,我很乐意澄清/提供更多信息。

1 个答案:

答案 0 :(得分:0)

我不知道这是否正确使用GEE模型而又不了解更多关于数据结构的内容,但我不知道任何有关GEE模型的特殊情况会阻止使用速率变量预测因子(与普通最小二乘回归模型相比)。如果模型在没有mero_rate预测器的情况下没问题,那么它也可能没问题。也许需要注意的是,它与其他预测因素不太相关。

就解释而言,我认为你已经得到了它。对于每单位时间的x + 1事件的mero_rate值,发病率的对数增加β单位,与每单位时间的x事件的mero_rate值相比,所有其他条件相等。