18 



UNITED STATES NATIONAL MUSEUM BULLETIN 285 



SUMMUL(LL) = SUMMUL(LL) + DIFMUL 



IF (LGT)26,2?,26 

 26 IF (TX(I,1))20, Ik, 28 

 28 IF (TX(I,2))20, 14,13 

 2? IF (X(l,Jl))l2,l4,12 



12 IF (X(I,J2))13, 14,13 



13 SUMSQ2(l) = SUMSQ2(1) + 

 SUMSQ2( 2) = SUMSQ2( 2) 



14 CONTINUE 



DIFSQR(I,l) 



+ DIFSQR(I, 2) 



BVRMA 75 

 BVRMA 76 

 BVRMA 77 

 BVRMA 78 

 BVRMA 79 

 BVRMA 80 

 BVRMA 8l 

 BVRMA 82 

 BVRMA 83 



.5 

 .5 



STDDVY(LL) = (SUMSQ2(l) / ( PAIRS (LL) - 1.)) , 



STDDVX(LL) = (SUMSQ2(2) / (PAIRS(LL) - 1.)) * 



SY2 = STDDVY(LL) ** 2. 



SX2 = STDDVX(LL) ** 2. 



R(LL) = SUMMUL(LL) / (SUMSQ2( l) * SUMSQ2( 2)) ** 



A(LL) = STDDVY(LL) / STDDVX(LL) 



RSQ(LL) = R(LL) ** 2. 



STERRA(LL) = A(LL) * ((1. - RSQ(LL)) / PAIRS(LL)) 



B(LL) = AMEANY(LL) - AMEANX(LL) * A(LL) 



(l.-R(LL)) * 



* .5 

 A(LL) + B(LL) 



BVRMA 84 

 BVRMA 85 

 BVRMA 86 

 BVRMA 87 

 BVRMA 88 

 BVRMA 89 

 BVRMA 90 

 .5 BVRMA 91 



BVRMA 92 

 (SX2+SY2) / ((AMEANX(LL)**2.) BVRMA 93 



BVRMA 94 

 BVRMA 95 

 BVRMA 96 

 BVRMA 97 

 BVRMA 98 

 BVRMA 99 

 BVRMAIOO 

 BVRMA101 



DRMA(LL) = 100. * (2. 

 1 + (AMEANY(LL)**2.))) 

 YMAXl(LL) = SPMAX(J2) 

 IF (LGT)49,48,49 



48 TOY = SPMAX(Jl) / 49. 

 GO TO 50 



49 TOY = (SPMAX(Jl) - SLGMIN(Jl)) / 49. 



50 TOPl(LL) = (YMAXl(LL) - SPMAX(Jl)) / TOY 

 BASEl(LL) = B(LL) / TOY 



DO 17 1=1, K BVRMA102 



IF (LGT)15,16,15 BVRMA103 



15 X(I,J1) = TX(I,1) BVRMA104 

 17 X(I,J2) = TX(I 2) BVRMA105 



SPMAX(Jl) = TSPMAX BVRMA106 



SPMAX(J2) = TSPMXX BVRMA107 



16 CONTINUE BVRMA108 

 5^ NN = NBIVAR BVRMA109 



IF (M+MM)36,31,36 BVRMA110 



31 IF (89-K)45,46,46 BVRMA111 



45 WRITE OUTPUT TAPE 6,19 BVRMA112 



46 WRITE OUTPUT TAPE 6 , l8 BVRMA113 

 GO TO 39 BVRMA114 



36 IF (M)43,42,43 BVRMA115 



42 IF (70-K)44,46,46 BVRMA116 

 44 WRITE OUTPUT TAPE 6,19 BVRMA117 



43 WRITE OUTPUT TAPE 6 ,40, (LGCOL( J) , J=l , IM) BVRMA118 

 39 WRITE OUTPUT TAPE 6 , FMT8 , (R( J) , J=2 , NN) , (A( J) , J=2 , NN) , (STERRA( J) , J=BVRMA119 



12,NN),(B(J),J=2,NN), (DRMA( J) , J=2 , NN) , (AMEANY( J) , J=2 , NN) , (AMEANX( J)BVRMA120 



2,J=2,NN),(YMAX1(J),J=2,NN),(T0P1(J),J=2,NN),(BASE1(J),J=2,NN) BVRMA121 



DO 53 J=2,NN 



BVRMA122 



