//JOB(ERCC07,STEPHENS,T=200,L=10000)
//OPTIONS(QUOTES,OPT)
//ALGOL
'COMMENT'
                         ***************
                         *             *
                         *   ALGG055   *
                         *             *
                         ***************
;
'BEGIN'
'COMMENT'    YC02   ;
'COMMENT'  ANALYSIS OF VARIANCE TABLE FOR MG,NA,K,P.;
'COMMENT'   A0,A6,A14;
'INTEGER' SETS,N,R,C,X,Y,I,J,Z;
'REAL'  K,GM,TB1,TB2,TB3,TB4,T1,T2,T3;
'REAL'  'ARRAY'  DATA[1:7,1:11],  AV[1:17,1:6],AGN[1:3,1:7],
             ALS[1:3,1:3], NG[1:3,1:4],GLSN[1:3,1:7],SN[1:5,1:7]
;
'COMMENT'   SET TO ZERO;
SETS:=N:=I:=J:=0;
T1:=T2:=T3:=0;
'FOR' I:= 1 'STEP' 1 'UNTIL' 17 'DO'
'BEGIN'
AV[I,1]:=0 'END';
'COMMENT' STANDARDS  AND DEGRES OF FREEDOM;
SETS:=READ;
T1:= 12.706;
T2:= 3.182;
T3:= 2.086;
AV[1,1]:=47;
AV[2,1]:=7;
AV[3,1]:=3;
AV[4,1]:=AV[5,1]:=AV[6,1]:=AV[7,1]:=AV[8,1]:=AV[11,1]:=1;
AV[12,1]:=AV[13,1]:=2;
AV[9,1]:=3;
AV[10,1]:=AV[14,1]:=AV[15,1]:=AV[16,1]:=5;
AV[17,1]:=20;
'COMMENT'   START OF LOOP BEGINING WITH SET TO ZERO;
     START:
R:=C:=X:=Y:=Z:=0;   K:=0.0;
       K:=TB1:=TB2:=TB3:=TB4:=0;
'FOR' R:=1  'STEP' 1 'UNTIL' 7  'DO'  'FOR' C:=1 'STEP' 1 'UNTIL'  11 'D
O'
             DATA[R,C]:=0;
'FOR' X:= 1 'STEP' 1 'UNTIL' 17 'DO'        'FOR' Y:=2 'STEP' 1 'UNTIL'
6 'DO'
            AV[X,Y]:=0;
'FOR'I:=1,2,3 'DO' 'FOR'J:=1,2,3 'DO'
'BEGIN'
      ALS[I,J]:=0 'END';
'FOR' I:= 1,2,3 'DO'       'FOR' J:= 1 'STEP' 1 'UNTIL' 4 'DO'
'BEGIN'
NG[I,J]:= 0 'END';
'FOR' I:=1,2,3 'DO'        'FOR'J:= 1 'STEP' 1 'UNTIL' 7 'DO'
'BEGIN'
     GLSN[I,J]:= 0 'END';
'FOR' I:= 1,2,3 'DO'   'FOR'  J:= 1 'STEP' 1 'UNTIL' 7 'DO'
'BEGIN'
      AGN[I,J]:=0 'END';
'FOR'I:=1 'STEP' 1'UNTIL' 5'DO'  'FOR'J:=1 'STEP' 1'UNTIL' 7 'DO'
'BEGIN'
     SN[I,J]:=0 'END';
'COMMENT'  WRITE TITLE AND READ IN DATA;
NEWLINES(1);SPACES(3);WRITE TEXT('('TITLE')');SPACES(5);
COPY TEXT ('('<>')');
'FOR' R:= 1 'STEP' 1 'UNTIL' 6 'DO' 'FOR' C:= 1 'STEP' 1 'UNTIL' 8 'DO'
         DATA[R,C]:=READ;
'COMMENT'  SUB%TOTAL AND SUB TABLES;
'FOR' R:= 1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
      DATA[R,9]:= DATA[R,1]+DATA[R,3]+DATA[R,5]+DATA[R,7];
      DATA[R,10]:=DATA[R,2]+DATA[R,4]+DATA[R,6]+DATA[R,8];
      DATA[R,11]:=DATA[R,9]+DATA[R,10]  'END';
