System/3 simulator V3.11-0 5703-XM1 COPYRIGHT HJS CORP. 1970 ENTER CONFIGURE COMMAND OR PRESS PROG START KEY ENTER DATE - MM/DD/YY >> 01/06/22 READY >> LOGON HENK/VOLF1 READY >> LISTCAT RRRR TEST DEMO2 TT DEMO3 READY >> EDIT RRRR RRRR COPIED TO WORK FILE BASIC PROGRAM FILE 068 LINES, 012 DISK UNITS IN FILE. DATE LAST MODIFIED - 01/01/41 READY >> LIST 0010 REM * SIMPLE LINEAR REGRESSION AND * 0011 REM * CORRELATION ANALYSIS * 0012 REM 0020 REM - DATA LOCATION L, SUMX X, SUMY Y, 0021 REM - SUMX SQ X2, SUMY SQ Y2P SUMXY 7, 0022 REM - NUMBER IN SAMPLE N, XBAR X1, YBAR Y1, 0023 REM - SLOPE B, INTERCEPT A, STD ERROR S1, 0024 REM - CORRELATION COEF R 0025 REM 0050 DIM L(2,50) 0070 X, Y, X2, Y2, Z = 0.0 0090 PRINT 0095 PRINT 'ENTER NUMBER IN SAMPLE' 0100 INPUT N1 0105 PRINT 0106 PRINT 'ENTER X,Y DATA LOCATION' 0110 FOR N = 1 TO N1 0120 INPUT L(1,N), L(2,N) 0130 NEXT N 0200 FOR N = 1 TO N1 0201 REM - DERIVE SUMS 0210 X = X + L(1,N) 0220 Y = Y + L(2,N) 0230 X2 = X2 + L(1,N) ^ 2 0240 Y2 = Y2 + L(2,N) ^ 2 0250 Z = Z + L(1,N) * L(2,N) 0260 NEXT N 0299 REM - DERIVE CONVERGENCE POINT (DATA MEAN) 0300 X1 = X / N 0310 Y1 = Y / N 0319 REM - DERIVE SLOPE AND INTERCEPT 0320 B = (N * Z - X * Y) / (N * X2 - X ^ 2) 0330 A = Y1 - B * X1 0399 REM - DERIVE STANDARD ERROR 0400 S1 = SQR ((Y2 - A * Y - B * Z) / (N - 2.0)) 0410 X0 = SQR (N * X2 - X ^ 2) 0420 Y0 = SQR (N * Y2 - Y ^ 2) 0430 REM - DERIVE CORRELATION COEFFICIENT 0440 R = (N * Z - X * Y) / (X0 * Y0) 0490 PRINT 0491 PRINT 0500 PRINT 0501 PRINT , '- SCATTER-GRAM ANALYSIS BY LINEAR -' 0502 PRINT , ' REGRESSION AND CORRELATION METHOD ' 0508 PRINT 0509 PRINT 0510 PRINT 'X VARIABLE','Y VARIABLE' 0520 FOR N = 1 TO N1 0521 PRINT L(1,N), L(2,N) 0522 NEXT N 0529 PRINT 0530 PRINT 'CHARACTERISTICS' 0535 PRINT 0536 PRINT 'NUMBER IN SAMPLE', N 0539 PRINT 0540 PRINT 'SUMX = ' X, 'SUMY = ' Y, 'SUMXY = ' Z 0549 PRINT 0550 PRINT 'SUMX SQ = ' X2, 'SUMY SQ = ' Y2 0559 PRINT 0560 PRINT 'XMEAN = ' X1, 'YMEAN = ' Y1 0569 PRINT 0570 PRINT 'SLOPE B = ' B, 'INTERCEPT A = ' A 0571 PRINT 0575 PRINT 'STANDARD ERROR = ' S1 0579 PRINT 0580 PRINT 'R = ' R 0899 STOP 0900 END READY >> DISABLE 100 READY >> OFF READY