L.4.2.1 BIOACC.SAS program statements (Cont.)

SET C3;

DO PCTDIFF=10,25,50,100,200,300;

SEDCONC=MEANCONC+((PCTDIFF/100)*MEANCONC);

D=SEDCONC-MEANCONC;

TBETA=D*SQRT(N/(2*MSE))-TALPHA;

POWER=PROBT(TBETA,DF);

OUTPUT;

END;

LABEL

PCTDIFF=' % INCREASE IN CONC. ABOVE REFERENCE'

SEDCONC='DREDGED SEDIMENT BIOACCUMULATION'

TBETA='T VALUE FOR (1-BETA,DF)'

POWER='POWER (1-BETA)';

PROC PRINT LABEL NOOBS; VAR PCTDIFF SEDCONC D TBETA POWER;

TITLE 'POWER OF LSD TO DETECT % INCREASE IN CONCENTRATION ABOVE REFERENCE';

TITLE2 'MEAN CONTAMINANT CONCENTRATION GIVEN N, MSE AND DF SHOWN ABOVE';

DATA C5;

SET C3;

DO POWER=.5,.6,.7,.8,.9,.95,.99;

TBETA=TINV(POWER,DF);

D=((TBETA+TALPHA)*SQRT(2*MSE))/SQRT(N);

SEDCONC=MEANCONC+D;

PCTDIFF=(D*100)/MEANCONC;

OUTPUT;

END;

LABEL

SEDCONC='DREDGED SEDIMENT BIOACCUMULATION'

PCTDIFF=' % INCREASE IN CONC. ABOVE REFERENCE'

TBETA='T VALUE FOR (1-BETA,DF)'

POWER='POWER (1-BETA)';

PROC PRINT LABEL NOOBS; VAR POWER D SEDCONC PCTDIFF TBETA;

TITLE 'MINIMUM DREDGED SEDIMENT BIOACCUMULATION THAT CAN BE DETECTED BY LSD';

TITLE2 'AS SIGNIFICANT GIVEN SPECIFIED POWER AND N, MSE, AND DF SHOWN ABOVE';

/* Calculation of upper confidence limits (UCL) for comparison of mean */

/* dredged sediment bioaccumulation with an action level. */

DATA D;

MERGE C1 Y; BY MERGEVAR;

IF TRT=1 THEN DELETE;

TALPHA1=TINV(.95,DF);

TALPHA2=TINV(.95,N-1);

UCL1=MEANCONC+TALPHA1*(SQRT(MSE/N));

UCL2=MEANCONC+TALPHA2*(SQRT(S2/N));

DMIN1=TALPHA1*SQRT(MSE/N);

DMIN2=TALPHA2*SQRT(S2/N);

LABEL

UCL1='UCL (EQUAL VARIANCES)'

UCL2='UCL (UNEQUAL VARIANCES)'

TALPHA1='T VALUE FOR (1-ALPHA=.95,DF)'

TALPHA2='T VALUE FOR (1-ALPHA=.95,N-1)'

DMIN1='MINIMUM SIGNIFICANT DIFFERENCE'

DMIN2='MINIMUM SIGNIFICANT DIFFERENCE'

MSE='MEAN SQUARE ERROR'

S2='VARIANCE'

MEANCONC='MEAN BIOACCUMULATION';

TITLE 'COMPARISON OF MEAN DREDGED SEDIMENT BIOACCUMULATION WITH ACTION

LEVEL:';

PROC PRINT LABEL NOOBS; VAR TRT MEANCONC UCL1 MSE TALPHA1 DF DMIN1;

TITLE2 'UPPER CONFIDENCE LIMITS (UCL) WHEN VARIANCES ARE EQUAL';

PROC PRINT LABEL NOOBS; VAR TRT MEANCONC UCL2 S2 TALPHA2 N DMIN2;

TITLE2 'UPPER CONFIDENCE LIMITS (UCL) WHEN VARIANCES ARE UNEQUAL';

/* Calculate power of dredged sediment-action level comparisons using */

L52

Appendix L Selected Resource Documents