U.S. ENVIRONMENTAL PROTECTION AGENCY
                NATIONAL EUTROPHICATION SURVEY
                          WORKING PAPER SERIES
                                       STATISTICAL MANIPULATION OF
                                      NATIONAL EUTRDPHICATION SURVEY
                                       WATER QUALITY DATA IN STORET
                                      FORMAT FOR ACQUIRING
                                   ANALYSIS CAPABILITIES OF
                                        WORKING PAPER No, 472
           PACIFIC NORTHWEST ENVIRONMENTAL RESEARCH LABORATORY
                          An Associate Laboratory of the
              NATIONAL ENVIRONMENTAL RESEARCH CENTER - CORVALLIS, OREGON
                                    and
        NATIONAL ENVIRONMENTAL RESEARCH CENTER - LAS VEGAS, NEVADA
•&GPO	697-O32

-------
                STATISTICAL MANIPULATION OF
               NATIONAL EUTRDPHICATION SURVEY
                WATER QUALITY DATA  IN STORET
               FORMAT FOR ACQUIRING RAPID.DATA
            ANALYSIS CAPABILITIES OF STORET DATA
                  WORKING PAPER No, 472
NATIONAL BJTROPHICATION SURVEY
         OCTOBER 1975

-------
          STATISTICAL MANIPULATION OF

         NATIONAL EUTROPHICATION SURVEY

          WATER QUALITY DATA IN STORET
        Format for Acquiring Rapid Data
      Analysis Capabilities of STORET Data
                       by
                 James D.  Bliss
         Water and Land Quality Branch
         Monitoring Operations Division

     Michael  J. Friedland  and Julia Hodson
              Data Services Branch
    Office of Program Management and Support

Environmental Monitoring and Support Laboratory
               Las Vegas,  Nevada
             Working Paper No. 472
       OFFICE OF RESEARCH AND DEVELOPMENT
      U.S. ENVIRONMENTAL PROTECTION AGENCY
                  October 1975

-------
     Effective June 29, 1975, the  National Environmental Research
Center-Las Vegas  was designated the Environmental Monitoring  and
Support  Laboratory-Las Vegas,  and  the  National   Environmental
Research Center-Corvallis was designated the Con/all is Environmen-
tal Research Laboratory.

     The mention  of trade  names or commercial products  does not
constitute U.S. Environmental  Protection  Agency  endorsement  or
recommendation for use.

-------
                        n





                     CONTENTS





                                             Page



FOREWORD	  iii



CONCLUSIONS 	    1



RECOMMENDATIONS 	    1



INTRODUCTION 	    2



CONCEPT OF DATA ANALYSIS 	    4



"STORET.TO.WYLBUR" 	    6



A WORD ABOUT SAS 	    9



REFERENCES AND BIBLIOGRAPHY 	   10



APPENDIX 	   11

-------
                                  m
                               FOREWORD
    The National  Eutrophication  Survey  was   initiated  in   1972   in
response to an Administration commitment to investigate the nationwide
threat   of   accelerated   eutrophication  to  freshwater lakes   and
reservoirs.

    The Survey was designed to  develop,  in  conjunction  with  State
environmental    agencies,    information   on    nutrient   sources,
concentrations, and impact on selected freshwater lakes as a  basis  for
formulating comprehensive  and  coordinated national,  regional,   and
State   management   practices  relating  to  point  source  discharge
reduction and nonpoint source pollution in lake watersheds.

    The mathematical   and  statistical  procedures  selected   for   the
Survey's eutrophication analyses are based on  related concepts that:

    (1)  A  generalized  representation  or model  relating   sources,
         concentrations and impacts can be constructed.

    (2)  By applying  measurements of  relevant  parameters associated
         with   lake    degradation,   the  generalized  model   can   be
         transformed  into an operational representation of a  lake,  its
         drainage basin, and related nutrients.

    (3)  With such a  transformation, an assessment  of  the  potential
         for eutrophication control can be made.

    An individual report is prepared for each  study lake.  In the  lake
report,  the  first  stage  of  evaluation  of lake and watershed  data
collected from the study lake and its drainage  basin  is  documented.
The  report  is designed  to provide State environmental agencies  with
specific   information    for    basin    planning    water    quality
criteria/standards  review,  clean lakes, and  water quality monitoring
activities  mandated   by  the  Federal  Water   Pollution  Control   Act
Amendments   of   1972.   Each  report  is  published  as  a   National
Eutrophication Survey Working Paper.

    Beyond  the  single  lake  analysis,  broader  based  correlations
