Document ID: EPA-HQ-OW-2002-0033-0216
Agency: epa
Document Type: Supporting & Related Material
Title: 
Posted Date: 2003-04-14T04:00Z

130
"
Lead
Benefits
(
Sub
12).
sas"
Using
data
from
the
lead
specific
data
set,
changes
in
blood
lead
levels
are
estimated
and
benefits
calculated
for
reductions
in
cases
of
hypertension,
CHD,
brain
infarctions,
and
mortality
for
recreational
fishermen
and
their
families.
131
/*
Program:
Lead
Benefits
(
Sub
12).
sas
Author:
Stuart
Smith,
Abt
Associates,
Inc.
Date:
07/
05/
2000
Updated:
11/
08/
2002
This
program
estimates
benefits
from
reduced
consumption
of
fish
tissue
containing
lead.
Reductions
in
MP&
M
facility
lead
loadings
result
in
lower
probabilities
of
lead
related
illnesses.
*/

title2
'
Lead
Related
Benefits
Recreational
(
Sub
12)';

data
lead;
set
rundat.
lead;

/*
Combine
the
exposed
populations
recreational
populations,
fresh
and
saltwater
*/

exppop
=
fwsptpop
+
swsptpop;

/*
Calculate
the
blood
lead
levels
under
the
baseline
and
post
compliance
scenarios
*/

PbB0m
=
4.55;
/*
Typical
PbB
conc.
in
absence
of
fish
cons.
for
men
*/
PbB0f
=
3.45;
/*
Typical
PbB
conc.
in
absence
of
fish
cons.
for
women
*/
bcf
=
49;
/*
Bioconcentration
factor
for
lead
L/
KG
*/
bksf
=
0.4;
/*
Biokinetic
slope
factor
for
lead
*/
ingm
=
17.5;
/*
Average
daily
fish
consumption
(
g/
day)
males
age
15­
44
*/
ingf
=
17.5;
/*
Average
daily
fish
consumption
(
g/
day)
females
age
15­
44
*/
ingm2
=
17.5;
/*
Average
daily
fish
consumption
(
g/
day)
males
over
age
45
*/
ingf2
=
17.5;
/*
Average
daily
fish
consumption
(
g/
day)
females
over
age
45
*/
af
=
0.03;
/*
Absorbtion
factor
*/
ef
=
365;
/*
Exposure
Frequency
*/
at
=
365;
/*
Average
timing,
or
duration
of
exposure
in
days
*/
cf
=
0.001;
/*
Conversion
factor
kg/
g
*/

/*
Uses
ingestion
for
males
ages
15­
44
*/
PbB_
bm
=
PbB0m
+
((
cnc_
base*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);
PbB_
1m
=
PbB0m
+
((
cnc_
trt1*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);
PbB_
2m
=
PbB0m
+
((
cnc_
trt2*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);

/*
Uses
ingestion
for
males
ages
45
and
older
*/
PbB_
bme
=
PbB0m
+
((
cnc_
base*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);
PbB_
1me
=
PbB0m
+
((
cnc_
trt1*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);
PbB_
2me
=
PbB0m
+
((
cnc_
trt2*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);

/*
Uses
ingestion
for
females
ages
15­
44
*/
PbB_
bf
=
PbB0f
+
((
cnc_
base*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);
PbB_
1f
=
PbB0f
+
((
cnc_
trt1*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);
PbB_
2f
=
PbB0f
+
((
cnc_
trt2*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);

/*
Uses
ingestion
for
females
ages
45
and
older
*/
PbB_
bfe
=
PbB0f
+
((
cnc_
base*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);
132
PbB_
1fe
=
PbB0f
+
((
cnc_
trt1*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);
PbB_
2fe
=
PbB0f
+
((
cnc_
trt2*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);

drop
fwsptpop
fwsubpop
swsptpop
swsubpop;
run;

data
men;
set
lead;

/*
Set
up
some
constants
*/
/*
All
percentages
based
on
the
1996
Census
numbers
*/

pctmale
=
0.319518;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
20­
74
*/
pctm4054
=
0.098405;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
40­
54
*/
pctm4059
=
0.118945;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
40­
59
*/
pctm4574
=
0.12933;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
45­
74
*/
pctm5564
=
0.03829;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
55­
64
*/
pctm6064
=
0.01775;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
60­
64
*/
pctm6574
=
0.031427;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
65­
74
*/

bbpm4054
=
79.9;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
40­
54
*/
bbpm4059
=
80.0;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
40­
59
*/
bbpm4574
=
78.8;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
45­
74
*/
bbpm5564
=
79.0;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
55­
64
*/
bbpm6064
=
77.8;
/*
Mean
diastolic
bp
from
NHANES
III,
males
age
60­
64
*/
bbpm6574
=
76.4;
/*
Mean
diastolic
bp
from
NHANES
III,
males
age
65­
74
*/

/**************************************************************************/
/*
Begin
Benefits
for
Males
*/
/**************************************************************************/

/*
Calculate
the
change
in
the
probability
of
hypertension
*/

pr_
hyp1
=
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
bm))))
­
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
1m))));

pr_
hyp2
=
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
bm))))
­
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
2m))));

ben_
hyp1
=
pr_
hyp1
*
exppop
*
pctmale
*
1141;
ben_
hyp2
=
pr_
hyp2
*
exppop
*
pctmale
*
1141;
hyppop
=
exppop
*
pctmale;
hyp_
red1
=
pr_
hyp1
*
exppop
*
pctmale;
hyp_
red2
=
pr_
hyp2
*
exppop
*
pctmale;

/*
Calculate
Benefits
from
Changes
in
CHD,
non­
fatal
cases
only,
which
are
2/
3
of
the
cases
*/

dbpm1
=
1.4
*
log(
PbB_
bm/
PbB_
1m);
dbpm2
=
1.4
*
log(
PbB_
bm/
PbB_
2m);

dbpm1e
=
1.4
*
log(
PbB_
bme/
PbB_
1me);
dbpm2e
=
1.4
*
log(
PbB_
bme/
PbB_
2me);
133
/*
Ages
40­
59,
change
in
the
10
year
prob
of
occ.
*/
d1m4059
=
bbpm4059
­
dbpm1;
d2m4059
=
bbpm4059
­
dbpm2;

pr1_
4059
=
(
1/(
1
+
exp(
4.996
­
0.030365*
bbpm4059)))
­
(
1/(
1
+
exp(
4.996
­
0.030365*
d1m4059)));

pr2_
4059
=
(
1/(
1
+
exp(
4.996
­
0.030365*
bbpm4059)))
­
(
1/(
1
+
exp(
4.996
­
0.030365*
d2m4059)));

/*
Ages
60­
64,
change
in
the
2
year
prob
of
occ.
*/
d1m6064
=
bbpm6064
­
dbpm1e;
d2m6064
=
bbpm6064
­
dbpm2e;

pr1_
6064
=
(
1/(
1
+
exp(
5.19676
­
0.02351*
bbpm6064)))
­
(
1/(
1
+
exp(
5.19676
­
0.02351*
d1m6064)));

pr2_
6064
=
(
1/(
1
+
exp(
5.19676
­
0.02351*
bbpm6064)))
­
(
1/(
1
+
exp(
5.19676
­
0.02351*
d2m6064)));

/*
Ages
65­
74,
change
in
the
2
year
prob
of
occ.
*/
d1m6574
=
bbpm6574
­
dbpm1e;
d2m6574
=
bbpm6574
­
dbpm2e;

pr1_
6574
=
(
1/(
1
+
exp(
4.90723
­
0.02031*
bbpm6574)))
­
(
1/(
1
+
exp(
4.90723
­
0.02031*
d1m6574)));

pr2_
6574
=
(
1/(
1
+
exp(
4.90723
­
0.02031*
bbpm6574)))
­
(
1/(
1
+
exp(
4.90723
­
0.02031*
d2m6574)));

b14059
=
(
pr1_
4059
*
exppop
*
pctm4059
*
76347
*
0.66)/
10;
b24059
=
(
pr2_
4059
*
exppop
*
pctm4059
*
76347
*
0.66)/
10;
b16064
=
(
pr1_
6064
*
exppop
*
pctm6064
*
76347
*
0.66)/
2;
b26064
=
(
pr2_
6064
*
exppop
*
pctm6064
*
76347
*
0.66)/
2;
b16574
=
(
pr1_
6574
*
exppop
*
pctm6574
*
76347
*
0.66)/
2;
b26574
=
(
pr2_
6574
*
exppop
*
pctm6574
*
76347
*
0.66)/
2;

/*
Number
of
reduced
cases
Multiply
by
the
2/
3
(.
66)
and
divide
by
the
years
of
occurance
associated
with
the
probability
to
get
cases
for
one
year.
*/

r14059
=
(
pr1_
4059
*
exppop
*
pctm4059
*
0.66)/
10;
r24059
=
(
pr2_
4059
*
exppop
*
pctm4059
*
0.66)/
10;
r16064
=
(
pr1_
6064
*
exppop
*
pctm6064
*
0.66)/
2;
r26064
=
(
pr2_
6064
*
exppop
*
pctm6064
*
0.66)/
2;
r16574
=
(
pr1_
6574
*
exppop
*
pctm6574
*
0.66)/
2;
r26574
=
(
pr2_
6574
*
exppop
*
pctm6574
*
0.66)/
2;

/*
Calculate
the
exposed
populations
at
each
age
group
*/
popm4059
=
exppop
*
pctm4059;
popm6064
=
exppop
*
pctm6064;
134
popm6574
=
exppop
*
pctm6574;

/*
Calculate
Benefits
from
Reduced
Cerebrovascular
Accidents
and
Brain
Infarctions,
affects
men
aged
45­
74
Non­
fatal
only,
which
is
70%
of
all
cases
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/
d1m4574
=
bbpm4574
­
dbpm1e;
d2m4574
=
bbpm4574
­
dbpm2e;

pr1_
cba
=
(
1/(
1
+
exp(
8.58889
­
0.04066*
bbpm4574)))
­
(
1/(
1
+
exp(
8.58889
­
0.04066*
d1m4574)));

pr2_
cba
=
(
1/(
1
+
exp(
8.58889
­
0.04066*
bbpm4574)))
­
(
1/(
1
+
exp(
8.58889
­
0.04066*
d2m4574)));

pr1_
bi
=
(
1/(
1
+
exp(
9.9516
­
0.04840*
bbpm4574)))
­
(
1/(
1
+
exp(
9.9516
­
0.04840*
d1m4574)));

pr2_
bi
=
(
1/(
1
+
exp(
9.9516
­
0.04840*
bbpm4574)))
­
(
1/(
1
+
exp(
9.9516
­
0.04840*
d2m4574)));

/*
Number
of
reduced
cases
Multiply
by
the
70%
(.
7)
and
divide
by
the
years
of
occurance
associated
with
the
probability
to
get
cases
for
one
year.
*/

b1_
cba
=
(
pr1_
cba
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b2_
cba
=
(
pr2_
cba
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b1_
bi
=
(
pr1_
bi
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b2_
bi
=
(
pr2_
bi
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;

r1_
cba
=
(
pr1_
cba
*
exppop
*
pctm4574
*
0.7)/
2;
r2_
cba
=
(
pr2_
cba
*
exppop
*
pctm4574
*
0.7)/
2;
r1_
bi
=
(
pr1_
bi
*
exppop
*
pctm4574
*
0.7)/
2;
r2_
bi
=
(
pr2_
bi
*
exppop
*
pctm4574
*
0.7)/
2;

/*
Calculate
the
exposed
population*/
popm4574
=
exppop
*
pctm4574;

/*
Calculate
Benefits
from
Premature
Mortality
*/

/*
Ages
40­
54,
change
in
the
12
year
prob
of
death
*/
d1m4054
=
bbpm4054
­
dbpm1;
d2m4054
=
bbpm4054
­
dbpm2;

m1_
4054
=
(
1/(
1
+
exp(
5.3158
­
0.03516*
bbpm4054)))
­
(
1/(
1
+
exp(
5.3158
­
0.03516*
d1m4054)));
135
m2_
4054
=
(
1/(
1
+
exp(
5.3158
­
0.03516*
bbpm4054)))
­
(
1/(
1
+
exp(
5.3158
­
0.03516*
d2m4054)));

/*
Ages
55­
64,
change
in
the
2
year
prob
of
death
*/
d1m5564
=
bbpm5564
­
dbpm1e;
d2m5564
=
bbpm5564
­
dbpm2e;

m1_
5564
=
(
1/(
1
+
exp(
4.89528
­
0.01866*
bbpm5564)))
­
(
1/(
1
+
exp(
4.89528
­
0.01866*
d1m5564)));

m2_
5564
=
(
1/(
1
+
exp(
4.89528
­
0.01866*
bbpm5564)))
­
(
1/(
1
+
exp(
4.89528
­
0.01866*
d2m5564)));

/*
Ages
65­
74,
change
in
the
2
year
prob
of
death
*/
d1m6574
=
bbpm6574
­
dbpm1e;
d2m6574
=
bbpm6574
­
dbpm2e;

m1_
6574
=
(
1/(
1
+
exp(
3.05723
­
0.00547*
bbpm6574)))
­
(
1/(
1
+
exp(
3.05723
­
0.00547*
d1m6574)));

m2_
6574
=
(
1/(
1
+
exp(
3.05723
­
0.00547*
bbpm6574)))
­
(
1/(
1
+
exp(
3.05723
­
0.00547*
d2m6574)));

/*
Number
of
reduced
cases
and
Benefits
Divide
by
the
years
of
occurance
associated
with
the
probability
to
get
cases
for
one
year.
*/

m14054l
=
(
m1_
4054
*
exppop
*
pctm4054
*
2500000)/
12;
m24054l
=
(
m2_
4054
*
exppop
*
pctm4054
*
2500000)/
12;
m14054h
=
(
m1_
4054
*
exppop
*
pctm4054
*
13600000)/
12;
m24054h
=
(
m2_
4054
*
exppop
*
pctm4054
*
13600000)/
12;
mid40541
=
(
m1_
4054
*
exppop
*
pctm4054
*
6500000)/
12;
mid40542
=
(
m2_
4054
*
exppop
*
pctm4054
*
6500000)/
12;

m15564l
=
(
m1_
5564
*
exppop
*
pctm5564
*
2500000)/
2;
m25564l
=
(
m2_
5564
*
exppop
*
pctm5564
*
2500000)/
2;
m15564h
=
(
m1_
5564
*
exppop
*
pctm5564
*
13600000)/
2;
m25564h
=
(
m2_
5564
*
exppop
*
pctm5564
*
13600000)/
2;
mid55641
=
(
m1_
5564
*
exppop
*
pctm5564
*
6500000)/
2;
mid55642
=
(
m2_
5564
*
exppop
*
pctm5564
*
6500000)/
2;

m16574l
=
(
m1_
6574
*
exppop
*
pctm6574
*
2500000)/
2;
m26574l
=
(
m2_
6574
*
exppop
*
pctm6574
*
2500000)/
2;
m16574h
=
(
m1_
6574
*
exppop
*
pctm6574
*
13600000)/
2;
m26574h
=
(
m2_
6574
*
exppop
*
pctm6574
*
13600000)/
2;
mid65741
=
(
m1_
6574
*
exppop
*
pctm6574
*
6500000)/
2;
mid65742
=
(
m2_
6574
*
exppop
*
pctm6574
*
6500000)/
2;

r14054l
=
(
m1_
4054
*
exppop
*
pctm4054)/
12;
r24054l
=
(
m2_
4054
*
exppop
*
pctm4054)/
12;
r14054h
=
(
m1_
4054
*
exppop
*
pctm4054)/
12;
r24054h
=
(
m2_
4054
*
exppop
*
pctm4054)/
12;
136
r15564l
=
(
m1_
5564
*
exppop
*
pctm5564)/
2;
r25564l
=
(
m2_
5564
*
exppop
*
pctm5564)/
2;
r15564h
=
(
m1_
5564
*
exppop
*
pctm5564)/
2;
r25564h
=
(
m2_
5564
*
exppop
*
pctm5564)/
2;
r16574l
=
(
m1_
6574
*
exppop
*
pctm6574)/
2;
r26574l
=
(
m2_
6574
*
exppop
*
pctm6574)/
2;
r16574h
=
(
m1_
6574
*
exppop
*
pctm6574)/
2;
r26574h
=
(
m2_
6574
*
exppop
*
pctm6574)/
2;

/*
Calculate
the
exposed
populations
*/
popm4054
=
exppop
*
pctm4054;
popm5564
=
exppop
*
pctm5564;
run;

data
women;
set
lead;

/*
Set
up
some
constants
*/
/*
All
percentages
based
on
the
1996
Census
numbers
*/
pctfem
=
0.511735;
/*
Percentage
of
U.
S.
Pop.
who
are
female
*/
pctf4574
=
0.143282;
/*
Percentage
of
U.
S.
Pop.
who
are
female,
age
45­
74
*/

bbpf4574
=
74.8;
/*
Mean
diastolic
bp
from
NHANES
III,
females
age
45­
74
*/

/**************************************************************************/
/*
Begin
Benefits
for
Females
*/
/**************************************************************************/

/*
No
dose­
response
function
for
female
hypertension
*/

/*
Calculate
Benefits
from
Changes
in
CHD,
non­
fatal
cases
only
which
are
2/
3
of
the
cases
*/

dbpf1
=
(
0.6
*
1.4)
*
log(
PbB_
bfe/
PbB_
1fe);
dbpf2
=
(
0.6
*
1.4)
*
log(
PbB_
bfe/
PbB_
2fe);

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/
d1f4574
=
bbpf4574
­
dbpf1;
d2f4574
=
bbpf4574
­
dbpf2;

pr1_
4574
=
(
1/(
1
+
exp(
6.9401
­
0.03072*
bbpf4574)))
­
(
1/(
1
+
exp(
6.9401
­
0.03072*
d1f4574)));

pr2_
4574
=
(
1/(
1
+
exp(
6.9401
­
0.03072*
bbpf4574)))
­
(
1/(
1
+
exp(
6.9401
­
0.03072*
d2f4574)));

b1fchd
=
(
pr1_
4574
*
exppop
*
pctf4574
*
76347
*
0.66)/
2;
b2fchd
=
(
pr2_
4574
*
exppop
*
pctf4574
*
76347
*
0.66)/
2;

red_
chd1
=
(
pr1_
4574
*
exppop
*
pctf4574
*
0.66)/
2;
red_
chd2
=
(
pr2_
4574
*
exppop
*
pctf4574
*
0.66)/
2;
137
popf4574
=
exppop
*
pctf4574;

/*
Calculate
Benefits
from
Reduced
Cerebrovascular
Accidents
and
Brain
Infarctions,
affects
women
aged
45­
74
Non­
fatal
only
which
is
70%
of
all
cases
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/

pr1f_
bi
=
(
1/(
1
+
exp(
10.6716
­
0.0544*
bbpf4574)))
­
(
1/(
1
+
exp(
10.6716
­
0.0544*
d1f4574)));

pr2f_
bi
=
(
1/(
1
+
exp(
10.6716
­
0.0544*
bbpf4574)))
­
(
1/(
1
+
exp(
10.6716
­
0.0544*
d2f4574)));

pr1f_
cba
=
(
1/(
1
+
exp(
9.07737
­
0.04287*
bbpf4574)))
­
(
1/(
1
+
exp(
9.07737
­
0.04287*
d1f4574)));

pr2f_
cba
=
(
1/(
1
+
exp(
9.07737
­
0.04287*
bbpf4574)))
­
(
1/(
1
+
exp(
9.07737
­
0.04287*
d2f4574)));

b1f_
cba
=
(
pr1f_
cba
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b2f_
cba
=
(
pr2f_
cba
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b1f_
bi
=
(
pr1f_
bi
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b2f_
bi
=
(
pr2f_
bi
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;

r1f_
cba
=
(
pr1f_
cba
*
exppop
*
pctf4574
*
0.7)/
2;
r2f_
cba
=
(
pr2f_
cba
*
exppop
*
pctf4574
*
0.7)/
2;
r1f_
bi
=
(
pr1f_
bi
*
exppop
*
pctf4574
*
0.7)/
2;
r2f_
bi
=
(
pr2f_
bi
*
exppop
*
pctf4574
*
0.7)/
2;

/*
Calculate
Benefits
from
Premature
Mortality
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
death
*/

m1_
4574
=
(
1/(
1
+
exp(
5.40374
­
0.01511*
bbpf4574)))
­
(
1/(
1
+
exp(
5.40374
­
0.01511*
d1f4574)));

m2_
4574
=
(
1/(
1
+
exp(
5.40374
­
0.01511*
bbpf4574)))
­
(
1/(
1
+
exp(
5.40374
­
0.01511*
d2f4574)));

m14574l
=
(
m1_
4574
*
exppop
*
pctf4574
*
2500000)/
2;
m24574l
=
(
m2_
4574
*
exppop
*
pctf4574
*
2500000)/
2;
m14574h
=
(
m1_
4574
*
exppop
*
pctf4574
*
13600000)/
2;
m24574h
=
(
m2_
4574
*
exppop
*
pctf4574
*
13600000)/
2;
mid45741
=
(
m1_
4574
*
exppop
*
pctf4574
*
6500000)/
2;
mid45742
=
(
m2_
4574
*
exppop
*
pctf4574
*
6500000)/
2;

r14574l
=
(
m1_
4574
*
exppop
*
pctf4574)/
2;
r24574l
=
(
m2_
4574
*
exppop
*
pctf4574)/
2;
r14574h
=
(
m1_
4574
*
exppop
*
pctf4574)/
2;
r24574h
=
(
m2_
4574
*
exppop
*
pctf4574)/
2;
run;
138
/*
Sum
benefits
and
cases
and
print
out
results
*/

proc
means
sum
data=
men;
var
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
output
out
=
rundat.
msptsamp(
drop
=
_
type_
_
freq_)
sum
=
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
title3
'
Mens
Lead
Related
Benefits
(
Sample
Basis)';
quit;

proc
means
sum
data=
men;
weight
event_
wt;
var
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
output
out
=
rundat.
msptwght(
drop
=
_
type_
_
freq_)
sum
=
139
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
title3
'
Mens
Lead
Related
Benefits
(
National
Basis)';
quit;

%
macro
maleleadtable(
tbl_
in=,
tbl_
out=);
data
rundat.&
tbl_
out
(
keep=
gender
categ
agecat
type
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h);
set
rundat.&
tbl_
in;

length
gender
$
50
categ
$
50
agecat
$
50
type
$
50;

array
_
hyper[
4]
ben_
hyp1
ben_
hyp2
hyp_
red1
hyp_
red2;
array
_
chd[
12]
b14059
b24059
b16064
b26064
b16574
b26574
r14059
r24059
r16064
r26064
r16574
r26574;
array
_
cba[
4]
b1_
cba
b2_
cba
r1_
cba
r2_
cba;
array
_
bi[
4]
b1_
bi
b2_
bi
r1_
bi
r2_
bi;
array
_
mort[
24]
m14054l
mid40541
m14054h
m24054l
mid40542
m24054h
m15564l
mid55641
m15564h
m25564l
mid55642
m25564h
m16574l
mid65741
m16574h
m26574l
mid65742
m26574h
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
array
ben[
6]
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
array
cases[
2]
case1
case2;

do
i
=
1
to
9;
if
i
=
1
then
do;
ben[
1]
=
_
hyper[
1];
ben[
2]
=
_
hyper[
1];
ben[
3]
=
_
hyper[
1];
ben[
4]
=
_
hyper[
2];
ben[
5]
=
_
hyper[
2];
ben[
6]
=
_
hyper[
2];
cases[
1]
=
_
hyper[
3];
cases[
2]
=
_
hyper[
4];
categ
=
"
Hypertension";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
2
then
do;
ben[
1]
=
_
chd[
1];
ben[
2]
=
_
chd[
1];
ben[
3]
=
_
chd[
1];
ben[
4]
=
_
chd[
2];
ben[
5]
=
_
chd[
2];
ben[
6]
=
_
chd[
2];
cases[
1]
=
_
chd[
7];
cases[
2]
=
_
chd[
8];
140
categ
=
"
CHD";
agecat
=
"
Age
40­
59";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
3
then
do;
ben[
1]
=
_
chd[
3];
ben[
2]
=
_
chd[
3];
ben[
3]
=
_
chd[
3];
ben[
4]
=
_
chd[
4];
ben[
4]
=
_
chd[
4];
ben[
6]
=
_
chd[
4];
cases[
1]
=
_
chd[
9];
cases[
2]
=
_
chd[
10];
categ
=
"
CHD";
agecat
=
"
Age
60­
64";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
4
then
do;
ben[
1]
=
_
chd[
5];
ben[
2]
=
_
chd[
5];
ben[
3]
=
_
chd[
5];
ben[
4]
=
_
chd[
6];
ben[
4]
=
_
chd[
6];
ben[
6]
=
_
chd[
6];
cases[
1]
=
_
chd[
11];
cases[
2]
=
_
chd[
12];
categ
=
"
CHD";
agecat
=
"
Age
65­
74";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
5
then
do;
ben[
1]
=
_
cba[
1];
ben[
2]
=
_
cba[
1];
ben[
3]
=
_
cba[
1];
ben[
4]
=
_
cba[
2];
ben[
5]
=
_
cba[
2];
ben[
6]
=
_
cba[
2];
cases[
1]
=
_
cba[
3];
cases[
2]
=
_
cba[
4];
categ
=
"
CBA";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
6
then
do;
ben[
1]
=
_
bi[
1];
ben[
2]
=
_
bi[
1];
ben[
3]
=
_
bi[
1];
ben[
4]
=
_
bi[
2];
ben[
5]
=
_
bi[
2];
ben[
6]
=
_
bi[
2];
cases[
1]
=
_
bi[
3];
cases[
2]
=
_
bi[
4];
categ
=
"
BI";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
7
then
do;
ben[
1]
=
_
mort[
1];
ben[
2]
=
_
mort[
2];
ben[
3]
=
_
mort[
3];
ben[
4]
=
_
mort[
4];
ben[
5]
=
_
mort[
5];
ben[
6]
=
_
mort[
6];
cases[
1]
=
_
mort[
19];
cases[
2]
=
_
mort[
20];
categ
=
"
Mortality";
agecat
=
"
Age
40­
54";
type
=
"
Sport";
141
gender
=
"
Male";
output;
end;
if
i
=
8
then
do;
ben[
1]
=
_
mort[
7];
ben[
2]
=
_
mort[
8];
ben[
3]
=
_
mort[
9];
ben[
4]
=
_
mort[
10];
ben[
5]
=
_
mort[
11];
ben[
6]
=
_
mort[
12];
cases[
1]
=
_
mort[
21];
cases[
2]
=
_
mort[
22];
categ
=
"
Mortality";
agecat
=
"
Age
55­
64";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
if
i
=
9
then
do;
ben[
1]
=
_
mort[
13];
ben[
2]
=
_
mort[
14];
ben[
3]
=
_
mort[
15];
ben[
4]
=
_
mort[
16];
ben[
5]
=
_
mort[
17];
ben[
6]
=
_
mort[
18];
cases[
1]
=
_
mort[
23];
cases[
2]
=
_
mort[
24];
categ
=
"
Mortality";
agecat
=
"
Age
65­
74";
type
=
"
Sport";
gender
=
"
Male";
output;
end;
end;
run;
%
mend;

%
maleleadtable(
tbl_
in=
msptsamp,
tbl_
out=
tblmssmp);
%
maleleadtable(
tbl_
in=
msptwght,
tbl_
out=
tblmswgt);

proc
means
sum
data=
women;
var
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
output
out
=
rundat.
fsptsamp(
drop
=
_
type_
_
freq_)
sum
=
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
title3
'
Female
Lead
Related
Benefits
(
Sample
Basis)';
quit;

proc
means
sum
data=
women;
weight
event_
wt;
var
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
142
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
output
out
=
rundat.
fsptwght(
drop
=
_
type_
_
freq_)
sum
=
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
title3
'
Female
Lead
Related
Benefits
(
National
Basis)';
quit;

%
macro
femaleleadtable(
tbl_
in=,
tbl_
out=);
data
rundat.&
tbl_
out
(
keep=
gender
categ
agecat
type
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h);
set
rundat.&
tbl_
in;

length
gender
$
50
categ
$
50
agecat
$
50
type
$
50;

array
_
chd[
4]
b1fchd
b2fchd
red_
chd1
red_
chd2;
array
_
cba[
4]
b1f_
cba
b2f_
cba
r1f_
cba
r2f_
cba;
array
_
bi[
4]
b1f_
bi
b2f_
bi
r1f_
bi
r2f_
bi;
array
_
mort[
8]
m14574l
mid45741
m14574h
m24574l
mid45742
m24574h
r14574l
r24574l;
array
ben[
6]
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
array
cases[
2]
case1
case2;

do
i
=
1
to
4;
if
i
=
1
then
do;
ben[
1]
=
_
chd[
1];
ben[
2]
=
_
chd[
1];
ben[
3]
=
_
chd[
1];
ben[
4]
=
_
chd[
2];
ben[
5]
=
_
chd[
2];
ben[
6]
=
_
chd[
2];
cases[
1]
=
_
chd[
3];
cases[
2]
=
_
chd[
4];
categ
=
"
CHD";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Female";
output;
end;
if
i
=
2
then
do;
ben[
1]
=
_
cba[
1];
ben[
2]
=
_
cba[
1];
ben[
3]
=
_
cba[
1];
ben[
4]
=
_
cba[
2];
ben[
5]
=
_
cba[
2];
ben[
6]
=
_
cba[
2];
cases[
1]
=
_
cba[
3];
cases[
2]
=
_
cba[
4];
categ
=
"
CBA";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Female";
output;
end;
if
i
=
3
then
do;
ben[
1]
=
_
bi[
1];
ben[
2]
=
_
bi[
1];
ben[
3]
=
_
bi[
1];
ben[
4]
=
_
bi[
2];
ben[
5]
=
_
bi[
2];
ben[
6]
=
_
bi[
2];
cases[
1]
=
_
bi[
3];
cases[
2]
=
_
bi[
4];
143
categ
=
"
BI";
agecat
=
"
All";
type
=
"
Sport";
gender
=
"
Female";
output;
end;
if
i
=
4
then
do;
ben[
1]
=
_
mort[
1];
ben[
2]
=
_
mort[
2];
ben[
3]
=
_
mort[
3];
ben[
4]
=
_
mort[
4];
ben[
5]
=
_
mort[
5];
ben[
6]
=
_
mort[
6];
cases[
1]
=
_
mort[
7];
cases[
2]
=
_
mort[
8];
categ
=
"
Mortality";
agecat
=
"
Age
45­
74";
type
=
"
Sport";
gender
=
"
Female";
output;
end;
end;
run;
%
mend;

%
femaleleadtable(
tbl_
in=
fsptsamp,
tbl_
out=
tblfssmp);
%
femaleleadtable(
tbl_
in=
fsptwght,
tbl_
out=
tblfswgt);

data
rundat.
tblspt_
s;
set
rundat.
tblmssmp
rundat.
tblfssmp;
run;

data
rundat.
tblspt_
w;
set
rundat.
tblmswgt
rundat.
tblfswgt;
run;

proc
means
sum
data=
men;
var
event_
wt;
title3
'
Number
of
Affected
Reaches
Nationally';
quit;

proc
sql;

create
table
malepop
as
select
sum(
exppop)
as
tpops,
sum(
exppop*
pctmale)
as
pm2074s,
sum(
exppop*
pctm4574)
as
pm4574s,
sum(
exppop*
event_
wt)
as
tpopn,
sum(
exppop*
pctmale*
event_
wt)
as
pm2074n,
sum(
exppop*
pctm4574*
event_
wt)
as
pm4574n
from
men;

create
table
fempop
as
select
sum(
exppop)
as
tpops,
sum(
exppop*
pctf4574)
as
pf4574s,
sum(
exppop*
event_
wt)
as
tpopn,
sum(
exppop*
pctf4574*
event_
wt)
as
pf4574n
from
women;
144
quit;

proc
print
data=
malepop;
var
tpops
pm2074s
pm4574s
tpopn
pm2074n
pm4574n;
label
tpops
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Smpl)'
pm2074s
=
'
Exposed
Males
Age
20­
74
(
Smpl)'
pm4574s
=
'
Exposed
Males
Age
45­
74
(
Smpl)'
tpopn
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Natl)'
pm2074n
=
'
Exposed
Males
Age
20­
74
(
Natl)'
pm4574n
=
'
Exposed
Males
Age
45­
74
(
Natl)';
title3
'
Exposed
Population
­
Males';
quit;

proc
print
data=
fempop;
var
tpops
pf4574s
tpopn
pf4574n;
label
tpops
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Smpl)'
pf4574s
=
'
Exposed
Females
Age
45­
74
(
Smpl)'
tpopn
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Natl)'
pf4574n
=
'
Exposed
Females
Age
45­
74
(
Natl)';
title3
'
Exposed
Population
­
Females';
quit;

