Hi All   I'm reproducing a Kaplan-Meier stepplot in layout lattice.   Here is the result.         I'm wondering why in the right graph, the axis table with patients at risk is located a bit higher than the corresponding table on the left.   I attach the code I used, if of any help.   proc template; define statgraph Fig_2_KM; begingraph / collation=binary designheight=480 designwidth=1280; layout lattice / rows=1 columns=2 rowgutter=10 columngutter=10 ; /*Unmatched population*/ legenditem type=text name="log" / label="Logrank p=0.0007" lineattrs=(color=black); layout overlay / xaxisopts=(Label="Years of Follow-up" LabelAttrs=(Weight=normal) type=linear) y2axisopts=(labelFitPolicy=Split) yaxisopts=(Label=" " labelFitPolicy=Split LabelAttrs=(Weight=normal) type=linear linearopts=(viewmin=0) ) y2axisopts=(labelFitPolicy=Split); entry "in the Unmatched Population" / valign=top location=outside textattrs=(weight=bold); entry "Product-Limit Survival Estimates"/ valign=top location=outside textattrs=(weight=bold size=12); StepPlot X=Time Y=Survival / primary=true Group=Stratum Lineattrs=(Pattern=1 Thickness=1) LegendLabel="Survival Probability" NAME="STEP"; scatterplot x=time y=censored / markerattrs=(symbol=plus) name='c'; scatterplot x=time y=censored / markerattrs=(symbol=plus) GROUP=stratum; discretelegend 'c' 'log' / location=inside halign=right valign=top order=columnmajor valueattrs=(size=10) border=false; drawtext textattrs=(size=10) 'Survival Probability' / x=-8 y=60 anchor=bottom xspace=wallpercent yspace=wallpercent rotate=90 width=50; /*--Subjects at risk-Unmatched--*/ InnerMargin / align=bottom opaque=true; AxisTable Value=AtRisk X=tAtRisk / labelPosition=min Class=Stratum LabelAttrs=(Weight=normal) ValueAttrs=(Size=10) colorGroup=Stratum Display=(Label Values) title="Subjects at risk" titleattrs=(size=10); EndInnerMargin; DiscreteLegend "STEP"; endlayout; /*Matched population*/ legenditem type=text name="log2" / label="Logrank p=0.9417" lineattrs=(color=black); layout overlay / xaxisopts=(Label="Years of Follow-up" LabelAttrs=(Weight=normal) type=linear) y2axisopts=(labelFitPolicy=Split) yaxisopts=(Label=" " labelFitPolicy=Split LabelAttrs=(Weight=normal) type=linear linearopts=(viewmin=0) ) y2axisopts=(labelFitPolicy=Split); entry "in the Matched Population" / valign=top location=outside textattrs=(weight=bold); entry "Product-Limit Survival Estimates"/ valign=top location=outside textattrs=(weight=bold size=12); StepPlot X=Time2 Y=Survival2 / primary=true group=Stratum2 Lineattrs=(Pattern=1 Thickness=1) LegendLabel="Survival Probability" NAME="STEP2"; scatterplot x=time2 y=censored2 / markerattrs=(symbol=plus) name='c2' legendlabel="Censored"; scatterplot x=time2 y=censored2 / markerattrs=(symbol=plus) GROUP=stratum2; discretelegend 'c2' 'log2' / location=inside halign=right valign=top order=columnmajor valueattrs=(size=10) border=false; drawtext textattrs=(size=10) 'Survival Probability' / x=-8 y=60 anchor=bottom xspace=wallpercent yspace=wallpercent rotate=90 width=50; /*--Subjects at risk-Unmatched--*/ InnerMargin / align=bottom opaque=true; AxisTable Value=AtRisk2 X=tAtRisk2 / labelPosition=min class=Stratum2 LabelAttrs=(Weight=normal) ValueAttrs=(Size=10) colorGroup=Stratum2 Display=(Label Values) title="Subjects at risk" titleattrs=(size=10); EndInnerMargin; DiscreteLegend "STEP"; endlayout; endlayout; endgraph; end; run; Thanks in advance     Sincerely   A ... View more

Hi all!   I have a dataset like this     Group   A   B   C   X     1       5    3   1   12     1       6    4   1   9     1       7    6   3   10     2       4    2   8   7     2       3    7   5   8     2       5    2   4   13     3       2    10  5   9     3       4    9   3  15     3       7    6   7   14     For each group, I need to extract only one row per group at a specific value of X (namely the max value of X in each group)   My final dataset should be like this   Group   A   B   C   X     1       5    3   1   12     2       5    2   4   13     3       4    9   3  15     Tks in advance     ... View more

Hi all     I would like to automatically use a value, derived from an output database, into a new data step.   I explain with an example:   proc logistic data=sample plots=none; model Y(event=1)=X / outroc=rocX; ods output ParameterEstimates=param; run; /*param DBS has two columns, Variable (in this example named Intercept and X) and Estimate*/   I would like to automatically include the value of Intercept and X, taken from the dbs param, into the following DATA Step   data rocX2(keep=cutoff prob Sensitivity Specificity Youden); set rocX; logit=log(_prob_/(1-_prob_)); cutoff=(logit-/*Intercept*/)//*X*/; prob=_prob_; Sensitivity=_SENSIT_; Specificity=1-_1MSPEC_; Youden=_SENSIT_+ (1-_1MSPEC_)-1; run; tsk in advance ... View more

Hi all   is there a macro/code that mimics in SAS what pROC and bootLR do in R?   my aim is to calculate bootstrap threshold, sensitivity, specificity, LT+ and LR- along with 95% CI   Tks! ... View more