'FOR' C:= 1'STEP'  1 'UNTIL' 11 'DO'     'FOR' R:= 1 'STEP' 1 'UNTIL' 6
'DO'
      DATA[7,C]:= DATA[7,C]+DATA[R,C];
      TB1:= DATA[7,1]+DATA[7,2];
      TB2:=DATA[7,3]+DATA[7,4];
      TB3:=DATA[7,5]+DATA[7,6];
      TB4:=DATA[7,7]+DATA[7,8];
'FOR' R:= 1,2,3,4,5,6 'DO'
'BEGIN'
ALS[1,1]:=ALS[1,1]+DATA[R,1]+DATA[R,5];
ALS[1,2]:=ALS[1,2]+DATA[R,2]+DATA[R,6];
ALS[2,1]:=ALS[2,1]+DATA[R,3]+DATA[R,7];
ALS[2,2]:=ALS[2,2]+DATA[R,4]+DATA[R,8]'END';
'FOR' I:= 1,2,3 'DO'
'BEGIN'
           NG[1,I]:=DATA[I,11];
      NG[2,I]:=DATA[I+3,11] 'END';
'FOR'  R := 1 'DO'
'BEGIN'
GLSN[1,1]:=GLSN[1,1]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:= 2 'DO'
GLSN[1,2]:=GLSN[1,2]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:= 3 'DO'
GLSN[1,3]:=GLSN[1,3]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:= 4 'DO'
GLSN[1,4]:=GLSN[1,4]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:= 5 'DO'
GLSN[1,5]:=GLSN[1,5]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:=6 'DO'
GLSN[1,6]:=GLSN[1,6]+DATA[R,1]+DATA[R,2]+DATA[R,5]+DATA
[R,6];
'FOR' R:= 1 'DO'
GLSN[2,1]:=GLSN[2,1]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8];
'FOR' R:= 2 'DO'
GLSN[2,2]:=GLSN[2,2]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8];
'FOR' R:= 3 'DO'
GLSN[2,3]:=GLSN[2,3]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8];
'FOR' R:= 4 'DO'
GLSN[2,4]:=GLSN[2,4]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8];
'FOR' R:=5 'DO'
GLSN[2,5]:=GLSN[2,5]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8];
'FOR' R:= 6 'DO'
GLSN[2,6]:=GLSN[2,6]+DATA[R,3]+DATA[R,4]+DATA[R,7]+DATA
[R,8] 'END';
'FOR' R:=1 'DO'
'BEGIN'
     AGN[1,1]:=AGN[1,1]+DATA[R,9];
     AGN[1,2]:=AGN[1,2]+DATA[R+1,9];
     AGN[1,3]:=AGN[1,3]+DATA[R+2,9];
  AGN[1,4]:=AGN[1 ,4]+DATA[R+3,9];
     AGN[1,5]:=AGN[1,5]+DATA[R+4,9];
      AGN[1,6]:=AGN[1,6]+DATA[R+5,9];
      AGN[2,1]:=AGN[2,1]+DATA[R,10];
      AGN[2,2]:=AGN[2,2]+DATA[R+1,10];
      AGN[2,3]:=AGN[2,3]+DATA[R+2,10];
      AGN[2,4]:=AGN[2,4]+DATA[R+3 ,10];
      AGN[2,5]:=AGN[2,5]+DATA[R+4 ,10];
      AGN[2,6]:=AGN[2,6]+DATA[R+5,10]'END';
'FOR'R:=1 'DO'
'BEGIN'
      SN[1,1]:=SN[1,1]+DATA[R,1]+DATA[R,5];
      SN[1,4]:=SN[1,4]+DATA[R,3]+DATA[R,7];
   SN[3,1]:=SN[3,1]+DATA[R,2]+DATA[R,6];
      SN[3,4]:=SN[3,4]+DATA[R,4]+DATA[R,8];
'FOR'R:= 2 'DO'
    SN[1,2]:=SN[1,2]+DATA[R,1]+DATA[R,5];
       SN[1,5]:=SN[1,5]+DATA[R,3]+DATA[R,7];
       SN[3,2]:=SN[3,2]+DATA[R,2]+DATA[R,6];
       SN[3,5]:=SN[3,5]+DATA[R,4]+DATA[R,8];
