I need to output EPS vector files for publication, when I run this code I get the error: WARNING: SASREPORT13 destination does not support EPS images. Using the default static format.   ods graphics / reset outputfmt=eps noborder; ods listing device=eps dpi=300 file= 'AJE-00824-2018 Sheridan Figure 1a.eps' gpath="G:\CA Birth Cohort AP\Final Data\";   proc sgplot data=week noborder; /*code****/ run; ods listing close; ods graphics / reset; ... View more

I would like my yaxis values to have two decimal places, this code does not work:   proc sgplot data=week noborder; format HazardRatio 4.2; scatter x=week y=HazardRatio / yerrorlower=WaldLower datalabelattrs=(size=12pt color=black) yerrorupper=WaldUpper errorbarattrs=(color=black) markerattrs=hr (symbol=CircleFilled size=6 color=black); refline 1 / axis=y lineattrs=(pattern=shortdash) transparency=.5; yaxis label="Hazard Ratio " min=.96 max=1.06 type=log logbase=10 valueattrs=(size=14pt color=black family="Times New Roman") labelattrs=(size=14pt color=black family="Times New Roman"); xaxis label="Week of Gestation" min=1 max=37 valueattrs=(size=14pt color=black family="Times New Roman") labelattrs=(size=14pt color=black family="Times New Roman") minor; run;   Thanks! ... View more

I am trying to overlay two plots of count vs. predicted probability with confidence intervals from glimmix output. I would like to show the predicted probabilities by group (glimmix was run by group, only two groups). I attempted to do this using proc plm and effectplot, but I don't see a way to overlay plots when there is no interaction term present in the model.   Thank you! ... View more

I am using proc phreg to estimate hazard ratios, my model is as follows:    time*event(0) = week|avg education race; random zipcode;   I would like to get the hazard ratio for 'week'. I was originally doing this by hand by summing the estimate from week and the interaction term, but assume there is a way to do this using estimate/contrast?    Thank you! ... View more