between  nutrient  concentrations  (and loading) and trophic  condition
are  being  made  to   strengthen  the  rationale  and  data  base   for
refinement  of  nutrient  water  quality  criteria  for  the   Nation's
freshwater  lakes.   Likewise,  multivariate   evaluations   for   the
relationships   between   land   use,  nutrient  export,  and  trophic

-------
                                   IV
condition, by lake class or use, are being developed to assist in  the
formulation   of   planning   guidelines  and  policies  by  the  U.S.
Environmental Protection Agency and to augment plans implementation by
the States.

-------
                             CONCLUSIONS


    As  a principal participating staff of the National  Eutrophication
Survey*,  the  Environmental  Monitoring  and  Support  Laboratory-Las
Vegas1  Water  and  Land  Quality  Branch of the Monitoring Operations
Division has collected water quality data from more than 800  selected
lakes  throughout  the  contiguous United States.  To dates the Survey
has yielded over two million data points.

    Analyses of these water quality data address  problems  unique  to
lakes  and  not  amenable  to those computer programs available in the
U.S. Environmental Protection Agency's STORE! (STOrage and  RETrieval)
system.

    Consequently,  a job control language packages "STORET„TO.WYLBUR9"
was developed to convert raw data in STORET to a  standardized  format
for statistical evaluation and manipulation of Survey data.

    "STORET.TO.WYLBUR"  has  proven  to  be a very effective means for
acquiring rapid data analysis capabilities of STORET data and  may  be
implemented   with   facility   by  both  experienced  and  relatively
inexperienced computer programmers.
                           RECOMMENDATIONS
    The "STORET.TO.WYLBUR" program has satisfied the  Water  and  Land
Quality  Branch's  additional  requirements and is easily used by most
computer personnel.  This data analysis scheme may well be  ideal  for
generalized  use.   It  is  recommended  that  the program be added as
another program option internal to the Agency's STORET system and  the
supporting  documentation  be generated to assist those wishing to use
this computer program.
*The  Special  Studies  Branch9  Criteria  and Assessment  Division9  Corvallis
  Environmental  Research  Laboratory,  Corvallis9 Oregon is also a  principal
  participating  staff  involved with this  pr9Ject.  The Corvallis  laboratory
  was  primarily  responsible for the collection of municipal sewage treatment
  plant and tributary  water quality data  (USEPA 1975).

-------
                             INTRODUCTION
    The purpose of this report is to demonstrate a system for a STORET
user to easily acquire rapid and flexible data  analysis  capabilities
with  minimal  effort.  The data management scheme outlined is for the
use of organizations which haves or will have, a data base  in  STORET
and require programming options beyond those available in STORET.

    The  following  definitions of terms may be helpful to enhance the
reader's comprehension:

    STORET:  a central computer-oriented system used by the Agency  to
    categorize, store, and retrieve water quality data.

    PGM=RET:   a  STORET  program  command  calling  for  a variety of
    retrieval programs internal to STORET.

    MORE File:   a  disk  file  created  by  PGM=RET  retrieval  which
    contains   condensed   IBM   (International   Business   Machines)
    hexadecimal output  (STORET  User's  Manual,  Vol. 2S  Chapter 10,
    Section M, 1973).

    WYLBUR:  a text "editor" and remote job entry facility.

    EDIT:   a  format  compatible  with the WYLBUR text-editing system
    available at OSI  (Optimum Systems Incorporated).

    SAS  (Statistical  Analysis   System)   and   BIOMED   (Biomedical
    Statistical Computer Programs):  two software packages which allow
    convenient implementation of various statistical techniques.
    The  procedures involved in implementing this scheme were designed
to enable users with varying levels of systems expertise to use  them.
However, it is expected that a user of this system have some knowledge

-------
of the mechanics of conducting  STORE!  retrievals, particularly PGM=RE!
and  associated  options   and the more common WYLBUR editing commands.
The ability to tailor the raw data  retrievals from  STORE!  will  help
appreciably in preparing  data for further analysis.
HISTORICAL BACKGROUND

    In  early  1974,   the only computer operation within the Water and
