   IN THIS EXAMPLE, STORET IS USED TO SELECT DATA FROM WELLS BOTH       00000400
   UP- AND DOWN-GRADIENT FROM A LAND-FILL SITE.                         00000500
   SAS IS INVOKED TO PERFORM VARIOUS STATISTICAL ANALYSES RELATED       00000600
   TO EPA RCRA REGULATIONS.                                             00000700
                                                                        00000800
   A ONE-SIDED COCHRAN'S APPROXIMATION TTEST IS PERFORMED.              00000900
                                                                        00001000
   STORET IS GRATEFUL TO ANDY OBRACHTA, CORPS OF ENGINEERS, TULSA       00001100
   OKLAHOMA FOR THIS EXAMPLE.                                           00001200
                                                                        00001210
   ANDY CAN BE REACHED AT (918)581-6170 OR FTS 745-6170                 00001220
   (PHONE NUMBERS CURRENT AS OF SEPTEMBER 1985)                         00001230
                                                                        00001300
                                                                        00005100
        ************************                                        00005200
        *    LANDFILL.TTEST    *                                        00005300
        ************************                                        00005400
                                                                        00005401
LANDFILL.TTEST USES PROC MEANS TO PRODUCE A T STATISTIC.                00005410
IT IS DIFFERENT FROM PROC TTEST IN THAT IT USES THE COCHRANS            00005411
APPROXIMATION AND IS A ONE-TAILED TEST (A ONE-TAILED TEST CHECKS FOR    00005412
INCREASES, WHILE A TWO-TAILED TEST CHECKS FOR EITHER INCREASES OR       00005413
DECREASES.)                                                             00005420
                                                                        00005500
PGM=RET,MORE=SAS,PRT=NO,                                                00005600
A=12OKSWT,                                                              00005700
S=177,S=178,S=179,S=180,S=181,S=182,S=183,S=184,S=185,S=186,S=187,S=188,00005800
P=1000,                                                                 00005900
NOECHO,                                                                 00006000
SASPARMS=BEGIN,                                                         00006100
DATA LANDFILL;                                                          00006300
*                                                                       00006400
INFILE FCF;                                                             00006500
INPUT @26 YY $2. @; IF YY^='99';                                        00006501
INPUT @1   STATION    $23.                                              00006502
      @36  ARSENIC    RB4.                                              00006503
      @236 RMK        $1.                                               00006504
         ;                                                              00006505
*  TREAT REMARK "K" AS HALFWAY BETWEEN ZERO AND VALUE GIVEN     ;       00006510
IF RMK='K' THEN ARSENIC=ARSENIC/2;                                      00006600
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00006610
*         IDENTIFY UP-GRADIENT WELLS AS SUCH.                     ;     00006620
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00006630
IF STATION='12OKSWT  183' THEN GRADIENT='UP';                           00006700
IF STATION='12OKSWT  185' THEN GRADIENT='UP';                           00006800
IF STATION='12OKSWT  177' THEN GRADIENT='UP';                           00006900
IF STATION='12OKSWT  178' THEN GRADIENT='UP';                           00007000
IF STATION='12OKSWT  179' THEN GRADIENT='UP';                           00007100
IF STATION='12OKSWT  180' THEN GRADIENT='UP';                           00007200
IF STATION='12OKSWT  181' THEN GRADIENT='UP';                           00007300
IF STATION='12OKSWT  182' THEN GRADIENT='UP';                           00007400
IF STATION='12OKSWT  184' THEN GRADIENT='UP';                           00007500
IF STATION='12OKSWT  186' THEN GRADIENT='UP';                           00007600
IF STATION='12OKSWT  187' THEN GRADIENT='UP';                           00007700
IF STATION='12OKSWT  188' THEN GRADIENT='UP';                           00007800
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00007810
*         IDENTIFY DOWN-GRADIENT WELLS AS SUCH.                   ;     00007820
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00007830
IF STATION='12OKSWT  177' THEN GRADIENT='DN';                           00007900
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00007910
PROC SORT; BY GRADIENT;                                                 00008000
PROC MEANS NOPRINT;                                                     00008100
VAR ARSENIC;                                                            00008200
BY GRADIENT;                                                            00008300
OUTPUT OUT=STATS MEAN=XBAR VAR=S N=N;                                   00008400
DATA _NULL_;                                                            00008500
RETAIN X1 X2 S1 S2 N1 N2;                                               00008600
SET STATS END=EOF;                                                      00008700
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00008710
*         NOTE: DUE TO SORT ORDER BY WELL, X1==>D-G               ;     00008720
*                                          X2==>U-G               ;     00008721
*  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  -  ;     00008730
IF _N_=1 THEN DO;                                                       00008800
X1=XBAR;S1=S;N1=N;                                                      00008900
END;                                                                    00009000
ELSE IF _N_=2 THEN DO;                                                  00009100
X2=XBAR;S2=S;N2=N;                                                      00009200
END;                                                                    00009300
ELSE DO;                                                                00009400
FILE PRINT;                                                             00009500
PUT 'ONLY TWO CLASS LEVELS ALLOWED. NO T TESTS WILL BE DONE.';          00009600
STOP;END;                                                               00009700
IF EOF THEN DO;                                                         00009800
VAR1=S1/N1;VAR2=S2/N2;                                                  00009900
SUMVAR=VAR1+VAR2;                                                       00010000
TPRIME=(X1-X2)/SQRT(SUMVAR);                                            00010100
CCT05=(VAR1*TINV(.95,N1-1)+VAR2*TINV(.95,N2-1))/SUMVAR;                 00010200
FILE PRINT;                                                             00010300
PUT 'T''=' TPRIME/                                                      00010400
    '.05 CRITICAL VALUE IS ' CCT05;                                     00010500