'FOR'R:= 3 'DO'
       SN[1,3]:=SN[1,3]+DATA[R,1]+DATA[R,5];
     SN[1,6]:=SN[1,6]+DATA[R,3]+DATA[R,7];
        SN[3,3]:=SN[3,3]+DATA[R,2]+DATA[R,6];
        SN[3,6]:=SN[R,6]+DATA[R,4]+DATA[R,8];
'FOR'R:= 4 'DO'
        SN[2,1]:=SN[2,1]+DATA[R,1]+DATA[R,5];
        SN[2,4]:=SN[2,4]+DATA[R,3]+DATA[R,7];
        SN[4,1]:=SN[4,1]+DATA[R,2]+DATA[R,6];
    SN[4,4]:=SN[4,4]+DATA[R,4]+DATA[R,8];
'FOR'R:= 5 'DO'
      SN[2,2]:=SN[2,2]+DATA[R,1]+DATA[R,5];
    SN[2,5]:=SN[2,5]+DATA[R,3]+DATA[R,7];
       SN[4,2]:=SN[4,2]+DATA[R,2]+DATA[R,6];
       SN[4,5]:=SN[4,5]+DATA[R,4]+DATA[R,8];
'FOR'R:= 6 'DO'
       SN[2,3]:=SN[2,3]+DATA[R,1]+DATA[R,5];
       SN[2,6]:=SN[2,6]+DATA[R,3]+DATA[R,7];
       SN[4,3]:=SN[4,3]+DATA[R,2]+DATA[R,6];
       SN[4,6]:=SN[4,6]+DATA[R,4]+DATA[R,8]'END';
'FOR'I:=1,2 'DO'  'FOR'J:=1,2 'DO'
'BEGIN'
      ALS[I,3]:=ALS[I,3]+ALS[I,J] 'END';
'FOR' I:= 1 ,2 'DO'         'FOR' J:= 1,2,3 'DO'
'BEGIN'
     NG[I,4]:=NG[I,4]+NG[I,J]'END';
'FOR'J:=1,2,3 'DO'
  'BEGIN'
      ALS[3,J]:=ALS[1,J]+ALS[2,J] 'END';
'FOR' J:= 1,2, 3,4'DO'
'BEGIN'
NG[3,J]:=NG[1,J]+NG[2,J]'END';
'FOR' I:=1,2 'DO'   'FOR' J:= 1, 2,3,4,5,6 'DO'
'BEGIN'
      AGN[I,7]:=AGN[I,7]+AGN[I,J];
'FOR' J:= 1 'STEP' 1 'UNTIL' 7 'DO'   'FOR' I:= 1,2'DO'
      AGN[3,J]:= AGN[3,J]+AGN[I,J]  'END';
'FOR' I:=1,2 'DO'            'FOR' J:=1,2,3,4,5,6  'DO'
'BEGIN'
      GLSN[I,7]:=GLSN[I,7]+GLSN[I,J]'END';
'FOR' J:= 1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
       GLSN[3,J]:=GLSN[1,J]+GLSN[2,J]'END' ;
'FOR' I:= 1,2,3,4 'DO'
'FOR' J:= 1,2,3,4,5,6 'DO'
'BEGIN'
     SN[I,7]:=SN[I,7]+SN[I,J];
 'END';
