I am working on an analysis of RERI from a PHREG output. I am first trying to generate a lsmeans output that shows the four absolute estimates for defined binary 'Delivery Ischemic Event' (1,0) and 'C-Section' (1,0). When I run this, a control level is automatically selected. Is there a way to remove this in this statement so there is data populated in all four rows, eliminating the control/reference group in the lsmeans output? In this attached case; a control level was set at Delivery Ischemic Event=0, C-Section=0. My code and the output I am trying to modify is below.
ods exclude none;
Title "PHREG Delivery Index Event: Any Ischemic Event C/S RERI";
PROC PHREG DATA=DELIVDATA;
strata NRD_STRATUM;
weight DISCWT;
class DELIVANYISCHEMIC (ref='0') DELIVCS (ref='0') /param=glm;
model days_frombirth*EXPOSUREDX(0) = DELIVANYISCHEMIC DELIVCS DELIVANYISCHEMIC*DELIVCS / type3;
hazardratio DELIVANYISCHEMIC/diff=all;
lsmeans DELIVANYISCHEMIC*DELIVCS /cl means diff exp plots=none;
ods output lsmeans=ph_lsmeans;
run;
quit;
Title "LS Means Output";
PROC PRINT DATA=work.ph_lsmeans;
run;
Many Thanks, let me know if any additional information is required
