Risk Ratio proc genmod predictor变量有多个级别

Question

我正在使用proc genmod来估计风险比率。我的两个预测变量有两个以上的水平。我的gender变量是‘male’, ‘female’, ‘other’，race是‘white’，’non-hispanic black’，‘hispanic’。下面是我建立模型的方法。我只得到1个风险比。是否有可能为预测变量的每个级别确定风险比？我以前没有这样做过，所以任何建议都是值得感谢的。

谢谢。

Proc genmod data=rr_genmod;
Class gender(ref=’male’) race (ref=’white’);
Model outcomeA(ref=’1’)=gender race/ dist=binomial link=log;
ESTIMATE 'BETA' GENDER 1-1/EXP;
ESTIMATE 'BETA' RACE 1-1/EXP;
RUN;

Answer 1

我会从LSMEANS开始，然后可能从LSMESTIMATE开始进行更复杂的估计。语法应该是这样的。

LSMEANS GENDER RACE / DIFF CL EXP;

还有其他可能有用的选项，但您应该参考文档以了解这些详细信息。