GM:=DATA[7,11]/48;
K:=DATA[7,11]**2/48;
NEWLINES(3);SPACES(3);WRITE TEXT('('GENERAL MEAN ')');
SPACES(3);
PRINT(GM,6,5);
NEWLINES(2);
'COMMENT' WRITE MEANS;
NEWLINES(3);SPACES(3);
WRITE TEXT('('SUB.BLOCK%MEANS')');
NEWLINES(3);SPACES(3);WRITE TEXT('('BL1YR.1')');SPACES(2);
WRITE TEXT('('BL1YR.2')');
SPACES(2);
WRITE TEXT('('BL1YR.3')');
SPACES(2);
WRITE TEXT('('BL2YR.1')');
SPACES(2);
WRITE TEXT('('BL2YR.2')');
SPACES(2);
WRITE TEXT('('BL2YR.3')');
SPACES(2);
WRITE TEXT('('BL3YR.1')');
SPACES(2);
WRITE TEXT('('BL3YR.2')');
SPACES(2);
WRITE TEXT('('BL3YR.3')');
SPACES(2);
WRITE TEXT('('BL4YR.1')');
SPACES(2);
WRITE TEXT('('BL4YR.2')');
SPACES(2);
WRITE TEXT('('BL4YR.3')');
SPACES(2);
'FOR'C:=1'STEP'1'UNTIL' 8 'DO'
PRINT(DATA[R,C]/6,4,3);
NEWLINES(3);SPACES(3);
WRITE TEXT('('BLOCK %MEANS')');
NEWLINES(4);SPACES(8);
WRITE TEXT('('BL1')');
SPACES(13);
WRITE TEXT('('BL2')');
SPACES(13);
WRITE TEXT('('BL3')');
SPACES(13);
WRITE TEXT('('BL4')');
NEWLINES(2);
SPACES(3);
PRINT(TB1/12,6,5);
SPACES(3);
PRINT(TB2/12,6,5);
SPACES(3);
PRINT(TB3/12,6,5);
SPACES(3);
PRINT(TB4/12,6,5);
NEWLINES(4);SPACES(3);
WRITE TEXT ('(' MEAN%OF%G%GC,N* TREAMENTS,AND%GVGC%
TIMES%N*TREATMENTS*AT%MEAN%YEAR.')');
NEWLINES(3);SPACES(38);
WRITE TEXT('('NITROGEN')');
SPACES(2);
WRITE TEXT('('RATE')');
SPACES(16);
WRITE TEXT('('MEAN%OF%GVGC')');
NEWLINES(1);SPACES(25);
WRITE TEXT('('N1')');
SPACES(15);
WRITE TEXT('('N2')');
SPACES(14);
WRITE TEXT('('N3')');
SPACES(6);
WRITE TEXT('('%TREATMENTS')');
NEWLINES(2);SPACES(13);
WRITE TEXT('('G')');
SPACES(2);
'FOR'J:=1,2,3 'DO'
'BEGIN'
SPACES(3);
PRINT(NG[1,J]/8,6,5);
'END';
SPACES(3);
PRINT(NG[1,4]/24,6,5);
NEWLINES(1);SPACES(1);
WRITE TEXT ('('%GRASS%TR.')');
NEWLINES(1);SPACES(13);
WRITE TEXT('('GC')');SPACES(1);
'FOR'J:=1,2,3 'DO'
'BEGIN'
SPACES(3);
PRINT(NG[2,J]/8,6,5);
'END';
SPACES(3);
PRINT(NG[2,4]/24,6,5);
NEWLINES(2);
WRITE TEXT ('('N*TREATMENT*MEAN')');
'FOR'J:=1,2,3 'DO'
'BEGIN'
SPACES(3);
PRINT(NG[3,J]/16,6,5);
'END';
SPACES(3);
PRINT(NG[3,4]/48,6,5);
NEWLINES(3);SPACES(20);
WRITE TEXT ('('MEANS*OF*LS,GC,GCN,')');
NEWLINES(4);SPACES(44);
WRITE TEXT('('GRASS')');
SPACES(2);
WRITE TEXT('('TREATMENTS')');
SPACES(45);
WRITE TEXT('('MEANS*OF*LS')');
NEWLINES(2);SPACES(25);
WRITE TEXT('('GN1')');
SPACES(10);
WRITE TEXT('('GN2')');
SPACES(10);
WRITE TEXT('('GN3')');
SPACES(10);
WRITE TEXT('('GCN1')');
SPACES(10);
WRITE TEXT('('GCN2')');
SPACES(10);
WRITE TEXT('('GCN3')');
SPACES(10);
WRITE TEXT('('TREATMENTS')');
NEWLINES(2);SPACES(1);
WRITE TEXT('('METHOD*OF')');
SPACES(3);
WRITE TEXT('('L')');
'FOR'J:=1 'STEP' 1'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(GLSN[1,J]/4,6,5);
'END';
SPACES(3);
PRINT(GLSN[1,7]/24,6,5);
NEWLINES(1);SPACES(1);
WRITE TEXT ('('N*APLICATION')');
NEWLINES(1);SPACES(13);
WRITE TEXT('('S')');
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(GLSN[2,J]/4,6,5);
'END';
SPACES(3);
PRINT(GLSN[2,7]/24,6,5);
NEWLINES(2);SPACES(2);
WRITE TEXT('('GR.*TR.MEANS')');
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(GLSN[3,J]/8,6,5);
'END';
SPACES(3);
PRINT(GLSN[3,7]/48,6,5);
NEWLINES(8);SPACES(14);
WRITE TEXT ('(' MEANS*OF*GN,GCN,AGE*OF*LEY')');
NEWLINES(7);SPACES(44);
WRITE TEXT('('GRASS')');
SPACES(2);
WRITE TEXT('('TREATMENTS')');
SPACES(45);
WRITE TEXT('('MEANS*OF*AGE')');
NEWLINES(3);SPACES(22);
WRITE TEXT('('GN1')');
SPACES(12);
WRITE TEXT('('GN2')');
SPACES(11);
WRITE TEXT('('GN3')');
SPACES(12);
WRITE TEXT('('GCN1')');
SPACES(12);
WRITE TEXT('('GCN2')');
SPACES(10);
WRITE TEXT('('GCN3')');
SPACES(10);
WRITE TEXT('('TREATMENTS')');
NEWLINES(2);SPACES(3);
WRITE TEXT('('AGE*OF')');
SPACES(2);
WRITE TEXT('(' 1 ')');
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(AGN[1,J]/4,6,5);
'END';
SPACES(3);
PRINT(AGN[1,7]/24,6,5);
NEWLINES(1);SPACES(5);
WRITE TEXT ('('LEYS')');
SPACES(2);
WRITE TEXT('(' 2')');
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(AGN[2,J]/4,6,5);
'END';
SPACES(3);
PRINT(AGN[2,7]/24,6,5);
NEWLINES(2);
WRITE TEXT('('GR.*TR.MEANS')');
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(AGN[3,J]/3,6,5);
'END';
SPACES(3);
PRINT(AGN[3,7]/48,6,5);
NEWLINES(10);SPACES(13);
WRITE TEXT ('('MEANS *OF G,GC,AGE*OF*LEY,N*TREATMENTS')');
NEWLINES(6);SPACES(35);
WRITE TEXT('(' L ')');
SPACES(35);
WRITE TEXT('(' S ')');
SPACES(33);
WRITE TEXT('('MEANS*OF*')');
NEWLINES(3);SPACES(20);
WRITE TEXT('('N1')');
SPACES(11);
WRITE TEXT('('N2')');
SPACES(12);
WRITE TEXT('('N3')');
SPACES(14);
WRITE TEXT('('N1')');
SPACES(12);
WRITE TEXT('('N2')');
SPACES(13);
WRITE TEXT('('N3')');
SPACES(10);
WRITE TEXT('('GRASS*TREATMENTS')');
NEWLINES(2);SPACES(5);
WRITE TEXT('('G')');
SPACES(3);
'FOR' J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(SN[1,J]/2,6,5);
'END';
SPACES(3);
PRINT(SN[1,7]/12,6,5);
NEWLINES(1);SPACES(2);WRITE TEXT('(' 1')');SPACES(2);
NEWLINES(1);SPACES(5);WRITE TEXT('('GC')');SPACES(2);
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(SN[2,J]/2,6,5);
'END';
SPACES(3);
PRINT(SN[2,7]/12,6,5);
NEWLINES(2);SPACES(5);WRITE TEXT('('G')');SPACES(3);
'FOR'J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(SN[3,J]/2,6,5);
'END';
SPACES(3);
PRINT(SN[3,7]/12,6,5);
NEWLINES(1);SPACES(2);WRITE TEXT('(' 2 ')');SPACES(2);
NEWLINES(1);SPACES(5);WRITE TEXT('('GC')');SPACES(2);
'FOR' J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
SPACES(3);
PRINT(SN[4,J]/2,6,5);
'END';
SPACES(3);
PRINT(SN[4,7]/12,6,5);
NEWLINES(2);
WRITE TEXT('('MEAN*OF*N')');
'FOR' J:=1 'STEP' 1 'UNTIL' 6 'DO'
'BEGIN'
PRINT(SN[5,J]/8,6,5);
SPACES(3);
'END';
SPACES(3);
PRINT(SN[5,7]/48,6,5);
NEWLINES(10);SPACES(3);
WRITE TEXT('('CORRECTION*FACTOR*K')');
SPACES(3);
PRINT(K,6,5);
NEWLINES(2);
NEWLINES(2);SPACES(3);
WRITE TEXT('('GENERAL MEAN ')');
SPACES(3);
PRINT(GM,6,5);
NEWLINES(2);
 'COMMENT' SUMSOF SQUARES;