/*
End
Lead
Benefits
(
Sub
12).
sas
*/
145
Program
Output
146
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
79
Lead
Related
Benefits
Recreational
(
Sub
12)
Mens
Lead
Related
Benefits
(
National
Basis)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Sum
                        
ben_
hyp1
44511.55
ben_
hyp2
64576.82
hyppop
133090.73
b14059
1787.60
b24059
2594.12
b16064
502.3133758
b26064
728.9498374
b16574
777.0629889
b26574
1127.67
popm4059
49544.87
popm6064
7393.51
popm6574
13090.48
b1_
cba
4156.82
b2_
cba
6032.22
b1_
bi
2342.05
b2_
bi
3398.68
popm4574
53870.59
m14054l
74537.83
m14054h
405485.81
m24054l
108167.53
m24054h
588431.38
m15564l
40590.06
m15564h
220809.92
m25564l
58903.91
m25564h
320437.27
m16574l
19827.98
m16574h
107864.19
m26574l
28774.46
m26574h
156533.09
mid40541
193798.36
mid40542
281235.59
mid55641
105534.15
mid55642
153150.17
mid65741
51552.74
mid65742
74813.61
popm4054
40989.22
popm5564
15949.16
hyp_
red1
39.0110017
hyp_
red2
56.5966874
r14059
0.0234141
r24059
0.0339781
r16064
0.0065793
r26064
0.0095479
r16574
0.0101780
147
r26574
0.0147703
r1_
cba
0.0124034
r2_
cba
0.0179994
r1_
bi
0.0069884
r2_
bi
0.0101412
r14054l
0.0298151
r24054l
0.0432670
r15564l
0.0162360
r25564l
0.0235616
r16574l
0.0079312
r26574l
0.0115098
                        
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
81
Lead
Related
Benefits
Recreational
(
Sub
12)
Female
Lead
Related
Benefits
(
National
Basis)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Sum
                        
b1fchd
1229.46
b2fchd
1784.50
popf4574
59682.10
b1f_
cba
1777.69
b2f_
cba
2580.20
b1f_
bi
1088.31
b2f_
bi
1579.60
m14574l
43015.20
m14574h
234002.69
m24574l
62434.83
m24574h
339645.45
mid45741
111839.52
mid45742
162330.55
red_
chd1
0.0161036
red_
chd2
0.0233735
r1f_
cba
0.0070726
r2f_
cba
0.0102653
r1f_
bi
0.0043298
r2f_
bi
0.0062844
r14574l
0.0172061
r14574h
0.0172061
r24574l
0.0249739
r24574h
0.0249739
                        
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
82
Lead
Related
Benefits
Recreational
(
Sub
12)
Number
of
Affected
Reaches
Nationally
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Analysis
Variable
:
EVENT_
WT
Event
Weight
148
Sum
            
23.7361201
            
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
83
Lead
Related
Benefits
Recreational
(
Sub
12)
Exposed
Population
­
Males
15:
45
Thursday,
February
6,
2003
Obs
tpops
pm2074s
pm4574s
tpopn
pm2074n
pm4574n
1
82069
26222.52
10613.98
416535.94
133090.73
53870.59
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
84
Lead
Related
Benefits
Recreational
(
Sub
12)
Exposed
Population
­
Females
15:
45
Thursday,
February
6,
2003
Obs
tpops
pf4574s
tpopn
pf4574n
1
82069
11759.01
416535.94
59682.10
149
"
Subsistence
Lead
Benefits
(
Sub
17).
sas"
Using
data
from
the
lead
specific
data
set,
changes
in
blood
lead
levels
are
estimated
and
benefits
calculated
for
reductions
in
cases
of
hypertension,
CHD,
brain
infarctions,
and
mortality
for
subsistence
fishermen
and
their
families.
150
/*
Program:
Subsistence
Lead
Benefits
(
Sub
17).
sas
Author:
Stuart
Smith,
Abt
Associates,
Inc.
Date:
07/
11/
2000
Updated:
11/
08/
2002
This
program
estimates
benefits
from
reduced
consumption
of
fish
tissue
containing
lead.
Reductions
in
MP&
M
facility
lead
loadings
result
in
lower
probabilties
of
lead
related
illnesses.
*/

%
let
tbl_
path
=
%
sysfunc(
pathname(
rundat));

title2
'
Subsistence
Lead
Related
Benefits
(
Sub
17)';

data
lead;
set
rundat.
lead;

/*
Combine
the
exposed
populations
for
subsistence;
fresh
and
saltwater.
*/

exppop
=
fwsubpop
+
swsubpop;

/*
Calculate
the
blood
lead
levels
under
the
baseline
and
post
compliance
scenarios
*/

PbB0m
=
5.40;
/*
Typical
PbB
conc.
in
absence
of
fish
cons.
for
men
*/
PbB0f
=
4.12;
/*
Typical
PbB
conc.
in
absence
of
fish
cons.
for
women
*/
bcf
=
49;
/*
Bioconcentration
factor
for
lead
L/
KG
*/
bksf
=
0.4;
/*
Biokinetic
slope
factor
for
lead
*/
ingm
=
142.4;
/*
Average
daily
fish
consumption
(
g/
day)
males
age
15­
44
*/
ingf
=
142.4;
/*
Average
daily
fish
consumption
(
g/
day)
females
age
15­
44
*/
ingm2
=
142.4;
/*
Average
daily
fish
consumption
(
g/
day)
males
over
age
45
*/
ingf2
=
142.4;
/*
Average
daily
fish
consumption
(
g/
day)
females
over
age
45
*/
af
=
0.03;
/*
Absorbtion
factor
*/
ef
=
365;
/*
Exposure
Frequency
*/
at
=
365;
/*
Average
timing,
or
duration
of
exposure
in
days
*/
cf
=
0.001;
/*
Conversion
factor
kg/
g
*/

/*
Uses
ingestion
for
males
ages
15­
44
*/
PbB_
bm
=
PbB0m
+
((
cnc_
base*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);
PbB_
1m
=
PbB0m
+
((
cnc_
trt1*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);
PbB_
2m
=
PbB0m
+
((
cnc_
trt2*
bcf*
ingm*
af*
bksf*
ef*
cf)/
at);

/*
Uses
ingestion
for
males
over
age
45
*/
PbB_
bme
=
PbB0m
+
((
cnc_
base*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);
PbB_
1me
=
PbB0m
+
((
cnc_
trt1*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);
PbB_
2me
=
PbB0m
+
((
cnc_
trt2*
bcf*
ingm2*
af*
bksf*
ef*
cf)/
at);

/*
Uses
ingestion
for
females
ages
15­
44
*/
PbB_
bf
=
PbB0f
+
((
cnc_
base*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);
PbB_
1f
=
PbB0f
+
((
cnc_
trt1*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);
PbB_
2f
=
PbB0f
+
((
cnc_
trt2*
bcf*
ingf*
af*
bksf*
ef*
cf)/
at);
151
/*
Uses
ingestion
for
females
over
age
45
*/
PbB_
bfe
=
PbB0f
+
((
cnc_
base*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);
PbB_
1fe
=
PbB0f
+
((
cnc_
trt1*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);
PbB_
2fe
=
PbB0f
+
((
cnc_
trt2*
bcf*
ingf2*
af*
bksf*
ef*
cf)/
at);

drop
fwsptpop
fwsubpop
swsptpop
swsubpop;
run;

/*
Print
reaches
in
the
lead
dataset
*/
proc
print;
var
reachnmb
PbB_
bm
PbB_
1m
PbB_
2m
cnc_
base
cnc_
trt1
cnc_
trt2
exppop
event_
wt;
quit;

data
men;
set
lead;

/*
Set
up
some
constants
*/
/*
All
percentages
based
on
the
1996
Census
numbers
*/

pctmale
=
0.319518;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
20­
74
*/
pctm4054
=
0.098405;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
40­
54
*/
pctm4059
=
0.118945;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
40­
59
*/
pctm4574
=
0.12933;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
45­
74
*/
pctm5564
=
0.03829;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
55­
64
*/
pctm6064
=
0.01775;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
60­
64
*/
pctm6574
=
0.031427;
/*
Percentage
of
U.
S.
Pop.
who
are
male,
age
65­
74
*/

bbpm4054
=
81.9;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
40­
54
*/
bbpm4059
=
81.8;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
40­
59
*/
bbpm4574
=
81.1;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
45­
74
*/
bbpm5564
=
80.6;
/*
Mean
diastolic
bp
from
NHANES
III,
males,
age
55­
64
*/
bbpm6064
=
79.5;
/*
Mean
diastolic
bp
from
NHANES
III,
males
age
60­
64
*/
bbpm6574
=
79.5;
/*
Mean
diastolic
bp
from
NHANES
III,
males
age
65­
74
*/

/**************************************************************************/
/*
Begin
Benefits
for
Males
*/
/**************************************************************************/

/*
Calculate
the
change
in
the
probability
of
hypertension
*/

pr_
hyp1
=
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
bm))))
­
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
1m))));

pr_
hyp2
=
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
bm))))
­
(
1/(
1
+
exp(
2.744
­
0.793*
log(
PbB_
2m))));

ben_
hyp1
=
pr_
hyp1
*
exppop
*
pctmale
*
1141;
ben_
hyp2
=
pr_
hyp2
*
exppop
*
pctmale
*
1141;
hyppop
=
exppop
*
pctmale;
hyp_
red1
=
pr_
hyp1
*
exppop
*
pctmale;
hyp_
red2
=
pr_
hyp2
*
exppop
*
pctmale;
152
/*
Calculate
Benefits
from
Changes
in
CHD,
non­
fatal
cases
only
which
are
2/
3
of
the
cases
*/

dbpm1
=
1.4
*
log(
PbB_
bm/
PbB_
1m);
dbpm2
=
1.4
*
log(
PbB_
bm/
PbB_
2m);

dbpm1e
=
1.4
*
log(
PbB_
bme/
PbB_
1me);
dbpm2e
=
1.4
*
log(
PbB_
bme/
PbB_
2me);

/*
Ages
40­
59,
change
in
the
10
year
prob
of
occ.
*/
d1m4059
=
bbpm4059
­
dbpm1;
d2m4059
=
bbpm4059
­
dbpm2;

pr1_
4059
=
(
1/(
1
+
exp(
4.996
­
0.030365*
bbpm4059)))
­
(
1/(
1
+
exp(
4.996
­
0.030365*
d1m4059)));

pr2_
4059
=
(
1/(
1
+
exp(
4.996
­
0.030365*
bbpm4059)))
­
(
1/(
1
+
exp(
4.996
­
0.030365*
d2m4059)));

/*
Ages
60­
64,
change
in
the
2
year
prob
of
occ.
*/
d1m6064
=
bbpm6064
­
dbpm1e;
d2m6064
=
bbpm6064
­
dbpm2e;

pr1_
6064
=
(
1/(
1
+
exp(
5.19676
­
0.02351*
bbpm6064)))
­
(
1/(
1
+
exp(
5.19676
­
0.02351*
d1m6064)));

pr2_
6064
=
(
1/(
1
+
exp(
5.19676
­
0.02351*
bbpm6064)))
­
(
1/(
1
+
exp(
5.19676
­
0.02351*
d2m6064)));

/*
Ages
65­
74,
change
in
the
2
year
prob
of
occ.
*/
d1m6574
=
bbpm6574
­
dbpm1e;
d2m6574
=
bbpm6574
­
dbpm2e;

pr1_
6574
=
(
1/(
1
+
exp(
4.90723
­
0.02031*
bbpm6574)))
­
(
1/(
1
+
exp(
4.90723
­
0.02031*
d1m6574)));

pr2_
6574
=
(
1/(
1
+
exp(
4.90723
­
0.02031*
bbpm6574)))
­
(
1/(
1
+
exp(
4.90723
­
0.02031*
d2m6574)));

b14059
=
(
pr1_
4059
*
exppop
*
pctm4059
*
76347
*
0.66)/
10;
b24059
=
(
pr2_
4059
*
exppop
*
pctm4059
*
76347
*
0.66)/
10;
b16064
=
(
pr1_
6064
*
exppop
*
pctm6064
*
76347
*
0.66)/
2;
b26064
=
(
pr2_
6064
*
exppop
*
pctm6064
*
76347
*
0.66)/
2;
b16574
=
(
pr1_
6574
*
exppop
*
pctm6574
*
76347
*
0.66)/
2;
b26574
=
(
pr2_
6574
*
exppop
*
pctm6574
*
76347
*
0.66)/
2;

/*
Number
of
reduced
cases
*/
r14059
=
(
pr1_
4059
*
exppop
*
pctm4059
*
0.66)/
10;
r24059
=
(
pr2_
4059
*
exppop
*
pctm4059
*
0.66)/
10;
r16064
=
(
pr1_
6064
*
exppop
*
pctm6064
*
0.66)/
2;
r26064
=
(
pr2_
6064
*
exppop
*
pctm6064
*
0.66)/
2;
r16574
=
(
pr1_
6574
*
exppop
*
pctm6574
*
0.66)/
2;
r26574
=
(
pr2_
6574
*
exppop
*
pctm6574
*
0.66)/
2;
153
popm4059
=
exppop
*
pctm4059;
popm6064
=
exppop
*
pctm6064;
popm6574
=
exppop
*
pctm6574;

/*
Calculate
Benefits
from
Reduced
Cerebrovascular
Accidents
and
Brain
Infarctions,
affects
men
aged
45­
74
Non­
fatal
only
which
is
70%
of
all
cases
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/
d1m4574
=
bbpm4574
­
dbpm1e;
d2m4574
=
bbpm4574
­
dbpm2e;

pr1_
cba
=
(
1/(
1
+
exp(
8.58889
­
0.04066*
bbpm4574)))
­
(
1/(
1
+
exp(
8.58889
­
0.04066*
d1m4574)));

pr2_
cba
=
(
1/(
1
+
exp(
8.58889
­
0.04066*
bbpm4574)))
­
(
1/(
1
+
exp(
8.58889
­
0.04066*
d2m4574)));

pr1_
bi
=
(
1/(
1
+
exp(
9.9516
­
0.04840*
bbpm4574)))
­
(
1/(
1
+
exp(
9.9516
­
0.04840*
d1m4574)));

pr2_
bi
=
(
1/(
1
+
exp(
9.9516
­
0.04840*
bbpm4574)))
­
(
1/(
1
+
exp(
9.9516
­
0.04840*
d2m4574)));

b1_
cba
=
(
pr1_
cba
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b2_
cba
=
(
pr2_
cba
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b1_
bi
=
(
pr1_
bi
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;
b2_
bi
=
(
pr2_
bi
*
exppop
*
pctm4574
*
335135
*
0.7)/
2;

r1_
cba
=
(
pr1_
cba
*
exppop
*
pctm4574
*
0.7)/
2;
r2_
cba
=
(
pr2_
cba
*
exppop
*
pctm4574
*
0.7)/
2;
r1_
bi
=
(
pr1_
bi
*
exppop
*
pctm4574
*
0.7)/
2;
r2_
bi
=
(
pr2_
bi
*
exppop
*
pctm4574
*
0.7)/
2;

popm4574
=
exppop
*
pctm4574;

/*
Calculate
Benefits
from
Premature
Mortality
*/

/*
Ages
40­
54,
change
in
the
12
year
prob
of
death
*/
d1m4054
=
bbpm4054
­
dbpm1;
d2m4054
=
bbpm4054
­
dbpm2;

m1_
4054
=
(
1/(
1
+
exp(
5.3158
­
0.03516*
bbpm4054)))
­
(
1/(
1
+
exp(
5.3158
­
0.03516*
d1m4054)));

m2_
4054
=
(
1/(
1
+
exp(
5.3158
­
0.03516*
bbpm4054)))
­
(
1/(
1
+
exp(
5.3158
­
0.03516*
d2m4054)));

/*
Ages
55­
64,
change
in
the
2
year
prob
of
death
*/
d1m5564
=
bbpm5564
­
dbpm1e;
d2m5564
=
bbpm5564
­
dbpm2e;
154
m1_
5564
=
(
1/(
1
+
exp(
4.89528
­
0.01866*
bbpm5564)))
­
(
1/(
1
+
exp(
4.89528
­
0.01866*
d1m5564)));

m2_
5564
=
(
1/(
1
+
exp(
4.89528
­
0.01866*
bbpm5564)))
­
(
1/(
1
+
exp(
4.89528
­
0.01866*
d2m5564)));

/*
Ages
65­
74,
change
in
the
2
year
prob
of
death
*/
d1m6574
=
bbpm6574
­
dbpm1e;
d2m6574
=
bbpm6574
­
dbpm2e;

m1_
6574
=
(
1/(
1
+
exp(
3.05723
­
0.00547*
bbpm6574)))
­
(
1/(
1
+
exp(
3.05723
­
0.00547*
d1m6574)));

m2_
6574
=
(
1/(
1
+
exp(
3.05723
­
0.00547*
bbpm6574)))
­
(
1/(
1
+
exp(
3.05723
­
0.00547*
d2m6574)));

m14054l
=
(
m1_
4054
*
exppop
*
pctm4054
*
2500000)/
12;
m24054l
=
(
m2_
4054
*
exppop
*
pctm4054
*
2500000)/
12;
m14054h
=
(
m1_
4054
*
exppop
*
pctm4054
*
13600000)/
12;
m24054h
=
(
m2_
4054
*
exppop
*
pctm4054
*
13600000)/
12;
mid40541
=
(
m1_
4054
*
exppop
*
pctm4054
*
6500000)/
12;
mid40542
=
(
m2_
4054
*
exppop
*
pctm4054
*
6500000)/
12;

m15564l
=
(
m1_
5564
*
exppop
*
pctm5564
*
2500000)/
2;
m25564l
=
(
m2_
5564
*
exppop
*
pctm5564
*
2500000)/
2;
m15564h
=
(
m1_
5564
*
exppop
*
pctm5564
*
13600000)/
2;
m25564h
=
(
m2_
5564
*
exppop
*
pctm5564
*
13600000)/
2;
mid55641
=
(
m1_
5564
*
exppop
*
pctm5564
*
6500000)/
2;
mid55642
=
(
m2_
5564
*
exppop
*
pctm5564
*
6500000)/
2;

m16574l
=
(
m1_
6574
*
exppop
*
pctm6574
*
2500000)/
2;
m26574l
=
(
m2_
6574
*
exppop
*
pctm6574
*
2500000)/
2;
m16574h
=
(
m1_
6574
*
exppop
*
pctm6574
*
13600000)/
2;
m26574h
=
(
m2_
6574
*
exppop
*
pctm6574
*
13600000)/
2;
mid65741
=
(
m1_
6574
*
exppop
*
pctm6574
*
6500000)/
2;
mid65742
=
(
m2_
6574
*
exppop
*
pctm6574
*
6500000)/
2;

r14054l
=
(
m1_
4054
*
exppop
*
pctm4054)/
12;
r24054l
=
(
m2_
4054
*
exppop
*
pctm4054)/
12;
r14054h
=
(
m1_
4054
*
exppop
*
pctm4054)/
12;
r24054h
=
(
m2_
4054
*
exppop
*
pctm4054)/
12;
r15564l
=
(
m1_
5564
*
exppop
*
pctm5564)/
2;
r25564l
=
(
m2_
5564
*
exppop
*
pctm5564)/
2;
r15564h
=
(
m1_
5564
*
exppop
*
pctm5564)/
2;
r25564h
=
(
m2_
5564
*
exppop
*
pctm5564)/
2;
r16574l
=
(
m1_
6574
*
exppop
*
pctm6574)/
2;
r26574l
=
(
m2_
6574
*
exppop
*
pctm6574)/
2;
r16574h
=
(
m1_
6574
*
exppop
*
pctm6574)/
2;
r26574h
=
(
m2_
6574
*
exppop
*
pctm6574)/
2;

popm4054
=
exppop
*
pctm4054;
popm5564
=
exppop
*
pctm5564;
155
run;

data
women;
set
lead;

/*
Set
up
some
constants
*/
/*
All
percentages
based
on
the
1996
Census
numbers
*/
pctfem
=
0.511735;
/*
Percentage
of
U.
S.
Pop.
who
are
female
*/
pctf4574
=
0.143282;
/*
Percentage
of
U.
S.
Pop.
who
are
female,
age
45­
74
*/
bbpf4574
=
76.5;
/*
Mean
diastolic
bp
from
NHANES
III,
females
age
45­
74
*/

/**************************************************************************/
/*
Begin
Benefits
for
Females
*/
/**************************************************************************/

/*
No
dose­
response
function
for
female
hypertension
*/

/*
Calculate
Benefits
from
Changes
in
CHD,
non­
fatal
cases
only
which
are
2/
3
of
the
cases
*/

dbpf1
=
(
0.6
*
1.4)
*
log(
PbB_
bfe/
PbB_
1fe);
dbpf2
=
(
0.6
*
1.4)
*
log(
PbB_
bfe/
PbB_
2fe);

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/
d1f4574
=
bbpf4574
­
dbpf1;
d2f4574
=
bbpf4574
­
dbpf2;

pr1_
4574
=
(
1/(
1
+
exp(
6.9401
­
0.03072*
bbpf4574)))
­
(
1/(
1
+
exp(
6.9401
­
0.03072*
d1f4574)));

pr2_
4574
=
(
1/(
1
+
exp(
6.9401
­
0.03072*
bbpf4574)))
­
(
1/(
1
+
exp(
6.9401
­
0.03072*
d2f4574)));

b1fchd
=
(
pr1_
4574
*
exppop
*
pctf4574
*
76347
*
0.66)/
2;
b2fchd
=
(
pr2_
4574
*
exppop
*
pctf4574
*
76347
*
0.66)/
2;

red_
chd1
=
(
pr1_
4574
*
exppop
*
pctf4574
*
0.66)/
2;
red_
chd2
=
(
pr2_
4574
*
exppop
*
pctf4574
*
0.66)/
2;

popf4574
=
exppop
*
pctf4574;

/*
Calculate
Benefits
from
Reduced
Cerebrovascular
Accidents
and
Brain
Infarctions,
affects
women
aged
45­
74
Non­
fatal
only
which
is
70%
of
all
cases
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
occ.
*/

pr1f_
bi
=
(
1/(
1
+
exp(
10.6716
­
0.0544*
bbpf4574)))
­
(
1/(
1
+
exp(
10.6716
­
0.0544*
d1f4574)));

pr2f_
bi
=
(
1/(
1
+
exp(
10.6716
­
0.0544*
bbpf4574)))
­
156
(
1/(
1
+
exp(
10.6716
­
0.0544*
d2f4574)));

pr1f_
cba
=
(
1/(
1
+
exp(
9.07737
­
0.04287*
bbpf4574)))
­
(
1/(
1
+
exp(
9.07737
­
0.04287*
d1f4574)));

pr2f_
cba
=
(
1/(
1
+
exp(
9.07737
­
0.04287*
bbpf4574)))
­
(
1/(
1
+
exp(
9.07737
­
0.04287*
d2f4574)));

