Ok, I am extremely new to SAS. In fact, this is the 1st time I've used it since I took the class. We use EG 6.1 I have to create the attached map or something close to it with the # of hospitals in each county here in Ohio. I've managed to find examples of code to create the Ohio Map. I've also found code to shade it by the number of hospitals. The issue that i am having is being able to shade the map and label it at the same time. I even managed to import a field into a dataset with the number of hospitals. When i try to combine the 2, i get an error because the variable i assigned is not working (I guess). I'm using the CNTYNAME set and then i added the a field hospitals, which had the number of hospitals for that county. I also added another field called CNTYHOSP that combined the name & # of hospitals in the county together. Then i imported that back into SAS and called it OHCNTYNAME. All of our SAS experts are out of the office until Monday, so...any advice would be greatly appreciated.
This code gave me the county name labels
/****************************************************************/
/* This sample program annotates the county names on a map */
/* of Ohio. It uses the %MAPLABEL annotate macro to create */
/* the annotate data set used with the GMAP procedure. */
/* This annotate macro is available beginning with Version 9 of */
/* the SAS System. */
/****************************************************************/
/* Set the graphics environment */
goptions reset=all cback=white border htitle=12pt htext=10pt;
/* Create a projected map of Ohio with county boundaries */
proc gproject data=maps.counties out=ohio;
/* State FIPS for Ohio */
where state=39;
id state county;
run;
quit;
/* Create a data set containing the names of the counties */
/* in Ohio. */
data ohnames;
set WORK.OHCNTYNAME2;
/* State FIPS code for Ohio */
where state=39;
run;
/* Make the annotate macros available */
%annomac;
/* Use the MAPLABEL annotate macro to create an annotate data set */
/* to label the counties. */
/* The macro parameters are: */
/* Input Map data set. */
/* Attribute data set that contains variable for the label. */
/* Resulting output data set for use in ANNO= in GMAP. */
/* Variable for label on. Can be text or numeric. */
/* Space-separated list of IDs that the map and attribute */
/* data sets are sorted on. */
/* FONT= Font for the label. */
/* COLOR= Color of the label. Default is BLACK. */
/* SIZE= Size of the label. Default is 1. */
/* HSYS= UNIT system for SIZE=. Default is 3 (PERCENT). */
%maplabel(ohio,ohnames,anno,cntyhosp,state county,font='Albany AMT/',
color=black,size=1.1);
title1 "County Map of Ohio";
pattern1 v=me c=black r=15;
/* Generate a county map of Ohio with county labels */
proc gmap data=ohio map=ohio;
id state county;
choro county / anno=anno nolegend coutline=grayaa;
run;
quit;
THIS ONE ALLOWED ME TO GET THE LEGEND AND COLOR IT
value ohiohos_
1 = '1'
2 = '2'
3 - 5 = '3 - 5'
6 - 21 = '6-21'
;
run;
* 2000 census population, New York State counties;
data numofhos0;
input county ohiohos @@;
datalines;
001 1 003 5 005 1 007 3
009 2 011 1 013 3 015 0
017 5 019 0 021 1 023 3
025 1 027 1 029 2 031 1
033 2 035 21 037 1 039 3
041 1 043 1 045 2 047 1
049 17 051 1 053 1 055 2
057 2 059 2 061 16 063 2
065 1 067 1 069 1 071 2
073 1 075 1 077 2 079 1
081 3 083 1 085 2 087 1
089 1 091 1 093 5 095 13
097 1 099 7 101 1 103 2
105 0 107 1 109 1 111 0
113 12 115 0 117 1 119 2
121 0 123 1 125 1 127 0
129 1 131 1 133 1 135 0
137 0 139 3 141 1 143 3
145 2 147 1 149 1 151 6
153 12 155 4 157 3 159 1
161 2 163 0 165 1 167 2
169 2 171 2 173 1 175 1
;
run;
* fill patterns for the map areas (gray-scale fills);
pattern1 v=s c=grayff;
pattern2 v=s c=grayda;
pattern3 v=s c=grayaa;
pattern4 v=s c=gray68;
proc gmap
data=numofhos0
map=ohio;
id county;
choro ohiohos / discrete;
format ohiohos ohiohos_.;
title1 'NUMBER OF OHIO HOSPITALS';
run;
quit;
I just wanted to say thank you so much. I did make a couple of changes. I added the gproject code because it seemed to produce a straighter map of Ohio. We switched out the countynm varible for ohiohos to display the number of hospitals and added some border lines since some of the counties had the same number of hospitals. We will work on the color scheme, but thank you all so much. Hopefully i can get to the point where i am helping others. I added the pics of the map from your code and my code.
proc gproject data=maps.counties out=ohio;
/* State FIPS for Ohio */
where state=39;
id state county;
run;
quit;
proc format library=work;
value ohiohos_
1 = '1'
2 = '2'
3 - 5 = '3 - 5'
6 - 21 = '6-21'
;
run;
data work.numofhos0;
input county ohiohos @@;
HospText = put(ohiohos,ohiohos_.);
datalines;
001 1 003 5 005 1 007 3
009 2 011 1 013 3 015 0
017 5 019 0 021 1 023 3
025 1 027 1 029 2 031 1
033 2 035 21 037 1 039 3
041 1 043 1 045 2 047 1
049 17 051 1 053 1 055 2
057 2 059 2 061 16 063 2
065 1 067 1 069 1 071 2
073 1 075 1 077 2 079 1
081 3 083 1 085 2 087 1
089 1 091 1 093 5 095 13
097 1 099 7 101 1 103 2
105 0 107 1 109 1 111 0
113 12 115 0 117 1 119 2
121 0 123 1 125 1 127 0
129 1 131 1 133 1 135 0
137 0 139 3 141 1 143 3
145 2 147 1 149 1 151 6
153 12 155 4 157 3 159 1
161 2 163 0 165 1 167 2
169 2 171 2 173 1 175 1
;
run;
proc sql;
create table work.ohnames as
select distinct county, countynm label='County Name'
from sashelp.zipcode
where state=39;
quit;
%annomac;
%maplabel(work.ohio,work.ohnames,work.nameanno,countynm, county,font='Albany AMT/',
color=black,size=1.5);
%maplabel(work.ohio,work.numofhos0,work.Hospanno,ohiohos, county,font='Albany AMT/',
color=black,size=1.7);
data work.combinedanno;
set
work.nameanno (in=Name)
work.Hospanno (in=Hosp);
if Name then position = '2';
Else if Hosp then position = '8';
/* this is where you could put if statements to adjust
the x and or y coordinates to display name or value
*/
run;
title1 "Map of Ohio Hospitals By County";
/* Generate a county map of Ohio with county labels */
pattern1 v=s c=grayff;
pattern2 v=s c=grayda;
pattern3 v=s c=grayaa;
pattern4 v=s c=gray68;
pattern5 v=s c=bibg;
proc gmap data=work.numofhos0 map=work.ohio;
id county;
choro ohiohos / anno=work.combinedanno
discrete
coutline=black
;
format ohiohos ohiohos_.;
label ohiohos= "Number of Hospitals";
run;
quit; title;
proc gproject data=maps.counties out=ohio;
/* State FIPS for Ohio */
where state=39;
id state county;
run;
quit;
proc sort data=mapsgfk.us_counties out=centroids;
by state county segment;
run;
data centroids_temp;
retain yi yj xi xj a cx cy x0 y0 0;
set centroids(keep=state county segment lat long rename=(lat=yj long=xj));
by state county segment;
if(first.segment) then do;
cx = 0;
cy = 0;
a = 0;
x0 = xj;
y0 = yj;
end;
else if(not first.segment) then do;
ta = (xi*yj - xj*yi);
cx + ((xi+xj)*ta);
cy + ((yi+yj)*ta);
a + ta;
end;
if(last.segment) then do;
ta = (xj*y0 - x0*yj);
cx + ((xj+x0)*ta);
cy + ((yj+y0)*ta);
a = ta + a * 0.5;
cx = cx / (6*a);
cy = cy / (6*a);
output;
end;
xi = xj;
yi = yj;
run;
proc sql;
create table centroid_weight as
select
state, county, sum(a) as sum
from centroids_temp
group by state, county;
quit;
proc sql;
create table centroids as
select a.state, a.county,
sum(cx*(a / sum)) as lat,
sum(cy*(a / sum)) as long
from centroids_temp a
inner join centroid_weight b
on (a.state = b.state and a.county = b.county)
group by a.state, a.county;
quit;
proc sql;
drop table centroids_temp;
drop table centroid_weight;
quit;
PROC SQL;
CREATE TABLE WORK.FILTER_FOR_CENTROIDS AS
SELECT t1.STATE,
t1.COUNTY,
t1.lat,
t1.long
FROM WORK.CENTROIDS t1
WHERE t1.STATE = 39;
QUIT;
proc format library=work;
value ohiohos_
1 = '1'
2 = '2'
3 - 5 = '3 - 5'
6 - 7 = '6 - 7'
7 - 21 = '7 or more'
;
run;
data work.numofhos0;
input county ohiohos @@;
HospText = put(ohiohos,ohiohos_.);
datalines;
001 1 003 5 005 1 007 3
009 2 011 1 013 3 015 0
017 5 019 0 021 1 023 3
025 1 027 1 029 2 031 1
033 2 035 21 037 1 039 3
041 1 043 1 045 2 047 1
049 17 051 1 053 1 055 2
057 2 059 2 061 16 063 2
065 1 067 1 069 1 071 2
073 1 075 1 077 2 079 1
081 3 083 1 085 2 087 1
089 1 091 1 093 5 095 13
097 1 099 7 101 1 103 2
105 0 107 1 109 1 111 0
113 12 115 0 117 1 119 2
121 0 123 1 125 1 127 0
129 1 131 1 133 1 135 0
137 0 139 3 141 1 143 3
145 2 147 1 149 1 151 6
153 12 155 4 157 3 159 1
161 2 163 0 165 1 167 2
169 2 171 2 173 1 175 1
;
run;
proc sql;
create table work.ohnames as
select distinct county, countynm label='County Name'
from sashelp.zipcode
where state=39;
quit;
%annomac;
%maplabel(work.ohio,work.ohnames,work.nameanno,countynm, county,font='Albany AMT/''boldt',
color=black,size=1.5);
%maplabel(work.ohio,work.numofhos0,work.Hospanno,ohiohos, county,font='Albany AMT/''boldt',
color=black,size=1.7);
data work.combinedanno;
set
WORK.CENTROIDS
work.nameanno (in=Name)
WORK.CENTROIDS
work.Hospanno (in=Hosp);
if Name then position = '2';
Else if Hosp then position = '8';
/* this is where you could put if statements to adjust
the x and or y coordinates to display name or value
*/
run;
title1 "Map of Ohio Hospitals By County";
/* Generate a county map of Ohio with county labels */
pattern1 v=s color=azure;
pattern2 v=s c=grayda;
pattern3 v=s color=biy;
pattern4 v=s color=bisque;
pattern5 v=s c=bibg;
pattern6 v=s color=cyan;
proc gmap data=work.numofhos0 map=work.ohio;
id county;
choro ohiohos / anno=work.combinedanno
discrete
coutline=black
;
format ohiohos ohiohos_.;
label ohiohos= "Number of Hospitals";
run;
quit; title;
I've been off for a couple of days... anyway, added some code that I found about centering the county labels, which is the centroid bit. Now the issue remains as previous, that the lake counties are not formatting proprely. Is there a way to scale the output to make it larger? I'm thinking that i could get the labels to fit if the output size of the map was bigger, maybe? Still working on color scheme
Hello everyone again. Hopefully this will be my last post on this subject. My map looks fairly good thanks to the help of all of you out there. I just have one problem and i can't seem to figure it out. It's the labels for the Ottawa (123) & Erie (43). See how messed they are. I tried playing with the postion code and it managed to the county & number close to each other. I someone can help me figure out those 2 counties. I tried if then statements but i
proc gproject data=maps.counties out=ohio;
/* State FIPS for Ohio */
where state=39;
id state county;
run;
quit;
proc sort data=mapsgfk.us_counties out=centroids;
by state county segment;
run;
data centroids_temp;
retain yi yj xi xj a cx cy x0 y0 0;
set centroids(keep=state county segment lat long rename=(lat=yj long=xj));
by state county segment;
if(first.segment) then do;
cx = 0;
cy = 0;
a = 0;
x0 = xj;
y0 = yj;
end;
else if(not first.segment) then do;
ta = (xi*yj - xj*yi);
cx + ((xi+xj)*ta);
cy + ((yi+yj)*ta);
a + ta;
end;
if(last.segment) then do;
ta = (xj*y0 - x0*yj);
cx + ((xj+x0)*ta);
cy + ((yj+y0)*ta);
a = ta + a * 0.5;
cx = cx / (6*a);
cy = cy / (6*a);
output;
end;
xi = xj;
yi = yj;
run;
proc sql;
create table centroid_weight as
select
state, county, sum(a) as sum
from centroids_temp
group by state, county;
quit;
proc sql;
create table centroids as
select a.state, a.county,
sum(cx*(a / sum)) as lat,
sum(cy*(a / sum)) as long
from centroids_temp a
inner join centroid_weight b
on (a.state = b.state and a.county = b.county)
group by a.state, a.county;
quit;
proc sql;
drop table centroids_temp;
drop table centroid_weight;
quit;
PROC SQL;
CREATE TABLE WORK.FILTER_FOR_CENTROIDS AS
SELECT t1.STATE,
t1.COUNTY,
t1.lat,
t1.long
FROM WORK.CENTROIDS t1
WHERE t1.STATE = 39;
QUIT;
proc format library=work;
value ohiohos_
1 = '1'
2 = '2'
3 - 5 = '3 - 5'
6 - 7 = '6 - 7'
7 - 21 = '7 or more'
;
run;
data work.numofhos0;
input county ohiohos @@;
HospText = put(ohiohos,ohiohos_.);
datalines;
001 1 003 5 005 1 007 3
009 2 011 1 013 3 015 0
017 5 019 0 021 1 023 3
025 1 027 1 029 2 031 1
033 2 035 21 037 1 039 3
041 1 043 1 045 2 047 1
049 17 051 1 053 1 055 2
057 2 059 2 061 16 063 2
065 1 067 1 069 1 071 2
073 1 075 1 077 2 079 1
081 3 083 1 085 2 087 1
089 1 091 1 093 5 095 13
097 1 099 7 101 1 103 2
105 0 107 1 109 1 111 0
113 12 115 0 117 1 119 2
121 0 123 1 125 1 127 0
129 1 131 1 133 1 135 0
137 0 139 3 141 1 143 3
145 2 147 1 149 1 151 6
153 12 155 4 157 3 159 1
161 2 163 0 165 1 167 2
169 2 171 2 173 1 175 1
;
run;
proc sql;
create table work.ohnames as
select distinct county, countynm label='County Name'
from sashelp.zipcode
where state=39;
quit;
%annomac;
%maplabel(work.ohio,work.ohnames,work.nameanno,countynm, county,font='Albany AMT/''boldt',
color=black,size=1.7);
%maplabel(work.ohio,work.numofhos0,work.Hospanno,ohiohos, county,font='Albany AMT/''boldt',
color=black,size=1.8);
data work.combinedanno;
set
WORK.CENTROIDS
work.nameanno (in=Name)
WORK.CENTROIDS
work.Hospanno (in=Hosp);
if Name then position = '2.1';
Else if Hosp then position = '5';
/* this is where you could put if statements to adjust
the x and or y coordinates to display name or value
*/
;
run;
title1 "Map of Ohio Hospitals By County";
/* Generate a county map of Ohio with county labels */
pattern1 v=s color=azure;
pattern2 v=s c=grayda;
pattern3 v=s color=biy;
pattern4 v=s color=bisque;
pattern5 v=s c=bibg;
pattern6 v=s color=cyan;
proc gmap data=work.numofhos0 map=work.ohio;
id county;
choro ohiohos / anno=work.combinedanno
discrete
coutline=black
;
format ohiohos ohiohos_.;
label ohiohos= "Number of Hospitals";
run;
quit; title;
guess i can;t figure out which variable i need to alter. Thanks again everyone, take care
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
-
2760 views
-
0 likes
-
-
