/*AHCA*/ libname AHCA odbc dsn=AHCA schema=DBO;
proc format;
	value $ county /*PCCOUNTY*/ 
'01' = 'Alachua'
'02' = 'Baker'
'03' = 'Bay'
'04' = 'Bradford'
'05' = 'Brevard'
'06' = 'Broward'
'07' = 'Calhoun'
'08' = 'Charlotte'
'09' = 'Citrus'
'10' = 'Clay'
'11' = 'Collier'
'12' = 'Columbia'
'13' = 'Miami-Dade'
'14' = 'Desoto'
'15' = 'Dixie'
'16' = 'Duval'
'17' = 'Escambia'
'18' = 'Flagler'
'19' = 'Franklin'
'20' = 'Gadsden'
'21' = 'Gilchrist'
'22' = 'Glades'
'23' = 'Gulf'
'24' = 'Hamilton'
'25' = 'Hardee'
'26' = 'Hendry'
'27' = 'Hernando'
'28' = 'Highlands'
'29' = 'Hillsborough'
'30' = 'Holmes'
'31' = 'Indian River'
'32' = 'Jackson'
'33' = 'Jefferson'
'34' = 'Lafayette'
'35' = 'Lake'
'36' = 'Lee'
'37' = 'Leon'
'38' = 'Levy'
'39' = 'Liberty'
'40' = 'Madison'
'41' = 'Manatee'
'42' = 'Marion'
'43' = 'Martin'
'44' = 'Monroe'
'45' = 'Nassau'
'46' = 'Okaloosa'
'47' = 'Okeechobee'
'48' = 'Orange'
'49' = 'Osceola'
'50' = 'Palm Beach'
'51' = 'Pasco'
'52' = 'Pinellas'
'53' = 'Polk'
'54' = 'Putnam'
'55' = 'Saint Johns'
'56' = 'Saint Lucie'
'57' = 'Santa Rosa'
'58' = 'Sarasota'
'59' = 'Seminole'
'60' = 'Sumter'
'61' = 'Suwannee'
'62' = 'Taylor'
'63' = 'Union'
'64' = 'Volusia'
'65' = 'Wakulla'
'66' = 'Walton'
'67' = 'Washington'
;
run;


/***Gather one year of records from AHCA dataset***/
data work.inpt (KEEP = AGE admsrc DISCHSTAT FACLNBR FAC_COUNTY OTHDIAG1-OTHDIAG30 POA1-POA30 PRINDIAG POA_PRIN_D SYS_RECID DIAG01-DIAG31 LOSDAYS);

set AHCA.HospitalInpatient_Details;
where (REPORTYEAR in (2019) AND TYPE_SERV='1') AND MOD_CODE IN ('CL01' 'CL02' 'CL07' 'CL10'); **class 1 hosp, class2, class 1 rural, class 3 special medical**;

/***Zero out any diagnoses that are not traumatic injuries***/
array diagnoses {31} PRINDIAG OTHDIAG1-OTHDIAG30;
array present {31} POA_PRIN_D POA1-POA30;
array new_diag {31} $ DIAG01-DIAG31;
do i=1 to 31; /*Zero out the diagnosis value if it is not an an injury.*/
	if 	(substr(diagnoses{i},1,1)="S" AND substr(diagnoses{i},3,1) NE "0" AND substr(diagnoses{i},7,1) IN ("A","B","C"))
		OR substr(diagnoses{i},1,3) IN ("T07","T14","T30","T31","T32")
		OR (substr(diagnoses{i},1,3) IN ("T20","T21","T22","T23","T24","T25","T26","T27","T28") AND substr(diagnoses{i},7,1)="A")
		OR (substr(diagnoses{i},1,5) IN ("T79A1","T79A2","T79A3","T79A4","T79A5","T79A6","T79A7","T79A8","T79A9") AND substr(diagnoses{i},7,1)="A")
		THEN new_diag{i}=diagnoses{i}; else new_diag{i}=" ";
end;
do i=1 to 31; /*Zero out the diagnosis value if it is not an an injury.*/
	if 	(substr(diagnoses{i},1,1)="S" AND substr(diagnoses{i},3,1) NE "0" AND substr(diagnoses{i},7,1) IN ("A","B","C"))
		OR substr(diagnoses{i},1,3) IN ("T07","T14","T30","T31","T32")
		OR (substr(diagnoses{i},1,3) IN ("T20","T21","T22","T23","T24","T25","T26","T27","T28") AND substr(diagnoses{i},7,1)="A")
		OR (substr(diagnoses{i},1,5) IN ("T79A1","T79A2","T79A3","T79A4","T79A5","T79A6","T79A7","T79A8","T79A9") AND substr(diagnoses{i},7,1)="A")
		THEN diagnoses{i}=1; else diagnoses{i}=0;
