I'm attempting to automate graphing of outcome variable mean values & error bars by embedding proc sgplot inside a macro do loop.   Unfortunately my code is only plotting the first variable over and over when it loops rather than plotting all 28 variables.   Here's a sample of my code & data for 2 variables that replicates the issue.   Thoughts what I'm doing wrong? data macrovars3; input trmt $ prez1 prez2 postz1 postz2; cards; No 0 1 0 0 No 1 1 0 1 No 1 0 1 0 No 1 0 1 0 No 0 0 0 0 No 0 0 0 0 No 0 0 0 0 No 0 0 0 0 No 0 0 0 0 Y 1 0 1 0 Y 1 1 1 0 Y 1 1 0 1 Y 1 1 1 1 Y 0 0 0 0 Y 0 0 1 0 Y 1 0 0 0 Y 1 1 0 1 Y 1 1 0 1 Y 1 0 1 0 Y 1 0 1 0 ; run; /* Macro variables for medical diagnosis categories */ %let z1 = CHF; %let z2 = VALVE; /* Macro graphing loop */ %macro RunPlot(DSName, NumVars); %do i=1 %to &NumVars; /* Stack datasets for pre and post outcomes and diagnosis, create period indicator. */ data cohort_match_real_data_ge0_pre (keep=trmt period diagnosis); set &DSName; length period $4.; period = "Pre"; diagnosis = prez&i; run; data cohort_match_real_data_ge0_post (keep=trmt period diagnosis); set &DSName; length period $4.; period = "Post"; diagnosis = postz&i; run; data cohort_match_prepost; set cohort_match_real_data_ge0_pre cohort_match_real_data_ge0_post; run; /* Calculate the means and the standard errors. */ proc sort data=cohort_match_prepost out=cohort_match_prepost_2; by trmt period; run; ods trace on; proc means data=cohort_match_prepost_2; by trmt period; var diagnosis; ods output summary=meansout; run; /* Calculate the upper and lower error bar values. */ data meansout; set meansout; lowerCI_diagnosis = diagnosis_Mean - 1.96*diagnosis_StdDev/diagnosis_N**.5; upperCI_diagnosis = diagnosis_Mean + 1.96*diagnosis_StdDev/diagnosis_N**.5; run; /* Plot mean values of pre and post diagnoses */ ods listing close; ods html file='scatterplot.html' path='.' style=styles.statistical; ods graphics / reset attrpriority=color width=600px height=400px imagefmt=gif; proc sgplot data=meansout; scatter x=period y=diagnosis_Mean / group=trmt yerrorlower=lowerCI_diagnosis yerrorupper=upperCI_diagnosis markerattrs=(symbol=CircleFilled); series x=period y=diagnosis_Mean / group=trmt lineattrs=(thickness=2); yaxis grid gridattrs=(color=white); xaxis grid gridattrs=(color=white pattern=shortdash) offsetmin=0.2 offsetmax=0.2; styleattrs backcolor=CXFFFFFF wallcolor=CXE0E0E0; /* keylegend "pre" "post";*/ title1 "Plot pre_&&z&i and post_&&z&i Means with Standard Error Bars"; run; ods html close; ods listing; %end; %mend; %RunPlot(macrovars3, 2) ... View more

I'm running a large number of logistic regressions and would like to specify particular output statistics from the proc logistic ods output dataset.   For example, using the SAS Remission dataset from here:    data Remission; input remiss cell smear infil li blast temp; label remiss='Complete Remission'; datalines; 1 .8 .83 .66 1.9 1.1 .996 1 .9 .36 .32 1.4 .74 .992 0 .8 .88 .7 .8 .176 .982 0 1 .87 .87 .7 1.053 .986 1 .9 .75 .68 1.3 .519 .98 0 1 .65 .65 .6 .519 .982 1 .95 .97 .92 1 1.23 .992 0 .95 .87 .83 1.9 1.354 1.02 0 1 .45 .45 .8 .322 .999 0 .95 .36 .34 .5 0 1.038 0 .85 .39 .33 .7 .279 .988 0 .7 .76 .53 1.2 .146 .982 0 .8 .46 .37 .4 .38 1.006 0 .2 .39 .08 .8 .114 .99 0 1 .9 .9 1.1 1.037 .99 1 1 .84 .84 1.9 2.064 1.02 0 .65 .42 .27 .5 .114 1.014 0 1 .75 .75 1 1.322 1.004 0 .5 .44 .22 .6 .114 .99 1 1 .63 .63 1.1 1.072 .986 0 1 .33 .33 .4 .176 1.01 0 .9 .93 .84 .6 1.591 1.02 1 1 .58 .58 1 .531 1.002 0 .95 .32 .3 1.6 .886 .988 1 1 .6 .6 1.7 .964 .99 1 1 .69 .69 .9 .398 .986 0 1 .73 .73 .7 .398 .986 ;  I'd like to run the following ods output, but SAS doesn't like me combining two conditions:   ods trace on; proc logistic data=Remission; model remiss(event='1') = cell smear infil li blast temp; ods output ParameterEstimates = logit_model (where=(Variable=temp)(keep=Variable Estimate ProbChiSq)); run;  Is there a way to make this work? I've tried a few variations with & without parentheses but no success.   ... View more

I have a list of 28 variable names I'd like to convert to x1-x28 in order to use a macro do loop for running proc genmod.   I've tried using an array as follows to assign values to x1-x28, but this doesn't work:   data macrovars (keep=y x1-x28); set health_data; array x[28] pre_CHF pre_VALVE pre_PULMCIRC pre_PERIVASC pre_PARA pre_NEURO pre_CHRNLUNG pre_DM pre_DMCX pre_HTN_C pre_HYPOTHY pre_RENLFAIL pre_LIVER pre_ULCER pre_AIDS pre_LYMPH pre_METS pre_TUMOR pre_ARTH pre_COAG pre_OBESE pre_WGHTLOSS pre_BLDLOSS pre_ANEMDEF pre_ALCOHOL pre_DRUG pre_PSYCH pre_DEPRESS; output x1-x28; run; Is there a way to do this using arrays? ... View more