'FOR' C:= 1 'STEP' 1 'UNTIL' 8 'DO' 'FOR' R:= 1 'STEP' 1'UNTIL' 6 'DO'
AV[1,2]:=AV[1,2]+DATA[R,C]**2-K/48;
'FOR' C:= 1 'STEP' 1'UNTIL' 8 'DO'
AV[2,2]:=AV[2,2]+DATA[7,C]**2-K/8;
      AV[3,2]:=(TB1**2+TB2**2+TB3**2+TB4**2)/12-K;
      AV[4,2]:=((TB1+TB3)**2 + (TB2+TB4)**2)/24-K;
      AV[5,2]:=((TB1+TB2)**2 + (TB3+TB4)**2)/24-K;
       AV[6,2]:=AV[3,2]-AV[4,2]-AV[5,2];
'FOR' R:= 1 'STEP' 1 'UNTIL' 6 'DO'
      AV[7,2]:=(DATA[7,9]**2 + DATA[7,10]**2)/24-K;
      AV[8,2]:=AV[8,2]+(ALS[1,1]**2+ALS[1,2]**2+ALS[2,1]**2
      + ALS[2,2]**2)/12-K-AV[4,2]-AV[7,2];
        AV[9,2]:=AV[2,2]-AV[3,2]-AV[7,2]-AV[8,2];