b1f_
cba
=
(
pr1f_
cba
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b2f_
cba
=
(
pr2f_
cba
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b1f_
bi
=
(
pr1f_
bi
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;
b2f_
bi
=
(
pr2f_
bi
*
exppop
*
pctf4574
*
251351
*
0.7)/
2;

r1f_
cba
=
(
pr1f_
cba
*
exppop
*
pctf4574
*
0.7)/
2;
r2f_
cba
=
(
pr2f_
cba
*
exppop
*
pctf4574
*
0.7)/
2;
r1f_
bi
=
(
pr1f_
bi
*
exppop
*
pctf4574
*
0.7)/
2;
r2f_
bi
=
(
pr2f_
bi
*
exppop
*
pctf4574
*
0.7)/
2;

/*
Calculate
Benefits
from
Premature
Mortality
*/

/*
Ages
45­
74,
change
in
the
2
year
prob
of
death
*/

m1_
4574
=
(
1/(
1
+
exp(
5.40374
­
0.01511*
bbpf4574)))
­
(
1/(
1
+
exp(
5.40374
­
0.01511*
d1f4574)));

m2_
4574
=
(
1/(
1
+
exp(
5.40374
­
0.01511*
bbpf4574)))
­
(
1/(
1
+
exp(
5.40374
­
0.01511*
d2f4574)));

m14574l
=
(
m1_
4574
*
exppop
*
pctf4574
*
2500000)/
2;
m24574l
=
(
m2_
4574
*
exppop
*
pctf4574
*
2500000)/
2;
m14574h
=
(
m1_
4574
*
exppop
*
pctf4574
*
13600000)/
2;
m24574h
=
(
m2_
4574
*
exppop
*
pctf4574
*
13600000)/
2;
mid45741
=
(
m1_
4574
*
exppop
*
pctf4574
*
6500000)/
2;
mid45742
=
(
m2_
4574
*
exppop
*
pctf4574
*
6500000)/
2;

r14574l
=
(
m1_
4574
*
exppop
*
pctf4574)/
2;
r24574l
=
(
m2_
4574
*
exppop
*
pctf4574)/
2;
r14574h
=
(
m1_
4574
*
exppop
*
pctf4574)/
2;
r24574h
=
(
m2_
4574
*
exppop
*
pctf4574)/
2;
run;

/*
Sum
up
benefits
and
reduced
cases
and
print
results
*/

proc
means
sum
data=
men;
var
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
157
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
output
out
=
rundat.
msubsamp(
drop
=
_
type_
_
freq_)
sum
=
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
title3
'
Mens
Lead
Related
Benefits
(
Sample
Basis)';
quit;

proc
means
sum
data=
men;
weight
event_
wt;
var
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
output
out
=
rundat.
msubwght(
drop
=
_
type_
_
freq_)
sum
=
ben_
hyp1
ben_
hyp2
hyppop
b14059
b24059
b16064
b26064
b16574
b26574
popm4059
popm6064
popm6574
b1_
cba
b2_
cba
b1_
bi
b2_
bi
popm4574
m14054l
m14054h
m24054l
m24054h
m15564l
m15564h
m25564l
m25564h
m16574l
m16574h
m26574l
m26574h
mid40541
mid40542
mid55641
mid55642
mid65741
mid65742
popm4054
popm5564
hyp_
red1
hyp_
red2
r14059
r24059
r16064
r26064
r16574
r26574
r1_
cba
r2_
cba
r1_
bi
r2_
bi
r14054l
r24054l
158
r15564l
r25564l
r16574l
r26574l;
title3
'
Mens
Lead
Related
Benefits
(
National
Basis)';
quit;

%
macro
maleleadtable(
tbl_
in=,
tbl_
out=);
data
rundat.&
tbl_
out
(
keep=
gender
categ
agecat
type
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h);
set
rundat.&
tbl_
in;

length
gender
$
50
categ
$
50
agecat
$
50
type
$
50;

array
_
hyper[
4]
ben_
hyp1
ben_
hyp2
hyp_
red1
hyp_
red2;
array
_
chd[
12]
b14059
b24059
b16064
b26064
b16574
b26574
r14059
r24059
r16064
r26064
r16574
r26574;
array
_
cba[
4]
b1_
cba
b2_
cba
r1_
cba
r2_
cba;
array
_
bi[
4]
b1_
bi
b2_
bi
r1_
bi
r2_
bi;
array
_
mort[
24]
m14054l
mid40541
m14054h
m24054l
mid40542
m24054h
m15564l
mid55641
m15564h
m25564l
mid55642
m25564h
m16574l
mid65741
m16574h
m26574l
mid65742
m26574h
r14054l
r24054l
r15564l
r25564l
r16574l
r26574l;
array
ben[
6]
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
array
cases[
2]
case1
case2;

do
i
=
1
to
9;
if
i
=
1
then
do;
ben[
1]
=
_
hyper[
1];
ben[
2]
=
_
hyper[
1];
ben[
3]
=
_
hyper[
1];
ben[
4]
=
_
hyper[
2];
ben[
5]
=
_
hyper[
2];
ben[
6]
=
_
hyper[
2];
cases[
1]
=
_
hyper[
3];
cases[
2]
=
_
hyper[
4];
categ
=
"
Hypertension";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
2
then
do;
ben[
1]
=
_
chd[
1];
ben[
2]
=
_
chd[
1];
ben[
3]
=
_
chd[
1];
ben[
4]
=
_
chd[
2];
ben[
5]
=
_
chd[
2];
ben[
6]
=
_
chd[
2];
cases[
1]
=
_
chd[
7];
cases[
2]
=
_
chd[
8];
categ
=
"
CHD";
agecat
=
"
Age
40­
59";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
3
then
do;
ben[
1]
=
_
chd[
3];
ben[
2]
=
_
chd[
3];
ben[
3]
=
_
chd[
3];
ben[
4]
=
_
chd[
4];
ben[
4]
=
_
chd[
4];
ben[
6]
=
_
chd[
4];
cases[
1]
=
_
chd[
9];
cases[
2]
=
_
chd[
10];
categ
=
"
CHD";
agecat
=
"
Age
60­
64";
159
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
4
then
do;
ben[
1]
=
_
chd[
5];
ben[
2]
=
_
chd[
5];
ben[
3]
=
_
chd[
5];
ben[
4]
=
_
chd[
6];
ben[
4]
=
_
chd[
6];
ben[
6]
=
_
chd[
6];
cases[
1]
=
_
chd[
11];
cases[
2]
=
_
chd[
12];
categ
=
"
CHD";
agecat
=
"
Age
65­
74";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
5
then
do;
ben[
1]
=
_
cba[
1];
ben[
2]
=
_
cba[
1];
ben[
3]
=
_
cba[
1];
ben[
4]
=
_
cba[
2];
ben[
5]
=
_
cba[
2];
ben[
6]
=
_
cba[
2];
cases[
1]
=
_
cba[
3];
cases[
2]
=
_
cba[
4];
categ
=
"
CBA";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
6
then
do;
ben[
1]
=
_
bi[
1];
ben[
2]
=
_
bi[
1];
ben[
3]
=
_
bi[
1];
ben[
4]
=
_
bi[
2];
ben[
5]
=
_
bi[
2];
ben[
6]
=
_
bi[
2];
cases[
1]
=
_
bi[
3];
cases[
2]
=
_
bi[
4];
categ
=
"
BI";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
7
then
do;
ben[
1]
=
_
mort[
1];
ben[
2]
=
_
mort[
2];
ben[
3]
=
_
mort[
3];
ben[
4]
=
_
mort[
4];
ben[
5]
=
_
mort[
5];
ben[
6]
=
_
mort[
6];
cases[
1]
=
_
mort[
19];
cases[
2]
=
_
mort[
20];
categ
=
"
Mortality";
agecat
=
"
Age
40­
54";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
if
i
=
8
then
do;
ben[
1]
=
_
mort[
7];
ben[
2]
=
_
mort[
8];
ben[
3]
=
_
mort[
9];
ben[
4]
=
_
mort[
10];
ben[
5]
=
_
mort[
11];
ben[
6]
=
_
mort[
12];
cases[
1]
=
_
mort[
21];
cases[
2]
=
_
mort[
22];
categ
=
"
Mortality";
agecat
=
"
Age
55­
64";
type
=
"
Sub";
gender
=
"
Male";
output;
160
end;
if
i
=
9
then
do;
ben[
1]
=
_
mort[
13];
ben[
2]
=
_
mort[
14];
ben[
3]
=
_
mort[
15];
ben[
4]
=
_
mort[
16];
ben[
5]
=
_
mort[
17];
ben[
6]
=
_
mort[
18];
cases[
1]
=
_
mort[
23];
cases[
2]
=
_
mort[
24];
categ
=
"
Mortality";
agecat
=
"
Age
65­
74";
type
=
"
Sub";
gender
=
"
Male";
output;
end;
end;
run;
%
mend;

%
maleleadtable(
tbl_
in=
msubsamp,
tbl_
out=
tblmbsmp);
%
maleleadtable(
tbl_
in=
msubwght,
tbl_
out=
tblmbwgt);

proc
means
sum
data=
women;
var
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
output
out
=
rundat.
fsubsamp(
drop
=
_
type_
_
freq_)
sum
=
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
title3
'
Female
Lead
Related
Benefits
(
Sample
Basis)';
quit;

proc
means
sum
data=
women;
weight
event_
wt;
var
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
output
out
=
rundat.
fsubwght(
drop
=
_
type_
_
freq_)
sum
=
b1fchd
b2fchd
popf4574
b1f_
cba
b2f_
cba
b1f_
bi
b2f_
bi
m14574l
m14574h
m24574l
m24574h
mid45741
mid45742
red_
chd1
red_
chd2
r1f_
cba
r2f_
cba
r1f_
bi
r2f_
bi
r14574l
r14574h
r24574l
r24574h;
title3
'
Female
Lead
Related
Benefits
(
National
Basis)';
quit;
161
%
macro
femaleleadtable(
tbl_
in=,
tbl_
out=);
data
rundat.&
tbl_
out
(
keep=
gender
categ
agecat
type
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h);
set
rundat.&
tbl_
in;

length
gender
$
50
categ
$
50
agecat
$
50
type
$
50;

array
_
chd[
4]
b1fchd
b2fchd
red_
chd1
red_
chd2;
array
_
cba[
4]
b1f_
cba
b2f_
cba
r1f_
cba
r2f_
cba;
array
_
bi[
4]
b1f_
bi
b2f_
bi
r1f_
bi
r2f_
bi;
array
_
mort[
8]
m14574l
mid45741
m14574h
m24574l
mid45742
m24574h
r14574l
r24574l;

array
ben[
6]
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
array
cases[
2]
case1
case2;

do
i
=
1
to
4;
if
i
=
1
then
do;
ben[
1]
=
_
chd[
1];
ben[
2]
=
_
chd[
1];
ben[
3]
=
_
chd[
1];
ben[
4]
=
_
chd[
2];
ben[
5]
=
_
chd[
2];
ben[
6]
=
_
chd[
2];
cases[
1]
=
_
chd[
3];
cases[
2]
=
_
chd[
4];
categ
=
"
CHD";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Female";
output;
end;
if
i
=
2
then
do;
ben[
1]
=
_
cba[
1];
ben[
2]
=
_
cba[
1];
ben[
3]
=
_
cba[
1];
ben[
4]
=
_
cba[
2];
ben[
5]
=
_
cba[
2];
ben[
6]
=
_
cba[
2];
cases[
1]
=
_
cba[
3];
cases[
2]
=
_
cba[
4];
categ
=
"
CBA";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Female";
output;
end;
if
i
=
3
then
do;
ben[
1]
=
_
bi[
1];
ben[
2]
=
_
bi[
1];
ben[
3]
=
_
bi[
1];
ben[
4]
=
_
bi[
2];
ben[
5]
=
_
bi[
2];
ben[
6]
=
_
bi[
2];
cases[
1]
=
_
bi[
3];
cases[
2]
=
_
bi[
4];
categ
=
"
BI";
agecat
=
"
All";
type
=
"
Sub";
gender
=
"
Female";
output;
end;
if
i
=
4
then
do;
ben[
1]
=
_
mort[
1];
ben[
2]
=
_
mort[
2];
ben[
3]
=
_
mort[
3];
ben[
4]
=
_
mort[
4];
ben[
5]
=
_
mort[
5];
ben[
6]
=
_
mort[
6];
cases[
1]
=
_
mort[
7];
cases[
2]
=
_
mort[
8];
categ
=
"
Mortality";
agecat
=
"
Age
45­
74";
type
=
"
Sub";
162
gender
=
"
Female";
output;
end;
end;
run;
%
mend;

%
femaleleadtable(
tbl_
in=
fsubsamp,
tbl_
out=
tblfbsmp);
%
femaleleadtable(
tbl_
in=
fsubwght,
tbl_
out=
tblfbwgt);

data
rundat.
tblsub_
s;
set
rundat.
tblmbsmp
rundat.
tblfbsmp;
run;

data
rundat.
tblsub_
w;
set
rundat.
tblmbwgt
rundat.
tblfbwgt;
run;

proc
means
sum
data=
men;
var
event_
wt;
title3
'
Number
of
Affected
Reaches
Nationally';
quit;

proc
sql;

create
table
malepop
as
select
sum(
exppop)
as
tpops,
sum(
exppop*
pctmale)
as
pm2074s,
sum(
exppop*
pctm4574)
as
pm4574s,
sum(
exppop*
event_
wt)
as
tpopn,
sum(
exppop*
pctmale*
event_
wt)
as
pm2074n,
sum(
exppop*
pctm4574*
event_
wt)
as
pm4574n
from
men;

create
table
fempop
as
select
sum(
exppop)
as
tpops,
sum(
exppop*
pctf4574)
as
pf4574s,
sum(
exppop*
event_
wt)
as
tpopn,
sum(
exppop*
pctf4574*
event_
wt)
as
pf4574n
from
women;
quit;

proc
print
data=
malepop;
var
tpops
pm2074s
pm4574s
tpopn
pm2074n
pm4574n;
label
tpops
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Smpl)'
pm2074s
=
'
Exposed
Males
Age
20­
74
(
Smpl)'
pm4574s
=
'
Exposed
Males
Age
45­
74
(
Smpl)'
tpopn
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Natl)'
pm2074n
=
'
Exposed
Males
Age
20­
74
(
Natl)'
pm4574n
=
'
Exposed
Males
Age
45­
74
(
Natl)';
title3
'
Exposed
Population
­
Males';
quit;
163
proc
print
data=
fempop;
var
tpops
pf4574s
tpopn
pf4574n;
label
tpops
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Smpl)'
pf4574s
=
'
Exposed
Females
Age
45­
74
(
Smpl)'
tpopn
=
'
Total
Exposed
Pop
(
male
+
fem)
(
Natl)'
pf4574n
=
'
Exposed
Females
Age
45­
74
(
Natl)';
title3
'
Exposed
Population
­
Females';
quit;

data
rundat.
leadtbls;
if
_
n_
=
1
then
do;
do
i
=
1
to
13;
set
rundat.
tblspt_
s
point
=
i;
output;
set
rundat.
tblsub_
s
point
=
i;
output;
end;
end;
run;

data
rundat.
leadtblw;
if
_
n_
=
1
then
do;
do
i
=
1
to
13;
set
rundat.
tblspt_
w
point
=
i;
output;
set
rundat.
tblsub_
w
point
=
i;
output;
end;
end;
run;

proc
summary
data=
rundat.
leadtbls
nway;
class
gender
categ;
var
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h;
output
out
=
rundat.
leadsums
(
drop=
_
type_
_
freq_)
sum
=
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h;
run;

proc
summary
data=
rundat.
leadtblw
nway;
class
gender
categ;
var
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h;
output
out
=
rundat.
leadsumw
(
drop=
_
type_
_
freq_)
sum
=
case1
ben1l
ben1m
ben1h
case2
ben2l
ben2m
ben2h;
run;

proc
summary
data=
rundat.
leadtbls;
var
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
output
out
=
rundat.
ldalls
(
drop=
_
type_
_
freq_)
sum
=
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
run;

proc
summary
data=
rundat.
leadtblw;
164
var
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
output
out
=
rundat.
ldallw
(
drop=
_
type_
_
freq_)
sum
=
ben1l
ben1m
ben1h
ben2l
ben2m
ben2h;
run;

