PHREG and LSMEANS: Eliminating Control Group / Generating Absolute Risk Estimate

Lynn1 — Thu, 11 Feb 2021 20:42:38 GMT

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

topic PHREG and LSMEANS: Eliminating Control Group / Generating Absolute Risk Estimate in Statistical Procedures

PHREG and LSMEANS: Eliminating Control Group / Generating Absolute Risk Estimate