I used FIRTH option in proc logistic and proc phreg as there is a zero event  in one of the two groups and more than 60% events in another group.  I got extreme values of OR and HR (greater than 100) and very wide CI and p-value <0.0001. Can I interpret the results as the same as for regular proc logistic and proc phreg? Any input/opinion based on your experience and knowledge would be appreciated. ... View more

 anyone used proc MI for longitudinal data in long format before?   data da;  input id$ visitn trt val; data lines; 001 1 A 32 001 2 A 29 001 3 A 25 001 4 A 30 001 5 A 22 002 1 A . 002 2 A 20 002 3 A . 002 4 A 18 002 5 A 28 003 1 A  15 003 2 A 10 003 3 A . 003 4 A 18 003 5 A 18 004 1 B 15 004 2 B 17 004 3 B 20 004 4 B 18 004 5 B . 005 1 B 16 005 2 B 14 005 3 B 20 005 4 B . 005 5 B . 006 1 B  19 006 2 B 17 006 3 B 20 006 4 B . 006 5 B 18 ; run;     proc mi data=da nimpute=50 out=da1; mcmc chain=multiple; var visitn val; by id; run;     ... View more

I tried to use thee following code to add two reference lines in thee plot, why the label for second line doesn't show in the plot?     proc sgplot data=indata; series x=visit y=result/group=usubjid datalabelattrs=(size=12pt) markers markerattrs=(size=12pt symbol=trianglefilled) lineattrs=(thickness=3 pattern=solid) name="series_name"; keylegend "series_name" / valueattrs=(size=14); refline 0 / axis=y lineattrs=(thickness=3 pattern=dash) label=("Upper Level of Healthy"); refline LLOQn / axis=y lineattrs=(thickness=3 pattern=dash) label=("LLOQ") labelloc=inside;   xaxis label="Visit" labelattrs=(size=16pt) VALUEATTRS=(Size=12         Weight=Bold); yaxis label="sample data" labelattrs=(size=16pt) VALUEATTRS=( Size=12          Weight=Bold); run; ... View more