end;

/*How POA indicators are coded: Y - Yes Present at time that order for inpatient admission occurs.
N - No Not present at the time that the order for inpatient admission occurs.
U - Unknown Documentation is insufficient to determine if condition is present on 
admission. W - Clinically Undetermined Provider is unable to clinically determine whether the 
condition was present on admission or not. E - Exempt The condition is exempt from POA reporting.*/
do i=1 to 31;	/*Zero out the diagnosis value if it is not present on admission.*/
	if present{i} IN ('N', 'E') THEN diagnoses{i}=0;
	if present{i} IN ('N', 'E') THEN new_diag{i}=" ";
end;

/***Elimiate records for patients with no injuries***/
injury = PRINDIAG+OTHDIAG1+OTHDIAG2+OTHDIAG3+OTHDIAG4+OTHDIAG5+OTHDIAG6+OTHDIAG7+OTHDIAG8+OTHDIAG9+OTHDIAG10+OTHDIAG11+
OTHDIAG12+OTHDIAG13+OTHDIAG14+OTHDIAG15+OTHDIAG16+OTHDIAG17+OTHDIAG18+OTHDIAG19+OTHDIAG20+OTHDIAG21+OTHDIAG22+OTHDIAG23+OTHDIAG24+
OTHDIAG25+OTHDIAG26+OTHDIAG27+OTHDIAG28+OTHDIAG29+OTHDIAG30;
if injury=0 THEN DELETE;

RUN;

/***Import survial risk ratios***/
data WORK.SRR    ;
%let _EFIERR_ = 0; /* set the ERROR detection macro variable */
infile 'C:\Users\cardks\Desktop\SRR_adult_2016to2019.csv' delimiter =',' MISSOVER DSD lrecl=32767 firstobs=2 ;
informat DIAGNOSIS $8. ;
informat Incidence best32. ;
informat Mortality best32. ;
informat SRR best32. ;
format DIAGNOSIS $8. ;
format Incidence best12. ;
format Mortality best12. ;
format SRR best12. ;
input DIAGNOSIS  $ Incidence Mortality SRR;
if _ERROR_ then call symputx('_EFIERR_',1);  /* set ERROR detection macro variable */
run;
data srr; set srr; DIAGNOSIS=compress(DIAGNOSIS,'.');run;