Land Quality Branch consisted of STORE!  (lake) data retrievals  via  a
remote  terminal  (WYLBUR).  Because  the data  processing and analytical
requirements were not satisfied by the program options internal to the
STORE! system9 a  serious need to convert  the  raw  data  in  STORE!  to
some other accessible form was recognized.
    !he  Branch  alerted  the  Laboratory's  data  services staff
Services Branch, Office of Program Management   and Support)  to  this
problem  and  requested  their assistance.   As  a results a job control
language (JCL) package, "STORE!.TO.WYLBUR,"  was written.  This program
contains the appropriate IBM JCL  and FORTRAN (FORmula  and TRANslatlon,
a computer language) programming  (Friedland  1974)   for conversion  of
the  PGM=RET  output,  or  raw data  retrieval  of STORE!,  into  an
acceptable input form for use with so-called packaged   programs (e.g.,
SAS,  BIOMED) or user-created programs.   After  a short training course
on SAS9 data evaluation was effected.    Within  about   3  months,  the
Branch  personnel could handle most analytical  requirements with ease.
Using "STORE!.TO.WYLBUR" in conjunction  with SAS has proven   to  be   a
powerful and highly flexible tool for data analysis.
PROGRAM SUMMARY

    Initially,  STORE! retrieval  data are transferred  into  a MORE File
by PGM=RE!.  "STORE!.TO.WYLBUR"   accesses  this   and   converts   it  to
standard  card-images  and  then   into  EDI!  format.    Data  are then
accessible to the WYLBUR user for  input  into  mathematical  packages
such  as  SASS  BIOMED, etc., and subsequently are analyzed easily  (or
merged  with  other  data  files)  in  accordance  with   the    unique
requirements  of  a  given  research  program.  Both novice and  expert
users are able to use this data  handling method  successfully.

-------
                       CONCEPT OF DATA ANALYSIS
    Generally, data analysis operations within the Branch are  divided
into  three  primary areas of involvement as shown in Figure 1,  Th©se
areas are:

    (1)  STORET - used primarily as a data storage  facility.   Except
         for transfer of raw data in the "STORET.TO.WYLBUR" operation,
         minimal use is made of the programs internal to STORET.

    (2)  "STORET.TO.WYLBUR" - converts data Into EDIT format from  the
         MORE  File  generated  by the STORET raw data retrieval.  The
         stylized  output   format   is   compatible   with   packaged
         programming or user-created programming.

    (3)  SAS - used for  its  extensive  data  massaging  capabilities
         (equivalent  to  FORTRAN)  as  well  as  numerous statistical
         procedures  which  generate   plots,   implement   regression
         analyses, etc.
    Other  operations  include merging or concatenation of other files
(Block 4, Figure 1) with the STORET data   (Block ])„   SAS  can  also
generate data cards to be used In other programs (Block 50 Figure 1).

-------
   FIGURE 1. Output flow diagram demonstrating data processing logic
 1.
   STORE!
  OUTPUT
