8999 REM REVISED JUNE, 1972   ANDY ROTH C.C.G.
9000 REM ***** LNREG1 ***** STATISTICS
9001 REM ***** VERSION 1 ***** 7/31/69 *****
9002 REM STATISTICS FOR LINEAR REGRESSION
9003 LET X1,Y1,X2,Y2,Z=0
9004 INPUT"NUMBER OF POINTS";N
9005 FOR I=1 TO N
9006 PRINT"POINT";I;"(X,Y)";:INPUT X,Y
9007 LET X1=X1+X
9008 LET Y1=Y1+Y
9009 LET X2=X2+X*X
9010 LET Y2=Y2+Y*Y
9011 LET Z=Z+X*Y
9012 NEXT I
9013 LET S1=N*X2-X1*X1
9014 LET S2=N*Z-X1*Y1
9015 LET B=S2/S1
9016 LET Y3=Y1/N
9017 LET X3=X1/N
9018 LET B1=Y3-B*X3
9019 LET N1=N-1
9020 LET N2=N1-1
9021 LET S3=(Y2-Y1*Y3-B*S2/N)
9022 LET S4=S3/N2
9023 PRINT "NUMBER = "N,"SLOPE = "B
9024 PRINT "MEAN OF X = "X1/N,"OF Y = "Y1/N
9025 PRINT "Y-INTERCEPT = "B1
9026 PRINT "SUM-OF-SQUARES","TOTAL",Y2
9027 PRINT "    ","MEAN",Y3*Y1
9028 PRINT "    ","SLOPE",B*S2/N
9029 PRINT "    ","RESIDUAL",S3
9030 PRINT "STANDARD DEVIATIONS"
9031 PRINT "    ","X  ",SQR(S1/N/N1)
9032 PRINT "    ","Y  ",SQR((Y2-Y1*Y3)/N1)
9033 PRINT "    ","ERROR",SQR(S4)
9034 PRINT "    ","Y-BAR",SQR(S4/N)
9035 PRINT "    ","SLOPE",SQR(S4/S1*N)
9036 PRINT "F-RATIO FOR SLOPE = "B*S2/N/S4
9900 DATA 5,480,56,500,61,520,78,540,71,560,82
9999 END