/***Assign SRRs to specific injuries***/
DATA work.SRR; SET work.SRR; DIAG01=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG01; RUN; 
PROC SORT data=work.inpt; BY DIAG01; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG01; IF InSet=1; SRR01=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG02=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG02; RUN; 
PROC SORT data=work.inpt; BY DIAG02; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG02; IF InSet=1; SRR02=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG03=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG03; RUN; 
PROC SORT data=work.inpt; BY DIAG03; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG03; IF InSet=1; SRR03=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG04=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG04; RUN; 
PROC SORT data=work.inpt; BY DIAG04; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG04; IF InSet=1; SRR04=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG05=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG05; RUN; 
PROC SORT data=work.inpt; BY DIAG05; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG05; IF InSet=1; SRR05=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG06=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG06; RUN; 
PROC SORT data=work.inpt; BY DIAG06; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG06; IF InSet=1; SRR06=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG07=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG07; RUN; 
PROC SORT data=work.inpt; BY DIAG07; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG07; IF InSet=1; SRR07=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG08=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG08; RUN; 
PROC SORT data=work.inpt; BY DIAG08; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG08; IF InSet=1; SRR08=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG09=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG09; RUN; 
PROC SORT data=work.inpt; BY DIAG09; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG09; IF InSet=1; SRR09=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG10=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG10; RUN; 
PROC SORT data=work.inpt; BY DIAG10; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG10; IF InSet=1; SRR10=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG11=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG11; RUN; 
PROC SORT data=work.inpt; BY DIAG11; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG11; IF InSet=1; SRR11=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG12=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG12; RUN; 
PROC SORT data=work.inpt; BY DIAG12; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG12; IF InSet=1; SRR12=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG13=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG13; RUN; 
PROC SORT data=work.inpt; BY DIAG13; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG13; IF InSet=1; SRR13=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG14=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG14; RUN; 
PROC SORT data=work.inpt; BY DIAG14; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG14; IF InSet=1; SRR14=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG15=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG15; RUN; 
PROC SORT data=work.inpt; BY DIAG15; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG15; IF InSet=1; SRR15=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG16=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG16; RUN; 
PROC SORT data=work.inpt; BY DIAG16; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG16; IF InSet=1; SRR16=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG17=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG17; RUN; 
PROC SORT data=work.inpt; BY DIAG17; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG17; IF InSet=1; SRR17=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG18=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG18; RUN; 
PROC SORT data=work.inpt; BY DIAG18; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG18; IF InSet=1; SRR18=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG19=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG19; RUN; 
PROC SORT data=work.inpt; BY DIAG19; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG19; IF InSet=1; SRR19=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG20=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG20; RUN; 
PROC SORT data=work.inpt; BY DIAG20; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG20; IF InSet=1; SRR20=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG21=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG21; RUN; 
PROC SORT data=work.inpt; BY DIAG21; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG21; IF InSet=1; SRR21=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG22=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG22; RUN; 
PROC SORT data=work.inpt; BY DIAG22; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG22; IF InSet=1; SRR22=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG23=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG23; RUN; 
PROC SORT data=work.inpt; BY DIAG23; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG23; IF InSet=1; SRR23=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG24=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG24; RUN; 
PROC SORT data=work.inpt; BY DIAG24; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG24; IF InSet=1; SRR24=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG25=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG25; RUN; 
PROC SORT data=work.inpt; BY DIAG25; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG25; IF InSet=1; SRR25=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG26=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG26; RUN; 
PROC SORT data=work.inpt; BY DIAG26; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG26; IF InSet=1; SRR26=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG27=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG27; RUN; 
PROC SORT data=work.inpt; BY DIAG27; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG27; IF InSet=1; SRR27=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG28=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG28; RUN; 
PROC SORT data=work.inpt; BY DIAG28; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG28; IF InSet=1; SRR28=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG29=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG29; RUN; 
PROC SORT data=work.inpt; BY DIAG29; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG29; IF InSet=1; SRR29=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG30=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG30; RUN; 
PROC SORT data=work.inpt; BY DIAG30; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG30; IF InSet=1; SRR30=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

DATA work.SRR; SET work.SRR; DIAG31=DIAGNOSIS; RUN; 
PROC SORT data=work.SRR; BY DIAG31; RUN; 
PROC SORT data=work.inpt; BY DIAG31; RUN;
DATA work.inpt; MERGE work.inpt (IN=InSet) work.SRR; BY DIAG31; IF InSet=1; SRR31=SRR; DROP DIAGNOSIS Incidence Mortality; RUN;

/***Calculate ICISS score***/
DATA work.inpt; SET work.inpt;
array new_diag {31} DIAG01-DIAG31;
array survival {31} SRR01-SRR31;
ICISS=1;
do i=1 to 31;
	IF new_diag{i}=' ' THEN survival{i}=1;
	IF survival{i}=. THEN survival{i}=1;
	ICISS=ROUND(ICISS*survival{i}, 0.001);
end;

drop diag01-diag31 prindiag OTHDIAG1-OTHDIAG30 POA_PRIN_D POA1-POA30 SRR01-SRR31 i;
RUN;
DATA work.inpt; SET work.inpt;
length facility_county $2. TSA $2.;

/***Segment records into severely injured and moderate/mild***/
IF ICISS GE 0.85 THEN Severity='Moderate/Mild'; ELSE IF ICISS LT 0.85 THEN Severity='Severe';