(STANDARD PGM\  (MORE FILE)
  OPTIONS)     \       [4-
 2.
  STORET.TO.WYLBUR
 OTHER FILES -
BIOLOGICAL, ETC.
                    (EDIT)
      (MERGING/
   'CONCATENATION)

    •(EDIT)
                3.
                    SAS
                (CARDS)
5.
  FORTRAN, BIO MED,
   ETC. PROGRAMS
  OUTPUT
        1
DATA MANIPULATION
  SORTING
  PRINTING
  STATISTICAL
       ANALYSIS

-------
                           "STORET.TO.WYLBUR"
 USING "STORET.TO. WYLBUR"
     The   necessary  changes   for  setting up a job are located 1n three
 areas  of the  JCL  package.  On Line  1  in Figure 2 is the job card which
 must be  adjusted  to  comply with the initials and account of each  user
 (see  page  2-1,   OSI   User's   Manual  1973).   The second segment to be
 changed  is  Lines  4-8,  the STORET  retrieval coding for  the  data  that
 are  to   be  converted.   If  "PRT=NO"  (a STORET program command meaning
 "do  not  print")  is not included,  STORET raw data printout will also be
 generated.  On Line  16,  the  number  of parameters found in  the  STORET
 coding is given.   Line 17 identifies  the disk file where the converted
 data  are  written  and  the  file  name.  Disk and data set name must be
 adjusted to fit .Individual requirements as shown below.
    17 //EXEC EDSCARDS, NAME='CNAAA.BBB.CCC',DISK=TSOXXX,TYPE=OLD

        . ACCOUNT NO. - AAA

         INITIAL - BBB

       ,  DATA SET NAME - CCC

         VOLUME USED - XXX
Since the file type is old, a dummy file with the file  name  assigned
by the user should be saved on the disk pack Indicated.  A copy of the
above  program  entitled  "&CNA805.RWT.STORET.TO.WYBLUR"  is currently
available to interested persons on TS0001 (one  of  16*  on-line  disk
packs available through OSI).
*Number of on-line disk packs subject to change.

-------
               FIGURE 2.  Listing of "STORET.TO.WYLBUR"
 1.     //IIIWY JOB (AAAA,BIN,01,01,01),'  STR-TO-WYL',MSGLEVEL=(1,1)
 2.     // EXEC WQDIST
 3.   .  7/DIST.CARDFD DD *
 4.     PGM RET,PURP=104/EPA,A 11EPALES,MORE=4, RWT 35,B P,
 5.     S=380000,5=380499,
 6.     P=671,P=665,
 7.     1=77777702,
 8.     PRT=NO,
 9.     // EXEC FORTGLG,PARM='NOMAP'
10.     //LKED.SYSIN DD DSN=CNA805.RNS.RDSTOR,
11.   •••//•  UNIT=3330,VOL=SER=TS0004,DISP=SHR
12.     //GO.FT15F001 DD DSN=&FCF,DISP=(OLD,DELETE)
13.     //GO.FT12F001 DD DSN=&LAKFIL,DISP=(NEW,PAS£),UNIT=SYSDA,
14;     //  SPACE=(TRK,(25,10)),DCB=(RECFM=FB,LRECL=80,BLKSIZE=1680)
15.    -7/GO.SYSIN DD *
16.     02  NUMBER OF PARAMETERS RETRIEVED IN 1-2 (1-50)
17.     // EXEC EDSCARDS,NAME*'CNA805.RWT.DATA.SET1,DISK=TS0001,TYPE=OLD
18;     //CARDS DD DSN-&LAKFIL,DISP*OLD
 "STORET.TO.WYLBUR" OUTPUT

     The  output of "STORET.TO.WYLBUR" written onto a disk pack follows
 the  format  shown  in  Figure 3.   Because  only  a  maximum  of  six
 parameters  can be written per line, a sequence number is given at the
 end of each line indicating whether it is  the  first,  second,  etc.,
 line  of information for a given station.  Each parameter result field
 is 10-columns wide.  At the end of each  line  is  given  the  station
 STORET  number  in Columns 61-66; year, 67-68; month, 69-70; date, 71-
 72; time of day, 73-76; depth, 77-79; and sequence in  Column 80.   As
 missing  values  are  represented  by  "-123.000,"  it is necessary to
 replace these numbers with a string of eight blanks using  the  WYLBUR
 "CHANGE"   command   before  making  use  of  these  data  sets.   The
 programming allows up to 50 parameters to be converted to EDIT  format
 in any one given run.

-------
                                  8
 FIGURE 3.   The data set created by "STORET.TO.WYLBUR" on TS0002
                                                          c              u
                                                          o              c
                                                         •r-            JC 0)
                                                         •«->     0)
                                                         +*.    (O   -i-   0) 
-------
ADVANTAGES
     (1)  Anyone who can set up a STORE! raw data retrieval  will  have
         minimal problems in converting it into EDIT format on the OSI
         computer  system.

     (2)  The  data  output  of  "STORET.TO.WYLBUR"  has  a  formalized
         structure  which provides a data file on which other programs
         (;user-written, BIOMED, SAS, etc.) can operate.

     (3)  Knowledge of IBM  JCL  is  minimal.   The  "STORET.TO.WYLBUR"
         program  is  ideal  for users who have neither time nor long-
         term requirements to learn the intricacies of IBM JCL.
DISADVANTAGES

     (1)  Parameter order must be clearly determined on the  data  file
         created, or gross errors  may occur.

     (2)  Large   "STORET.TO.WYLBUR"  data  files   are   expensive   to
         manipulate using the WYLBUR editor.
                           A WORD ABOUT SAS
    Since  SAS  is  the principal software package used by the Branch,
its capabilities can be clearly  stated.   First  of  all,  SAS  is  a
unified system of data manipulation, editing, and statistical analysis
(Service 1972).  Programming and assignment statements are designed to
enable the user to effectively modify data prior to final analysis and
display.   Statistical  analysis, plotting, and sorting can usually be
implemented  using  only  several  lines  of  coding  associated  with
"PROCEDURE"  statements.   Subsetting,  concatenation,  and  data  set
merging have proven  to  be  invaluable  tools  in  data  manipulation
available  in  SAS.   In  fact, SAS has been used extensively to merge
water  quality  information  originating  in  STORET  with  biological
information  to  create  needed  data  files.  SAS may also be used to
process information with the subsequent output passed from the  system
into  card  format  for  input  into  other software packages or user-
written programming.  A sample SAS program built around  a  data  file
created  by "STORET.TO.WYLBUR" is shown in the Appendix. The output of
the program is also included.

-------
                                  10


                              REFERENCES
Friedland, M. J.  1974.  A Method for Improving User Access to STORET.
    Proceedings No. 1, Office of Research and Development ADP Workshop,
    Bethany College, October 2-4.  pp. 37-45.

Service, Jolayne.  1972.  A .User's Guide to the Statistical Analysis
    System.  (Based on Statistical Analysis System Manual prepared by
    A- J. Barr and J. H. Goodnight).  Student Supply Stores, North
    Carolina State University, Raleigh, North Carolina.

U.S. Environmental Protection Agency.  1974.  STORET Handbook, Vol. 2,
    Chapter 10.

U.S. Environmental Protection Agency.  1975.  National Eutrophication
    Survey Methods 1973-1976.  Working Paper No. 175.  Environmental
    Monitoring and Support Laboratory, Las Vegas, Nevada, and
    Corvallis Environmental Research Laboratory, Corvallis, Oregon.
                             BIBLIOGRAPHY

Optimum Systems Incorporated.  1973.  User's Manual.  Bethesda,
    Maryland.

Perkins, Carroll G.  1974.  A Guide to the Supplementary Procedure
    Library for the Statistical Analysis System.  North Carolina
    State University, Raleigh, North Carolina.

U.S. Environmental Protection Agency.  1974.  Handbook - Water Quality
    Control Information System - STORET.  Volumes 1, 2, and 3.
    Washington, D. C.

           1975.  User's Manual for EPA Scientific Applications
    Software.  Washington, D. C.

-------
                                11

                            APPENDIX
The  program  (&CNA^5.Rfcn\ST^                             .—


     //RWTWY  JOB  (A805.BIN,01,01,01),' STR-TO-WYL1,MS6LEVEL«(1,1)
     // EXEC  WQDIST
     //DIST.CARDFD DD *
     PGM RET,PURP=104/EPA,A 11EPALES,MORE=4,RWT 35.B P,
     5=380000,5=380499,
     P=671,P=665,
     1=77777702,
     PRT=NO,
     // EXEC  FORTGLG,PARM='NOMAP'
     //LKED.SYSIN DD DSN=CNA805.RNS.RDSTOR,
     //   UNIT=3330,VOL=SER=TS0004,DISP=SHR
     //GO.FT15F001. DD DSN=&FCF,DISP=(OLD,DELETE)
     //GO.FT12F001 DD DSN=&LAKFIL,DISP=(NEW,PASS),UNIT=SYSDA,
     //  SPACE=(TRK,(25,10)),DCB=(RECFM=FB,LRECL=80,BLKSIZE=1680)
     //GO.SYSIN DD *
     02  NUMBER OF PARAMETERS RETRIEVED  IN 1-2  (1-50)
     // EXEC  EDSCARDS,NAME='CNA805.RHT.DATA.SET',DISK»TS0002.TYPE«OLD
     //CARDS  DD DSN»&LAKFIL,DISP«OLD

Part of data set written  (note that the data set Is a standard UYLBUR
File and Is  accessed Into working memory by the "USE" command):

Data Set Name:  DATA.SET                 Location:  TS0002


0.170     0.260  -  0.000     0.000     0.000     0.00038010174 430 920  01
0.161     0.245     0.000     0.000     0.000     0.00038010174 430 920  51
0.160     0.243     0.000     0.000     0.000     0.00038010174 430 920 151
0.155     0.239     0.000     0.000     0.000     0.00038010174 430 920 231
0.279     0.283     0.000     0.000     0.000     0.00038010174 7171120  01
0.277     0.292     0.000     0.000     0.000     0.00038010174 7171120  51
0.307     0.326     0.000     0.000     0.000     0.00038010174 7171120 151
0.274     0.318     0.000     0.000     0.000     0.00038010174 7171120 201
0.220     0.275     0.000     0.000     0.000     0.00038010174 9171410  01
0.204     0.245     0.000     0.000     0.000     0.00038010174 9171410 151
0.193     0.245     0.000     0.000     0.000     0.00038010174 9171410 271
0.147     0.261     0.000     0.000     0.000     0.00038010274 430 945  01
0.138     0.247     0.000     0.000     0.000     0.00038010274 430 945  81
0.273     0.377     0.000     0.000     0.000     0.00038010274 7171030  01
0.280     0.314     0.000     0.000     0.000     0.00038010274 7171030  51
0.281     0.319     0.000     0.000     0.000     0.00038010274 7171030 171
0.088     0.224     0.000     0.000     0.000     0.00038010274 9171455  01


                                  ETC.

-------
                                     12
    The SAS program below 1s built around the preceding data set.   Check the
    external copy command 1n the OSI Manual for rapid combining methods of
    pre-existing programming and data files:
//RWT JOB  (A805,45,1,1),'SAS RUN1,MSGLEVEL«(0,0)
// EXEC SAS
//SAS.SYSIN DD *
TITLE 'ORTHOPHOSPHATE VS TOTAL PHOSPHORUS FOR SELECTED
DATA RAW;
INPUT CHLA 1-10 PT 11-20;
CARDS;
                                                  LAKES'
     0.170
     0.161
     0.160
     0.155
     0.279
     0.277
     0.307
     0.
     0.
     0.
  ,193
  ,147
  ,138
0.273
0.280
0.281
0.088
          0.260
          0.245
          0.243
          0.239
          0,283
          0.292
          0.326
0.245
0.261
0.247
0.377
0.314
0.319
0.224
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
          0.000     0.000     0.000     0.00038010174
0,
0,
0,
0,
0.
0.
0.
000
000
000
000
000
000
000
                                 ETC
PROCEDURE PRINT;
PROCEDURE PLOT;
PROCEDURE CORR;
PROCEDURE REGR; MODEL PROTHO=PT;
/*
                                                     01
                                                     51
                                          430 920
                                          430 920
                                          430 920 151
                                          430 920 231
                                          7171120
                                          7171120
                                                     01
                                                     51
0.000     0.000     0.00038010174
0.000     0.000     0.00038010274
0.000     0.000     0.00038010274
0.000     0.000     0.00038010274
0.000     0.000     0.00038010274
0.000     0.000     0.00038010274
0.000     0.000     0.00038010274
                                            7171120 151
9171410 271
430 945  01
430 945  81
7171030
7171030
7171030
9171455
                                                                              01
                                                                              51
                                                                             171
                                                                              01
    Following is the output generated by the preceding program.   The results
    can be fetched using "FETCH XXX DDN=FT03F001" where XXX is the job number.
    This output demonstrates the capabilities of the SAS "PROCEDURE" statements
    and was a modified program example (Service  1972):

-------
ORTHOPHOSPHATE  vs  TOTAL  PHOSPHORUS FOR  SELECTED LAKES
19t 1« THURSDAYt OCTOBER 9t 1975
                     TITLE  »0«THOPHOSPHATE VS  TOTAL PHOSPHOHUS FOR SELECTED LAKES'I
                     DATA RAWI
                     INPUT  PORTHO  1-10 PT  11-20«
                     CARDS
                      106  OBSERVATIONS  IN  DATA  SET  RAto            2   VARIABLES


                     PROCEDURE PRINT!


                                            ORThOPHOSPHATE VS TOTAL  PHOSPHORUS FOR SELECTED LAKES
                                                       OHS
                                                              POHTHO
                                                                          PT
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
16
19
20
21
22
23
24
26
26
27
26
29
30
31
32
33
34
35
36
37
38
39
.170
.161
.160
.155
.279
.277
.307
.274
.220
.204
.193
.147
.136
.273
.280
.281
.088
.097
.119
.117
.119
.304
.317
.299
.OSA
.066
.074
.117
.273
.273
.286
.278
.117
.117
.112
.023
.Olb
.016
.013
.260
.245
.243
.239
.283
.292
.326
.318
.275
.245
.245
.261
.247
.377
.314
.319
.224
.216
.237
.237
.236
.477
.412
.395
.202
.196
.197
.247
.303
.323
.298
.306
.212
.213
.286
.130
.131
.129
.059
                                                                                                                                     co
                                                                   ETC.

-------
                    PROCEDURE PLOT I
    0.40000000
    0.32000000
    0.24000000
      PORTHO
    0.16000000
    0.00000000
    0.00000000
                                           ORTMOPHOSPHATE VS TOTAL PHOSPHORUS FOR SELECTED LAKES
                                                       PLOT OF PORTHO vs PT
                                                                                          A   A
                                                                                       A A AAAAA
                        A A

                  A
                  A   A AA

                   A
                     A
                  A   A
                 AA
                      A
                   A
                                                                                                      A   A
                                                                         B    C A
                                               BAA
                                                A
                                                        AAA
                                         A AAABA        A
                                       ABE A ABO  AACB A AA  A
                                       AA    A  BA A
                    -0.02500000          0.07500000
               LEGEND I A « 1 OBS • B • 2 OBS » ETC.
OHTHOPHOSPHATE vs TOTAL PHOSPHORUS FOR SELECTED LAKES
0.17500000
           PT
0.27300000          0.37500000          0.47500000     I

                  19118 THURSDAY* OCTOBER 9* 1975

-------
                    PROCEDURE CORRI
ORTMO
T
- 106
ORTNO
T
N
106
106

PORTHO
1.000000
0.0000
0.921752
0.0001
SUM
10.99700000
19.64900000
C
PT
0.921752
0.0001
1.000000
0.0000
                                          OHTHOPHOSPHATE VS TOTAL PHOSPHORUS FOB SELECTED LAKES
                                               MEAN              MIN VALUE          MAX VALUE
                                             0.10374528         0.00600000         0.31700000
                                             0.18536798         0.04100000         0.47700000
                                          OHTHOPMOSPHATE VS TOTAL PHOSPHORUS FOR SELECTED LAKES
                                   CORRELATION COEFFICIENTS   /   PROB >  I * I UNDER HO: RHO=0
 CORRECTED SS
   1.08141012
   1.27634065
STANDARD DEV
  0.10148470
  0.11033890
 RTHOPHOSPHAU VS TOTAL PHOSPHORUS FOR SELECTED LAKES
      THURSDAY. OCTOBER 9. 1979
                    PROCEDURE RE6RI MODEL POKTHO=PTJ


PROC REGR
DATA SET
: ORTHOPHOSPHATE VS TOTAL PHOSPHORUS FOR SELECTED LAKES
: RAW NUMBER
VARIABLES : PORTHO PT
ORTHOPHOSPHATE VS TOTAL
ANALYSIS OF VARIANCE TABLE . REGRESSION COEFFICIENTS
SOURCE
REGRESSION
ERROR
CORRECTED TOTAL
SOURCE
PT
SOURCE
INTERCEPT
PT
OF
1
104
105
OF
1
B VALUES
-0.05340699
0.84778571
SUM OF SQUARES
0.91879534
0.16261479
1.08141012
SEQUENTIAL SS
0.91879534
T FOR HO:B«0
-7.08762
24.24075
OF VARIABLES » 2

PHOSPHORUS FOR SELECTED
* AND STATISTICS OF FIT
MEAN SQUARE F
0.91879534 587.
0.00156360
f VALUE PROB > F
587.61393 0.0001
PROB > IT)
0.0001
0.0001
NUMBER OF CLASSES * 0

LAKES
FOR DEPENDENT VARIABLE
VALUE PROB > F
61393 0.0001
PARTIAL SS
0.91879534
STD ERR B
0.00753525
0.03497358

PORTHO
H-SOUARE C.V.
0.84962709 38.11492 %
STD DEV PORTHO MEAN
0.03954243 0.10375
F VALUE PROB > F
587.61393 0.0001
STO B VALUES
0.0
0.9217521B
                                                                                                                                          IT.
ORTHOPHOSPHATE VS TOTAL PHOSPHORUS FOR SELECTED LAKES
19118 THURSDAY* OCTOBER 9* 1975

-------