'FOR' J:= 1'STEP' 1 'UNTIL' 6 'DO'
AV[10,2]:=AV[10,2]+GLSN[3,J]**2/8-K/6;
  'FOR' I:= 1,2 'DO'
     AV[11,2]:=AV[11,2] + NG[I,4]**2/24-K/2;
'FOR' J:=1,2,3 'DO'
     AV[12,2]:=AV[12,2] + NG[3,J]**2/16-K/3;
AV[13,2]:=AV[10,2]-AV[11,2]-AV[12,2];
'FOR' I:=1,2 'DO'  'FOR' J:=1,2,3,4,5,6 'DO'
     AV[14,2]:=AV[14,2] + AGN[I,J]**2/4-K/12 -
     AV[10,2]/12-AV[7,2]/12;
'FOR' I:=1,2 'DO'  'FOR' J:=1 'STEP' 1 'UNTIL' 6 'DO'
     AV[15,2]:=AV[15,2] + GLSN[I,J]**2/4-K/12 -
     AV[10,2]/12-AV[4,2]/12;
'FOR'I:=1,2,3,4 'DO'   'FOR'J:=1,2,3,4,5,6 'DO'
     AV[16,2]:=AV[16,2] + SN[I,J]**2/2-K/24-AV[10,2]/24
     -AV[7,2]/24 -
