/*This file maps Kansas Congressional District Three townships by 2000 Population */
/* If you have questions, please contact:
	Xanthippe Stevens
	Information Specialist
	The University of Kansas
	Policy Research Institute
	1541 Lilac Lane, Suite 607
	Lawrence, KS 66044
	(785) 864-9111
	xan@ku.edu */

libname tiger 'd:\data\tiger\';
libname data 'd:\data\census\2000\pl94171\';

data sumlvl60;
	set data.sumlvl60; /*Data on township level*/
		where county='045' or county='091' or county='121' or county='209';
	keep cousubr name arealand totpop ;
	length cousubr 8;
	cousubr=cousub;
	totpop=P0010001;
run;
data dist3;
	set tiger.cousub20045 tiger.cousub20091 tiger.cousub20121 tiger.cousub20209;
run;
/* Project data */
proc gproject data=dist3 out=task4 degree asis eastlong;
  id line;
run;
proc sort data=task4;
	by cousubr segment;
run;
data testanno;
	length function color $8 cousub 8;
	set task4 (keep=x y line type tlid geo cousubr segment);
		by cousubr segment;
	cousub=cousubr;
	if first.segment then do
		row=line;
		line=1;
		xsys = "2";
		ysys = "2";
		function="poly";
		color="cxFFFF99";
		style="msolid";
		output;
		end;
	else do
		row=line;
		function = "polycont";
		xsys = "2";
		ysys = "2";
		line = 1;
		style="msolid";
		color="cxFF9999";
		output;
		end;
	drop cousubr;
run;

data tanno3;
	set testanno(rename=(cousub=cousubr));
	output;
run;

pattern1 value=msolid;
goptions reset=global colors=(cxEEEEFF cxCCCCFF cx9999FF  cx6666DD 
		cx3333FF cx000099  cx2222FF)
	gunit=pct border ftext=swissb ctext=black htitle=6 htext=3;

proc gmap map=tanno3 data=sumlvl60;
	id cousubr;
	choro totpop/levels=6
				nolegend
				coutline=black;
run;

quit;