proc contents data = bact1; run; data big_data.comp2_noentero; length group $30; set comp2; group = 'No Non-tidal Enterococci'; run; data big_data.comp_entero; length group $30; set big_data.comp_entero; group = 'Non-tidal Enterococci Included'; run; data entero_noentero; set big_data.comp_entero big_data.comp2_noentero; run; proc sql; select distinct group from entero_noentero; quit; proc sql; create table plot_big as select x.date1, x.big_all as big_ent, x.nobig as nobig_ent, y.big_all as big_noent, y.nobig as nobig_noent from big_data.comp_entero as x left join big_data.comp2_noentero as y on x.date1 = y.date1; quit; ods graphics on; ods listing sge = on; footnote; Title "BIG Moving Geomean With and Without Non-Tidal Enterococci Data"; proc sgplot data=plot_big noautolegend; yaxis label='Ratio of Geometric Mean to PCR Standard' min=2 max=8 values= ( 3 to 9 by 1); xaxis label = 'Year'; y2axis label = ' ' min=3 max = 9 ; pbspline x=date1 y=big_ent / smooth = 0.05 nomarkers lineattrs= (color=red thickness=2) legendlabel = 'BIG Area- Non-tidal Enterococci Included' name = 'w'; pbspline x=date1 y=big_noent / smooth = 0.05 nomarkers lineattrs= (pattern = shortdash color=green thickness=2) legendlabel = 'BIG Area - No Non-tidal Enterococci' name = 'o'; keylegend 'w' 'o' / location = inside position = topright; run; title " Moving Geomean Outside BIG Area With and Without Non-Tidal Enterococci Data"; proc sgplot data=plot_big noautolegend; yaxis label='Ratio of Geometric Mean to PCR Standard' min=2 max=8 values= ( 1 to 3 by 0.5); xaxis label = 'Year'; y2axis label = ' ' min = 1 max = 3 ; pbspline x=date1 y=nobig_ent / smooth = 0.05 nomarkers lineattrs= (color=red thickness=2) legendlabel = 'Outside BIG Area- Non-tidal Enterococci Included' name = 'w'; pbspline x=date1 y=nobig_noent / smooth = 0.05 nomarkers lineattrs= (pattern = shortdash color=green thickness=2) legendlabel = 'Outside BIG Area - No Non-tidal Enterococci' name = 'o'; keylegend 'w' 'o' / location = inside position = topright; run;