AV[4,2]/24-AV[14,2]/24-AV[15,2]/24-AV[8,2]/24;
AV[17,2]:=AV[1,2]-AV[2,2]-AV[10,2]-AV[14,2]-AV[15,2]-AV[16,
2];
'COMMENT' MEAN SQUARES;
'FOR' I:=2 'STEP' 1 'UNTIL' 17 'DO'
AV[I,3]:=AV[I,2]/AV[I,1];
'COMMENT' VARIANCE RATIOS;
'FOR' I:=2,3,4,5 'DO'
AV[I,4]:=AV[I,3]/AV[6,3];
'FOR' I:=7,8 'DO'
AV[I,4]:=AV[I,3]/AV[9,3];
'FOR' I:=10 'STEP' 1 'UNTIL' 16 'DO'
AV[I,4]:=AV[I,3]/AV[17,3];
'COMMENT' LEVEL OF SIGNIFICANCE;
'IF' AV[2,4] >=  5.981          'THEN' AV[2,5]:=1.0 'ELSE'
'IF' AV[2,4] >= 0.39          'THEN' AV[2,5]:=5.0 'ELSE'
AV[2,5]:=-1.0;
'IF' AV[3,4]>=     5.403      'THEN' AV[3,5]:=1.0 'ELSE'
'IF' AV[3,4]  >=   216        'THEN' AV[3,5]:=5.0 'ELSE'
AV[3,5]:=-1.0;
'IF' AV[4,4]  >=   4.052       'THEN' AV[4,5]:=1.0 'ELSE'
'IF' AV[4,4]  >=  0.161           'THEN' AV[4,5]:=5.0    'ELSE'
AV[4,5]:=-1.0;
'IF' AV[5,4]  >=  4.052          'THEN' AV[5,5]:=1.0'ELSE'
'IF' AV[5,4]  >=  0.161           'THEN' AV[5,5]:=5.0    'ELSE'
AV[5,5]:=-1.0;
'FOR' I:=7,8 'DO'
'BEGIN'
'IF' AV[I,4]  >= 0.81         'THEN' AV[I,5]:=1.0    'ELSE'
'IF' AV[I,4]  >= .55             'THEN' AV[I,5]:=5.0    'ELSE'
AV[I,5]:=-1.0 'END';
'FOR' I:=10,14,15,16 'DO'
'BEGIN'
'IF' AV[I,4] >= 6.46         'THEN' AV[I,5]:=0.1 'ELSE'
'IF' AV[I,4] >= 4.10                'THEN' AV[I,5]:=1.0    'ELSE'
'IF' AV[I,4] >= 6.46             'THEN' AV[I,5]:=5.0 'ELSE'
AV[I,5]:=-1.0 'END';
'IF' AV[11,4] >= 14.82             'THEN' AV[11,5]:=0.1    'ELSE'
'IF' AV[11,4] >= 8.10             'THEN' AV[11,5]:=1.0    'ELSE'
'IF' AV[11,4] >= 4.35             'THEN' AV[11,5]:=5.0    'ELSE'
AV[11,5]:=-1.0;
'FOR' I:=12,13 'DO'
'BEGIN'
'IF' AV[I,4] >= 9.95            'THEN' AV[I,5]:=0.1    'ELSE'
'IF' AV[I,4] >= 5.85             'THEN' AV[I,5]:=1.0    'ELSE'
'IF' AV[I,4] >= 3.49             'THEN' AV[I,5]:=5.0    'ELSE'
AV[I,5]:=-1.0 'END';
'COMMENT' LEAST  SIGNIFICANT DIFFERENCES;
'IF' AV[2,5] >= 0.1 'THEN' AV[2,6]:=SQRT(AV[6,3]/3)*T1
'ELSE' AV[2,6]:=-1.0;
'IF' AV[3,5] >= 0.1 'THEN' AV[3,6]:=SQRT(AV[6,3]/6)*T1
'ELSE' AV[3,6]:=-1.0;
'IF' AV[4,5] >= 0.1 'THEN' AV[4,6]:=SQRT(AV[6,3]/12)*T1
'ELSE' AV[4,6]:=-1.0;
'IF' AV[5,5] >= 0.1 'THEN' AV [5,6]:=SQRT(AV[6,3]/12)*T1
'ELSE' AV[5,6]:=-1.0;
'IF' AV[7,5] >= 0.1 'THEN' AV[7,6]:=SQRT(AV[9,3]/12)*T2
'ELSE' AV[7,6]:=-1.0;
'IF' AV[8,5] >= 0.1 'THEN' AV[8,6]:=SQRT(AV[9,3]/3)*T2
'ELSE' AV[8,6]:=-1.0;
'IF' AV[10,5] >= 0.1 'THEN' AV[10,6]:=SQRT(AV[17,3]/4)*T3
'ELSE' AV[10,6]:=-1.0;
'IF' AV[11,5] >=  0.1 'THEN' AV[11,6]:=SQRT(AV[17,3]/12)*T3
'ELSE' AV[11,6]:=-1.0;
'IF' AV[12,5] >= 0.1 'THEN' AV[12,6]:=SQRT(AV[17,3]/8)*T3
'ELSE' AV[12,6]:=-1.0;
'IF' AV[13,5] >= 0.1 'THEN' AV[13,6]:= SQRT(AV[17,3]/2)*T3
'ELSE' AV[13,6]:=-1.0;
'IF' AV[14,5] >= 0.1 'THEN' AV[14,6]:=SQRT(AV[17,3])*T3
'ELSE' AV[14,6]:=-1.0;
'IF' AV[15,5] >= 0.1 'THEN' AV[15,6]:=SQRT(AV[17,3]/2)*T3
'ELSE' AV[15,6]:=-1.0;
'IF' AV[16,5] >= 0.1 'THEN' AV[16,6]:= SQRT(AV[17,3]*4)*T3
'ELSE' AV[16,6]:=-1.0;
'COMMENT' WRITE ANALYSIS OF VARIANCE TABLE;
NEWLINES(1);SPACES(3);
WRITE TEXT ('(' ANALYSIS*OF*VARIANCE*TABLE.')');
NEWLINES(1);SPACES(3);
WRITE TEXT('('*BREAKDOWN*AS*FOLLOVS.
(1)TOTAL.(2)SUB.BLOCKS.(3)BLOCKS.(4)LVS.(5)REPLICATION.(6)ERROR.(A).
')');
NEWLINES(2);SPACES(1);
WRITE TEXT('('(7)AGE OF LEY.(8) AGE*LVS.(9)ERROR.(B).(10)GRASS')');
WRITE TEXT('('TREATMENTS')');
NEWLINES(2);SPACES(3);
WRITE TEXT('('(11)GVGC
.(12)NITROGEN.(13)INTERACTION.(14)GRASS
TR,*AGE.(15)GRASSTR,*LVS.')');
NEWLINES(2);SPACES(3);
WRITE TEXT('('(16)GRASS TR,*AGE*LVS.
(17)ERROR.(C).
')');
NEWLINES(3);SPACES(7);
WRITE TEXT('(' D*OF*F')');
SPACES(10);
WRITE TEXT('(' S*S*D')');
SPACES(12);
WRITE TEXT('('*M*S')');
SPACES(11);
WRITE TEXT('(' V*R')');
SPACES(12);
WRITE TEXT('(' L*OF
*S')');
SPACES(51);
WRITE TEXT('(' L*S*D')');
NEWLINES(2);
'FOR' I:= 1 'STEP' 1 'UNTIL' 17 'DO'
'BEGIN'
'FOR'J:= 1'STEP' 1 'UNTIL' 5'DO'
'BEGIN'
     'IF' AV[I,J]<0.0000 'THEN'
     'BEGIN'
      SPACES(8);
      WRITE TEXT ('('  N.*S')');
     'END'