data
exlleadw
(
keep=
titlecol
col1­
col10
label="
Adult
Lead
Benefits
­
Weighted");
set
rundat.
leadsumw
rundat.
ldallw;
length
titlecol
$
50
col1
­
col10
$
128;

if
_
n_
=
1
then
do;
titlecol
=
"
Table:
Adult
Lead
Benefits,
Weighted
Estimates";
output;
titlecol
=
"
Run
on
&
sysdate";
output;
titlecol
=
"
";
output;
col3
=
"
Option
1";
col7
=
"
Option
2";
output;
col3
=
"
Reduced
Cases";
col4
=
"
Low";
col5
=
"
Mid";
col6
=
"
High";
col7
=
"
Reduced
Cases";
col8
=
"
Low";
col9
=
"
Mid";
col10
=
"
High";
output;
end;

if
_
n_
=
1
or
_
n_
=
6
then
do;
col1
=
gender;
end;

if
_
n_
=
10
then
do;
col1
=
"
Total";
end;

col2
=
trim(
categ);
col3
=
trim(
case1);
col4
=
trim(
ben1l);
col5
=
trim(
ben1m);
col6
=
trim(
ben1h);
col7
=
trim(
case2);
col8
=
trim(
ben2l);
col9
=
trim(
ben2m);
col10
=
trim(
ben2h);
output;
run;

proc
export
data
=
exlleadw
outfile
=
"
C:\
Output
Files\
Adult
Lead.
xls"
dbms=
excel2000
replace;
run;

/*
End
Subsistence
Lead
Benefits
(
Sub
17).
sas
*/
165
Program
Output
166
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
87
Subsistence
Lead
Related
Benefits
(
Sub
17)
Mens
Lead
Related
Benefits
(
National
Basis)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Sum
                        
ben_
hyp1
16492.49
ben_
hyp2
23968.72
hyppop
6653.82
b14059
632.9234048
b24059
921.4129639
b16064
176.1924563
b26064
256.5073159
b16574
278.4498246
b26574
405.3839462
popm4059
2476.98
popm6064
369.6358339
popm6574
654.4532593
b1_
cba
1540.77
b2_
cba
2242.89
b1_
bi
883.9719789
b2_
bi
1286.73
popm4574
2693.24
m14054l
26727.52
m14054h
145397.72
m24054l
38909.27
m24054h
211666.45
m15564l
14106.13
m15564h
76737.33
m25564l
20536.70
m25564h
111719.66
m16574l
6800.87
m16574h
36996.73
m26574l
9901.87
m26574h
53866.19
mid40541
69491.56
mid40542
101164.11
mid55641
36675.93
mid55642
53395.42
mid65741
17682.26
mid65742
25744.87
popm4054
2049.24
popm5564
797.3721736
hyp_
red1
14.4544170
hyp_
red2
21.0067648
r14059
0.0082901
r24059
0.0120688
r16064
0.0023078
r26064
0.0033598
r16574
0.0036472
167
r26574
0.0053098
r1_
cba
0.0045975
r2_
cba
0.0066925
r1_
bi
0.0026377
r2_
bi
0.0038395
r14054l
0.0106910
r24054l
0.0155637
r15564l
0.0056425
r25564l
0.0082147
r16574l
0.0027203
r26574l
0.0039607
                        
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
89
Subsistence
Lead
Related
Benefits
(
Sub
17)
Female
Lead
Related
Benefits
(
National
Basis)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Sum
                        
b1fchd
432.7430596
b2fchd
630.7469027
popf4574
2983.78
b1f_
cba
639.0835090
b2f_
cba
931.4497257
b1f_
bi
399.0219674
b2f_
bi
581.5365850
m14574l
14750.25
m14574h
80241.36
m24574l
21500.72
m24574h
116963.92
mid45741
38350.65
mid45742
55901.87
red_
chd1
0.0056681
red_
chd2
0.0082616
r1f_
cba
0.0025426
r2f_
cba
0.0037058
r1f_
bi
0.0015875
r2f_
bi
0.0023136
r14574l
0.0059001
r14574h
0.0059001
r24574l
0.0086003
r24574h
0.0086003
                        
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
90
Subsistence
Lead
Related
Benefits
(
Sub
17)
Number
of
Affected
Reaches
Nationally
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Analysis
Variable
:
EVENT_
WT
Event
Weight
168
Sum
            
23.7361201
            
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
91
Subsistence
Lead
Related
Benefits
(
Sub
17)
Exposed
Population
­
Males
15:
45
Thursday,
February
6,
2003
Obs
tpops
pm2074s
pm4574s
tpopn
pm2074n
pm4574n
1
4103
1310.98
530.641
20824.55
6653.82
2693.24
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
92
Subsistence
Lead
Related
Benefits
(
Sub
17)
Exposed
Population
­
Females
15:
45
Thursday,
February
6,
2003
Obs
tpops
pf4574s
tpopn
pf4574n
1
4103
587.886
20824.55
2983.78
169
Lead
Benefits
for
Children
170
"
Cancer­
fish
(
Sub
4).
sas"
From
the
data
set
of
pollutant
concentrations
for
each
RF1
reach,
observations
for
in­
stream
lead
concentrations
are
outputed
into
a
separate
data
set.
See
Cancer
from
Fish
Consumption
for
program
code
and
output.
171
"
NeoNatal
Mortality
Benefits
(
Sub
15).
sas"
Using
data
from
the
lead
specific
data
set,
changes
in
blood
lead
levels
for
women
of
childbearing
age
are
estimated.
Reduced
risk
of
neonatal
mortality
is
calculated
and
monetized.
172
/*
Program:
NeoNatal
Mortality
Benefits
(
Sub
15).
sas
Author:
Stuart
Smith,
Abt
Associates,
Inc.
Date:
07/
11/
2000
Updated:
11/
08/
2002
This
program
estimates
benefits
from
reduced
consumption
of
fish
tissue
containing
lead.
Reductions
in
MP&
M
facility
lead
loadings
result
in
lower
probabilties
of
lead
related
illnesses.
*/

title2
'
NeoNatal
Mortality
Lead
Benefits
(
Sub
15)';

data
lead;
set
rundat.
lead;

/*
Combine
the
exposed
recreational
and
subsistence
populations;
fresh
and
salt
water
*/

sptpop
=
fwsptpop
+
swsptpop;
subpop
=
fwsubpop
+
swsubpop;

/*
Calculate
the
blood
lead
levels
under
the
baseline
and
post
compliance
scenarios
*/

PbB0f
=
2.05;
/*
Typical
PbB
conc.
in
absence
of
fish
cons
for
women
Aged
15­
45
*/
bcf
=
49;
/*
Bioconcentration
factor
for
lead
L/
KG
*/
bksf
=
0.4;
/*
Biokinetic
slope
factor
for
lead
*/
ing
=
17.5;
/*
Average
daily
fish
consumption
(
g/
day)
recreational
*/
ing2
=
142.4;
/*
Average
daily
fish
consumption
(
g/
day)
subsistence
*/
af
=
0.03;
/*
Absorbtion
factor
*/
ef
=
365;
/*
Exposure
Frequency
*/
at
=
365;
/*
Average
timing,
or
duration
of
exposure
in
days
*/
cf
=
0.001;
/*
Conversion
factor
kg/
g
*/

/*
Recreational
*/
PbB_
bfr
=
(
PbB0f
+
(
cnc_
base*
bcf*
ing*
af*
bksf*
ef*
cf))/
at;
PbB_
1fr
=
(
PbB0f
+
(
cnc_
trt1*
bcf*
ing*
af*
bksf*
ef*
cf))/
at;
PbB_
2fr
=
(
PbB0f
+
(
cnc_
trt2*
bcf*
ing*
af*
bksf*
ef*
cf))/
at;

/*
Subsistence
*/
PbB_
bfs
=
(
PbB0f
+
(
cnc_
base*
bcf*
ing2*
af*
bksf*
ef*
cf))/
at;
PbB_
1fs
=
(
PbB0f
+
(
cnc_
trt1*
bcf*
ing2*
af*
bksf*
ef*
cf))/
at;
PbB_
2fs
=
(
PbB0f
+
(
cnc_
trt2*
bcf*
ing2*
af*
bksf*
ef*
cf))/
at;

/*
Changes
in
blood
lead
levels
*/
diff1r
=
PbB_
bfr
­
PbB_
1fr;
diff2r
=
PbB_
bfr
­
PbB_
2fr;

diff1s
=
PbB_
bfs
­
PbB_
1fs;
diff2s
=
PbB_
bfs
­
PbB_
2fs;

drop
fwsptpop
fwsubpop
swsptpop
swsubpop;
173
run;

proc
print;
var
reachnmb
PbB_
bfr
PbB_
1fr
PbB_
2fr
PbB_
bfs
PbB_
1fs
PbB_
2fs
diff1r
diff2r
diff1s
diff2s
cnc_
base
cnc_
trt1
cnc_
trt2
sptpop
subpop
event_
wt;
quit;

data
neo;
set
lead;

pctfem
=
0.22531;
/*
Percentage
of
Women
Aged
15­
45
(
Childbearing
Age)
*/

expspt
=
sptpop*
pctfem;
expsub
=
subpop*
pctfem;

red1r
=
diff1r
*
sptpop
*
pctfem
*
0.0001;
red2r
=
diff2r
*
sptpop
*
pctfem
*
0.0001;

red1s
=
diff1s
*
subpop
*
pctfem
*
0.0001;
red2s
=
diff2s
*
subpop
*
pctfem
*
0.0001;

ben1_
lr
=
red1r
*
2500000;
ben1_
hr
=
red1r
*
13600000;
ben2_
lr
=
red2r
*
2500000;
ben2_
hr
=
red2r
*
13600000;
mid1r
=
red1r
*
6500000;
mid2r
=
red2r
*
6500000;

ben1_
ls
=
red1s
*
2500000;
ben1_
hs
=
red1s
*
13600000;
ben2_
ls
=
red2s
*
2500000;
ben2_
hs
=
red2s
*
13600000;
mid1s
=
red1s
*
6500000;
mid2s
=
red2s
*
6500000;

tot1_
l
=
ben1_
lr
+
ben1_
ls;
tot1_
h
=
ben1_
hr
+
ben1_
hs;
tot2_
l
=
ben2_
lr
+
ben2_
ls;
tot2_
h
=
ben2_
hr
+
ben2_
hs;
mid1
=
mid1r
+
mid1s;
mid2
=
mid2r
+
mid2s;
run;

proc
means
sum
data=
neo;
var
red1r
red2r
red1s
red2s
ben1_
lr
ben1_
ls
ben1_
hr
ben1_
hs
ben2_
lr
ben2_
ls
ben2_
hr
ben2_
hs
tot1_
l
tot1_
h
tot2_
l
tot2_
h
mid1r
mid1s
mid2r
mid2s
mid1
mid2
expspt
expsub;
label
red1r
=
'
Number
of
Reductions
rec
(
Opt
1)'
red2r
=
'
Number
of
Reductions
rec
(
Opt
2)'
red1s
=
'
Number
of
Reductions
sub
(
Opt
1)'
red2s
=
'
Number
of
Reductions
sub
(
Opt
2)'
174
ben1_
lr
=
'
Sport
Low
(
Opt1)'
ben1_
ls
=
'
Sub
Low
(
Opt1)'
ben1_
hr
=
'
Sport
High
(
Opt1)'
ben1_
hs
=
'
Sub
High
(
Opt1)'
ben2_
lr
=
'
Sport
Low
(
Opt2)'
ben2_
ls
=
'
Sub
Low
(
Opt2)'
ben2_
hr
=
'
Sport
High
(
Opt2)'
ben2_
hs
=
'
Sub
High
(
Opt2)'
tot1_
l
=
'
Total
Benefits
­
Low
(
Opt
1)'
tot1_
h
=
'
Total
Benefits
­
High
(
Opt
1)'
tot2_
l
=
'
Total
Benefits
­
Low
(
Opt
2)'
tot2_
h
=
'
Total
Benefits
­
High
(
Opt
2)'
mid1
=
'
Mid
Benefits
(
Opt1)'
mid2
=
'
Mid
Benefits
(
Opt2)'
expspt
=
'
Exposed
Sport
Pop
(
15­
45yrs)'
expsub
=
'
Exposed
Sub
Pop
(
15­
45yrs)';
title3
'
Sample
Basis';
quit;

proc
means
sum
data=
neo;
weight
event_
wt;
var
red1r
red2r
red1s
red2s
ben1_
lr
ben1_
ls
ben1_
hr
ben1_
hs
ben2_
lr
ben2_
ls
ben2_
hr
ben2_
hs
tot1_
l
tot1_
h
tot2_
l
tot2_
h
mid1r
mid1s
mid2r
mid2s
mid1
mid2
expspt
expsub;
label
red1r
=
'
Number
of
Reductions
rec
(
Opt
1)'
red2r
=
'
Number
of
Reductions
rec
(
Opt
2)'
red1s
=
'
Number
of
Reductions
sub
(
Opt
1)'
red2s
=
'
Number
of
Reductions
sub
(
Opt
2)'
ben1_
lr
=
'
Sport
Low
(
Opt1)'
ben1_
ls
=
'
Sub
Low
(
Opt1)'
ben1_
hr
=
'
Sport
High
(
Opt1)'
ben1_
hs
=
'
Sub
High
(
Opt1)'
ben2_
lr
=
'
Sport
Low
(
Opt2)'
ben2_
ls
=
'
Sub
Low
(
Opt2)'
ben2_
hr
=
'
Sport
High
(
Opt2)'
ben2_
hs
=
'
Sub
High
(
Opt2)'
tot1_
l
=
'
Total
Benefits
­
Low
(
Opt
1)'
tot1_
h
=
'
Total
Benefits
­
High
(
Opt
1)'
tot2_
l
=
'
Total
Benefits
­
Low
(
Opt
2)'
tot2_
h
=
'
Total
Benefits
­
High
(
Opt
2)'
mid1
=
'
Mid
Benefits
(
Opt1)'
mid2
=
'
Mid
Benefits
(
Opt2)'
expspt
=
'
Exposed
Sport
Pop
(
15­
45yrs)'
expsub
=
'
Exposed
Sub
Pop
(
15­
45yrs)';
output
out
=
rundat.
chld_
neo
(
drop=
_
type_
_
freq_)
sum
=
red1r
red2r
red1s
red2s
ben1_
lr
ben1_
ls
ben1_
hr
ben1_
hs
ben2_
lr
ben2_
ls
ben2_
hr
ben2_
hs
tot1_
l
tot1_
h
tot2_
l
tot2_
h
mid1r
mid1s
mid2r
mid2s
mid1
mid2
expspt
expsub;
title3
'
National
Basis';
quit;

data
rundat.
neo_
tbl
(
keep=
gender
categ
type
exppop
case1
case2
benlow1
benlow2
benmid1
benmid2
benhigh1
benhigh2);
175
set
rundat.
chld_
neo;

length
gender
$
50
categ
$
50
type
$
50;

gender
=
"
Children";
categ
=
"
Neo­
Natal
Mortality";
type
=
"
Sport";
exppop
=
expspt;
case1
=
red1r;
case2
=
red2r;
benlow1
=
ben1_
lr;
benlow2
=
ben2_
lr;
benmid1
=
mid1r;
benmid2
=
mid2r;
benhigh1
=
ben1_
hr;
benhigh2
=
ben2_
hr;
output;

gender
=
"
Children";
categ
=
"
Neo­
Natal
Mortality";
type
=
"
Sub";
exppop
=
expsub;
case1
=
red1s;
case2
=
red2s;
benlow1
=
ben1_
ls;
benlow2
=
ben2_
ls;
benmid1
=
mid1s;
benmid2
=
mid2s;
benhigh1
=
ben1_
hs;
benhigh2
=
ben2_
hs;
output;
run;

/*
End
NeoNatal
Mortality
Benefits
(
Sub
15).
sas
*/
176
Program
Output
177
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
95
NeoNatal
Mortality
Lead
Benefits
(
Sub
15)
National
Basis
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Label
Sum
                                                            
red1r
Number
of
Reductions
rec
(
Opt
1)
0.1088800
red2r
Number
of
Reductions
rec
(
Opt
2)
0.1579182
red1s
Number
of
Reductions
sub
(
Opt
1)
0.0442938
red2s
Number
of
Reductions
sub
(
Opt
2)
0.0642432
ben1_
lr
Sport
Low
(
Opt1)
272200.00
ben1_
ls
Sub
Low
(
Opt1)
110734.43
ben1_
hr
Sport
High
(
Opt1)
1480767.99
ben1_
hs
Sub
High
(
Opt1)
602395.29
ben2_
lr
Sport
Low
(
Opt2)
394795.39
ben2_
ls
Sub
Low
(
Opt2)
160608.04
ben2_
hr
Sport
High
(
Opt2)
2147686.91
ben2_
hs
Sub
High
(
Opt2)
873707.71
tot1_
l
Total
Benefits
­
Low
(
Opt
1)
382934.43
tot1_
h
Total
Benefits
­
High
(
Opt
1)
2083163.28
tot2_
l
Total
Benefits
­
Low
(
Opt
2)
555403.42
tot2_
h
Total
Benefits
­
High
(
Opt
2)
3021394.62
mid1r
707719.99
mid1s
287909.51
mid2r
1026468.01
mid2s
417580.89
mid1
Mid
Benefits
(
Opt1)
995629.51
mid2
Mid
Benefits
(
Opt2)
1444048.90
expspt
Exposed
Sport
Pop
(
15­
45yrs)
93849.71
expsub
Exposed
Sub
Pop
(
15­
45yrs)
4691.98
                                                            
178
Integrated
Exposure
Uptake
Biokinetic
Model
(
IEUBK)
The
model
used
a
range
of
lead
concentrations
in
fish
tissue
to
estimate
the
corresponding
blood
lead
level
associated
with
a
given
lead
concentration
in
the
fish
tissue.
179
"
Gen4.
txt"
Output
data
set
from
the
IEUBK
model
run.
See
DCN
37746
for
data
set.
180
"
Create
Lookup
Table.
sas"
Takes
the
output
file
from
the
IEUBK
model
run
and
generates
another
SAS
program
which
when
run
creates
a
format
used
to
transform
lead
concentrations
in
fish
tissue
into
blood
lead
levels.
181
/*
Program:
Create
Lookup
Table.
sas
Author:
Stuart
Smith,
Abt
Associates
Inc.
Date:
05/
18/
2000
Revised:
03/
28/
2002
This
program
reads
in
the
output
file
from
the
IEUBK
model
and
generates
a
SAS
data
set
with
concentrations
and
predicted
blood
lead
levels.
A
second
file
is
produced
which
contains
row
numbers
mapped
to
concentrations
values
for
lookup
by
the
benefits
program.
*/

libname
mpm2
'
H:\
ENVR\
Triton05
MP&
M\
Benefits\
MP&
M
Phase
2
Analysis\
Data\
SAS
Sets
for
Analysis';
/*
Static
Data
Sets
*/
filename
ubk
'
H:\
ENVR\
IEUBK\
Gen4.
txt';
/*
IEUBK
output
file
*/

data
readin;
infile
ubk
firstobs=
11;
input
id
1­
8
conc
56­
61
bloodPb
72­
76;
run;

data
mpm2.
ubklook4;
set
readin;
run;

/*
Generate
the
format
sas
file
*/
data
frmt;
set
readin;
file
'
H:\
ENVR\
IEUBK\
cnc_
fmt3.
sas';
put
conc
'
=
'
id;
run;

/*
End
Create
Lookup
Table.
sas
*/
182
"
Ublook4.
sd2"
Lookup
table
data
set.
See
DCN
37669
for
data
set.
183
"
Cnc_
fmt3.
sas"
Creates
a
library
format
for
transforming
lead
concentrations
in
fish
tissue
into
blood
lead
levels.
See
Appendix
B:
Programs
and
Estimation
Results
for
the
Proposed
Option
National
Benefits
Analysis,
DCN
20479
for
program
code.
184
"
Child
Lead
Benefits
(
Sub
16).
sas"
Using
data
from
the
lead
specific
data
set
and
results
of
the
IEUBK
lead
model,
child
blood
lead
levels
are
estimated
based
on
fish
consumption
in
recreational
and
subsistence
fishing
families.
185
/*
Program:
Child
Lead
Benefits
(
Sub
16).
sas
Author:
Stuart
Smith,
Abt
Associates,
Inc.
Date:
06/
12/
2000
Updated:
11/
27/
2002
This
program
estimates
benefits
from
reduced
consumption
of
fish
tissue
containing
lead.
Reductions
in
MP&
M
facility
lead
loadings
result
in
lower
probabilties
of
lead
related
illnesses.
*/

libname
library
'
H:\
ENVR\
Triton05
MP&
M\
Benefits\
MP&
M
Phase
2
Analysis\
Data\
SAS
Sets
for
Analysis\';

%
let
tbl_
path
=
%
sysfunc(
pathname(
rundat));

title2
'
Child
Lead
Related
Benefits
(
Sub
16)';

data
lead;
set
rundat.
lead;

/*
Combine
the
exposed
populations,
recreational
and
subsistence;
fresh
and
salt
water
*/

expspt
=
fwsptpop
+
swsptpop;
expsub
=
fwsubpop
+
swsubpop;

/*
Calculate
the
lead
concentration
in
fish
tissue
available
for
intake
*/
/*
=
instream
conc
(
ug/
L)*
biocon
factor
(
kg/
L)*
ingestion
rate(
kg/
day)
*/

bcncspt
=
cnc_
base*
49*
0.00603;
t1cncspt
=
cnc_
trt1*
49*
0.00603;
t2cncspt
=
cnc_
trt2*
49*
0.00603;

bcncsub
=
cnc_
base*
49*
0.03033;
t1cncsub
=
cnc_
trt1*
49*
0.03033;
t2cncsub
=
cnc_
trt2*
49*
0.03033;

keep
reachnmb
expspt
expsub
bcncspt
t1cncspt
t2cncspt
bcncsub
t1cncsub
t2cncsub
event_
wt;
run;

proc
print;
var
reachnmb
expspt
expsub
bcncspt
t1cncspt
t2cncspt
bcncsub
t1cncsub
t2cncsub
event_
wt;
title3
'
Exposed
Pops
and
Pb
conc
in
fish
tissue';
quit;

data
lead2;
set
lead;

/*
Determine
sport
blood
levels
*/
/*
cnc_
fmt.
is
a
special
format
which
translates
the
concentration
in
the
fish
tissue
into
an
index
value
for
a
lookup
table
of
corresponding
blood
lead
levels
*/
186
/*
Get
base
pbp
level
*/
lookup
=
input(
put(
round(
bcncspt,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
bspt
=
bloodPb;

/*
Get
t1
pbp
level
*/
lookup
=
input(
put(
round(
t1cncspt,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
1spt
=
bloodPb;

/*
Get
t2
pbp
level
*/
lookup
=
input(
put(
round(
t2cncspt,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
2spt
=
bloodPb;

/*
Determine
subsistence
blood
levels
*/
/*
Get
base
pbp
level
*/
lookup
=
input(
put(
round(
bcncsub,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
bsub
=
bloodPb;

/*
Get
t1
pbp
level
*/
lookup
=
input(
put(
round(
t1cncsub,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
1sub
=
bloodPb;

/*
Get
t2
pbp
level
*/
lookup
=
input(
put(
round(
t2cncsub,.
1),
cnc_
fmt.),
4.0);
if
lookup
=
9999
then
lookup
=
1;
set
mpm2.
ubklook4
point
=
lookup;
pb_
2sub
=
bloodPb;

run;

proc
print;
var
reachnmb
pb_
bspt
pb_
1spt
pb_
2spt
pb_
bsub
pb_
1sub
pb_
2sub;
title3
'
Blood
lead
levels';
quit;

data
lead3;
set
lead2;

kidspct
=
0.1177;
/*
Percent
of
kids
aged
0­
7,
U.
S.
Statistical
Abstract
1997
*/

popspt
=
expspt*
kidspct;
popsub
=
expsub*
kidspct;

GM1_
spt
=
pb_
bspt
­
pb_
1spt;
187
GM2_
spt
=
pb_
bspt
­
pb_
2spt;

GM1_
sub
=
pb_
bsub
­
pb_
1sub;
GM2_
sub
=
pb_
bsub
­
pb_
2sub;

/*
Avoided
lost
IQ
points
*/
avIQ1spt
=
(
GM1_
spt)*
0.25*(
popspt/
7);
avIQ2spt
=
(
GM2_
spt)*
0.25*(
popspt/
7);

avIQ1sub
=
(
GM1_
sub)*
0.25*(
popsub/
7);
avIQ2sub
=
(
GM2_
sub)*
0.25*(
popsub/
7);

IQ1_
spt
=
avIQ1spt
*
9419;
IQ2_
spt
=
avIQ2spt
*
9419;

IQ1_
sub
=
avIQ1sub
*
9419;
IQ2_
sub
=
avIQ2sub
*
9419;

/*
Fewer
children
with
IQ
below
70
*/

if
GM1_
spt
>
0
then
do;
num1spt
=
(
70
­
(
100
+
0.25*
GM1_
spt))/
15;
sv701spt
=
0.02275
­
cdf('
NORMAL',
num1spt);
b70_
1spt
=
sv701spt
*
(
popspt/
7)
*
58012;
end;
else
do;
sv701spt
=
0;
b70_
1spt
=
0;
end;

if
GM2_
spt
>
0
then
do;
num2spt
=
(
70
­
(
100
+
0.25*
GM2_
spt))/
15;
sv702spt
=
0.02275
­
cdf('
NORMAL',
num2spt);
b70_
2spt
=
sv702spt
*
(
popspt/
7)
*
58012;
end;
else
do;
sv702spt
=
0;
b70_
2spt
=
0;
end;

if
GM1_
sub
>
0
then
do;
num1sub
=
(
70
­
(
100
+
0.25*
GM1_
sub))/
15;
sv701sub
=
0.02275
­
cdf('
NORMAL',
num1sub);
b70_
1sub
=
sv701sub
*
(
popsub/
7)
*
58012;
end;
else
do;
sv701sub
=
0;
b70_
1sub
=
0;
end;

if
GM2_
sub
>
0
then
do;
num2sub
=
(
70
­
(
100
+
0.25*
GM2_
sub))/
15;
sv702sub
=
0.02275
­
cdf('
NORMAL',
num2sub);
188
b70_
2sub
=
sv702sub
*
(
popsub/
7)
*
58012;
end;
else
do;
sv702sub
=
0;
b70_
2sub
=
0;
end;

/*
Number
of
cases
of
IQ<
70
saved
*/

c70_
1spt
=
sv701spt
*
(
popspt/
7);
c70_
2spt
=
sv702spt
*
(
popspt/
7);
c70_
1sub
=
sv701sub
*
(
popsub/
7);
c70_
2sub
=
sv702sub
*
(
popsub/
7);

check1
=
pdf('
LOGNORMAL',
20,
pb_
bspt,
2.09);
check2
=
pdf('
LOGNORMAL',
20,
pb_
bsub,
2.09);

spt1_
20=(
check1
­
pdf('
LOGNORMAL',
20,
pb_
1spt,
2.09))*
0.2*(
popspt/
7)*
16485;
spt2_
20=(
check1
­
pdf('
LOGNORMAL',
20,
pb_
2spt,
2.09))*
0.2*(
popspt/
7)*
16485;
sub1_
20=(
check2
­
pdf('
LOGNORMAL',
20,
pb_
1sub,
2.09))*
0.2*(
popsub/
7)*
16485;
sub2_
20=(
check2
­
pdf('
LOGNORMAL',
20,
pb_
2sub,
2.09))*
0.2*(
popsub/
7)*
16485;

red201sp=(
check1
­
pdf('
LOGNORMAL',
20,
pb_
1spt,
2.09))*
0.2*(
popspt/
7);
red202sp=(
check1
­
pdf('
LOGNORMAL',
20,
pb_
2spt,
2.09))*
0.2*(
popspt/
7);
red201sb=(
check2
­
pdf('
LOGNORMAL',
20,
pb_
1sub,
2.09))*
0.2*(
popsub/
7);
red202sb=(
check2
­
pdf('
LOGNORMAL',
20,
pb_
2sub,
2.09))*
0.2*(
popsub/
7);
run;

proc
means
sum
data=
lead3;
var
avIQ1spt
avIQ1sub
avIQ2spt
avIQ2sub
IQ1_
spt
IQ1_
sub
IQ2_
spt
IQ2_
sub
b70_
1spt
b70_
1sub
b70_
2spt
b70_
2sub
sv701spt
sv702spt
sv701sub
sv702sub
c70_
1spt
c70_
2spt
c70_
1sub
c70_
2sub
spt1_
20
spt2_
20
sub1_
20
sub2_
20
red201sp
red202sp
red201sb
red202sb
popspt
popsub;
title3
'
Benefits
(
sample)';
quit;

proc
means
sum
data=
lead3;
weight
event_
wt;
var
avIQ1spt
avIQ1sub
avIQ2spt
avIQ2sub
IQ1_
spt
IQ1_
sub
IQ2_
spt
IQ2_
sub
b70_
1spt
b70_
1sub
b70_
2spt
b70_
2sub
sv701spt
sv702spt
sv701sub
sv702sub
c70_
1spt
c70_
2spt
c70_
1sub
c70_
2sub
spt1_
20
spt2_
20
sub1_
20
sub2_
20
red201sp
red202sp
red201sb
red202sb
popspt
popsub;
output
out
=
rundat.
chld_
ld
sum
=
avIQ1spt
avIQ1sub
avIQ2spt
avIQ2sub
IQ1_
spt
IQ1_
sub
IQ2_
spt
IQ2_
sub
189
b70_
1spt
b70_
1sub
b70_
2spt
b70_
2sub
sv701spt
sv702spt
sv701sub
sv702sub
c70_
1spt
c70_
2spt
c70_
1sub
c70_
2sub
spt1_
20
spt2_
20
sub1_
20
sub2_
20
red201sp
red202sp
red201sb
red202sb
popspt
popsub;
title3
'
Benefits
(
national)';
quit;

data
rundat.
chld_
tbl
(
keep=
gender
categ
type
exppop
case1
case2
benlow1
benlow2
benmid1
benmid2
benhigh1
benhigh2);
set
rundat.
chld_
ld;

length
gender
$
50
categ
$
50
type
$
50;

array
IQ[
8]
IQ1_
spt
avIQ1spt
IQ2_
spt
avIQ2spt
IQ1_
sub
avIQ1sub
IQ2_
sub
avIQ2sub;
array
RD[
8]
b70_
1spt
c70_
1spt
b70_
2spt
c70_
2spt
b70_
1sub
c70_
1sub
b70_
2sub
c70_
2sub;
array
PB[
8]
spt1_
20
red201sp
spt2_
20
red202sp
sub1_
20
red201sb
sub2_
20
red202sb;

gender
=
"
Children";

if
_
n_
=
1
then
do;
do
i
=
1
to
6;
if
i
=
1
then
do;
categ
=
"
Avoided
IQ
Loss";
type
=
"
Sport";
case1
=
IQ[
2];
case2
=
IQ[
4];
benlow1
=
IQ[
1];
benlow2
=
IQ[
3];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
if
i
=
2
then
do;
categ
=
"
Avoided
IQ
Loss";
type
=
"
Sub";
case1
=
IQ[
6];
case2
=
IQ[
8];
benlow1
=
IQ[
5];
benlow2
=
IQ[
7];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
if
i
=
3
then
do;
categ
=
"
Reduced
IQ
<
70";
type
=
"
Sport";
case1
=
RD[
2];
case2
=
RD[
4];
benlow1
=
RD[
1];
benlow2
=
RD[
3];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
if
i
=
4
then
do;
categ
=
"
Reduced
IQ
<
70";
type
=
"
Sub";
190
case1
=
RD[
6];
case2
=
RD[
8];
benlow1
=
RD[
5];
benlow2
=
RD[
7];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
if
i
=
5
then
do;
categ
=
"
Reduced
Pb
>
20";
type
=
"
Sport";
case1
=
PB[
2];
case2
=
PB[
4];
benlow1
=
PB[
1];
benlow2
=
PB[
3];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
if
i
=
6
then
do;
categ
=
"
Reduced
Pb
>
20";
type
=
"
Sub";
case1
=
PB[
6];
case2
=
PB[
8];
benlow1
=
PB[
5];
benlow2
=
PB[
7];
benmid1
=
benlow1;
benmid2
=
benlow2;
benhigh1
=
benlow1;
benhigh2
=
benlow2;
output;
end;
end;
end;
run;

data
rundat.
chld_
all;
set
rundat.
neo_
tbl
rundat.
chld_
tbl;
run;

proc
summary
data=
rundat.
chld_
all;
var
benlow1
benmid1
benhigh1
benlow2
benmid2
benhigh2;
output
out
=
rundat.
chld_
tot
(
drop=
_
type_
_
freq_)
sum
=
benlow1
benmid1
benhigh1
benlow2
benmid2
benhigh2;
run;

data
exlchld
(
keep=
titlecol
col1­
col11
label="
Child
Lead
Benefits
­
Weighted");
set
rundat.
chld_
all
rundat.
chld_
tot;
length
titlecol
$
50
col1
­
col11
$
128;

if
_
n_
=
1
then
do;
titlecol
=
"
Table:
Child
Lead
Benefits,
Weighted
Estimates";
output;
titlecol
=
"
Run
on
&
sysdate";
output;
titlecol
=
"
";
output;
col4
=
"
Option
1";
col8
=
"
Option
2";
output;
col2
=
"
Category";
col3
=
"
Type";
col4
=
"
Reduced
Cases";
col5
=
"
Low";
col6
=
"
Mid";
191
col7
=
"
High";
col8
=
"
Reduced
Cases";
col9
=
"
Low";
col10
=
"
Mid";
col11
=
"
High";
output;
end;

if
_
n_
=
9
then
col1
=
"
Total";
else
col1
=
gender;

col2
=
trim(
categ);
col3
=
trim(
type);
col4
=
trim(
case1);
col5
=
trim(
benlow1);
col6
=
trim(
benmid1);
col7
=
trim(
benhigh1);
col8
=
trim(
case2);
col9
=
trim(
benlow2);
col10
=
trim(
benmid2);
col11
=
trim(
benhigh2);
output;
run;

proc
export
data
=
exlchld
outfile
=
"
C:\
Output
Files\
Child
Lead.
xls"
dbms=
excel2000
replace;
run;

/*
End
Child
Lead
Benefits
(
Sub
16).
sas
*/
192
Program
Output
193
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
99
Child
Lead
Related
Benefits
(
Sub
16)
Benefits
(
national)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Sum
                        
avIQ1spt
0
avIQ1sub
31.9909293
avIQ2spt
0
avIQ2sub
36.8381039
IQ1_
spt
0
IQ1_
sub
301322.56
IQ2_
spt
0
IQ2_
sub
346978.10
b70_
1spt
0
b70_
1sub
6636.58
b70_
2spt
0
b70_
2sub
7635.65
sv701spt
0
sv702spt
0
sv701sub
0.0069909
sv702sub
0.0090608
c70_
1spt
0
c70_
2spt
0
c70_
1sub
0.1144002
c70_
2sub
0.1316220
spt1_
20
0
spt2_
20
0
sub1_
20
­
447.4330014
sub2_
20
­
488.6185253
red201sp
0
red202sp
0
red201sb
­
0.0271418
red202sb
­
0.0296402
popspt
49026.28
popsub
2451.05
                        
194
Recreational
Benefits
195
Recreational
Fishing
196
"
Human
Health
AWQC
Exceedences
(
Sub
2).
sas"
The
number
of
human
health
based
AWQC
exceedences
for
each
RF1
reach
are
estimated
and
output
to
a
data
set.
See
Human
Health
Benefits
Cancer
from
Fish
Consumption
for
program
code
and
output.
197
"
Aquatic
Life
AWQC
Exceedences
(
Sub
3).
sas"
The
number
of
aquatic
life
based
AWQC
exceedences
for
each
RF1
reach
are
estimated
and
output
to
a
data
set.
198
/*
Program:
Aquatic
Life
AWQC
(
Sub
3).
sas
Author:
Stuart
Smith
(
based
on
Elena
Besedin,
Rob
Sartain)
,
Abt.
Associates,
Inc.
Date
Modified:
10/
31/
2002
(
1)
IDENTIFIES
MP&
M
DISCHARGE
REACHES
IN
WHICH
POLLUTANT
CONCENTRATIONS
EXCEED
THE
AWQC
VALUES
FOR
AQUATIC
LIFE;
(
2)
EXTRAPOLATES
THE
RESULTS
TO
THE
POPULATION
LEVEL
BY
USING
DIFFERENTIAL
WEIGHTING.

THE
TOTAL
NUMBER
OF
EXCEEDENCE
EVENTS
IS
THEN
ESTIMATED
BY
SUMMING
OVER
THE
RESULTS
FROM
THE
FIRST
AND
SECOND
RUN.
*/

title2
'
Output
from
Aquatic
Life
(
Sub
3)';

/*
DISCHARGE
EVENTS,
ASSOCIATED
LOADINGS
AND
WEIGHTS
HAVE
BEEN
CALCULATED
BY
ANOTHER
PROGRAM.
ALL
UNIT
CONVERSIONS
HAVE
OCCURRED.
SEE
"
Generate
Discharge
Events
and
Loadings
(
Sub
1).
sas"
FOR
SPECIFIC
PROGRAM
INFORMATION.
*/

/*
PULL
TOGETHER
THE
EVENTS
AND
CHEMICAL
SPECIFIC
PARAMETERS
*/

Proc
sort
data
=
mpm2.
chemdata;
by
casno;
quit;

Proc
sql;
create
table
stream
as
select
reachnmb,
event_
ld.
casno,
evnt_
flo,
base,
treat1,
treat2,
event_
wt,
segl,
mnflo
as
mean_
flo,
svten
as
flo7q10,
acute
as
adf,
chronic
as
cdf,
dcp,
wbtype,
biodegre
as
decay,
surface,
drainage,
volume
from
rundat.
event_
ld,
mpm2.
chemdata
where
event_
ld.
casno
=
chemdata.
casno
order
by
reachnmb,
casno;
quit;

/******************************************************************************/
/*
*/
/*
FRESHWATER
REACHES
­
CONCENTRATIONS
AND
AWQC
EXCEEDENCES
*/
/*
*/
/******************************************************************************/

/*
Calculate
chemical
concentrations
in
lakes
*/

data
Lake;
set
stream;
if
wbtype
=
'
lake';
if
decay;

keep
reachnmb
casno
base
treat1
treat2
adf
cdf
restime
decay
event_
wt;

/*
drainage
area
and
surface
area
available
*/
if
surface
^=
.
and
drainage
^=
.
then
do;
num
=
drainage;
/*
drainage
area
in
sq.
miles
*/
199
denom
=
surface*
0.0015625;
/*
convert
surface
area
from
acres
to
sq.
miles
*/
logdays
=
4.077
­
(
1.177*
log10(
num/
denom));
restime
=
10**
logdays;
time1
=
restime*
24;
/*
CONVERT
HYDRAULIC
TIME
FROM
DAYS
INTO
HOURS
*/
cnc_
aqcb
=
(
base/
evnt_
flo)/(
1+
time1*
decay);
cnc_
aqc1
=
(
treat1/
evnt_
flo)/(
1+
time1*
decay);
cnc_
aqc2
=
(
treat2/
evnt_
flo)/(
1+
time1*
decay);
good
=
1;
end;

/*
missing
area
measure(
s),
do
end­
of­
pipe
value
=
loadings/
flow
*/
else
do;
cnc_
aqcb
=
(
base/
evnt_
flo);
cnc_
aqc1
=
(
treat1/
evnt_
flo);
cnc_
aqc2
=
(
treat2/
evnt_
flo);
end;

cnc_
aqab
=
cnc_
aqcb;
cnc_
aqa1
=
cnc_
aqc1;
cnc_
aqa2
=
cnc_
aqc2;

keep
reachnmb
casno
cnc_
aqcb
cnc_
aqc1
cnc_
aqc2
cnc_
aqab
cnc_
aqa1
cnc_
aqa2
event_
wt;
run;

/*
Great
Lakes
reaches
*/

data
BigLake;
set
stream;
if
wbtype
=
'
glake';

if
cdf
^=.
then
do;
cnc_
aqcb
=
base/(
evnt_
flo*
cdf);
cnc_
aqc1
=
treat1/(
evnt_
flo*
cdf);
cnc_
aqc2
=
treat2/(
evnt_
flo*
cdf);
end;
else
if
dcp^=.
then
do;
benchld=
10000*
907.2;
/*
benchmark
load
is
10000
per
year;
907.2
is
coversion
factor
from
kg
to
ton*/
cnc_
aqcb=((
base/
10**
9)*
365*
dcp)/
benchld;
cnc_
aqc1=((
treat1/
10**
9)*
365*
dcp)/
benchld;
cnc_
aqc2=((
treat2/
10**
9)*
365*
dcp)/
benchld;
end;
else
do;
cnc_
aqcb
=
.;
cnc_
aqc1
=
.;
cnc_
aqc2
=
.;
end;
cnc_
aqab
=
cnc_
aqcb;
cnc_
aqa1
=
cnc_
aqc1;
cnc_
aqa2
=
cnc_
aqc2;

if
cnc_
aqcb=.
then
delete;
200
keep
reachnmb
casno
cnc_
aqcb
cnc_
aqc1
cnc_
aqc2
cnc_
aqab
cnc_
aqa1
cnc_
aqa2
event_
wt;
run;

/*
Stream
reaches
*/
/*
Calculate
1Q10.
Calculate
in­
stream
concentrations
based
on
7Q10
and
1Q10*/

data
stream2;
set
stream;
if
wbtype
=
'
stream';

/*
Account
for
missing
7q10
data,
averge
*/
if
flo7q10
=
.
or
flo7q10
=
0
then
flo7q10
=
909.7836;
if
mean_
flo
=
.
or
mean_
flo
=
0
then
mean_
flo
=
6711.2411;

/*
Convert
flo7q10
from
ft3/
sec
to
l/
day
using
conversion
factor
2446576
*/

svnq10
=
flo7q10*
2446576;

/*
Calculate
1q10
flow
*/

flo1q10=
0.843*(
svnq10**
0.993);

/*
all
flows
in
l/
day
*/

cnc_
aqcb
=
base/(
evnt_
flo+
svnq10);
cnc_
aqc1
=
treat1/(
evnt_
flo+
svnq10);
cnc_
aqc2
=
treat2/(
evnt_
flo+
svnq10);
cnc_
aqab
=
base/(
evnt_
flo+
flo1q10);
cnc_
aqa1
=
treat1/(
evnt_
flo+
flo1q10);
cnc_
aqa2
=
treat2/(
evnt_
flo+
flo1q10);

keep
reachnmb
casno
cnc_
aqcb
cnc_
aqc1
cnc_
aqc2
cnc_
aqab
cnc_
aqa1
cnc_
aqa2
event_
wt;
run;

proc
sort
data=
stream2;
by
reachnmb
casno;
quit;

/*
Combine
small
and
large
lakes
with
streams
*/

data
freshwat;
set
stream2
BigLake
Lake;
run;

proc
sort
data=
freshwat;
by
reachnmb
casno;
quit;

/*
Compare
in­
stream
pollutant
concentrations
with
freshwater
AWQCs
for
aquatic
life
*/
proc
sql;

/*
Merge
in
freshwater
aquatic
life
limits.
*/
201
create
table
awqc1
as
select
reachnmb,
freshwat.
casno,
chemical,
cnc_
aqcb,
cnc_
aqc1,
cnc_
aqc2,
cnc_
aqab,
cnc_
aqa1,
cnc_
aqa2,
faval__
i
as
faval,
fcval__
i
as
fcval,
event_
wt
from
freshwat,
mpm2.
chemdata
where
freshwat.
casno
=
chemdata.
casno
order
by
reachnmb,
casno;
quit;

/*
Compare
concentrations
to
awqc
limits
and
determine
exceedences
*/

data
stream4;
set
awqc1;

if
cnc_
aqcb
>
fcval
and
fcval
>
0
then
chronicb
=
1;
if
cnc_
aqc1
>
fcval
and
fcval
>
0
then
chronic1
=
1;
if
cnc_
aqc2
>
fcval
and
fcval
>
0
then
chronic2
=
1;
if
cnc_
aqcb
>
faval
and
faval
>
0
then
acuteb
=
1;
if
cnc_
aqc1
>
faval
and
faval
>
0
then
acute1
=
1;
if
cnc_
aqc2
>
faval
and
faval
>
0
then
acute2
=
1;

keep
reachnmb
casno
chemical
chronicb
chronic1
chronic2
acuteb
acute1
acute2
event_
wt;
run;

/***********************************************************************/
/*
END
FRESHWATER
REACHES
*/
/***********************************************************************/

/******************************************************************************/
/*
*/
/*
SALTWATER
REACHES
­
CONCENTRATIONS
AND
AWQC
EXCEEDENCES
*/
/*
*/
/******************************************************************************/

/*
Calculate
pollutant
concentrations
in
Oceans
*/

data
marine;
set
stream;
if
wbtype
=
'
bay';

if
cdf
^=.
then
do;
cnc_
aqcb
=
base/(
evnt_
flo*
cdf);
cnc_
aqc1
=
treat1/(
evnt_
flo*
cdf);
cnc_
aqc2
=
treat2/(
evnt_
flo*
cdf);
end;
else
if
dcp
^=.
then
do;
benchld=
10000*
907.2;
/*
benchmark
load
is
10000
per
year;
907.2
is
coversion
factor
from
kg
to
ton
*/

/*
facility
loadings
were
given
in
ug/
day
*/

cnc_
aqcb=((
base/
10**
9)*
365*
dcp)/
benchld;
202
cnc_
aqc1=((
treat1/
10**
9)*
365*
dcp)/
benchld;
cnc_
aqc2=((
treat2/
10**
9)*
365*
dcp)/
benchld;
end;
else
do;
cnc_
aqcb
=
.;
cnc_
aqc1
=
.;
cnc_
aqc2
=
.;
end;

cnc_
aqab
=
cnc_
aqcb;
cnc_
aqa1
=
cnc_
aqc1;
cnc_
aqa2
=
cnc_
aqc2;

if
cnc_
aqcb=.
then
delete;
keep
reachnmb
casno
cnc_
aqcb
cnc_
aqc1
cnc_
aqc2
cnc_
aqab
cnc_
aqa1
cnc_
aqa2
event_
wt;
run;

/*
Get
AWQC
limits
for
marine
environment
*/

proc
sql;

/*
Merge
in
Aquatic
life
saltwater
awqc
limits
*/

create
table
marine1
as
select
reachnmb,
marine.
casno,
chemical,
cnc_
aqcb,
cnc_
aqc1,
cnc_
aqc2,
cnc_
aqab,
cnc_
aqa1,
cnc_
aqa2,
event_
wt,
maval__
i
as
maval,
mcval__
i
as
mcval,
count(
distinct
reachnmb)
as
reaches
from
marine,
mpm2.
chemdata
where
marine.
casno
=
chemdata.
casno
order
by
reachnmb,
casno;
quit;

/*
Compare
concentrations
to
awqc
limits
and
determine
exceedences
*/

data
marine2;
set
marine1;
if
cnc_
aqcb
>
mcval
and
mcval
>
0
then
chronicb
=
1;
if
cnc_
aqc1
>
mcval
and
mcval
>
0
then
chronic1
=
1;
if
cnc_
aqc2
>
mcval
and
mcval
>
0
then
chronic2
=
1;
if
cnc_
aqcb
>
maval
and
maval
>
0
then
acuteb
=
1;
if
cnc_
aqc1
>
maval
and
maval
>
0
then
acute1
=
1;
if
cnc_
aqc2
>
maval
and
maval
>
0
then
acute2
=
1;

keep
reachnmb
casno
chemical
chronicb
chronic1
chronic2
acuteb
acute1
acute2
reaches
event_
wt;
run;

/*
Combine
exceedences
from
fresh
and
saltwater
*/
data
allreach;
set
marine2
stream4;

keep
reachnmb
casno
chemical
chronicb
chronic1
chronic2
acuteb
acute1
acute2
reaches
event_
wt;
run;
203
proc
sql;

/*
Estimate
AWQCs
exceedance
by
chemicals
on
a
sample
and
national
basis*/
/*
awt
replaced
with
1
*/
create
table
chm_
awqc
as
select
chemical,
sum(
chronicb)
as
crb,
sum(
chronic1)
as
cr1,
sum(
chronic2)
as
cr2,
sum(
acuteb)
as
acb,
sum(
acute1)
as
ac1,
sum(
acute2)
as
ac2,
sum(
chronicb*
event_
wt*
1)
as
ntlbw,
sum(
chronic1*
event_
wt*
1)
as
ntl1w,
sum(
chronic2*
event_
wt*
1)
as
ntl2w,
sum(
acuteb*
event_
wt*
1)
as
ntlbo,
sum(
acute1*
event_
wt*
1)
as
ntl1o,
sum(
acute2*
event_
wt*
1)
as
ntl2o,
case
when
sum(
chronicb)
>=
1
then
1
else
0
end
as
chemcrb,
case
when
sum(
chronic1)
>=
1
then
1
else
0
end
as
chemcr1,
case
when
sum(
chronic2)
>=
1
then
1
else
0
end
as
chemcr2,
case
when
sum(
acuteb)
>=
1
then
1
else
0
end
as
chemacb,
case
when
sum(
acute1)
>=
1
then
1
else
0
end
as
chemac1,
case
when
sum(
acute2)
>=
1
then
1
else
0
end
as
chemac2
from
allreach
where
chronicb>
0
or
acuteb>
0
group
by
chemical
order
by
chemical;

/*
Estimate
aquatic
life
AWQC
Exceedences
by
option/
reach
*/

create
table
aq_
life2
as
select
allreach.
reachnmb,
sum(
chronicb)
as
test_
cb,
sum(
chronic1)
as
test_
c1,
sum(
chronic2)
as
test_
c2,
sum
(
acuteb)
as
test_
ab,
sum(
acute1)
as
test_
a1,
sum
(
acute2)
as
test_
a2,
sum(
chronicb*
event_
wt*
1)
as
ntlb1,
sum(
chronic1*
event_
wt*
1)
as
ntnl11,
sum(
chronic2*
event_
wt*
1)
as
ntnl21,
sum(
acuteb
*
event_
wt*
1)
as
ntlb2,
sum(
acute1*
event_
wt*
1)
as
ntnl12,
sum(
acute2*
event_
wt*
1)
as
ntnl22,
case
when
sum(
chronicb)
>=
1
then
1
else
0
end
as
rchcrb,
case
when
sum(
chronic1)
>=
1
then
1
else
0
end
as
rchcr1,
case
when
sum(
chronic2)
>=
1
then
1
else
0
end
as
rchcr2,
case
when
sum(
acuteb)
>=
1
then
1
else
0
end
as
rchacb,
case
when
sum(
acute1)
>=
1
then
1
else
0
end
as
rchac1,
case
when
sum(
acute2)
>=
1
then
1
else
0
end
as
rchac2,
sum(
event_
wt)
as
reachwt
from
allreach
where
chronicb>
0
or
acuteb>
0
group
by
reachnmb
order
by
reachnmb;

/*
CALCULATE
THE
TOTAL
NUMBER
OF
AFFECTED
REACHES
ON
A
NATIONAL
BASIS
awt
replaced
with
a
1
*/
create
table
rch1
as
select
204
sum(
event_
wt)
as
wght1,
reachnmb,
chemical
from
allreach
where
chronicb
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
chrnb
as
select
max(
wght1)
as
chrnc_
b,
reachnmb
from
rch1
group
by
reachnmb
order
by
reachnmb;

create
table
rch2
as
select
sum(
event_
wt)
as
wght2,
reachnmb,
chemical
from
allreach
where
acuteb
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
actb
as
select
max(
wght2)
as
acutr_
b,
reachnmb
from
rch2
group
by
reachnmb
order
by
reachnmb;

create
table
rch3
as
select
sum(
event_
wt)
as
wght3,
reachnmb,
chemical
from
allreach
where
chronic1
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
chrn1
as
select
max(
wght3)
as
chrnc_
1,
reachnmb
from
rch3
group
by
reachnmb
order
by
reachnmb;

create
table
rch4
as
select
sum(
event_
wt)
as
wght4,
reachnmb,
chemical
from
allreach
where
acute1
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
act1
as
select
max(
wght4)
as
acutr_
1,
reachnmb
from
rch4
group
by
reachnmb
order
by
reachnmb;
205
create
table
rch5
as
select
sum(
event_
wt)
as
wght5,
reachnmb,
chemical
from
allreach
where
chronic2
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
chrn2
as
select
max(
wght5)
as
chrnc_
2,
reachnmb
from
rch5
group
by
reachnmb
order
by
reachnmb;

create
table
rch6
as
select
sum(
event_
wt)
as
wght6,
reachnmb,
chemical
from
allreach
where
acute2
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
act2
as
select
max(
wght6)
as
acutr_
2,
reachnmb
from
rch6
group
by
reachnmb
order
by
reachnmb;

create
table
rch7
as
select
sum(
event_
wt)
as
wght7,
reachnmb,
chemical
from
allreach
where
chronicb
>
0
or
acuteb
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
multi_
b
as
select
max(
wght7)
as
multi_
b,
reachnmb
from
rch7
group
by
reachnmb
order
by
reachnmb;

create
table
rch8
as
select
sum(
event_
wt)
as
wght8,
reachnmb,
chemical
from
allreach
where
chronic1
>
0
or
acute1
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
multi_
1
as
select
max(
wght8)
as
multi_
1,
reachnmb
from
rch8
206
group
by
reachnmb
order
by
reachnmb;

create
table
rch9
as
select
sum(
event_
wt)
as
wght9,
reachnmb,
chemical
from
allreach
where
chronic2
>
0
or
acute2
>
0
group
by
reachnmb,
chemical
order
by
reachnmb,
chemical;

create
table
multi_
2
as
select
max(
wght9)
as
multi_
2,
reachnmb
from
rch9
group
by
reachnmb
order
by
reachnmb;
quit;

data
rchwts;
merge
chrnb
actb
chrn1
act1
chrn2
act2
multi_
b
multi_
1
multi_
2;
by
reachnmb;

if
chrnc_
b=.
then
chrnc_
b=
0;
if
chrnc_
1=.
then
chrnc_
1=
0;
if
chrnc_
2=.
then
chrnc_
2=
0;
if
acutr_
b=.
then
acutr_
b=
0;
if
acutr_
1=.
then
acutr_
1=
0;
if
acutr_
2=.
then
acutr_
2=
0;
if
multi_
b=.
then
multi_
b=
0;
if
multi_
1=.
then
multi_
1=
0;
if
multi_
2=.
then
multi_
2=
0;
run;

quit;

/*
Output
the
affected
reaches
for
the
recreational
benefits
programs
*/

data
rundat.
rec_
al;
set
allreach;
if
chronicb
>
0
or
acuteb
>
0;
drop
reaches;
run;

/*
Print
out
results
*/

proc
means
sum
data=
aq_
life2;
var
test_
cb
test_
c1
test_
c2
test_
ab
test_
a1
test_
a2
ntlb1
ntnl11
ntnl21
ntlb2
ntnl12
ntnl22
rchcrb
rchcr1
rchcr2
rchacb
rchac1
rchac2;
label
reachnmb
='
Reach
ID'
test_
cb
=
'
Chronic
AWQC
Exeedences
(
Smpl/
BL)'
test_
c1
=
'
Chronic
AWQC
Exeedences
(
Smpl/
Opt
1)'
test_
c2
=
'
Chronic
AWQC
Exeedences
(
Smpl/
Opt
2)'
test_
ab
=
'
Acute
AWQC
Exeedences
(
Smpl/
BL)'
207
test_
a1
=
'
Acute
AWQC
Exeedences
(
Smpl/
Opt
1)'
test_
a2
=
'
Acute
AWQC
Exeedences
(
Smpl/
Opt
2)'
ntlb1
=
'
Chronic
AWQC
Exeedences
(
Ntl/
BL)'
ntnl11
=
'
Chronic
AWQC
Exeedences
(
Ntl/
OPt
1)'
ntnl21
=
'
Chronic
AWQC
Exeedences
(
Ntl/
OPt
2)'
ntlb2
=
'
Acute
AWQC
Exeedences
(
Ntl/
Bl)'
ntnl12
=
'
Acute
AWQC
Exeedences
(
Ntl/
Opt
1)'
ntnl22
=
'
Acute
AWQC
Exeedences
(
Ntl/
Opt
2)'
rchcrb
=
'
Number
of
Chronic
reaches
(
BL)'
rchcr1
=
'
Number
of
Chronic
reaches
(
Opt
1)'
rchcr2
=
'
Number
of
Chronic
reaches
(
Opt
2)'
rchacb
=
'
Number
of
Acute
reaches
(
BL)'
rchac1
=
'
Number
of
Acute
reaches
(
Opt
1)'
rchac2
=
'
Number
of
Acute
reaches
(
Opt
2)';

title3
'
Aquatic
Life
AWQC
Exceedences
by
Option
(
Sample/
National)';

proc
means
sum
data=
chm_
awqc;
var
crb
cr1
cr2
acb
ac1
ac2
ntlbw
ntl1w
ntl2w
ntlbo
ntl1o
ntl2o
chemcrb
chemcr1
chemcr2
chemacb
chemac1
chemac2;
label
crb
=
'
Chronic
AWQC
Exeedences
(
Smpl/
BL)'
cr1
=
'
Chronic
AWQC
Exeedences
(
Smpl/
Opt
1)'
cr2
=
'
Chronic
AWQC
Exeedences
(
Smpl/
Opt
2)'
acb
=
'
Acute
AWQC
Exeedences
(
Smpl/
BL)'
ac1
=
'
Acute
AWQC
Exeedences
(
Smpl/
Opt
1)'
ac2
=
'
Acute
AWQC
Exeedences
(
Smpl/
Opt
2)'
ntlbw
=
'
Chronic
AWQC
Exeedences
(
Ntl/
BL)'
ntl1w
=
'
Chronic
AWQC
Exeedences
(
Ntl/
OPt
1)'
ntl2w
=
'
Chronic
AWQC
Exeedences
(
Ntl/
OPt
2)'
ntlbo
=
'
Acute
AWQC
Exeedences
(
Ntl/
Bl)'
ntl1o
=
'
Acute
AWQC
Exeedences
(
Ntl/
Opt
1)'
ntl2o
=
'
Acute
AWQC
Exeedences
(
Ntl/
Opt
2)'
chemcrb
=
'
Number
of
Chronic
pol.
(
BL)'
chemcr1
=
'
Number
of
Chronic
pol.
(
Opt
1)'
chemcr2
=
'
Number
of
Chronic
pol.
(
Opt
2)'
chemacb
=
'
Number
of
Acute
pol.
(
BL)'
chemac1
=
'
Number
of
Acute
pol.
(
Opt
1)'
chemac2
=
'
Number
of
Acute
pol.
(
Opt
2)';

title3
'
AWQC
Exceedances
by
pollutant/
option
(
Sample
/
National)';

proc
means
sum
data=
rchwts;
var
chrnc_
b
chrnc_
1
chrnc_
2
acutr_
b
acutr_
1
acutr_
2
multi_
b
multi_
1
multi_
2;
label
chrnc_
b
=
'#
of
Reaches
Exceeding
Chronic
AWQC
(
Bl)
'
chrnc_
1
=
'#
of
Reaches
Exceeding
Chronic
AWQC
(
Opt
1)'
chrnc_
2
=
'#
of
Reaches
Exceeding
Chronic
AWQC
(
Opt
2)'
acutr_
b
=
'#
of
Reaches
Exceeding
Acute
AWQC
(
Bl)'
acutr_
1
=
'#
of
Reaches
Exceeding
Acute
AWQC
(
Opt
1)'
acutr_
2
=
'#
of
Reaches
Exceeding
Acute
AWQC
(
Opt
2)'
multi_
b
=
'
Either
(
BL)'
multi_
1
=
'
Either
(
Opt
1)'
multi_
2
=
'
Either
(
Opt
2)';
title3
'
IL
Reaches
Exceeding
AWQC
Limits
for
Aquatic
Life
(
National)';
208
run;

/*
Sum
the
number
of
exceedence
events
by
reach
and
event
weight
*/
proc
sql;
create
table
preme
as
select
reachnmb,
event_
wt,
sum(
chronicb)
as
testb1,
sum(
chronic1)
as
test1_
1,
sum(
chronic2)
as
test2_
1,
sum(
acuteb)
as
testb2,
sum(
acute1)
as
test1_
2,
sum(
acute2)
as
test2_
2
from
allreach
group
by
reachnmb,
event_
wt
order
by
reachnmb,
event_
wt;
quit;

/*
Determine
the
number
of
reaches
which
eliminate
all
exceedences
and
the
number
of
reaches
which
only
partially
reduce
*/

data
reduc;
set
preme;
if
testb1=.
then
testb1=
0;
if
test1_
1=.
then
test1_
1=
0;
if
test2_
1=.
then
test2_
1=
0;
if
testb2=.
then
testb2=
0;
if
test1_
2=.
then
test1_
2=
0;
if
test2_
2=.
then
test2_
2=
0;

if
testb1
>=
1
and
test1_
1
=
0
then
elim1_
1
=
1;
else
elim1_
1=
0;
if
testb1
>=
1
and
test2_
1
=
0
then
elim2_
1
=
1;
else
elim2_
1=
0;
if
testb2
>=
1
and
test1_
2
=
0
then
elim1_
2
=
1;
else
elim1_
2=
0;
if
testb2
>=
1
and
test2_
2
=
0
then
elim2_
2
=
1;
else
elim2_
2=
0;

if
testb1
>=
1
and
test1_
1
<
testb1
and
elim1_
1
=
0
then
red1_
1
=
1;
else
red1_
1=
0;
if
testb1
>=
1
and
test2_
1
<
testb1
and
elim2_
1
=
0
then
red2_
1
=
1;
else
red2_
1=
0;
if
testb2
>=
1
and
test1_
2
<
testb2
and
elim1_
2
=
0
then
red1_
2
=
1;
else
red1_
2=
0;
if
testb2
>=
1
and
test2_
2
<
testb2
and
elim2_
2
=
0
then
red2_
2
=
1;
else
red2_
2=
0;

if
elim1_
1=
1
or
elim1_
2=
1
then
elim1
=
1;
else
elim1
=
0;
if
elim2_
1=
1
or
elim2_
2=
1
then
elim2
=
1;
else
elim2
=
0;
if
red1_
1
=
1
or
red1_
2
=
1
then
reduc1
=
1;
else
reduc1
=
0;
if
red2_
1
=
1
or
red2_
2
=
1
then
reduc2
=
1;
else
reduc2
=
0;
run;

/*
Print
results
*/

proc
means
sum
data=
reduc;
var
elim1_
1
elim2_
1
elim1_
2
elim2_
2
red1_
1
red2_
1
red1_
2
red2_
2
reduc1
reduc2;
label
elim1_
1
=
'
Elim.
Chronic
(
Opt1)'
elim2_
1
=
'
Elim.
Chronic
(
Opt2)'
209
elim1_
2
=
'
Elim.
Acute
(
Opt1)'
elim2_
2
=
'
Elim.
Acute
(
Opt2)'
red1_
1
=
'
Redu.
Chronic
(
Opt1)'
red2_
1
=
'
Redu.
Chronic
(
Opt2)'
red1_
2
=
'
Redu.
Acute
(
Opt1)'
red2_
2
=
'
Redu.
Acute
(
Opt2)'
reduc1
=
'
Total
Reaches
Reducing
(
Opt1)'
reduc2
=
'
Total
Reaches
Reducing
(
Opt2)';
title3
'
AWQC
Eliminations
and
Reductions
(
Sample)';
title4
'
Number
of
Sample
Reaches';
quit;

proc
means
sum
data=
reduc;
weight
event_
wt;
var
elim1_
1
elim2_
1
elim1_
2
elim2_
2
red1_
1
red2_
1
red1_
2
red2_
2
reduc1
reduc2;
label
elim1_
1
=
'
Elim.
Chronic
(
Opt1)'
elim2_
1
=
'
Elim.
Chronic
(
Opt2)'
elim1_
2
=
'
Elim.
Acute
(
Opt1)'
elim2_
2
=
'
Elim.
Acute
(
Opt2)'
red1_
1
=
'
Redu.
Chronic
(
Opt1)'
red2_
1
=
'
Redu.
Chronic
(
Opt2)'
red1_
2
=
'
Redu.
Acute
(
Opt1)'
red2_
2
=
'
Redu.
Acute
(
Opt2)'
reduc1
=
'
Total
Reaches
Reducing
(
Opt1)'
reduc2
=
'
Total
Reaches
Reducing
(
Opt2)';
title3
'
AWQC
Eliminations
and
Reductions
(
National)';
title4
'
Number
of
National
Reaches';
quit;

proc
print
data=
aq_
life2;
quit;
proc
print
data=
chm_
awqc;
title3
'
Exceedences
by
Chemical';
quit;

/*
End
Aquatic
Life
AWQC
(
Sub
3).
sas
*/
210
Program
Output
211
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
9
Output
from
Aquatic
Life
(
Sub
3)
Aquatic
Life
AWQC
Exceedences
by
Option
(
Sample/
National)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Label
Sum
                                                                
test_
cb
Chronic
AWQC
Exeedences
(
Smpl/
BL)
23.0000000
test_
c1
Chronic
AWQC
Exeedences
(
Smpl/
Opt
1)
6.0000000
test_
c2
Chronic
AWQC
Exeedences
(
Smpl/
Opt
2)
.
test_
ab
Acute
AWQC
Exeedences
(
Smpl/
BL)
4.0000000
test_
a1
Acute
AWQC
Exeedences
(
Smpl/
Opt
1)
.
test_
a2
Acute
AWQC
Exeedences
(
Smpl/
Opt
2)
.
ntlb1
Chronic
AWQC
Exeedences
(
Ntl/
BL)
423.3144485
ntnl11
Chronic
AWQC
Exeedences
(
Ntl/
OPt
1)
52.8929261
ntnl21
Chronic
AWQC
Exeedences
(
Ntl/
OPt
2)
.
ntlb2
Acute
AWQC
Exeedences
(
Ntl/
Bl)
35.1274576
ntnl12
Acute
AWQC
Exeedences
(
Ntl/
Opt
1)
.
ntnl22
Acute
AWQC
Exeedences
(
Ntl/
Opt
2)
.
rchcrb
Number
of
Chronic
reaches
(
BL)
12.0000000
rchcr1
Number
of
Chronic
reaches
(
Opt
1)
6.0000000
rchcr2
Number
of
Chronic
reaches
(
Opt
2)
0
rchacb
Number
of
Acute
reaches
(
BL)
2.0000000
rchac1
Number
of
Acute
reaches
(
Opt
1)
0
rchac2
Number
of
Acute
reaches
(
Opt
2)
0
                                                                
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
10
Output
from
Aquatic
Life
(
Sub
3)
AWQC
Exceedances
by
pollutant/
option
(
Sample
/
National)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Label
Sum
                                                                
crb
Chronic
AWQC
Exeedences
(
Smpl/
BL)
23.0000000
cr1
Chronic
AWQC
Exeedences
(
Smpl/
Opt
1)
6.0000000
cr2
Chronic
AWQC
Exeedences
(
Smpl/
Opt
2)
.
acb
Acute
AWQC
Exeedences
(
Smpl/
BL)
4.0000000
ac1
Acute
AWQC
Exeedences
(
Smpl/
Opt
1)
.
ac2
Acute
AWQC
Exeedences
(
Smpl/
Opt
2)
.
ntlbw
Chronic
AWQC
Exeedences
(
Ntl/
BL)
423.3144485
ntl1w
Chronic
AWQC
Exeedences
(
Ntl/
OPt
1)
52.8929261
ntl2w
Chronic
AWQC
Exeedences
(
Ntl/
OPt
2)
.
ntlbo
Acute
AWQC
Exeedences
(
Ntl/
Bl)
35.1274576
ntl1o
Acute
AWQC
Exeedences
(
Ntl/
Opt
1)
.
ntl2o
Acute
AWQC
Exeedences
(
Ntl/
Opt
2)
.
chemcrb
Number
of
Chronic
pol.
(
BL)
9.0000000
chemcr1
Number
of
Chronic
pol.
(
Opt
1)
3.0000000
chemcr2
Number
of
Chronic
pol.
(
Opt
2)
0
chemacb
Number
of
Acute
pol.
(
BL)
4.0000000
212
chemac1
Number
of
Acute
pol.
(
Opt
1)
0
chemac2
Number
of
Acute
pol.
(
Opt
2)
0
                                                                
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
11
Output
from
Aquatic
Life
(
Sub
3)
IL
Reaches
Exceeding
AWQC
Limits
for
Aquatic
Life
(
National)
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Label
Sum
                                                                       
chrnc_
b
#
of
Reaches
Exceeding
Chronic
AWQC
(
Bl)
353.0595333
chrnc_
1
#
of
Reaches
Exceeding
Chronic
AWQC
(
Opt
1)
52.8929261
chrnc_
2
#
of
Reaches
Exceeding
Chronic
AWQC
(
Opt
2)
0
acutr_
b
#
of
Reaches
Exceeding
Acute
AWQC
(
Bl)
17.5637288
acutr_
1
#
of
Reaches
Exceeding
Acute
AWQC
(
Opt
1)
0
acutr_
2
#
of
Reaches
Exceeding
Acute
AWQC
(
Opt
2)
0
multi_
b
Either
(
BL)
353.0595333
multi_
1
Either
(
Opt
1)
52.8929261
multi_
2
Either
(
Opt
2)
0
                                                                       
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
13
Output
from
Aquatic
Life
(
Sub
3)
AWQC
Eliminations
and
Reductions
(
National)
Number
of
National
Reaches
15:
45
Thursday,
February
6,
2003
The
MEANS
Procedure
Variable
Label
Sum
                                                         
elim1_
1
Elim.
Chronic
(
Opt1)
300.1666072
elim2_
1
Elim.
Chronic
(
Opt2)
353.0595333
elim1_
2
Elim.
Acute
(
Opt1)
17.5637288
elim2_
2
Elim.
Acute
(
Opt2)
17.5637288
red1_
1
Redu.
Chronic
(
Opt1)
8.7818644
red2_
1
Redu.
Chronic
(
Opt2)
0
red1_
2
Redu.
Acute
(
Opt1)
0
red2_
2
Redu.
Acute
(
Opt2)
0
reduc1
Total
Reaches
Reducing
(
Opt1)
8.7818644
reduc2
Total
Reaches
Reducing
(
Opt2)
0
                                                         
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
15
Output
from
Aquatic
Life
(
Sub
3)
Exceedences
by
Chemical
15:
45
Thursday,
February
6,
2003
C
H
c
c
c
c
c
c
E
h
h
h
h
h
h
M
n
n
n
n
n
n
e
e
e
e
e
e
I
t
t
t
t
t
t
m
m
m
m
m
m
O
C
c
c
c
a
a
a
l
l
l
l
l
l
c
c
c
a
a
a
b
A
r
r
r
c
c
c
b
1
2
b
1
2
r
r
r
c
c
c
213
s
L
b
1
2
b
1
2
w
w
w
o
o
o
b
1
2
b
1
2
1
Aluminum
1
.
.
.
.
.
8.782
.
.
.
.
.
1
0
0
0
0
0
2
Cadmium
1
.
.
1
.
.
8.782
.
.
8.78186
.
.
1
0
0
1
0
0
3
Copper
4
1
.
1
.
.
276.064
8.7819
.
8.78186
.
.
1
1
0
1
0
0
4
Lead
2
.
.
.
.
.
17.564
.
.
.
.
.
1
0
0
0
0
0
5
Manganese
1
1
.
.
.
.
21.382
21.3819
.
.
.
.
1
1
0
0
0
0
6
Mercury
1
.
.
.
.
.
8.782
.
.
.
.
.
1
0
0
0
0
0
7
Nickel
1
.
.
.
.
.
8.782
.
.
.
.
.
1
0
0
0
0
0
8
Silver
11
4
.
1
.
.
64.395
22.7291
.
8.78186
.
.
1
1
0
1
0
0
9
Zinc
1
.
.
1
.
.
8.782
.
.
8.78186
.
.
1
0
0
1
0
0
214
"
Fishing
Recreation
(
Sub
7).
sas"
Data
on
human
health
based
and
aquatic
life
based
AWQC
exceedences
for
each
RF1
reach
are
merged,
and
the
total
number
of
AWQC
exceedences
is
generated.
Based
on
the
reduction
in
AWQC
exceedences,
the
baseline
value
of
a
fishing
day,
and
a
range
of
percentage
increases
in
value
from
AWQC
reductions,
benefits
to
recreational
fishermen
are
estimated.
215
/*
Program:
Fishing
Recreation
(
Sub
7).
sas
Authors:
Abt
Associates
Inc;
Rob
Sartain,
Elena
Besedin,
Stuart
Smith
Date:
04/
21/
2000
Updated:
09/
26/
2002
This
program
estimates
recreation
benefits
to
sport
anglers
stemming
from
improvments
in
or
the
elimination
of
awqc
exceedences
for
human
health
or
aquatic
life
on
reaches.
*/

title2
'
Fishing
Recreation
(
Sub
7)';

/*
COUNTY
LEVEL
FISHING
LICENSE
DATA
*/

/*
Load
fishing
license
data
here.
Sum
over
the
types
by
county.
Sum
over
the
counties
containing
a
given
reach.

Final
data
set
contains
reaches
with
the
total
possible
number
of
affected
anglers.
This
number
will
be
multiplied
by
the
ratio
of
reach
length
to
total
reach
miles
in
affected
counties
to
obtain
an
estimate
of
the
number
of
affected
anglers.
*/

data
licens;
set
mpm2.
fish_
lic;

fw_
lic
=
sum(
fw_
res,
fw_
reso,
fw_
sres,
fw_
smres,
fw_
dres,
trout,
fw_
life,
fw_
hunt1,
fw_
2hunt,
fw_
hunts,
fw_
hunt2,
fw_
salt1,
fw_
trout);
sw_
lic
=
sum(
fw_
salt1,
sw_
res,
sw_
dres);

if
fw_
lic=.
then
fw_
lic=
0;
if
sw_
lic=.
then
sw_
lic=
0;
cntyfip
=
fips;

keep
cntyfip
fw_
lic
sw_
lic;
run;

proc
sort
nodupkey
data=
licens;
by
cntyfip;
quit;

/*
DETERMINE
COUNTIES
AFFECTED
BY
REACH
*/
data
facs;
set
mpm2.
fac_
info;
run;

proc
sort
data=
facs;
by
siteid;
quit;

data
cnty;
set
facs;

cntyfip1=
stco1;
cntyfip2=
stco2;
cntyfip3=
stco3;
cntyfip4=
stco4;
cntyfip5=
stco5;
cntyfip6=
stco6;
cntyfip7=
stco7;
cntyfip8=
stco8;
216
keep
reachnmb
cntyfip1­
cntyfip8;
run;

proc
sort
nodupkey
data=
cnty;
by
reachnmb;
quit;

proc
transpose
data=
cnty
out=
cnty2;
by
reachnmb;
var
cntyfip1­
cntyfip8;
quit;

data
cnty3;
set
cnty2;
cntyfip
=
col1/
1;
if
cntyfip=.
then
delete;
keep
reachnmb
cntyfip;
run;

proc
sort
data=
cnty3;
by
cntyfip;
quit;

/*
OBTAIN
REACH
MILES
BY
COUNTY
*/

data
rchmile;
set
mpm2.
us_
rchmi;
cntyfip=
fips;
keep
countymi
cntyfip;
run;

proc
sort
data=
rchmile;
by
cntyfip;

/*
PULL
TOGETHER
OUR
REACHES
WITH
ASSOCIATED
COUNTIES
*/

data
tot_
mile;
merge
cnty3
rchmile;
by
cntyfip;

if
reachnmb=.
then
delete;
run;

proc
sort
data=
tot_
mile;
by
reachnmb
cntyfip;
quit;
proc
sort
nodupkey
data=
facs;
by
reachnmb;
quit;

/*
SUM
REACH
MILES
OVER
COUNTIES
CONTAINING
THE
SAME
REACH
FINISHED
DATA
SET
CONTAINS
TOTAL
NUMBER
OF
REACH
MILES
IN
COUNTIES
CONTAINING
A
GIVEN
REACH.
BRING
IN
THE
FISHING
LICNESE
DATA
BY
COUNTY.
*/

proc
sql;

create
table
miles
as
select
tot_
mile.
reachnmb,
sum(
countymi)
as
denom,
min(
segl)
as
leng,
sum(
fw_
lic)
as
fwlictot,
sum(
sw_
lic)
as
swlictot
from
tot_
mile,
facs,
licens
where
tot_
mile.
reachnmb
=
facs.
reachnmb
and
tot_
mile.
cntyfip
=
licens.
cntyfip
217
group
by
tot_
mile.
reachnmb
order
by
tot_
mile.
reachnmb;
quit;

proc
sort
data=
miles;
by
reachnmb;
quit;

/*
ESTIMATE
SPORT
ANGLERS
*/

data
pops;
set
miles;
by
reachnmb;
if
first.
reachnmb;

/*
no
family
adjustment,
just
anglers
*/
/*
no
subsistence
anglers
*/

if
leng=.
then
fwspt
=
fwlictot*
0.06;
else
fwspt
=
fwlictot*(
leng/
denom);

if
leng=.
then
swspt
=
swlictot*
0.06;
else
swspt
=
swlictot*(
leng/
denom);

fwsptpop
=
round(
fwspt);
swsptpop
=
round(
swspt);

sportpop
=
fwsptpop
+
swsptpop;

keep
reachnmb
sportpop
fwsptpop
swsptpop
fwlictot
swlictot;
run;

proc
sort
data=
pops;
by
reachnmb;
quit;

/*
Combine
Population
with
average
fishing
days
*/

proc
sql;

create
table
pops2
as
select
pops.*,
fip
from
pops,
mpm2.
rchstate
where
pops.
reachnmb
=
rchstate.
reachnmb
order
by
reachnmb;

create
table
pops3
as
select
pops2.*,
avgfw,
avgsw
from
pops2,
mpm2.
avgdfish
where
pops2.
fip
=
avgdfish.
fip
order
by
reachnmb;
quit;

/*
Pull
in
the
reach
data
sets
with
exceedence
events
and
combine
by
reach,
casno,
event
weight
*/
218
proc
sort
data=
rundat.
rec_
hh;
by
reachnmb
casno
event_
wt;
quit;
proc
sort
data=
rundat.
rec_
al;
by
reachnmb
casno
event_
wt;
quit;

data
rchexc;
merge
rundat.
rec_
hh
rundat.
rec_
al;
by
reachnmb
casno
event_
wt;

/*
Calculate
total
exceedences
for
a
reach,
both
sample
and
national
Sum
all
exceedences
for
a
reach,
that's
acute
and
chronic
for
AL
and
h2o/
org
and
org
for
HH
*/

if
orgh2o_
b=.
then
orgh2o_
b=
0;
if
org_
b=.
then
org_
b=
0;
if
chronicb=.
then
chronicb=
0;
if
acuteb=.
then
acuteb=
0;
if
orgh2ot1=.
then
orgh2ot1=
0;
if
org_
t1=.
then
org_
t1=
0;
if
chronic1=.
then
chronic1=
0;
if
acute1=.
then
acute1=
0;
if
orgh2ot2=.
then
orgh2ot2=
0;
if
org_
t2=.
then
org_
t2=
0;
if
chronic2=.
then
chronic2=
0;
if
acute2=.
then
acute2=
0;

base
=
orgh2o_
b
+
org_
b
+
chronicb
+
acuteb;
trt1
=
orgh2ot1
+
org_
t1
+
chronic1
+
acute1;
trt2
=
orgh2ot2
+
org_
t2
+
chronic2
+
acute2;
run;

proc
sql;
create
table
evntsum
as
select
reachnmb,
event_
wt,
sum(
base)
as
exc_
b,
sum(
trt1)
as
exc_
1,
sum(
trt2)
as
exc_
2
from
rchexc
group
by
reachnmb,
event_
wt
order
by
reachnmb,
event_
wt;
quit;

proc
print
data=
evntsum;
title3
'
Summary
of
Exceedences
by
Reach';
quit;

/*
Calculate
benefit
ratios
Ratio
=
#
exceedences
eliminated
/
#
baseline
exceedences
If
all
exceedences
eliminated
then
Ratio
=
1
and
full
benefit
is
realized.
Otherwise
the
fraction
of
the
benefit
realized
is
equal
to
the
fraction
of
exceedences
eliminated.
*/

data
benrat;
set
evntsum;
219
bfrac_
1
=
(
exc_
b
­
exc_
1)/
exc_
b;
bfrac_
2
=
(
exc_
b
­
exc_
2)/
exc_
b;

weight
=
event_
wt;

if
weight=.
then
weight
=
1;

keep
reachnmb
weight
bfrac_
1
bfrac_
2
exc_
b
exc_
1
exc_
2;
run;

/*
Pull
together
the
benefit
fractions
and
fishing
days
and
populations
*/

proc
sql;

create
table
benefit
as
select
benrat.*,
sportpop,
fwsptpop,
swsptpop,
avgfw,
avgsw
from
benrat,
pops3
where
benrat.
reachnmb
=
pops3.
reachnmb
order
by
benrat.
reachnmb;
quit;

data
benefit2;
set
benefit;

/*
Base
values
for
Fishing
Days
Values
in
$
2001
Low
($
99)
28.11
Upper
($
99)
60.43
*/

bvalfwl
=
fwsptpop
*
avgfw
*
28.11;
bvalfwh
=
fwsptpop
*
avgfw
*
60.43;
bvalswl
=
swsptpop
*
avgsw
*
28.11;
bvalswh
=
swsptpop
*
avgsw
*
60.43;
midfw
=
fwsptpop
*
avgfw
*
42.12;
midsw
=
swsptpop
*
avgsw
*
42.12;

tripsfw
=
fwsptpop*
avgfw;
tripssw
=
swsptpop*
avgsw;

/*
Benefits
under
Average
Values
Increase
range
of
9.83%
to
14.71%
Midpoint
12.27%
*/

avg1_
lf
=
bvalfwl
*
bfrac_
1
*
0.0983;
avg1_
hf
=
bvalfwh
*
bfrac_
1
*
0.1471;
avg1_
ls
=
bvalswl
*
bfrac_
1
*
0.0983;
avg1_
hs
=
bvalswh
*
bfrac_
1
*
0.1471;
mid1_
f
=
midfw
*
bfrac_
1
*
0.1227;
mid1_
s
=
midsw
*
bfrac_
1
*
0.1227;
220
avg2_
lf
=
bvalfwl
*
bfrac_
2
*
0.0983;
avg2_
hf
=
bvalfwh
*
bfrac_
2
*
0.1471;
avg2_
ls
=
bvalswl
*
bfrac_
2
*
0.0983;
avg2_
hs
=
bvalswh
*
bfrac_
2
*
0.1471;
mid2_
f
=
midfw
*
bfrac_
2
*
0.1227;
mid2_
s
=
midsw
*
bfrac_
2
*
0.1227;
run;

proc
means
sum;
var
avg1_
lf
avg1_
hf
avg1_
ls
avg1_
hs
avg2_
lf
avg2_
hf
avg2_
ls
avg2_
hs
mid1_
f
mid1_
s
mid2_
f
mid2_
s
fwsptpop
swsptpop
tripsfw
tripssw
bvalfwl
bvalfwh
bvalswl
bvalswh;
title3
'
Sample
Benefits
(
includes
partial)';
quit;

proc
means
sum;
weight
weight;
var
avg1_
lf
avg1_
hf
avg1_
ls
avg1_
hs
avg2_
lf
avg2_
hf
avg2_
ls
avg2_
hs
mid1_
f
mid1_
s
mid2_
f
mid2_
s
fwsptpop
swsptpop
tripsfw
tripssw
bvalfwl
bvalfwh
bvalswl
bvalswh;
output
out
=
rundat.
rec_
fsh(
drop
=
_
type_
_
freq_)
sum
=
avg1_
lf
avg1_
hf
avg1_
ls
avg1_
hs
avg2_
lf
avg2_
hf
avg2_
ls
avg2_
hs
mid1_
f
mid1_
s
mid2_
f
mid2_
s
fwsptpop
swsptpop
tripsfw
tripssw
bvalfwl
bvalfwh
bvalswl
bvalswh;
title3
'
National
Benefits
(
includes
partial)';
quit;

data
rundat.
rec_
fish(
keep
=
activity
benl1
benm1
benh1
benl2
benm2
benh2
pop);
set
rundat.
rec_
fsh;
activity
=
"
Fishing";
pop
=
fwsptpop
+
swsptpop;
benl1
=
avg1_
lf
+
avg1_
ls;
benm1
=
mid1_
f
+
mid1_
s;
benh1
=
avg1_
hf
+
avg1_
hs;
benl2
=
avg2_
lf
+
avg2_
ls;
benm2
=
mid2_
f
+
mid2_
s;
benh2
=
avg2_
hf
+
avg2_
hs;
run;

data
reaches;
set
benrat;

if
exc_
b
>
0
then
baseexc
=
1;
else
baseexc
=
0;
if
exc_
1
>
0
then
post1
=
1;
else
post1=
0;
if
exc_
2
>
0
then
post2
=
1;
else
post2=
0;
if
bfrac_
1
=
1
then
elm_
1
=
1;
else
elm_
1
=
0;
if
bfrac_
1
<
1
and
bfrac_
1
>
0
then
red_
1
=
1;
else
red_
1
=
0;
if
bfrac_
2
=
1
then
elm_
2
=
1;
else
elm_
2
=
0;
if
bfrac_
2
<
1
and
bfrac_
2
>
0
then
red_
2
=
1;
else
red_
2
=
0;
run;

proc
means
sum
data=
reaches;
var
baseexc
post1
post2
elm_
1
elm_
2
red_
1
red_
2;
label
221
baseexc
=
'
Baseline
Reaches
Exceeding
AWQCs'
post1
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt1)'
post2
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt2)'
elm_
1
=
'
Reaches
Eliminating
(
Opt1)'
elm_
2
=
'
Reaches
Eliminating
(
Opt2)'
red_
1
=
'
Reaches
Reducing
(
Opt1)'
red_
2
=
'
Reaches
Reducing
(
Opt2)';
title3
'
Reach
Information';
title4
'
Sample
Level';
quit;

proc
means
sum
data=
reaches;
weight
weight;
var
baseexc
post1
post2
elm_
1
elm_
2
red_
1
red_
2;
label
baseexc
=
'
Baseline
Reaches
Exceeding
AWQCs'
post1
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt1)'
post2
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt2)'
elm_
1
=
'
Reaches
Eliminating
(
Opt1)'
elm_
2
=
'
Reaches
Eliminating
(
Opt2)'
red_
1
=
'
Reaches
Reducing
(
Opt1)'
red_
2
=
'
Reaches
Reducing
(
Opt2)';
title3
'
Reach
Information';
title4
'
National
Level';
quit;

data
check;
set
benefit2;

if
bfrac_
1
<=
0
then
delete;
run;

data
check2;
merge
check
pops;
by
reachnmb;
if
weight=.
then
delete;

keep
reachnmb
fwlictot
swlictot
fwsptpop
swsptpop
weight;
run;

/*
End
Fishing
Recreation
(
Sub
7).
sas
*/
222
Program
Output
223
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
51
Fishing
Recreation
(
Sub
7)
15:
45
Thursday,
February
6,
2003
National
Benefits
(
includes
partial)

The
MEANS
Procedure
Variable
Sum
                        
avg1_
lf
28659321.49
avg1_
hf
92197003.81
avg1_
ls
53551.53
avg1_
hs
172275.21
avg2_
lf
46788745.11
avg2_
hf
150519338
avg2_
ls
128706.69
avg2_
hs
414049.27
mid1_
f
53602428.48
mid1_
s
100159.11
mid2_
f
87510458.49
mid2_
s
240724.16
fwsptpop
978291.79
swsptpop
4140.69
tripsfw
16932731.97
tripssw
46578.64
bvalfwl
475979096
bvalfwh
1023244993
bvalswl
1309325.43
bvalswh
2814746.92
                        
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
53
Fishing
Recreation
(
Sub
7)
15:
45
Thursday,
February
6,
2003
Reach
Information
National
Level
The
MEANS
Procedure
Variable
Label
Sum
                                                               
baseexc
Baseline
Reaches
Exceeding
AWQCs
394.4913361
post1
Post
Reaches
Exceeding
AWQCs
(
Opt1)
109.2789846
post2
Post
Reaches
Exceeding
AWQCs
(
Opt2)
0
elm_
1
Reaches
Eliminating
(
Opt1)
285.2123515
elm_
2
Reaches
Eliminating
(
Opt2)
394.4913361
red_
1
Reaches
Reducing
(
Opt1)
23.7361201
red_
2
Reaches
Reducing
(
Opt2)
0
                                                               
224
Near­
water
Activity
and
Boating
225
"
Human
Health
AWQC
Exceedences
(
Sub
2).
sas"
The
number
of
human
health
based
AWQC
exceedences
for
each
RF1
reach
are
estimated
and
output
to
a
data
set.
See
Human
Health
Benefits
Cancer
from
Fish
Consumption
for
program
code
and
output.
226
"
Aquatic
Life
AWQC
Exceedences
(
Sub
3).
sas"
The
number
of
aquatic
life
based
AWQC
exceedences
for
each
RF1
reach
are
estimated
and
output
to
a
data
set.
See
Recreational
Fishing
for
program
code
and
output.
227
"
Other
Recreation
(
Sub
8).
sas"
Data
on
human
health
based
and
aquatic
life
based
AWQC
exceedences
for
each
RF1
reach
are
merged,
and
the
total
number
of
AWQC
exceedences
is
generated.
Based
on
the
reduction
in
AWQC
exceedences,
the
baseline
value
of
a
boating
or
near­
water
activity
day,
and
a
range
of
percentage
increases
in
value
from
AWQC
reductions,
benefits
to
boaters
and
near­
water
activity
participants
are
estimated.
228
/*
Program:
Other
Recreation
(
Sub
8).
sas
Authors:
Abt
Associates
Inc;
Rob
Sartain,
Elena
Besedin,
Stuart
Smith
Date:
04/
21/
2000
Updated:
09/
26/
2002
This
program
estimates
recreation
benefits
to
activity
participants
stemming
from
improvements
in
or
the
elimination
of
awqc
exceedences
for
human
health
or
aquatic
life
on
reaches.
Valuations
for
viewing
and
boating
are
estimated.
*/

%
let
tbl_
path
=
%
sysfunc(
pathname(
rundat));

title2
'
Other
Recreation
(
Sub
8)';

/*
DETERMINE
COUNTIES
AFFECTED
BY
REACH
*/
data
facs;
set
mpm2.
fac_
info;
run;

proc
sort
data=
facs;
by
siteid;
quit;

data
cnty;
set
facs;

cntyfip1=
stco1;
cntyfip2=
stco2;
cntyfip3=
stco3;
cntyfip4=
stco4;
cntyfip5=
stco5;
cntyfip6=
stco6;
cntyfip7=
stco7;
cntyfip8=
stco8;

keep
reachnmb
cntyfip1­
cntyfip8;
run;

proc
sort
nodupkey
data=
cnty;
by
reachnmb;
quit;

proc
transpose
data=
cnty
out=
cnty2;
by
reachnmb;
var
cntyfip1­
cntyfip8;
quit;

data
cnty3;
set
cnty2;
cntyfip
=
col1/
1;
if
cntyfip=.
then
delete;
keep
reachnmb
cntyfip;
run;

proc
sort
data=
cnty3;
by
cntyfip;
quit;

/*
OBTAIN
REACH
MILES
BY
COUNTY
*/

data
rchmile;
set
mpm2.
us_
rchmi;
cntyfip=
fips;
keep
countymi
cntyfip;
229
run;

proc
sort
data=
rchmile;
by
cntyfip;

/*
PULL
TOGETHER
OUR
REACHES
WITH
ASSOCIATED
COUNTIES
*/

data
tot_
mile;
merge
cnty3
rchmile;
by
cntyfip;

if
reachnmb=.
then
delete;
run;

proc
sort
data=
tot_
mile;
by
reachnmb
cntyfip;
quit;
proc
sort
nodupkey
data=
facs;
by
reachnmb;
quit;

/*
SUM
REACH
MILES
OVER
COUNTIES
CONTAINING
THE
SAME
REACH
FINISHED
DATA
SET
CONTAINS
TOTAL
NUMBER
OF
REACH
MILES
IN
COUNTIES
CONTAINING
A
GIVEN
REACH.
BRING
IN
THE
FISHING
LICNESE
DATA
BY
COUNTY.
*/

proc
sql;
create
table
miles
as
select
tot_
mile.
reachnmb,
sum(
countymi)
as
denom,
min(
segl)
as
leng
from
tot_
mile,
facs
where
tot_
mile.
reachnmb
=
facs.
reachnmb
group
by
tot_
mile.
reachnmb
order
by
tot_
mile.
reachnmb;
quit;

/*
Merge
in
the
county
populations
and
state
abbreviations
*/
proc
sql;
create
table
rchpop
as
select
cnty3.
reachnmb,
cnty3.
cntyfip,
pop2000
as
totpop,
stateabv
as
state,
denom,
leng
from
cnty3,
mpm2.
cntypop,
mpm2.
rchstate,
miles
where
cnty3.
cntyfip
=
cntypop.
fips
and
cnty3.
reachnmb
=
rchstate.
reachnmb
and
cnty3.
reachnmb
=
miles.
reachnmb
order
by
reachnmb;

/*
ESTIMATE
Number
of
Participants
by
Activity
*/

create
table
pops
as
select
reachnmb,
sum(
totpop)
as
pop,
min(
state)
as
stateabv,
avg(
denom)
as
den,
avg(
leng)
as
lengt
from
rchpop
group
by
reachnmb
order
by
reachnmb;

create
table
pops2
as
select
reachnmb,
stateabv,
viewtrp,
boattrp,
swimtrp,
pop,
pop*
pctview
as
totview,
pop*
pctboat
as
totboat,
pop*
pctswim
as
totswim,
case
230
when
lengt
is
missing
then
pop*
pctview*
0.06
else
pop*
pctview*(
lengt/
den)
end
as
viewers,
case
when
lengt
is
missing
then
pop*
pctboat*
0.06
else
pop*
pctboat*(
lengt/
den)
end
as
boaters,
case
when
lengt
is
missing
then
pop*
pctswim*
0.06
else
pop*
pctswim*(
lengt/
den)
end
as
swimers
from
pops,
mpm2.
recpcts
where
pops.
stateabv
=
recpcts.
state
order
by
reachnmb;
quit;

/*
Pull
in
the
reach
data
sets
with
exceedence
events
and
combine
by
reach,
casno,
event
weight
*/

proc
sort
data=
rundat.
rec_
hh;
by
reachnmb
casno
event_
wt;
quit;
proc
sort
data=
rundat.
rec_
al;
by
reachnmb
casno
event_
wt;
quit;

data
rchexc;
merge
rundat.
rec_
hh
rundat.
rec_
al;
by
reachnmb
casno
event_
wt;

/*
Calculate
total
exceedences
for
a
reach,
both
sample
and
national
Sum
all
exceedences
for
a
reach,
that's
acute
and
chronic
for
AL
and
h2o/
org
and
org
for
HH
*/

if
orgh2o_
b=.
then
orgh2o_
b=
0;
if
org_
b=.
then
org_
b=
0;
if
chronicb=.
then
chronicb=
0;
if
acuteb=.
then
acuteb=
0;
if
orgh2ot1=.
then
orgh2ot1=
0;
if
org_
t1=.
then
org_
t1=
0;
if
chronic1=.
then
chronic1=
0;
if
acute1=.
then
acute1=
0;
if
orgh2ot2=.
then
orgh2ot2=
0;
if
org_
t2=.
then
org_
t2=
0;
if
chronic2=.
then
chronic2=
0;
if
acute2=.
then
acute2=
0;

base
=
orgh2o_
b
+
org_
b
+
chronicb
+
acuteb;
trt1
=
orgh2ot1
+
org_
t1
+
chronic1
+
acute1;
trt2
=
orgh2ot2
+
org_
t2
+
chronic2
+
acute2;
run;

proc
sql;
create
table
evntsum
as
select
reachnmb,
event_
wt,
sum(
base)
as
exc_
b,
sum(
trt1)
as
exc_
1,
sum(
trt2)
as
exc_
2
from
rchexc
231
group
by
reachnmb,
event_
wt
order
by
reachnmb,
event_
wt;
quit;

proc
print
data=
evntsum;
quit;

/*
Calculate
benefit
ratios
Ratio
=
#
exceedences
eliminated
/
#
baseline
exceedences
If
all
exceedences
eliminated
then
Ratio
=
1
and
full
benefit
is
realized.
Otherwise
the
fraction
of
the
benefit
realized
is
equal
to
the
fraction
of
exceedences
eliminated.
*/

data
benrat;
set
evntsum;

bfrac_
1
=
(
exc_
b
­
exc_
1)/
exc_
b;
bfrac_
2
=
(
exc_
b
­
exc_
2)/
exc_
b;

weight
=
event_
wt;

if
weight=.
then
weight
=
1;

keep
reachnmb
weight
bfrac_
1
bfrac_
2
exc_
b
exc_
1
exc_
2;
run;

/*
Pull
together
the
benefit
fractions
and
fishing
days
and
populations
*/

proc
sql;

create
table
benefit
as
select
benrat.*,
viewers,
boaters,
viewtrp,
boattrp,
swimers,
swimtrp
from
benrat,
pops2
where
benrat.
reachnmb
=
pops2.
reachnmb
order
by
benrat.
reachnmb;
quit;

data
benefit2;
set
benefit;

/*
Base
values
for
Viewing
and
Boating
and
Swimming
Days
Viewing
Day
($
2001)
22.73
­
28.73
Boating
Day
($
2001)
37.30
­
59.26
midpoint
view
26.28
midpoint
boat
48.30
*/

bvlow
=
viewers
*
22.73
*
viewtrp;
bvhigh
=
viewers
*
28.73
*
viewtrp;
bblow
=
boaters
*
37.30
*
boattrp;
bbhigh
=
boaters
*
59.26
*
boattrp;
232
midview
=
viewers
*
26.28
*
viewtrp;
midboat
=
boaters
*
48.30
*
boattrp;

tripsv
=
viewers*
viewtrp;
tripsb
=
boaters*
boattrp;

/*
Benefits
from
Average
values
Viewing
Increase
of
14.94%
to
21.34%
Boating
Increase
of
6.51%
to
11.35%

mid
viewing
18.14%
mid
boating
8.93%
*/

av1_
low
=
bvlow
*
bfrac_
1
*
0.1494;
av1_
high
=
bvhigh
*
bfrac_
1
*
0.2134;
av2_
low
=
bvlow
*
bfrac_
2
*
0.1494;
av2_
high
=
bvhigh
*
bfrac_
2
*
0.2134;

ab1_
low
=
bblow
*
bfrac_
1
*
0.0651;
ab1_
high
=
bbhigh
*
bfrac_
1
*
0.1135;
ab2_
low
=
bblow
*
bfrac_
2
*
0.0651;
ab2_
high
=
bbhigh
*
bfrac_
2
*
0.1135;

midvw1
=
midview
*
bfrac_
1
*
0.1814;
midbt1
=
midboat
*
bfrac_
1
*
0.0893;
midvw2
=
midview
*
bfrac_
2
*
0.1814;
midbt2
=
midboat
*
bfrac_
2
*
0.0893;
run;

proc
means
sum;
var
av1_
low
av1_
high
av2_
low
av2_
high
ab1_
low
ab1_
high
ab2_
low
ab2_
high
midvw1
midvw2
midbt1
midbt2
viewers
boaters
tripsv
tripsb;
label
av1_
low
=
'
Avg
View
Benefits
­
Low
(
Opt
1)'
av1_
high
=
'
Avg
View
Benefits
­
High(
Opt
1)'
av2_
low
=
'
Avg
View
Benefits
­
Low
(
Opt
2)'
av2_
high
=
'
Avg
View
Benefits
­
High(
Opt
2)'
ab1_
low
=
'
Avg
Boat
Benefits
­
Low
(
Opt
1)'
ab1_
high
=
'
Avg
Boat
Benefits
­
High(
Opt
1)'
ab2_
low
=
'
Avg
Boat
Benefits
­
Low
(
Opt
2)'
ab2_
high
=
'
Avg
Boat
Benefits
­
High(
Opt
2)'
midvw1
=
'
Mid
View
Ben
(
Opt1)'
midvw2
=
'
Mid
View
Ben
(
Opt2)'
midbt1
=
'
Mid
Boat
Ben
(
Opt1)'
midbt2
=
'
Mid
Boat
Ben
(
Opt2)'
viewers
=
'
Number
of
Viewing
Participants'
boaters
=
'
Number
of
Boating
Participants'
tripsv
=
'
Number
of
Viewing
Days'
tripsb
=
'
Number
of
Boating
Days'
title3
'
Sample
Benefits
(
includes
partial)';
title4
'
Viewing,
Boating
Activities
';
233
quit;

proc
means
sum;
weight
weight;
var
av1_
low
av1_
high
av2_
low
av2_
high
ab1_
low
ab1_
high
ab2_
low
ab2_
high
midvw1
midvw2
midbt1
midbt2
viewers
boaters
tripsv
tripsb;
output
out
=
rundat.
rec_
oth(
drop
=
_
type_
_
freq_)
sum
=
av1_
low
av1_
high
av2_
low
av2_
high
ab1_
low
ab1_
high
ab2_
low
ab2_
high
midvw1
midvw2
midbt1
midbt2
viewers
boaters
tripsv
tripsb;
label
av1_
low
=
'
Avg
View
Benefits
­
Low
(
Opt
1)'
av1_
high
=
'
Avg
View
Benefits
­
High(
Opt
1)'
av2_
low
=
'
Avg
View
Benefits
­
Low
(
Opt
2)'
av2_
high
=
'
Avg
View
Benefits
­
High(
Opt
2)'
ab1_
low
=
'
Avg
Boat
Benefits
­
Low
(
Opt
1)'
ab1_
high
=
'
Avg
Boat
Benefits
­
High(
Opt
1)'
ab2_
low
=
'
Avg
Boat
Benefits
­
Low
(
Opt
2)'
ab2_
high
=
'
Avg
Boat
Benefits
­
High(
Opt
2)'
midvw1
=
'
Mid
View
Ben
(
Opt1)'
midvw2
=
'
Mid
View
Ben
(
Opt2)'
midbt1
=
'
Mid
Boat
Ben
(
Opt1)'
midbt2
=
'
Mid
Boat
Ben
(
Opt2)'
viewers
=
'
Number
of
Viewing
Participants'
boaters
=
'
Number
of
Boating
Participants'
tripsv
=
'
Number
of
Viewing
Days'
tripsb
=
'
Number
of
Boating
Days'
title3
'
National
Benefits
(
includes
partial)';
title4
'
Viewing,
Boating
Activities
';
quit;

data
rundat.
rec_
othr(
keep
=
activity
benl1
benm1
benh1
benl2
benm2
benh2
pop);
set
rundat.
rec_
oth;

if
_
n_
=
1
then
do;
activity
=
"
Boating";
benl1
=
ab1_
low;
benm1
=
midbt1;
benh1
=
ab1_
high;
benl2
=
ab2_
low;
benm2
=
midbt2;
benh2
=
ab2_
high;
pop
=
boaters;
output;
activity
=
"
Viewing";
benl1
=
av1_
low;
benm1
=
midvw1;
benh1
=
av1_
high;
benl2
=
av2_
low;
benm2
=
midvw2;
benh2
=
av2_
high;
pop
=
viewers;
output;
234
end;
run;

data
allrec;
set
rundat.
rec_
fish
rundat.
rec_
othr;
run;

proc
means
sum
data=
allrec;
var
benl1
benm1
benh1
benl2
benm2
benh2;
output
out
=
rectot(
drop
=
_
type_
_
freq_)
sum
=
benl1
benm1
benh1
benl2
benm2
benh2;
quit;

data
totls(
keep
=
activity
benl1
benm1
benh1
benl2
benm2
benh2);
set
rectot;
if
_
n_
=
1
then
do;
activity
=
"
Recreational
Total";
output;
activity
=
"
Non­
Use
Benefits";
benl1
=
benl1
*
0.50;
benm1
=
benm1
*
0.50;
benh1
=
benh1
*
0.50;
benl2
=
benl2
*
0.50;
benm2
=
benm2
*
0.50;
benh2
=
benh2
*
0.50;
output;
activity
=
"
Total
Rec
and
Nonuse";
benl1
=
benl1
+
(
benl1
*
0.50);
benm1
=
benm1
+
(
benm1
*
0.50);
benh1
=
benh1
+
(
benh1
*
0.50);
benl2
=
benl2
+
(
benl2
*
0.50);
benm2
=
benm2
+
(
benm2
*
0.50);
benh2
=
benh2
+
(
benh2
*
0.50);
end;
run;

data
exlrec
(
keep=
titlecol
col1­
col8
label="
Recreational
and
Nonuse
Benefits
­
Weighted");
set
allrec
totls;
length
titlecol
$
50
col1
­
col8
$
128;

if
_
n_
=
1
then
do;
titlecol
=
"
Table:
Recreational
and
Non­
Use
Benefits,
Weighted
Estimates";
output;
titlecol
=
"
Run
on
&
sysdate";
output;
titlecol
=
"
";
output;
col2
=
"
Option
1";
col5
=
"
Option
2";
output;
col1
=
"
Activity";
col2
=
"
Low";
col3
=
"
Mid";
col4
=
"
High";
col5
=
"
Low";
col6
=
"
Mid";
col7
=
"
High";
235
col8
=
"
Population";
output;
end;

col1
=
trim(
activity);
col2
=
trim(
benl1);
col3
=
trim(
benm1);
col4
=
trim(
benh1);
col5
=
trim(
benl2);
col6
=
trim(
benm2);
col7
=
trim(
benh2);
col8
=
pop;
output;
run;

proc
export
data
=
exlrec
outfile
=
"
C:\
Output
Files\
Recreational
and
Nonuse
Benefits.
xls"
dbms=
excel2000
replace;
run;

data
reaches;
set
benrat;

if
exc_
b
>
0
then
baseexc
=
1;
else
baseexc
=
0;
if
exc_
1
>
0
then
post1
=
1;
else
post1=
0;
if
exc_
2
>
0
then
post2
=
1;
else
post2=
0;
if
bfrac_
1
=
1
then
elm_
1
=
1;
else
elm_
1
=
0;
if
bfrac_
1
<
1
and
bfrac_
1
>
0
then
red_
1
=
1;
else
red_
1
=
0;
if
bfrac_
2
=
1
then
elm_
2
=
1;
else
elm_
2
=
0;
if
bfrac_
2
<
1
and
bfrac_
2
>
0
then
red_
2
=
1;
else
red_
2
=
0;
run;

proc
means
sum
data=
reaches;
var
baseexc
post1
post2
elm_
1
elm_
2
red_
1
red_
2;
label
baseexc
=
'
Baseline
Reaches
Exceeding
AWQCs'
post1
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt1)'
post2
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt2)'
elm_
1
=
'
Reaches
Eliminating
(
Opt1)'
elm_
2
=
'
Reaches
Eliminating
(
Opt2)'
red_
1
=
'
Reaches
Reducing
(
Opt1)'
red_
2
=
'
Reaches
Reducing
(
Opt2)';
title3
'
Reach
Information';
title4
'
Sample
Level';
quit;

proc
means
sum
data=
reaches;
weight
weight;
var
baseexc
post1
post2
elm_
1
elm_
2
red_
1
red_
2;
label
baseexc
=
'
Baseline
Reaches
Exceeding
AWQCs'
post1
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt1)'
post2
=
'
Post
Reaches
Exceeding
AWQCs
(
Opt2)'
elm_
1
=
'
Reaches
Eliminating
(
Opt1)'
elm_
2
=
'
Reaches
Eliminating
(
Opt2)'
236
red_
1
=
'
Reaches
Reducing
(
Opt1)'
red_
2
=
'
Reaches
Reducing
(
Opt2)';
title3
'
Reach
Information';
title4
'
National
Level';
quit;

/*
Determine
the
number
of
counties
affected
by
the
benefiting
reaches
*/

data
goodrch1;
set
benefit2;

if
bfrac_
1
<=
0
then
delete;
keep
reachnmb
weight;
run;

data
goodrch2;
set
benefit2;

if
bfrac_
2
<=
0
then
delete;

keep
reachnmb
weight;
run;

proc
sql;

create
table
cntyb1
as
select
cnty3.
reachnmb,
cntyfip,
weight
from
cnty3,
goodrch1
where
cnty3.
reachnmb
=
goodrch1.
reachnmb;

create
table
cntyb2
as
select
cnty3.
reachnmb,
cntyfip,
weight
from
cnty3,
goodrch2
where
cnty3.
reachnmb
=
goodrch2.
reachnmb;

create
table
cntyb1a
as
select
count(
distinct
cntyfip)
as
numcnty1
from
cntyb1;

create
table
cntyb2a
as
select
count(
distinct
cntyfip)
as
numcnty2
from
cntyb2;

create
table
cntyb1b
as
select
distinct
cntyfip,
avg(
weight)
as
wght
from
cntyb1
group
by
cntyfip;

create
table
cntyb1c
as
select
reachnmb,
count(
cntyfip)
as
cnties
from
cntyb1
group
by
reachnmb;
237
quit;

proc
print
data=
cntyb1a;
title3
'
Number
of
affected
counties
Option
1';
title4
'
Sample';
quit;

proc
print
data=
cntyb2a;
title3
'
Number
of
affected
counties
Option
2';
title4
'
Sample';
quit;

proc
means
sum
data=
cntyb1b;
var
wght;
title3
'
Sum
of
Counties
times
Average
weight';
quit;

proc
means
data=
cntyb1b;
var
wght;
title3
'
Mean
County
weight';
quit;

proc
means
data=
cntyb1c;
var
cnties;
title3
'
Average
counties
traversed
by
a
reach';
quit;

data
check;
set
benefit2;

if
bfrac_
1
<=
0
then
delete;
run;

data
check2;
merge
check
pops2;
by
reachnmb;
if
weight=.
then
delete;

keep
reachnmb
pop
totview
totboat
viewers
swimers
weight;
run;

/*
End
Other
Recreation
(
Sub
8).
sas
*/
238
Program
Output
239
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
56
Other
Recreation
(
Sub
8)
15:
45
Thursday,
February
6,
2003
National
Benefits
(
includes
partial)
Viewing,
Boating
Activities
The
MEANS
Procedure
Variable
Label
Sum
                                                           
av1_
low
Avg
View
Benefits
­
Low
(
Opt
1)
56583713.44
av1_
high
Avg
View
Benefits
­
High(
Opt
1)
102157782
av2_
low
Avg
View
Benefits
­
Low
(
Opt
2)
79885242.54
av2_
high
Avg
View
Benefits
­
High(
Opt
2)
144226999
ab1_
low
Avg
Boat
Benefits
­
Low
(
Opt
1)
36511488.63
ab1_
high
Avg
Boat
Benefits
­
High(
Opt
1)
101134010
ab2_
low
Avg
Boat
Benefits
­
Low
(
Opt
2)
51175193.84
ab2_
high
Avg
Boat
Benefits
­
High(
Opt
2)
141751344
midvw1
Mid
View
Ben
(
Opt1)
79433565.34
midvw2
Mid
View
Ben
(
Opt2)
112144807
midbt1
Mid
Boat
Ben
(
Opt1)
64854230.29
midbt2
Mid
Boat
Ben
(
Opt2)
90900917.26
viewers
Number
of
Viewing
Participants
3118919.91
boaters
Number
of
Boating
Participants
2528083.11
tripsv
Number
of
Viewing
Days
23524290.01
tripsb
Number
of
Boating
Days
21075101.55
                                                           
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
59
Other
Recreation
(
Sub
8)
15:
45
Thursday,
February
6,
2003
Reach
Information
National
Level
The
MEANS
Procedure
Variable
Label
Sum
                                                               
baseexc
Baseline
Reaches
Exceeding
AWQCs
394.4913361
post1
Post
Reaches
Exceeding
AWQCs
(
Opt1)
109.2789846
post2
Post
Reaches
Exceeding
AWQCs
(
Opt2)
0
elm_
1
Reaches
Eliminating
(
Opt1)
285.2123515
elm_
2
Reaches
Eliminating
(
Opt2)
394.4913361
red_
1
Reaches
Reducing
(
Opt1)
23.7361201
red_
2
Reaches
Reducing
(
Opt2)
0
                                                               
240
POTW
Benefits
241
POTW
Inhibitions
from
MP&
M
Facilities
242
"
Generate
Discharge
Events
(
Sub1).
sas"
Generates
a
data
set
of
pollutant
loading
events
and
their
corresponding
weights
for
each
RF1
reach
receiving
MP&
M
discharges.
A
special
set
of
event
loadings
discharged
to
POTWs
by
MP&
M
indirect
dischargers
is
created
without
employing
the
POTW
removal
rates.
See
Estimating
MP&
M
Discharge
Loadings
for
program
code.
243
"
POTW
Inhibitions
(
Sub
9).
sas"
Using
the
set
of
loadings
discharged
to
POTWs
by
MP&
M
facilities,
POTW
specific
pollutant
discharge
events
are
calculated.
These
pollutants
concentrations
are
then
compared
with
POTW
inhibition
values
and
exceedences
are
determined.
244
/*
Program:
POTW
Inhibitions
(
Sub
9).
sas
Author:
Abt
Associates
Inc;
Stuart
Smith,
based
on
INH_
FNL.
sas
by
Elena
Besedin,
Rob
Sartain
Date:
04/
21/
2000
Updated:
09/
02/
2002
1.
Calculates
MP&
M
pollutant
concentrations
in
influent
flow
to
the
POTWs
associated
with
the
MP&
M
sample
facilities;
2.
Compares
the
estimated
pollutant
concentrations
with
the
inhibition
values;
3.
Identifies
POTWs
in
which
pollutant
concentrations
exceed
inhibition
values;
4.
Calculates
the
total
number
of
exceedence
events
on
a
sample
and
national
basis;
*/

title2
'
POTW
Inhibitions
(
Sub
9)';

/*
Get
indirect
discharge
loads
*/

proc
sql;
create
table
potw
as
select
potw_
in.
npdes,
casno,
base,
treat1,
treat2,
poflow,
potw_
in.
weights,
potw_
in.
siteid,
size
from
rundat.
potw_
in,
mpm2.
facsize
where
potw_
in.
siteid
=
facsize.
siteid
order
by
npdes,
casno;
quit;

/*
Get
potw
flow
for
each
unique
npdes
*/
data
evntflo1;
set
potw;
run;
proc
sort
data=
evntflo1;
by
npdes
siteid;
quit;

data
evntflo2;
set
evntflo1;
by
npdes
siteid;

if
not
first.
npdes
then
delete;

flow
=
poflow;
keep
npdes
flow;
run;

/*
Pull
in
the
npdes
flow
*/

proc
sql;
create
table
potw1
as
select
potw.
npdes,
casno,
weights,
flow
as
potflow,
base,
treat1,
treat2,
siteid,
size,
0
as
old_
wght
from
potw,
evntflo2
where
potw.
npdes
=
evntflo2.
npdes
and
potw.
npdes
^=
'
'
and
casno
^=
'
'
order
by
npdes,
casno,
size,
weights;
quit;
245
/*
Keep
a
copy
of
initial
set
to
then
drop
observations
from
as
the
event
building
loops
*/

data
starter;
set
potw1;
run;

/*
MACRO
­
This
macro
creates
a
set
of
events
to
each
potw
for
each
chemical
*/
/*
based
on
the
facility
weight.
Loadings
are
summed
by
potw
by
chemical
and
*/
/*
by
mp&
m
facility
size.
Then
(
if
there
was
a
previous
build)
the
weight
*/
/*
is
adjusted
by
the
previous
weight
and
an
output
set
is
appended
to
the
*/
/*
previous
build
output
set.
Finally,
the
facilities
with
the
current
weight
*/
/*
are
removed
from
the
starter
data
set
and
the
build
loop
will
be
run
again
*/

%
macro
buildloop(
first=);
proc
sql;
/*
Sum
by
npdes,
casno,
and
mp&
m
facility
size.
Count
the
number
of
mp&
m
*/
/*
facilities
by
size
included
in
the
summed
loadings.
*/
create
table
summed
as
select
npdes,
casno,
size,
sum(
base)
as
base,
sum(
treat1)
as
treat1,
sum(
treat2)
as
treat2,
min(
potflow)
as
potflow,
count(
1)
as
facs
from
starter
group
by
npdes,
casno,
size
order
by
npdes,
casno,
size;

/*
Grab
the
smallest
weight
of
any
of
the
facilities
discharging
to
the
*/
/*
potw.
This
is
the
current
facility
weight.
Also
keep
the
weight
used
on
*/
/*
the
previous
build
loop.
*/
create
table
wght
as
select
npdes,
min(
weights)
as
weights,
min(
old_
wght)
as
old_
wght
from
starter
group
by
npdes;

/*
Grab
the
loadings
for
this
event
and
calculate
the
weight
of
the
event
*/
/*
Subtract
the
previous
build
weight
to
account
for
only
the
remaining
*/
/*
portion
of
the
differning
weights.
*/

create
table
outit
as
select
summed.*,(
weights
­
old_
wght)
as
weights
from
summed
left
join
wght
on
summed.
npdes
=
wght.
npdes;

/*
Rebuild
the
starter
set
by
dropping
the
facilities
with
the
current
*/
/*
weight.
*/
create
table
starter
as
select
starter.
npdes,
casno,
starter.
weights,
potflow,
base,
treat1,
treat2,
siteid,
size,
wght.
weights
as
old_
wght
from
starter
left
join
wght
on
starter.
npdes
=
wght.
npdes
where
starter.
weights
^=
wght.
weights;
quit;

/*
If
this
is
the
initial
build
loop,
create
the
output
data
set
*/
%
IF
&
first
=
1
%
THEN
%
DO;
data
outset;
set
outit;
246
run;
%
END;

/*
If
this
is
a
subsequent
build
loop,
append
to
the
output
data
set
*/
%
IF
&
first
=
0
%
THEN
%
DO;
data
outset;
set
outset
outit;
run;
%
END;
%
mend
buildloop;

/*
Determine
the
number
of
loops
required
to
build
all
the
needed
events.
*/
/*
The
number
of
loops
is
equal
to
the
maximum
number
of
different
weights
*/
/*
for
facilities
discharging
to
a
single
potw.
If
there
are
5
facilities
*/
/*
with
different
weights
discharging
to
a
potw
and
this
is
the
maximum
*/
/*
number
of
different
weights
to
any
potw,
then
the
loop
must
be
run
5
times
*/
/*
to
account
for
all
the
differential
weights.
Potws
with
fewer
different
*/
/*
weights
will
drop
out
when
their
maximum
has
been
reached
(
i.
e.
they
will
*/
/*
no
longer
appear
in
the
starter
data
set.
*/
proc
sql
noprint;
/*
Get
the
number
of
distinct
weights
across
facilities
for
a
given
potw
*/
create
table
maxw
as
select
distinct
npdes,
count(
distinct
weights)
as
numwghts
from
starter
group
by
npdes;

/*
Get
the
maximum
number
of
different
weights
among
the
potws
*/
create
table
maxwg
as
select
max(
numwghts)
as
maxnum
from
maxw;

/*
Save
the
maximum
number
of
weights
in
a
macro
variable.
This
is
the
*/
/*
number
of
build
loops
needed
to
exhaust
the
data.
*/
select
maxnum
into
:
max_
loop
from
maxwg;
quit;

/*
Run
the
initial
build
loop
*/
%
buildloop(
first=
1);

/*
Use
a
temporary
data
set
of
1
observation
to
allow
for
the
looping
of
*/
/*
of
the
build
events.
*/
data
_
temp;
set
maxwg;

/*
We
did
one
loop
already,
so
subtract
it
from
the
total
needed
*/
stopit
=
&
max_
loop
­
1;

/*
Call
the
macro
for
each
of
the
remaining
loops
needed
*/
do
i
=
1
to
stopit;
call
execute('%
buildloop(
first=
0)');
end;
run;
247
/*
Save
the
output
data
set
as
the
event
file
for
use
in
the
remainder
of
*/
/*
this
program
and
in
the
sludge
benefits
program
(
Sub
10).
*/
data
rundat.
potwevnt;
set
outset;

/*
rename
*/
evnt_
flo
=
potflow;
event_
wt
=
weights;
run;

/*
Calculate
influent
concentrations
*/
proc
sql;
create
table
potw2
as
select
npdes,
casno,
sum(
base)/
avg(
evnt_
flo)
as
concb,
sum(
treat1)/
avg(
evnt_
flo)
as
conct1,
sum(
treat2)/
avg(
evnt_
flo)
as
conct2,
event_
wt
from
rundat.
potwevnt
group
by
casno,
npdes,
event_
wt
order
by
casno,
npdes,
event_
wt;
quit;

/*
Get
Inhibition
Values
*/

data
limits;
set
mpm2.
chemdata;
if
inhibit;
keep
casno
inhibit;
run;

proc
sort;
by
casno;
quit;

/*
Compare
influent
concentrations
with
inhibition
criteria
*/

data
inhibit;
merge
potw2
limits;
by
casno;
if
inhibit;

keep
npdes
casno
inhibit
concb
conct1
conct2
event_
wt;

if
concb
<
inhibit
then
testb=
0;
else
testb=
1;
if
conct1
<
inhibit
then
testn1=
0;
else
testn1=
1;
if
conct2
<
inhibit
then
testn2=
0;
else
testn2=
1;

keep
npdes
casno
testb
testn1
testn2
event_
wt;
run;

proc
sort
data=
inhibit;
by
npdes
casno
event_
wt;
quit;

/*
Proc
SQL
calculates:

1.
Total
number
exceedence
events
by
pollutant;
2.
Number
of
exceedence
events
per
POTW/
option;
3.
Total
Number
of
affected
POTWs
on
a
sample/
national
basis
*/
248
proc
sql;

/*
Total
number
of
exceedence
events
by
pollutant
*/

create
table
inhbn0
as
select
inhibit.
casno,
min(
chemical)
as
chemname,
sum(
testb)
as
occurb,
sum(
testn1)
as
occurn1,
sum(
testn2)
as
occurn2
from
inhibit,
mpm2.
chemdata
where
testb>
0
and
inhibit.
casno=
chemdata.
casno
group
by
inhibit.
casno
order
by
inhibit.
casno;

/*
Total
number
of
exceedence
events
by
POTW/
option
*/

create
table
inhbn1
as
select
npdes,
sum(
testb)
as
excurb,
sum(
testn1)
as
excurn1,
sum(
testn2)
as
excurn2,
event_
wt,
case
when
sum(
testb)>
0
then
1
else
0
end
as
potwb,
case
when
sum(
testn1)>
0
then
1
else
0
end
as
potwn1,
case
when
sum(
testn2)>
0
then
1
else
0
end
as
potwn2
from
inhibit
where
testb>
0
group
by
npdes,
event_
wt
order
by
npdes,
event_
wt;

/*
Total
Number
of
Exceedence
Events
by
Chemical
on
a
National
Basis
*/

create
table
inh_
ntl
as
select
chemical,
sum(
testb*
event_
wt)
as
chemb,
sum(
testn1*
event_
wt)
as
chemn1,
sum(
testn2*
event_
wt)
as
chemn2
from
inhibit,
mpm2.
chemdata
where
testb>
0
and
inhibit.
casno=
chemdata.
casno
group
by
chemical
order
by
chemical;

/*
Total
Number
of
Exceedence
Events
on
a
National
Basis
*/

create
table
inh_
ntl2
as
select
sum(
excurb*
event_
wt)
label='
Number
of
Exceedence
Events
(
Baseline)',
sum(
potwb*
event_
wt)
label='
Number
of
Affected
POTWs
(
Baseline)',
sum
(
excurn1*
event_
wt)
label='
Number
of
Exceedence
Events
(
Opt
1)',
sum(
potwn1*
event_
wt)
label='
Number
of
Affected
POTWs
(
Opt
1)',
sum
(
excurn2*
event_
wt)
label='
Number
of
Exceedence
Events
(
Opt
2)',
sum(
potwn2*
event_
wt)
label='
Number
of
Affected
POTWs
(
Opt
2)'
from
inhbn1;

/*
Total
Number
of
Exceedence
Events
on
a
Sample
Basis
*/

create
table
inh_
smpl
as
select
excurb
label='
Number
of
Exceedence
Events
(
Baseline)',
potwb
label='
Number
of
Affected
POTWs
(
Baseline)',
excurn1
label='
Number
of
Exceedence
Events
(
Opt
1)',
potwn1
label='
Number
of
Affected
POTWs
(
Opt
1)',
249
excurn2
label='
Number
of
Exceedence
Events
(
Opt
2)',
potwn2
label='
Number
of
Affected
POTWs
(
Opt
2)'
from
inhbn1;
quit;

proc
print
data=
inh_
ntl2
label
;

title3
'
Total
Number
of
POTWs
Affected
by
Inhibition
Problems
caused
by';
title4
'
MP&
M
discharges
on
a
national
basis.';

proc
print
data=
inh_
ntl;
label
chemb
='
Number
of
Exceedence
Events
(
Baseline)'
chemn1
='
Number
of
Exceedence
Events
(
Opt
1)'
chemn2
='
Number
of
Exceedence
Events
(
Opt
2)';
title3
'
Total
number
of
exceedence
events
caused
by
MP&
M';
title4
'
discharges
by
chemical
on
a
national
basis';

proc
print
data=
inhbn0;
label
occurb
='
Number
of
Exceedence
Events
(
Baseline)'
occurn1
='
Number
of
Exceedence
Events
(
Opt
1)'
occurn2
='
Number
of
Exceedence
Events
(
Opt
2)';
title3
'
Total
number
of
excedence
events
caused
MP&
M
discharges';
title4
'
by
chemical
on
a
sample
basis';

proc
print
data=
inh_
smpl;
sum
excurb
excurn1
excurn2
potwb
potwn1
potwn2;
label
excurb
=
'
Number
of
Exceedence
Events
(
Baseline)'
excurn1
=
'
Number
of
Exceedence
Events
(
Opt
1)'
excurn2
=
'
Number
of
Exceedence
Events
(
Opt
2)'
potwb=
'
Number
of
Affected
POTWs
(
Baseline)'
potwn1=
'
Number
of
Affected
POTWs
(
Opt
1)'
potwn2=
'
Number
of
Affected
POTWs
(
Opt
2)';
title3
'
Total
Number
of
Exceedence
Events
by
NPDES
on
a
sample
basis';
quit;

data
reduc;
set
inhbn1;

if
excurn1
<
excurb
and
excurn1
>
0
then
red1
=
1;
else
red1
=
0;
if
excurn2
<
excurb
and
excurn2
>
0
then
red2
=
1;
else
red2
=
0;
run;

proc
means
sum
data=
reduc;
var
red1
red2;
label
red1
=
'#
POTWs
with
red.
poll.
inhibit.(
Opt1)'
red2
=
'#
POTWs
with
red.
poll.
inhibit.(
Opt2)';
title3
'
POTWS
with
Reduced
Inhibitions';
title4
'
SAMPLE';
quit;

proc
means
sum
data=
reduc;
weight
event_
wt;
250
var
red1
red2;
label
red1
=
'#
POTWs
with
red.
poll.
inhibit.(
Opt1)'
red2
=
'#
POTWs
with
red.
poll.
inhibit.(
Opt2)';
title3
'
POTWS
with
Reduced
Inhibitions';
title4
'
NATIONAL';
quit;

/*
End
POTW
Inhibitions
(
Sub
9).
sas
*/
251
Program
Output
252
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
65
POTW
Inhibitions
(
Sub
9)
15:
45
Thursday,
February
6,
2003
Total
Number
of
POTWs
Affected
by
Inhibition
Problems
caused
by
MP&
M
discharges
on
a
national
basis.

Number
of
Number
of
Number
of
Number
of
Number
of
Number
of
Exceedence
Affected
Exceedence
Affected
Exceedence
Affected
Events
POTWs
Events
POTWs
(
Opt
Events
POTWs
(
Opt
Obs
(
Baseline)
(
Baseline)
(
Opt
1)
1)
(
Opt
2)
2)

1
139.186
51.1227
115.188
30.4689
0
0
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
66
POTW
Inhibitions
(
Sub
9)
15:
45
Thursday,
February
6,
2003
Total
number
of
exceedence
events
caused
by
MP&
M
discharges
by
chemical
on
a
national
basis
Obs
CHEMICAL
chemb
chemn1
chemn2
1
Chromium
30.4689
27.1252
0
2
Copper
27.1252
27.1252
0
3
Lead
39.2507
30.4689
0
4
Silver
42.3408
30.4689
0
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
70
POTW
Inhibitions
(
Sub
9)
15:
45
Thursday,
February
6,
2003
POTWS
with
Reduced
Inhibitions
NATIONAL
The
MEANS
Procedure
Variable
Label
Sum
                                                                  
red1
#
POTWs
with
red.
poll.
inhibit.(
Opt1)
3.3436324
red2
#
POTWs
with
red.
poll.
inhibit.(
Opt2)
0
                                                                  
253
Sludge
Benefits
254
"
POTW
Inhibitions
(
Sub
9).
sas"
Using
the
set
of
loadings
discharged
to
POTWs
by
MP&
M
facilities,
POTW
specific
pollutant
events
are
calculated
The
data
set
of
these
events
is
outputed.
See
POTW
Inhibitions
from
MP&
M
Facilities
for
program
code
and
output.
255
"
Sldg
MP&
M
Based
(
Sub
10).
sas"
Using
the
data
set
of
POTW
pollutants
events,
MP&
M
pollutant
concentrations
in
sewage
sludge
are
estimated.
These
concentrations
are
compared
with
use/
disposal
criteria
for
given
disposal
practices.
Shifts
in
POTW
disposal
practices
from
reduced
MP&
M
loadings
are
determined
and
monetized.
256
/*
Program:
Sldg
MP&
M
Based
(
Sub
10).
sas
Authors:
Stuart
Smith,
Elena
Besedin,
Rob
Sartain
(
based
on
sldg.
sas
for
the
IL
project.)
Date
Modified:
08/
12/
2002
Modified
By:
Stuart
Smith,
Abt.
Associates,
Inc.

1.
Calculates
MP&
M
pollutant
concentrations
in
the
sewage
sludge
generated
by
POTWs
associated
with
the
MP&
M
sample
facilities;
2.
Compares
the
estimated
sludge
pollutant
concentrations
with
the
the
sewage
sludge
use/
disposal
criteria;
3.
Identifies
POTWs
in
which
pollutant
concentrations
exceed
sewage
sludge
use/
disposal
criteria
for
a
given
use/
disposal
practice;
4.
Estimates
shifts
in
sewage
sludge
use/
disposal
practices
due
to
the
reduced
MP&
M
discharges.
*/

title2
'
Sludge
MP&
M
Based
(
Sub
10)';

/*
Separate
potws
by
flow
volume
STRATUM
=
1;
small
size
(
daily
flow
rate
from
2
to
10
mg/
day)
STRATUM
=
2;
medium
size
(
daily
flow
rate
from
10
to
50
mg/
day)
STRATUM
=
3;
large
size
(
daily
flow
rate
>
50
mg/
day)

convert
the
event
flow
in
L/
day
to
MG/
day
dayflow
=
(
evnt_
flo*
0.2641721)/
1000000;

if
dayflow
<
10
then
strat=
1;
if
dayflow
>=
10
and
dayflow
<
50
then
strat=
2;
if
dayflow
>=
50
then
strat=
3;
*/

/*
Bring
in
the
potw
event
loadings
file
and
generate
stratums
*/
proc
sql;

create
table
potw_
0
as
select
npdes,
casno,
base
as
ldb,
treat1
as
ldt1,
treat2
as
ldt2,
event_
wt
as
smplwt,
evnt_
flo
as
potwflow,
facs,
size
as
facsze,
case
when
((
evnt_
flo*
0.2641721)/
1000000)
<
10
then
1
when
((
evnt_
flo*
0.2641721)/
1000000)
>=
50
then
3
else
2
end
as
stratum
from
rundat.
potwevnt
order
by
npdes,
casno,
smplwt;

/*
Calculate
Influent
Concentrations
and
obatin
sludge
criterion
*/

create
table
potw_
t
as
select
npdes,
stratum,
casno,
facs,(
ldb/
potwflow)
as
concb,(
ldt1/
potwflow)
as
conc1,(
ldt2/
potwflow)
as
conc2,
la_
h,
la_
l,
sludcrit,
surfdis,
257
smplwt,
twf_
remr.
name,
potw_
rem
as
potwrem,
facsze
as
size
from
potw_
0,
mpm2.
sldg_
cri,
mpm2.
twf_
remr
where
potw_
0.
casno
=
sldg_
cri.
casnum
and
potw_
0.
casno
=
twf_
remr.
cas
and
potw_
rem
is
not
missing
order
by
npdes,
casno,
smplwt;
quit;

/*
Calculate
Sludge
Concentrations
*/
/*
sewage
sludge
generation
factor
is
5.96
*/

data
conc_
sl;
set
potw_
t;

/*
Calculate
the
total
concentrations
*/
if
stratum
=
1
then
do;
if
size
=
'
s'
then
avgfac
=
32.81;
if
size
=
'
m'
then
avgfac
=
2.48;
if
size
=
'
l'
then
avgfac
=
1.21;

if
casno
=
'
7439921'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Lead*/
if
casno
=
'
7440382'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Arse*/
if
casno
=
'
7440020'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Nick*/
if
casno
=
'
7440439'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Cadm*/
if
casno
=
'
7440508'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Copp*/
if
casno
=
'
7440666'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Zinc*/
if
casno
=
'
7782492'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Seln*/
if
casno
=
'
7439976'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Merc*/

/*
Calculate
loadings
from
non­
MP&
M
facilities.
These
loadings
don't
change
under
the
post­
complaince
scenario.
*/
cnc1
=
(
conc1*
avgfac)/
facs;
cnc2
=
(
conc2*
avgfac)/
facs;
end;
if
stratum
=
2
then
do;
if
size
=
'
s'
then
avgfac
=
72.12;
if
size
=
'
m'
then
avgfac
=
7.99;
if
size
=
'
l'
then
avgfac
=
2.66;

if
casno
=
'
7439921'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Lead*/
if
casno
=
'
7440382'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Arse*/
if
casno
=
'
7440020'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Nick*/
if
casno
=
'
7440439'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Cadm*/
if
casno
=
'
7440508'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Copp*/
if
casno
=
'
7440666'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Zinc*/
if
casno
=
'
7782492'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Seln*/
if
casno
=
'
7439976'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Merc*/

/*
Calculate
loadings
from
non­
MP&
M
facilities.
These
loadings
don't
change
under
the
post­
complaince
scenario.
*/
cnc1
=
(
conc1*
avgfac)/
facs;
cnc2
=
(
conc2*
avgfac)/
facs;
end;
if
stratum
=
3
then
do;
258
if
size
=
'
s'
then
avgfac
=
147.68;
if
size
=
'
m'
then
avgfac
=
24.46;
if
size
=
'
l'
then
avgfac
=
10.43;

if
casno
=
'
7439921'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Lead*/
if
casno
=
'
7440382'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Arse*/
if
casno
=
'
7440020'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Nick*/
if
casno
=
'
7440439'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Cadm*/
if
casno
=
'
7440508'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Copp*/
if
casno
=
'
7440666'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Zinc*/
if
casno
=
'
7782492'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Seln*/
if
casno
=
'
7439976'
then
cncb=(
avgfac*
concb)/(
facs);
/*
Merc*/

/*
Calculate
loadings
from
non­
MP&
M
facilities.
These
loadings
don't
change
under
the
post­
complaince
scenario.
*/
cnc1
=
(
conc1*
avgfac)/
facs;
cnc2
=
(
conc2*
avgfac)/
facs;
end;
run;

/*
Sum
up
the
loadings
across
chemicals
and
POTWs
*/

/*
Separate
POTWS
by
size
*/
proc
sql;

create
table
small
as
select
distinct
npdes,
1
as
small
from
conc_
sl
where
size='
s'
order
by
npdes;

create
table
med
as
select
distinct
npdes,
1
as
med
from
conc_
sl
where
size='
m'
order
by
npdes;

create
table
large
as
select
distinct
npdes,
1
as
large
from
conc_
sl
where
size='
l'
order
by
npdes;

/*
Sum
over
checmicals
and
POTWs
*/
create
table
bypotw
as
select
npdes,
casno,
min(
stratum)
as
stratum,
sum(
cncb)
as
bas,
sum(
cnc1)
as
treat1,
sum(
cnc2)
as
treat2,
avg(
la_
h)
as
la_
h,
avg(
la_
l)
as
la_
l,
avg(
sludcrit)
as
sludcrit,
avg(
surfdis)
as
surfdis,
smplwt,
avg(
potwrem)
as
potwrem,
min(
name)
as
name
from
conc_
sl
group
by
npdes,
casno,
smplwt
order
by
npdes,
casno,
smplwt;
quit;
259
/*
Merge
sums
and
sizes
*/
data
concsl2;
merge
bypotw
small
med
large;
by
npdes;

if
small=.
then
small
=
0;
if
med=.
then
med
=
0;
if
large=.
then
large
=
0;

array
sm{*}
sm1
­
sm8;
array
md{*}
md1
­
md8;
array
lg{*}
lg1
­
lg8;

/*
Percentage
of
chemical
coming
from
MP&
M
facilities
based
on
National
Sludge
Survey
of
POTWS
*/
sm1=
0.138;
sm2=
0.074;
sm3=
0.251;
sm4=
0.161;
sm5=
0.189;
sm6=
0.202;
sm7=
0.072;
sm8=
0.079;

md1=
0.198;
md2=
0.140;
md3=
0.244;
md4=
0.234;
md5=
0.216;
md6=
0.160;
md7=
0.085;
md8=
0.208;

lg1=
0.103;
lg2=
0.070;
lg3=
0.158;
lg4=
0.128;
lg5=
0.109;
lg6=
0.082;
lg7=
0.033;
lg8=
0.060;

/*
Calculate
the
amount
of
sludge
coming
from
Non­
MP&
M
facilities
*/
if
casno
=
'
7439921'
then
other=
bas/(
sm1*
small
+
md1*
med
+
lg1*
large);/*
Lead*/
if
casno
=
'
7440382'
then
other=
bas/(
sm2*
small
+
md2*
med
+
lg2*
large);/*
Arse*/
if
casno
=
'
7440020'
then
other=
bas/(
sm3*
small
+
md3*
med
+
lg3*
large);/*
Nick*/
if
casno
=
'
7440439'
then
other=
bas/(
sm4*
small
+
md4*
med
+
lg4*
large);/*
Cadm*/
if
casno
=
'
7440508'
then
other=
bas/(
sm5*
small
+
md5*
med
+
lg5*
large);/*
Copp*/
if
casno
=
'
7440666'
then
other=
bas/(
sm6*
small
+
md6*
med
+
lg6*
large);/*
Zinc*/
if
casno
=
'
7782492'
then
other=
bas/(
sm7*
small
+
md7*
med
+
lg7*
large);/*
Seln*/
if
casno
=
'
7439976'
then
other=
bas/(
sm8*
small
+
md8*
med
+
lg8*
large);/*
Merc*/

/*
Total
concentration
is
the
sum
of
MP&
M
loadings
and
loadings
from
Non­
MP&
M
facilities
*/
cncb
=
bas
+
other;
cnc1
=
treat1
+
other;
cnc2
=
treat2
+
other;

/*
Apply
sludge
generation
factor
and
account
for
removal
percentage
*/

sludgeb
=
cncb
*
(
potwrem/
100)
*
5.96;
sludge1
=
cnc1
*
(
potwrem/
100)
*
5.96;
sludge2
=
cnc2
*
(
potwrem/
100)
*
5.96;

/*
Calculate
Pre
and
Post
Compliance
Concentrations
*/

base
=
sludgeb;
post1
=
sludge1;
post2
=
sludge2;

run;
260
/*
Compare
Calculated
Sludge
Concentrations
With
the
Sludge
Disposal
Criteria*/

data
sludge1;
set
concsl2;

/*
Compare
to
limits
for
land­
applied
high
*/
if
base
<
la_
h
then
critb_
h
=
0;
else
critb_
h
=
1;
if
post1
<
la_
h
then
crit1_
h
=
0;
else
crit1_
h
=
1;
if
post2
<
la_
h
then
crit2_
h
=
0;
else
crit2_
h
=
1;

/*
Compare
to
limits
for
land­
applied
low
*/
if
base
<
la_
l
then
critb_
l
=
0;
else
critb_
l
=
1;
if
post1
<
la_
l
then
crit1_
l
=
0;
else
crit1_
l
=
1;
if
post2
<
la_
l
then
crit2_
l
=
0;
else
crit2_
l
=
1;

/*
Assume
those
not
qualified
for
la_
l
use
mix
disposal
*/
if
critb_
l
=
1
then
critb_
m
=
1;
else
critb_
m
=
0;
if
crit1_
l
=
1
then
crit1_
m
=
1;
else
crit1_
m
=
0;
if
crit2_
l
=
1
then
crit2_
m
=
1;
else
crit2_
m
=
0;

run;

/*
Calculate
Sludge
Quantity
produced
by
the
POTWs
associated
with
MP&
M
facilities
in
millions
of
dry
metric
tons
(
DMT)
*/
proc
sql;
create
table
sldg_
wht
as
select
npdes,((
poflow*
0.2641721*
365)/
1000000)*
0.6363
as
sl_
wght
from
rundat.
potw_
in
where
npdes
^=
'
'
and
npdes
^=
'
Not
Found'
order
by
npdes;
quit;

/*
Keep
one
observation
per
unique
NPDES
*/
proc
sort
nodupkey
data=
sldg_
wht;
by
npdes;
quit;

/*
Estimate
Number
of
Exceedence
Events
by
Stratum/
Chemicals
on
a
Sample
Basis
*/

proc
sql;

create
table
sldg_
1
as
select
stratum,
casno,
min(
name)
as
polltnt,
sum(
critb_
h)
as
occ_
bh,
sum(
crit1_
h)
as
occ_
1h,
sum(
crit2_
h)
as
occ_
2h,
sum(
critb_
l)
as
occ_
bl,
sum(
crit1_
l)
as
occ_
1l,
sum(
crit2_
l)
as
occ_
2l,
sum
(
critb_
m)
as
occ_
bm,
sum(
crit1_
m)
as
occ_
1m,
sum(
crit2_
m)
as
occ_
2m,
sum(
sl_
wght)
as
weight,
case
when
sum(
crit1_
h)>
0
then
sum(
sl_
wght)
else
0
end
as
weight1,
case
when
sum(
crit2_
h)>
0
then
sum(
sl_
wght)
else
0
261
end
as
weight2
from
sludge1,
sldg_
wht
where
critb_
h>
0
and
sludge1.
npdes
=
sldg_
wht.
npdes
group
by
stratum,
casno
order
by
stratum,
casno;

create
table
chem_
nat
as
select
stratum,
casno,
min(
name)
as
polltnt,
sum(
critb_
h*
smplwt)
as
occ_
bh,
sum(
crit1_
h*
smplwt)
as
occ_
1h,
sum(
crit2_
h*
smplwt)
as
occ_
2h,
sum(
critb_
l*
smplwt)
as
occ_
bl,
sum(
crit1_
l*
smplwt)
as
occ_
1l,
sum(
crit2_
l*
smplwt)
as
occ_
2l,
sum
(
critb_
m*
smplwt)
as
occ_
bm,
sum(
crit1_
m*
smplwt)
as
occ_
1m,
sum(
crit2_
m*
smplwt)
as
occ_
2m,
sum(
sl_
wght*
smplwt)
as
weight,
case
when
sum(
crit1_
h)>
0
then
sum(
sl_
wght*
smplwt)
else
0
end
as
weight1,
case
when
sum(
crit2_
h)>
0
then
sum(
sl_
wght*
smplwt)
else
0
end
as
weight2
from
sludge1,
sldg_
wht
where
critb_
h>
0
and
sludge1.
npdes
=
sldg_
wht.
npdes
group
by
stratum,
casno
order
by
stratum,
casno;

/*
Estimate
Number
of
Exceedence
Events
by
POTW/
Option
on
a
Sample
Basis
*/

create
table
sldg_
2
as
select
sludge1.
npdes,
stratum,
sum(
critb_
h)
as
hla_
b,
sum(
crit1_
h)
as
hla_
1,
sum(
crit2_
h)
as
hla_
2,
sum(
critb_
l)
as
lla_
b,
sum(
crit1_
l)
as
lla_
1,
sum(
crit2_
l)
as
lla_
2,
sum(
critb_
m)
as
mixb,
sum(
crit1_
m)
as
mix1,
sum(
crit2_
m)
as
mix2,
smplwt,
case
when
sum(
critb_
h)>
0
then
avg(
sl_
wght)
else
0
end
as
weight,
case
when
sum(
crit1_
h)>
0
then
avg(
sl_
wght)
else
0
end
as
weight1,
case
when
sum(
crit2_
h)>
0
then
avg(
sl_
wght)
else
0
end
as
weight2,
case
when
sum(
critb_
h)>
0
then
1
else
0
end
as
potw_
hb,
case
when
sum(
crit1_
h)
>
0
then
1
else
0
262
end
as
potw_
h1,
case
when
sum(
crit2_
h)
>
0
then
1
else
0
end
as
potw_
h2,
case
when
sum(
critb_
l)>
0
then
1
else
0
end
as
potw_
lb,
case
when
sum(
crit1_
l)>
0
then
1
else
0
end
as
potw_
l1,
case
when
sum(
crit2_
l)>
0
then
1
else
0
end
as
potw_
l2,
case
when
sum
(
critb_
m)>
0
then
1
else
0
end
as
potw_
mb,
case
when
sum(
crit1_
m)>
0
then
1
else
0
end
as
potw_
m1,
case
when
sum(
crit2_
m)>
0
then
1
else
0
end
as
potw_
m2
from
sludge1,
sldg_
wht
where
critb_
h>
0
and
sludge1.
npdes
=
sldg_
wht.
npdes
group
by
stratum,
sludge1.
npdes,
smplwt
order
by
stratum,
npdes,
smplwt;
quit;

data
sludge2;
set
sldg_
2;
proc
sort
nodupkey;
by
stratum
npdes
smplwt;
quit;

proc
sql;

/*
Estimate
the
Total
Number
of
Exceedence
Events
on
a
National
Basis
by
Sewage
Sludge
Disposal
Practice/
Option
*/

create
table
ntl_
slg1
as
select
stratum,
sum(
hla_
b*
smplwt)
as
ex_
bh,
sum(
hla_
1*
smplwt)
as
ex_
1h,
sum(
hla_
2*
smplwt)
as
ex_
2h,
sum(
lla_
b*
smplwt)
as
ex_
bl,
sum(
lla_
1*
smplwt)
as
ex_
1l,
sum(
lla_
2*
smplwt)
as
ex_
2l,
sum(
mixb*
smplwt)
as
ex_
bm,
sum(
mix1*
smplwt)
as
ex_
1m,
sum(
mix2*
smplwt)
as
ex_
2m,
sum(
weight*
smplwt)
as
t_
wght,
sum(
weight1*
smplwt)
as
t_
1wght,
sum(
weight2*
smplwt)
as
t_
2wght
from
sludge2
where
hla_
b>
0
group
by
stratum
263
order
by
stratum;

/*
Estimate
the
Total
Number
of
POTWs
That
Don't
Meet
a
Given
Sludge
Disposal
Standard
*/

create
table
ntl_
slg2
as
select
stratum,
sum(
potw_
hb*
smplwt)
as
ex_
bh,
sum(
potw_
h1*
smplwt)
as
ex_
1h,
sum(
potw_
h2*
smplwt)
as
ex_
2h,
sum(
potw_
lb*
smplwt)
as
ex_
bl,
sum(
potw_
l1*
smplwt)
as
ex_
1l,
sum(
potw_
l2*
smplwt)
as
ex_
2l,
sum(
potw_
mb*
smplwt)
as
ex_
bm,
sum(
potw_
m1*
smplwt)
as
ex_
1m,
sum(
potw_
m2*
smplwt)
as
ex_
2m,
sum(
weight*
smplwt)
as
t_
wght,
sum(
weight1*
smplwt)
as
t_
1wght,
sum(
weight2*
smplwt)
as
t_
2wght
from
sludge2
where
hla_
b>
0
group
by
stratum
order
by
stratum;
quit;

/*
Indentify
Shifts
in
Sludge
Use/
Disposal
Practices
from
Land
Application
Low
to
Land
Application
High
on
a
Sample
Basis
*/
proc
sql;

create
table
shifts1
as
select
stratum,
npdes,
weight,
count(
distinct
npdes)
as
upgrade,
smplwt,
case
when
hla_
1=
0
then
'
lla_
llh'
when
hla_
1>
0
and
lla_
1=
0
then
'
nswitch'
else
'
downgrd'
end
as
switch1,
case
when
hla_
2=
0
then
'
lla_
llh'
when
hla_
2>
0
and
lla_
2=
0
then
'
nswitch'
else
'
downgrd'
end
as
switch2
from
sludge2
where
hla_
b
>
0
and
lla_
b
=
0
order
by
stratum,
npdes,
smplwt;
quit;

data
shift11;
set
shifts1;
proc
sort;
by
stratum
npdes
switch1
switch2;
quit;

/*
Indentify
Shifts
in
the
Sewage
Sludge
Disposal
Practices
from
mix
to
(
1)
Land
Application
High;
(
2)
Land
Application
Low
;
*/

proc
sql;
create
table
shifts3
as
select
stratum,
npdes,
weight,
count(
distinct
npdes)
as
upgrade,
smplwt,
case
264
when
hla_
1=
0
then
'
mix_
llh'
when
hla_
1>
0
and
lla_
1=
0
then
'
mix_
lla'
else
'
nswitch'
end
as
switch1,
case
when
hla_
2=
0
then
'
mix_
llh'
when
hla_
2>
0
and
lla_
2=
0
then
'
mix_
lla'
else
'
nswitch'
end
as
switch2
from
sludge2
where
mixb>
0
order
by
stratum,
npdes,
smplwt;
quit;

data
shift13;
set
shifts3;
proc
sort;
by
stratum
npdes
switch1
switch2;
quit;

/*
Calculate
Savings
From
Shifts
in
Sewage
Sludge
Disposal
Practices
Note
that
Savings
From
the
Reduced
Record
Keeping
is
Estimated
Separately
*/

data
savings;
set
shift11
shift13;
if
npdes
ne
'
';

keep
stratum
npdes
switch1
switch2
weight
smplwt;
run;

/*
Estimate
Cost
Savings
from
Shifts
to
Less
Costly
Sewage
Sludge
Disposal
Practices
on
a
Sample
Basis
The
amount
newly
qualified
(
weight)
sludge
must
be
adjusted
since
not
all
can
be
land
applied.
43.3%
of
the
benefitting
sludge
can
be
land
applied.
Of
the
remaining
56.7%
of
the
sludge,
11%
can
be
bagged
and
sold,
which
is
6.4%
of
the
total
newly
qualified.

So,
the
benefitting
sludge
is
really
the
43.3%
that
can
be
land
applied
AND
the
6.4%
which
is
bagged
and
sold.

Thus,
49.7%
of
the
total
newly
qualified
sludge
for
land
application
is
used
in
the
cost
savings
estimates.
*/

proc
sql;

create
table
svngs
as
select
stratum,
npdes,
switch1,
switch2,
weight,
smplwt,
case
when
switch1='
lla_
llh'
then
weight*
0.684
when
switch1='
mix_
llh'
then
0.497*
weight*
0.75*
0.39*
103.82
+
0.497*
weight*
0.75*
0.61*
100.44
+
0.497*
weight*
0.25*
0.39*
95.91
+
0.497*
weight*
0.25*
0.61*
69.96
when
switch1='
mix_
lla'
then
0.497*
weight*
0.39*
103.82
+
0.497*
weight*
0.61*
100.44
when
switch1='
nswitch'
then
0
265
else
0
end
as
savel1,
case
when
switch2='
lla_
llh'
then
weight*
0.684
when
switch2='
mix_
llh'
then
0.497*
weight*
0.75*
0.39*
103.82
+
0.497*
weight*
0.75*
0.61*
100.44
+
0.497*
weight*
0.25*
0.39*
95.91
+
0.497*
weight*
0.25*
0.61*
69.96
when
switch2='
mix_
lla'
then
0.497*
weight*
0.39*
103.82
+
0.497*
weight*
0.61*
100.44
when
switch2='
nswitch'
then
0
else
0
end
as
savel2,
case
when
switch1='
lla_
llh'
then
weight*
1.362
when
switch1='
mix_
llh'
then
0.497*
weight*
0.75*
0.39*
103.82
+
0.497*
weight*
0.75*
0.61*
100.44
+
0.497*
weight*
0.25*
0.39*
95.91
+
0.497*
weight*
0.25*
0.61*
69.96
when
switch1='
mix_
lla'
then
0.497*
weight*
0.39*
103.82
+
0.497*
weight*
0.61*
100.44
when
switch1='
nswitch'
then
0
else
0
end
as
saveh1,
case
when
switch2='
lla_
llh'
then
weight*
1.362
when
switch2='
mix_
llh'
then
0.497*
weight*
0.75*
0.39*
103.82
+
0.497*
weight*
0.75*
0.61*
100.44
+
0.497*
weight*
0.25*
0.39*
95.91
+
0.497*
weight*
0.25*
0.61*
69.96
when
switch2='
mix_
lla'
then
0.497*
weight*
0.39*
103.82
+
0.497*
weight*
0.61*
100.44
when
switch2='
nswitch'
then
0
else
0
end
as
saveh2
from
savings
order
by
stratum,
npdes,
switch1,
switch2;
quit;

/*
Estimate
Cost
Savings
from
Shifts
to
Less
Costly
Sewage
Sludge
Disposal
Practices
on
a
National
Basis
*/

proc
sql;
create
table
ntlshft1
as
select
stratum,
switch1,
sum(
smplwt)
as
potwntl,
sum(
weight*
smplwt)
as
totwght,
sum(
savel1*
smplwt)
as
lowsave1,
sum(
saveh1*
smplwt)
as
hghsave1
from
svngs
group
by
stratum,
switch1
order
by
stratum,
switch1;

create
table
ntlshft2
as
select
stratum,
switch2,
sum(
smplwt)
as
potwntl,
sum(
weight*
smplwt)
as
totwght,
sum(
savel2*
smplwt)
as
lowsave2,
sum(
saveh2*
smplwt)
as
hghsave2
from
svngs
group
by
stratum,
switch2
order
by
stratum,
switch2;
266
quit;

proc
print
data=
sldg_
1;
label
occ_
bh=
'
Execeedence
of
LA­
High
Value
(
Baseline)'
occ_
1h=
'
Execeedence
of
LA­
High
Value
(
Option
1)'
occ_
2h=
'
Execeedence
of
LA­
High
Value
(
Option
2)'
occ_
bl=
'
Execeedence
of
LA­
Low
Value
(
Baseline)'
occ_
1l=
'
Execeedence
of
LA­
Low
Value
(
Option
1)'
occ_
2l=
'
Execeedence
of
LA­
Low
Value
(
Option
2)'
occ_
bs=
'
Execeedence
of
Surf.
Disp.
Value
(
Baseline)'
occ_
1s=
'
Execeedence
of
Surf.
Disp.
Value
(
Option
1)'
occ_
2s=
'
Execeedence
of
Surf.
Disp.
Value
(
Option
2)'
weight
=
'
Sludge
Quantity
DMT
(
Baseline)'
weight1
=
'
Sludge
Quantity
DMT
(
Option
1)'
weight2
=
'
Sludge
Quantity
DMT
(
Option
2)';
title3
'
Sewage
Sludge
Contamination
Problems
Caused
by
';
title4
'
MP&
M
dischages
by
Pollutant
(
Sample)';
quit;

proc
print
data=
chem_
nat;
label
occ_
bh=
'
Execeedence
of
LA­
High
Value
(
Baseline)'
occ_
1h=
'
Execeedence
of
LA­
High
Value
(
Option
1)'
occ_
2h=
'
Execeedence
of
LA­
High
Value
(
Option
2)'
occ_
bl=
'
Execeedence
of
LA­
Low
Value
(
Baseline)'
occ_
1l=
'
Execeedence
of
LA­
Low
Value
(
Option
1)'
occ_
2l=
'
Execeedence
of
LA­
Low
Value
(
Option
2)'
occ_
bs=
'
Execeedence
of
Surf.
Disp.
Value
(
Baseline)'
occ_
1s=
'
Execeedence
of
Surf.
Disp.
Value
(
Option
1)'
occ_
2s=
'
Execeedence
of
Surf.
Disp.
Value
(
Option
2)'
weight
=
'
Sludge
Quantity
DMT
(
Baseline)'
weight1
=
'
Sludge
Quantity
DMT
(
Option
1)'
weight2
=
'
Sludge
Quantity
DMT
(
Option
2)';
title3
'
Sewage
Sludge
Contamination
Problems
Caused
by
';
title4
'
MP&
M
dischages
by
Pollutant
(
National)';
quit;

proc
means
sum
data=
sldg_
2;
var
hla_
b
hla_
1
hla_
2
lla_
b
lla_
1
lla_
2
mixb
mix1
mix2
weight
weight1
weight2
potw_
hb
potw_
h1
potw_
h2
potw_
lb
potw_
l1
potw_
l2
potw_
mb
potw_
m1
potw_
m2;
label
hla_
b='
Execeedences
of
LA­
High
Value
(
BL)'
hla_
1='
Execeedences
of
LA­
High
Value
(
Opt1)'
hla_
2='
Execeedences
of
LA­
High
Value
(
Opt2)'
lla_
b=
'
Execeedences
of
LA­
Low
Value
(
BL)'
lla_
1=
'
Execeedences
of
LA­
Low
Value
(
Opt1)'
lla_
2=
'
Execeedences
of
LA­
Low
Value
(
Opt2)'
mixb=
'
Execeedences
of
Mix
Value
(
BL)'
mix1=
'
Execeedences
of
Mix
Value
(
Opt1)'
mix2=
'
Execeedences
of
Mix
Value
(
Opt2)'
weight
=
'
Sludge
Quantity
DMT
(
Baseline)'
weight1
=
'
Sludge
Quantity
DMT
(
Option
1)'
weight2
=
'
Sludge
Quantity
DMT
(
Option
2)'
potw_
hb
=
'
Number
of
POTW
exceeding
LA­
High
(
BL)'
267
potw_
h1
=
'
Number
of
POTW
exceeding
LA­
High
(
Opt1)'
potw_
h2
=
'
Number
of
POTW
exceeding
LA­
High
(
Opt2)'
potw_
lb
=
'
Number
of
POTW
exceeding
LA­
Low
(
Bl)'
potw_
l1
=
'
Number
of
POTW
exceeding
LA­
Low
(
Opt1)'
potw_
l2
=
'
Number
of
POTW
exceeding
LA­
Low
(
Opt2)'
potw_
mb
=
'
Number
of
POTW
exceeding
Mix
(
Bl)'
potw_
m1
=
'
Number
of
POTW
exceeding
Mix
(
Opt1)'
potw_
m2
=
'
Number
of
POTW
exceeding
Mix
(
Opt2)';

title3
'
Total
Number
of
Exceedances
by
Disposal
Method/
Stratum
on
(
Sample)';
quit;

proc
means
sum
data=
svngs;
var
weight
savel1
saveh1
savel2
saveh2;
label
switch1='
Shift
Category
Opt
1'
switch2='
Shift
Category
Opt
2'
weight
=
'
Sudge
DMT
Newly
Qualified
for
a
Practice'
savel1
=
'
Low
Cost
Savings
Estimates
Opt
1'
saveh1
=
'
High
Cost
Savings
Estimates
Opt
1'
savel2
=
'
Low
Cost
Savings
Estimates
Opt
2'
saveh2
=
'
High
Cost
Savings
Estimates
Opt
2';
title3
'
Cost
Savings
Resulting
from
Shifts
to
Less
Costly';
title4
'
Sewage
Sludge
Disposal
Practices
on
a
Sample
Basis';
quit;

proc
means
sum
data=
ntl_
slg2;
var
ex_
bh
ex_
1h
ex_
2h
ex_
bl
ex_
1l
ex_
2l
ex_
bm
ex_
1m
ex_
2m
t_
wght
t_
1wght
t_
2wght;
label
ex_
bh
='
Execeedences
of
LA­
High
Value
(
BL)'
ex_
1h='
Execeedences
of
LA­
High
Value
(
Opt1)'
ex_
2h='
Execeedences
of
LA­
High
Value
(
Opt2)'
ex_
bl=
'
Execeedences
of
LA­
Low
Value
(
BL)'
ex_
1l=
'
Execeedences
of
LA­
Low
Value
(
Opt1)'
ex_
2l=
'
Execeedences
of
LA­
Low
Value
(
Opt2)'
ex_
bm=
'
Execeedences
of
Mix
Value(
BL)'
ex_
1m=
'
Execeedences
of
Mix
Value
(
Opt1)'
ex_
2m=
'
Execeedences
of
Mix
Value
(
Opt2)'
t_
wght
=
'
Sludge
Quantity
DMT
(
BL)'
t_
1wght
=
'
Sludge
Quantity
DMT
(
Opt1)'
t_
2wght
=
'
Sludge
Quantity
DMT
(
Opt2)';
title3
'
Total
Number
of
POTWs
exceeding
a
Sludge
Disposal
Method(
National)';
quit;

proc
print
data=
ntlshft1;
sum
lowsave1
hghsave1;
label
switch1='
Shift
Category'
potwntl=
'
Number
of
Upgrading
POTWs'
totwght
=
'
Sludge
(
DMT)
Newly
Qualified
for
a
Practice'
Lowsave1
=
'
Low
Cost
Savings
Estimates'
Hghsave1
=
'
High
Cost
Savings
Estimates';

title3
'
Cost
Savings
Resulting
From
Shifts
to
Less
Costly';
title4
'
Sludge
Disposal
Practices
on
a
National
Basis';
title5
'
Option
1';
268
quit;

proc
print
data=
ntlshft2;
sum
lowsave2
hghsave2;
label
switch2='
Shift
Category'
potwntl=
'
Number
of
Upgrading
POTWs'
totwght
=
'
Sludge
(
DMT)
Newly
Qualified
for
a
Practice'
Lowsave2
=
'
Low
Cost
Savings
Estimates'
Hghsave2
=
'
High
Cost
Savings
Estimates';

title3
'
Cost
Savings
Resulting
From
Shifts
to
Less
Costly';
title4
'
Sludge
Disposal
Practices
on
a
National
Basis';
title5
'
Option
2';
quit;

/*
End
Sldg
MP&
M
Based
(
Sub
10).
sas
*/
269
Program
Output
270
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
72
Sludge
MP&
M
Based
(
Sub
10)
15:
45
Thursday,
February
6,
2003
Sewage
Sludge
Contamination
Problems
Caused
by
MP&
M
dischages
by
Pollutant
(
National)

Obs
stratum
CASNO
polltnt
occ_
bh
occ_
1h
occ_
2h
occ_
bl
occ_
1l
1
1
7439921
Lead
220.171
220.171
204.956
93.271
93.271
2
1
7440020
Nickel
722.796
710.563
626.443
722.796
710.563
3
1
7440382
Arsenic
83.688
83.688
4.072
4.072
4.072
4
1
7440439
Cadmium
564.304
549.494
542.807
477.722
477.722
5
1
7440508
Copper
516.049
482.625
470.753
99.958
99.958
6
1
7440666
Zinc
216.127
200.911
185.779
65.691
65.691
7
1
7782492
Selenium
109.428
109.428
109.428
31.572
31.572
8
2
7439921
Lead
5.627
5.627
0.000
0.000
0.000
9
2
7440020
Nickel
32.138
32.138
32.138
32.138
32.138
10
2
7440439
Cadmium
114.705
114.705
92.115
65.326
65.326
11
2
7440508
Copper
9.569
9.569
5.627
0.000
0.000
12
3
7439921
Lead
8.217
8.217
8.217
0.000
0.000
13
3
7440020
Nickel
8.217
8.217
8.217
8.217
8.217
14
3
7440439
Cadmium
74.518
74.518
71.174
32.214
32.214
15
3
7440508
Copper
8.217
8.217
8.217
0.000
0.000
16
3
7440666
Zinc
8.217
8.217
8.217
0.000
0.000
Obs
occ_
2l
occ_
bm
occ_
1m
occ_
2m
weight
weight1
weight2
1
93.271
93.271
93.271
93.271
60462.62
60462.62
60462.62
2
626.443
722.796
710.563
626.443
148866.70
148866.70
148866.70
3
4.072
4.072
4.072
4.072
11123.45
11123.45
11123.45
4
453.978
477.722
477.722
453.978
240931.32
240931.32
240931.32
5
99.958
99.958
99.958
99.958
77686.19
77686.19
77686.19
6
62.347
65.691
65.691
62.347
56456.92
56456.92
56456.92
7
31.572
31.572
31.572
31.572
19949.57
19949.57
19949.57
8
0.000
0.000
0.000
0.000
13377.19
13377.19
0.00
9
32.138
32.138
32.138
32.138
190996.12
190996.12
190996.12
10
65.326
65.326
65.326
65.326
579175.21
579175.21
579175.21
11
0.000
0.000
0.000
0.000
33997.92
33997.92
33997.92
12
0.000
0.000
0.000
0.000
130888.39
130888.39
130888.39
13
8.217
8.217
8.217
8.217
130888.39
130888.39
130888.39
14
32.214
32.214
32.214
32.214
1824663.82
1824663.82
1824663.82
15
0.000
0.000
0.000
0.000
130888.39
130888.39
130888.39
16
0.000
0.000
0.000
0.000
130888.39
130888.39
130888.39
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
75
Sludge
MP&
M
Based
(
Sub
10)
15:
45
Thursday,
February
6,
2003
Total
Number
of
POTWs
exceeding
a
Sludge
Disposal
Method(
National)

The
MEANS
Procedure
Variable
Label
Sum
                                                                
ex_
bh
Execeedences
of
LA­
High
Value
(
BL)
1019.76
ex_
1h
Execeedences
of
LA­
High
Value
(
Opt1)
1004.95
ex_
2h
Execeedences
of
LA­
High
Value
(
Opt2)
975.6693539
271
ex_
bl
Execeedences
of
LA­
Low
Value
(
BL)
856.0269055
ex_
1l
Execeedences
of
LA­
Low
Value
(
Opt1)
856.0269055
ex_
2l
Execeedences
of
LA­
Low
Value
(
Opt2)
852.6832731
ex_
bm
Execeedences
of
Mix
Value(
BL)
856.0269055
ex_
1m
Execeedences
of
Mix
Value
(
Opt1)
856.0269055
ex_
2m
Execeedences
of
Mix
Value
(
Opt2)
852.6832731
t_
wght
Sludge
Quantity
DMT
(
BL)
2663743.61
t_
1wght
Sludge
Quantity
DMT
(
Opt1)
2647195.25
t_
2wght
Sludge
Quantity
DMT
(
Opt2)
2319026.34
                                                                
MP&
M
National
Benefits
Preferred(
Opt1=
Pref,
Opt2=
NULL)
76
Sludge
MP&
M
Based
(
Sub
10)
15:
45
Thursday,
February
6,
2003
Cost
Savings
Resulting
From
Shifts
to
Less
Costly
Sludge
Disposal
Practices
on
a
National
Basis
Option
1
Obs
stratum
switch1
potwntl
totwght
lowsave1
hghsave1
1
1
lla_
llh
14.810
16548.36
11319.08
22538.87
2
1
nswitch
815.724
243356.22
0.00
0.00
3
2
nswitch
114.705
579175.21
0.00
0.00
4
3
nswitch
74.518
1824663.82
0.00
0.00
========
========
11319.08
22538.87