L.4.1.1 WATTOX.SAS program statements

are also provided. WATTOX.SAS also includes calculation of LC50 using the

Probit procedure.

L.4.1.1 WATTOX.SAS program statements

LIBNAME Q 'C:\SAS\SASFILES';

OPTIONS LINESIZE=79 PAGESIZE=59 NODATE NONUMBER;

/* Identify the treatment codes. */

PROC FORMAT;

VALUE TRTFMT

0='DILUTION WATER '

1='100 percent ELUTRIATE '

2='50 percent ELUTRIATE '

3='25 percent ELUTRIATE '

4='12.5 percent ELUTRIATE';

/* Input the toxicity test data after the CARDS statement,

listing the

/* treatment code, replicate, and number of survivors. A

permanent SAS

/* data set is created in the directory specified in the

LIBNAME statement. */

DATA Q.WATCOL;

INPUT TRT REP SURV @@;

CARDS;

0 1 20 0 2 19 0 3 20 0 4 20 0 5 19

116127139145158

2 1 8 2 2 8 2 3 9 2 4 10 2 5 11

3 1 12 3 2 18 3 3 15 3 4 14 3 5 13

4 1 17 4 2 17 4 3 18 4 4 16 4 5 18

;

/* Input no. of organisms (M) per test container at start of

test. */

/* Calculate proportion of survivors (SURV/M) and take the

SQRT. */

/* Arcsine transform SQRT(SURV/M). */

/* Format, print, sort the data. Print no. of observations,

mean, and */

/* standard error for survival in each treatment. */

DATA A0;

SET Q.WATCOL;

M=20;

ARCSURV=ARSIN(SQRT(SURV/M));

LABEL TRT='TREATMENT GROUP'

REP='REPLICATE'

M='NO. OF ORGANISMS PER REPLICATE'

SURV='NUMBER OF SURVIVORS'

ARCSURV='ARCSINE TRANSFORMATION';

FORMAT TRT TRTFMT.;

TITLE 'WATER COLUMN TOXICITY DATA';

PROC PRINT LABEL; VAR TRT REP M SURV ARCSURV;

PROC SORT; BY TRT;

PROC MEANS NOPRINT; BY TRT; VAR SURV;

OUTPUT OUT=Y N=N SUM=TOTAL MEAN=MEANSURV STDERR=SE;

PROC PRINT LABEL; VAR TRT N MEANSURV SE;

LABEL MEANSURV='MEAN SURVIVAL';

L41

Appendix L Selected Resource Documents