IF TPRIME>CCT05 THEN PUT 'PARAMETER SHOWS AN INCREASE FROM UPGRADIENT TO00010600
 DOWNGRADIENT';                                                         00010700
ELSE PUT 'PARAMETER DOES NOT SHOW A SIGNIFICANT INCREASE FROM UPGRADIENT00010800
 TO DOWNGRADIENT';                                                      00010900
END;                                                                    00011000
TITLE1 83 LANDFILL;                                                     00011100
TITLE2 COCHRAN''S APPROXIMATION  BEHRENS-FISHER STUDENT''S T TEST;      00011200
TITLE3 ONE TAIL TEST;                                                   00011300
TITLE4  ARSENIC;                                                        00011400
STOPSAS,                                                                00011500
./III      JOB (AAAASTORP,MIII),STORET,TIME=(,10),                      00011600
./             MSGLEVEL=(1,1),PRTY=4                                    00011700
**ROUTE  PRINT HOLD                                                     00011800
**JOBPARM LINES=10                                                      00011900
                                                                        00012000
        ************************                                        00012100
        *    LANDFILL.MEANS    *                                        00012200
        ************************                                        00012300
                                                                        00012400
LANDFILL.MEANS COMPUTES MEANS AND OTHER SIMPLE STATISTICS ON EACH       00012500
OF THE WELLS AT THE SITE.                                               00012600
SAS IS USED TO TREAT REMARK CODE "K" ACCORDING TO THE AUTHOR'S          00012700
PREFERENCE.                                                             00012800
                                                                        00013000
PGM=RET,MORE=SAS,PRT=NO,                                                00013100
A=12OKSWT,                                                              00013200
S=177,S=178,S=179,S=180,S=181,S=182,S=183,S=184,S=185,S=186,S=187,S=188,00013300
P=1000,P=1005,P=1049,P=71890,P=82295,P=32732,P=950,P=930,               00013400
P=665,P=400,P=70304,P=95,P=680,P=81375,P=39370,P=39365,P=34301,         00013500
NOECHO,                                                                 00013600
SASPARMS=BEGIN,                                                         00013700
DATA LANDFILL;                                                          00013800
*                                                                       00013900
INCLUDE (FCFREAD)                                                       00014000
DO OVER P;                                                              00014100
IF R='K' THEN P=P/2;                                                    00014200
END;                                                                    00014300
DO I=1 TO 6;                                                            00014310
IF P^=. THEN P=P*.001; * UNITS CHANGE ON FIRST SIX PARAMETERS ;         00014400
END;                                                                    00014500
RENAME P1=ARSENIC P2=BARIUM P3=LEAD P4=MERCURY P5=CHLORIDE P6=PHENOL    00014700
       P7=FLUORIDE P8=SODIUM P9=PHOS P10=PH P11=TDS P12=CONDUCT         00014800
       P13=TOC P14=TOX P15=DDT P16=DDD P17=CHLORBNZ;                    00014900
LABEL P1='MG/L' P2='MG/L' P3='MG/L' P4='MG/L' P5='MG/L' P6='MG/L'       00015000
      P7='MG/L' P8='MG/L' P9='MG/L' P11='MG/L' P12='UMHOS/CM'           00015100
      P13='MG/L' P14='MG/L' P15='UG/L' P16='UG/L' P17='UG/L';           00015200
IF STATION='12OKSWT  177' THEN WELL='177';                              00015300
IF STATION='12OKSWT  178' THEN WELL='178';                              00015400
IF STATION='12OKSWT  179' THEN WELL='179';                              00015500
IF STATION='12OKSWT  180' THEN WELL='180';                              00015600
IF STATION='12OKSWT  181' THEN WELL='181';                              00015700
IF STATION='12OKSWT  182' THEN WELL='182';                              00015800
IF STATION='12OKSWT  183' THEN WELL='183';                              00015900
IF STATION='12OKSWT  184' THEN WELL='184';                              00016000
IF STATION='12OKSWT  185' THEN WELL='185';                              00016100
IF STATION='12OKSWT  186' THEN WELL='186';                              00016200
IF STATION='12OKSWT  187' THEN WELL='187';                              00016300
IF STATION='12OKSWT  188' THEN WELL='188';                              00016400
PROC MEANS MAXDEC=3 N MEAN CV STD MIN MAX;                              00016500
VAR ARSENIC BARIUM LEAD MERCURY CHLORIDE PHENOL FLUORIDE SODIUM PHOS    00016600
    PH TDS CONDUCT TOC TOX DDT DDD CHLORBNZ;                            00016700
