 |
||
|
|
||
|
Page Title:
L.4.2.1 BIOACC.SAS program statements |
||
| |||||||||||||||
|
|
 L.4.2.1 BIOACC.SAS program statements
LIBNAME Q 'C:\SAS\SASFILES';
OPTIONS LINESIZE=79 PAGESIZE=59 NODATE NONUMBER;
/* Identify the treatment codes. */
PROC FORMAT;
VALUE TRTFMT
1='REFERENCE '
2='SEDIMENT 1'
3='SEDIMENT 2'
4='SEDIMENT 3';
/*
Input the bioaccumulation data after the CARDS statement, listing the */
/*
treatment code, replicate, and contaminant concentration. A permanent */
/*
SAS data set is created in the directory specified in the LIBNAME */
/*
statement. */
DATA Q.BIOACC;
INPUT TRT REP CONC
@@;
CARDS;
1 1 .06 1 2 .05 1 3
.05
1
4
.08
1
5
.09
2 1 .16 2 2 .19 2 3
.18
2
4
.22
2
5
.31
3 1 .24 3 2 .10 3 3
.13
3
4
.18
3
5
.30
4 1 .13 4 2 .05 4 3
.17
4
4
.08
4
5
.22
;
/* Format, print, sort the data. Print no. of observations, mean, and */
/* standard error for concentration in each treatment for both */
/* untransformed and log10-transformed data. Calculate rankits. */
DATA A0;
SET Q.BIOACC;
LOGCONC=LOG10(CONC);
MERGEVAR=1;
LABEL TRT='TREATMENT GROUP'
REP='REPLICATE'
CONC='CONTAMINANT CONCENTRATION, ug/g'
LOGCONC='LOG10 CONCENTRATION';
FORMAT TRT TRTFMT.;
TITLE 'CONTAMINANT BIOACCUMULATION DATA';
PROC RANK NORMAL=BLOM OUT=A;
VAR CONC; RANKS RANKIT;
PROC PRINT LABEL; VAR TRT REP CONC LOGCONC RANKIT;
LABEL RANKIT='NORMALIZED RANK FOR CONCENTRATION';
PROC SORT; BY TRT;
PROC MEANS NOPRINT; BY TRT; VAR CONC LOGCONC; ID MERGEVAR;
OUTPUT OUT=Y N=N NLOG MEAN=MEANCONC MEANLOG VAR=S2 S2LOG STDERR=SE SELOG;
PROC PRINT LABEL; VAR TRT N MEANCONC S2 SE MEANLOG S2LOG SELOG;
LABEL
MEANCONC='MEAN CONTAMINANT CONC.'
S2='VARIANCE'
SE='STANDARD ERROR'
MEANLOG='MEAN LOG10 CONC.'
S2LOG='VARIANCE OF LOGS'
SELOG='STANDARD ERROR OF LOGS';
/* Test normality of residuals of untransformed and log-transformed data */
/* using Shapiro-Wilk's Test. */
L49
Appendix L Selected Resource Documents
|
|
Privacy Statement - Press Release - Copyright Information. - Contact Us - Support Integrated Publishing |