'ELSE'
'BEGIN'
SPACES(3);
PRINT(AV[I,J],6,5);
'END';
'END';
'IF' AV[I,6]<0.0 'THEN'
'BEGIN'
   SPACES(8);
 WRITE TEXT ('('
N.S.')');
'END'
 'ELSE' PRINT(AV[I,6],6,5); 'END';
NEWLINES(2);
'COMMENT' WRITE COEFFICIENTS OF VARIATION;
SPACES(3);
WRITE TEXT ('(' COEFFICIENT *OF *VARIATION*FOR*E1')');
SPACES(3);
PRINT((SQRT(AV[6,3])/GM)*100,6,5);
      NEWLINES(1);
SPACES(3);
WRITE TEXT ('(' COEFFICIENT*OF*VARIATION*FOR*E2')');
SPACES(3);
PRINT((SQRT(AV[9,3])/GM)*100,6,5);
      NEWLINES(1);
SPACES(3);
WRITE TEXT ('(' COEFFICIENT*OF*VARIATION*FOR*E3')');
SPACES(3);
PRINT((SQRT(ABS(AV[17,3]))/GM)*100,6,5);
      NEWLINES(1);
'COMMENT' LOOP MESSAGE;
     N:=N+1;
'IF' N=SETS 'THEN' 'GOTO' START;
FINISH:
               'END'
//RUN
-1;  DAVE'S TITLE  <>
5.9;0.4;-1.0;12.0;14.0;0.0;3.5;1.1;
0.2;2.1;3.0;0.4;8.0;2.2;4.0;3.12;
1.1;2.2;3.3;4.4;5.5;6.6;7.7;0.9;
13.0;2.1;0.20;5.8;13.4;21.1;1.0;0.001;
13.0;2.1;0.20;5.8;13.4;21.1;1.0;0.001;
1.1;2.2;3.3;4.4;5.5;6.6;7.7;0.9;
0.2;2.1;3.0;0.4;8.0;2.2;4.0;3.12;
5.9;0.4;-1.0;12.0;14.0;0.0;3.5;1.1;
13.0;2.1;0.20;5.8;13.4;21.1;1.0;0.001;0.3;