* goods movement, ssd - homepage | california air ... · pdf file/* goods movement, ssd cargo...
If you can't read please download the document
TRANSCRIPT
/* goods movement, ssdcargo handling equipmentto estimate health endpoints for POLA & POLB; results s/ be
divided by 0.70 to get SW.file cargo_mortality_20051005.sas 10/05/2005 h.t */
/* PART 1: 2004 Baseline */
PROC IMPORT OUT= WORK.conc DATAFILE= "X:\!! RD - Current Staff\HienTran\Goods Movt\Cargo Handling\Cargo-2004-baseline.dbf" DBMS=DBF REPLACE; GETDELETED=NO;RUN;
data conc;set conc;rename pla_plb_20=conc;run;
/* endpoint 1: mortality, central */data mort;
set conc;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do; brate=0.0102395962771332; /* orange */factor=0.5539;end;
beta=0.0046257;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total sum=pop effect;run;
data total1;set total;endpoint='mortality';run;
/* endpoint 1: mortality, low */data mort;
set conc;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do; brate=0.0102395962771332; /* orange */
factor=0.5539;end;
beta=0.002264684;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total1_low sum=pop effect;run;
/* endpoint 1: mortality, high */data mort;
set conc;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do; brate=0.0102395962771332; /* orange */factor=0.5539;end;
beta=0.006986716;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total1_hi sum=pop effect;run;
/* endpoint 2: asthma attacks, central */data asthma;
set conc;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.001440;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;run;
proc univariate data=asthma noprint;var population effect;output out=total2 sum=pop effect;run;
/* endpoint 2: asthma attacks, low */data asthma;
set conc;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.000350;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;
run;
proc univariate data=asthma noprint;var population effect;output out=total2_low sum=pop effect;run;
/* endpoint 2: asthma attacks, high */data asthma;
set conc;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.002530;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;run;
proc univariate data=asthma noprint;var population effect;output out=total2_hi sum=pop effect;run;
/* endpoint 3: work loss days, central */data wld;
set conc;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=2.3652;factor=0.6313; /* orange */end;
beta=0.004600;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total3 sum=pop effect;run;
/* endpoint 3: work loss days, low */data wld;
set conc;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=2.3652;factor=0.6313; /* orange */
end;beta=0.003894;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total3_low sum=pop effect;run;
/* endpoint 3: work loss days, high */data wld;
set conc;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=2.3652;factor=0.6313; /* orange */end;
beta=0.005306;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total3_hi sum=pop effect;run;
/* endpoint 4: mrad, central */data wld;
set conc;if county_id=19 then do;
brate=7.80005;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=7.80005;factor=0.6313; /* orange */end;
beta=0.007410;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total4 sum=pop effect;run;
/* endpoint 4: mrad, low */data wld;
set conc;if county_id=19 then do;
brate=7.80005;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=7.80005;factor=0.6313; /* orange */end;
beta=0.006038;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total4_low sum=pop effect;run;
/* endpoint 4: mrad, high */data wld;
set conc;if county_id=19 then do;
brate=7.80005;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=7.80005;factor=0.6313; /* orange */end;
beta=0.008782;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total4_hi sum=pop effect;run;
/* PART 2: 2007-2020 Regulation, Annual Average Effect */
PROC IMPORT OUT= WORK.conc_r DATAFILE= "X:\!! RD - Current Staff\HienTran\Goods Movt\Cargo Handling\Cargo-2007-2020-change.dbf" DBMS=DBF REPLACE; GETDELETED=NO;RUN;
data conc_r;set conc_r;conc=0;conc=adjusted_2;
run;
/* endpoint 1: mortality, central */data mort;
set conc_r;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do; brate=0.0102395962771332; /* orange */factor=0.5539;end;
beta=0.0046257;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total sum=pop effect;run;
data total1;set total;endpoint='mortality';run;
/* endpoint 1: mortality, low */data mort;
set conc_r;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do; brate=0.0102395962771332; /* orange */factor=0.5539;end;
beta=0.002264684;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total1_low sum=pop effect;run;
/* endpoint 1: mortality, high */data mort;
set conc_r;if county_id=19 then do;
brate=0.0121491919830609; /* los angeles */factor=0.5336 ;/* portion of age 30+ over all ages */end;
else if county_id=30 then do;
brate=0.0102395962771332; /* orange */factor=0.5539;end;
beta=0.006986716;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=mort noprint;var population effect;output out=total1_hi sum=pop effect;run;
/* endpoint 2: asthma attacks, central */data asthma;
set conc_r;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.001440;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;run;
proc univariate data=asthma noprint;var population effect;output out=total2 sum=pop effect;run;
/* endpoint 2: asthma attacks, low */data asthma;
set conc_r;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.000350;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;run;
proc univariate data=asthma noprint;var population effect;output out=total2_low sum=pop effect;run;
/* endpoint 2: asthma attacks, high */data asthma;
set conc_r;brate=0.027;factor=0.0561 ;/* portion of asthmatics over all ages */beta=0.002530;effect=(-1)*((brate /((1-brate)*exp(beta*conc)+brate)) - brate) *
population*factor*365;run;
proc univariate data=asthma noprint;var population effect;
output out=total2_hi sum=pop effect;run;
/* endpoint 3: work loss days, central */data wld;
set conc_r;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=2.3652;factor=0.6313; /* orange */end;
beta=0.004600;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total3 sum=pop effect;run;
/* endpoint 3: work loss days, low */data wld;
set conc_r;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles */end;
else if county_id=30 then do; brate=2.3652;factor=0.6313; /* orange */end;
beta=0.003894;effect= (-1) * (brate * (exp(-beta*conc) -1)) * population*factor;run;
proc univariate data=wld noprint;var population effect;output out=total3_low sum=pop effect;run;
/* endpoint 3: work loss days, high */data wld;
set conc_r;if county_id=19 then do;
brate=2.3652;factor=0.6218 ;/* portion of age 18-65 over all ages, los
angeles