/***Re-format facility county to character value***/
IF FAC_COUNTY=1 THEN facility_county='01'; 
else if FAC_COUNTY=2 THEN facility_county='02';
else if FAC_COUNTY=3 THEN facility_county='03';
else if FAC_COUNTY=4 THEN facility_county='04';
else if FAC_COUNTY=5 THEN facility_county='05';
else if FAC_COUNTY=6 THEN facility_county='06';
else if FAC_COUNTY=7 THEN facility_county='07';
else if FAC_COUNTY=8 THEN facility_county='08';
else if FAC_COUNTY=9 THEN facility_county='09';
else if FAC_COUNTY=10 THEN facility_county='10';
else if FAC_COUNTY=11 THEN facility_county='11';
else if FAC_COUNTY=12 THEN facility_county='12';
else if FAC_COUNTY=13 THEN facility_county='13';
else if FAC_COUNTY=14 THEN facility_county='14';
else if FAC_COUNTY=15 THEN facility_county='15';
else if FAC_COUNTY=16 THEN facility_county='16';
else if FAC_COUNTY=17 THEN facility_county='17';
else if FAC_COUNTY=18 THEN facility_county='18';
else if FAC_COUNTY=19 THEN facility_county='19';
else if FAC_COUNTY=20 THEN facility_county='20';
else if FAC_COUNTY=21 THEN facility_county='21';
else if FAC_COUNTY=22 THEN facility_county='22';
else if FAC_COUNTY=23 THEN facility_county='23';
else if FAC_COUNTY=24 THEN facility_county='24';
else if FAC_COUNTY=25 THEN facility_county='25';
else if FAC_COUNTY=26 THEN facility_county='26';
else if FAC_COUNTY=27 THEN facility_county='27';
else if FAC_COUNTY=28 THEN facility_county='28';
else if FAC_COUNTY=29 THEN facility_county='29';
else if FAC_COUNTY=30 THEN facility_county='30';
else if FAC_COUNTY=31 THEN facility_county='31';
else if FAC_COUNTY=32 THEN facility_county='32';
else if FAC_COUNTY=33 THEN facility_county='33';
else if FAC_COUNTY=34 THEN facility_county='34';
else if FAC_COUNTY=35 THEN facility_county='35';
else if FAC_COUNTY=36 THEN facility_county='36';
else if FAC_COUNTY=37 THEN facility_county='37';
else if FAC_COUNTY=38 THEN facility_county='38';
else if FAC_COUNTY=39 THEN facility_county='39';
else if FAC_COUNTY=40 THEN facility_county='40';
else if FAC_COUNTY=41 THEN facility_county='41';
else if FAC_COUNTY=42 THEN facility_county='42';
else if FAC_COUNTY=43 THEN facility_county='43';
else if FAC_COUNTY=44 THEN facility_county='44';
else if FAC_COUNTY=45 THEN facility_county='45';
else if FAC_COUNTY=46 THEN facility_county='46';
else if FAC_COUNTY=47 THEN facility_county='47';
else if FAC_COUNTY=48 THEN facility_county='48';
else if FAC_COUNTY=49 THEN facility_county='49';
else if FAC_COUNTY=50 THEN facility_county='50';
else if FAC_COUNTY=51 THEN facility_county='51';
else if FAC_COUNTY=52 THEN facility_county='52';
else if FAC_COUNTY=53 THEN facility_county='53';
else if FAC_COUNTY=54 THEN facility_county='54';
else if FAC_COUNTY=55 THEN facility_county='55';
else if FAC_COUNTY=56 THEN facility_county='56';
else if FAC_COUNTY=57 THEN facility_county='57';
else if FAC_COUNTY=58 THEN facility_county='58';
else if FAC_COUNTY=59 THEN facility_county='59';
else if FAC_COUNTY=60 THEN facility_county='60';
else if FAC_COUNTY=61 THEN facility_county='61';
else if FAC_COUNTY=62 THEN facility_county='62';
else if FAC_COUNTY=63 THEN facility_county='63';
else if FAC_COUNTY=64 THEN facility_county='64';
else if FAC_COUNTY=65 THEN facility_county='65';
else if FAC_COUNTY=66 THEN facility_county='66';
else if FAC_COUNTY=67 THEN facility_county='67';

/***Assign counties to TSA***/
IF facility_county IN ('17' '46' '57' '66') THEN TSA='01';
ELSE IF facility_county IN ('03' '23' '30' '67') THEN TSA='02';
ELSE IF facility_county IN ('07' '19' '20' '32' '33' '37' '39' '40' '62' '65') THEN TSA='03';
ELSE IF facility_county IN ('01' '04' '12' '15' '21' '24' '34' '38' '54' '61' '63') THEN TSA='04';
ELSE IF facility_county IN ('02' '10' '16' '45' '55') THEN TSA='05';
ELSE IF facility_county IN ('09' '27' '42') THEN TSA='06';
ELSE IF facility_county IN ('18' '64') THEN TSA='07';
ELSE IF facility_county IN ('35' '48' '49' '59' '60') THEN TSA='08';
ELSE IF facility_county IN ('51' '52') THEN TSA='09';
ELSE IF facility_county IN ('29') THEN TSA='10';
ELSE IF facility_county IN ('25' '28' '53') THEN TSA='11';
ELSE IF facility_county IN ('05' '31') THEN TSA='12';
ELSE IF facility_county IN ('14' '41' '58') THEN TSA='13';
ELSE IF facility_county IN ('43' '47' '56') THEN TSA='14';
ELSE IF facility_county IN ('08' '11' '22' '26' '36') THEN TSA='15';
ELSE IF facility_county IN ('50') THEN TSA='16';
ELSE IF facility_county IN ('06') THEN TSA='17';
ELSE IF facility_county IN ('13' '44') THEN TSA='18';

