Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-01 - 43,50110/stat10.sta
There are 2 other files named stat10.sta in the archive. Click here to see a list.
100'  NAME--STAT10
110'
120'  DESCRIPTION--COMPUTES THE SLOPE AND OTHER STATISTICS FOR A
130'  SIMPLE LINEAR REGRESSION WITH ONE INDEPENDENT VARIABLE
140'
150'  SOURCE--UNKNOWN
160'
170'  INSTRUCTIONS--PUT DATA IN LINE 900 AND FOLLOWING.
180'  FIRST DATA IS N, THE NUMBER OF POINTS, THEN THE DATA
190'  BY POINT(NOT BY SERIES).
200'  SAMPLE DATA ARE IN LINES 900 THROUGH 950.
210'
220'
230'  *  *  *  *  *  *   MAIN PROGRAM   *  *  *  *  *  *  *  *  *  *  *
240'
250  READ N
260  FOR I = 1 TO N
270    READ X,Y
280    LET X1 = X1 + X
290   LET Y1 = Y1 + Y
300   LET X2 = X2 + X*X
310   LET Y2 = Y2 + Y*Y
320   LET Z = Z + X*Y
330 NEXT I
340 LET S1 = N*X2 - X1*X1
350 LET S2 = N*Z - X1*Y1
360 LET B = S2/S1
370 LET Y3 = Y1/N
380 LET X3 = X1/N
390 LET B1 = Y3 - B*X3
400 LET N1 = N - 1
410 LET N2 = N1 - 1
420 LET S3 = (Y2 - Y1*Y3 - B*S2/N)
430 LET S4 = S3/N2
440 PRINT "NUMBER = "N, "SLOPE = "B
450 PRINT "MEAN OF X = "X1/N, "OF Y = "Y1/N
460 PRINT "Y-INTERCEPT = "B1
470 PRINT "SUM-OF-SQUARES","TOTAL", Y2
480 PRINT "    ", "MEAN", Y3*Y1
490 PRINT "    ", "SLOPE", B*S2/N
500 PRINT "    ", "RESIDUAL", S3
510 PRINT "STANDARD DEVIATIONS"
520 PRINT "    ", "X  ", SQR(S1/N/N1)
530 PRINT "    ", "Y  ", SQR((Y2 - Y1*Y3)/N1)
540 PRINT "    ", "ERROR", SQR(S4)
550 PRINT "    ", "Y-BAR", SQR(S4/N)
560 PRINT "    ", "SLOPE", SQR(S4/S1*N)
570 PRINT "F-RATIO FOR SLOPE = " B*S2/N/S4
580 STOP
900 DATA  5
910 DATA  480, 56
920 DATA  500, 61
930 DATA  520, 78
940 DATA  540, 71
950 DATA  560, 82
960END