我正在使用proc genmod
来估计风险比。我的两个预测变量具有 2 个以上的水平。 我的gender
变量是‘male’, ‘female’, ‘other’
,种族是‘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;
我会从LSMEANS开始,然后可能是LSMESTIMATE进行更复杂的估计。 语法是这样的。
LSMEANS GENDER RACE / DIFF CL EXP;
还有其他选项可能有用,但您应该查阅文档以了解这些详细信息。