RUN;

/***Assign trauma center status to hospitals***/
DATA work.inpt; SET work.inpt;

IF FACLNBR IN (
100131/*Aventura*/,
100093/*Baptist*/,
100026/*Bay*/,
100032/*Bayfront*/,
100213/*Blake*/,
100039/*Broward*/,
100086/*Broward North*/,
100161/*Central Regional*/,
100258/*Delray*/,
100223/*Ft. Walton Beach*/,
100017/*Halifax*/,
100019/*Holmes*/,
100022/*Jackson*/,
100208/*Jackson South*/,
100250/*All Childrens*/,
100209/*Kendall*/,
100157/*Lakeland*/,
100246/*Lawnwood*/,
100012/*Lee*/,
100179/*Memorial Jacksonville*/,
100038/*Memorial Regional*/,
110199/*Nicklaus Children*/,
100212/*Ocala*/,
100006/*Orlando*/,
100110/*Osceola*/,
100256/*Bayonet Point*/,
100113/*Shands*/,
100025/*Sacred Heart*/,
100087/*Sarasota*/,
100075/*St Joseph*/,
100010/*St Marys*/,
100135/*TMH*/,
100128/*Tampa*/
) THEN Trauma_Center="Yes"; ELSE Trauma_Center="No";

if faclnbr in (
100001,
100002,
100004,
100007,
100008,
100009,
100014,
100015,
100018,
100020,
100021,
100023,
100024,
100028,
100029,
100030,
100034,
100035,
100040,
100043,
100044,
100045,
100046,
100047,
100048,
100049,
100050,
100051,
100052,
100053,
100054,
100055,
100056,
100057,
100061,
100062,
100063,
100067,
100069,
100070,
100071,
100072,
100073,
100074,
100077,
100078,
100080,
100081,
100084,
100088,
100089,
100092,
100098,
100099,
100102,
100103,
100105,
100106,
100107,
100109,
100112,
100114,
100117,
100118,
100121,
100122,
100124,
100125,
100126,
100127,
100130,
100132,
100134,
100137,
100140,
100142,
100147,
100151,
100153,
100154,
100156,
100159,
100160,
100162,
100165,
100166,
100167,
100168,
100169,
100173,
100175,
100176,
100177,
100180,
100181,
100183,
100187,
100189,
100191,
100195,
100200,
100204,
100206,
100210,
100211,
100217,
100219,
100220,
100224,
100225,
100226,
100228,
100230,
100231,
100232,
100234,
100236,
100238,
100240,
100241,
100242,
100243,
100244,
100248,
100249,
100252,
100253,
100254,
100255,
100259,
100260,
100263,
100264,
100282,
104026,
110001,
110003,
110004,
110006,
110008,
110010,
110019,
110044,
110049,
110403,
111525,
111527,
120001,
120002,
120003,
120004,
120005,
120006,
120007,
120009,
120013,
23960017,
23960025,
23960032,
23960034,
23960039,
23960041,
23960046,
23960050,
23960051,
23960052,
23960057,
23960064,
23960084,
23960088,
23960092,
23960096,
23960098,
23960099,
23960100,
23960105,
23960108,
23960111,
23960117,
23960121,
23960123,
23960125,
23960157
)
then acute_care='Yes';
else acute_care='No';

RUN; 

/***Produce tables***/
ods rtf;
proc freq data=work.inpt; where Severity='Severe';
	tables TSA/list nopercent nocum;
title 'All severely injured patients by TSA 2019';
run;
quit;

title 'Severely Injured treated at acute care by TSA, 2019';
proc freq data=work.inpt;
	where Severity='Severe' and acute_care='Yes';
	tables TSA/list nopercent nocum;
run;
quit;
ods rtf close;