BY WELL;                                                                00016800
TITLE1 ***********************;                                         00016900
TITLE2 *     83 LANDFILL     *;                                         00017000
TITLE3 * STATISTICAL SUMMARY *;                                         00017100
TITLE4 ***********************;                                         00017200
STOPSAS,                                                                00017300
./III      JOB (A014STORP,MIII),STORET,TIME=(,10),                      00017400
./             MSGLEVEL=(1,1),PRTY=4                                    00017500
**ROUTE  PRINT HOLD                                                     00017600
**JOBPARM LINES=10                                                      00017700
                                                                        00017800
        ************************                                        00017900
        *    LANDFILL.BKGRND   *                                        00018000
        ************************                                        00018100
                                                                        00018200
LANDFILL.BKGRND PRODUCES GROSS SITE STATISTICS SIMILAR TO THOSE         00018300
IS THE STORET MEAN PROGRAM.                                             00018400
                                                                        00018500
PGM=RET,MORE=SAS,PRT=NO,                                                00018600
A=12OKSWT,                                                              00018700
S=177,S=178,S=179,S=180,S=181,S=182,S=183,S=184,S=185,S=186,S=187,S=188,00018800
P=1000,P=1005,P=1049,P=71890,P=82295,P=32732,P=950,P=930,               00018900
P=665,P=400,P=70304,P=95,P=680,P=81375,P=39370,P=39365,P=34301,         00019000
NOECHO,                                                                 00019100
SASPARMS=BEGIN,                                                         00019200
DATA LANDFILL;                                                          00019300
*                                                                       00019400
INCLUDE AXOCOTU.SAS.CONVERT                                             00019500
DO OVER R;                                                              00019600
IF R='K' THEN P=P/2;                                                    00019700
END;                                                                    00019800
DO I=1 TO 6;                                                            00019900
P=P/1000;                                                               00020000
END;                                                                    00020100
RENAME P1=ARSENIC P2=BARIUM P3=LEAD P4=MERCURY P5=CHLORIDE P6=PHENOL    00020200
       P7=FLUORIDE P8=SODIUM P9=PHOS P10=PH P11=TDS P12=CONDUCT         00020300
       P13=TOC P14=TOX P15=DDT P16=DDD P17=CHLORBNZ;                    00020400
LABEL P1='MG/L' P2='MG/L' P3='MG/L' P4='MG/L' P5='MG/L' P6='MG/L'       00020500
      P7='MG/L' P8='MG/L' P9='MG/L' P11='MG/L' P12='UMHOS/CM'           00020600
      P13='MG/L' P14='MG/L' P15='UG/L' P16='UG/L' P17='UG/L';           00020700
IF STATION='12OKSWT  177' THEN WELL='BKGRND';                           00020800
IF STATION='12OKSWT  178' THEN WELL='BKGRND';                           00020900
IF STATION='12OKSWT  179' THEN WELL='BKGRND';                           00021000
IF STATION='12OKSWT  180' THEN WELL='BKGRND';                           00021100
IF STATION='12OKSWT  181' THEN WELL='BKGRND';                           00021200
IF STATION='12OKSWT  182' THEN WELL='BKGRND';                           00021300
IF STATION='12OKSWT  183' THEN WELL='BKGRND';                           00021400
IF STATION='12OKSWT  184' THEN WELL='BKGRND';                           00021500
IF STATION='12OKSWT  185' THEN WELL='BKGRND';                           00021600
IF STATION='12OKSWT  186' THEN WELL='BKGRND';                           00021700
IF STATION='12OKSWT  187' THEN WELL='BKGRND';                           00021800
IF STATION='12OKSWT  188' THEN WELL='BKGRND';                           00021900
PROC MEANS MAXDEC=3 N MEAN CV STD MIN MAX;                              00022000
VAR ARSENIC BARIUM LEAD MERCURY CHLORIDE PHENOL FLUORIDE SODIUM PHOS    00022100
    PH TDS CONDUCT TOC TOX DDT DDD CHLORBNZ;                            00022200
BY WELL;                                                                00022300
TITLE1 ***********************;                                         00022400
TITLE2 *     83 LANDFILL     *;                                         00022500
TITLE3 * STATISTICAL SUMMARY *;                                         00022600
TITLE4 ***********************;                                         00022700
STOPSAS,                                                                00022800
./III      JOB (AAAASTORP,MIII),STORET,TIME=(,10),                      00022900
./             MSGLEVEL=(1,1),PRTY=4                                    00023000
**ROUTE  PRINT HOLD                                                     00023100
**JOBPARM LINES=10                                                      00023200
