home *** CD-ROM | disk | FTP | other *** search
/ PC-Blue - MS DOS Public Domain Library / PC-Blue MS-DOS Public Domain Library - NYACC.iso / vol030 / regress.bas < prev    next >
Encoding:
BASIC Source File  |  1987-01-11  |  1.5 KB  |  63 lines
  1. 100 CLS
  2. 110 A$=STRING$(80,205)
  3. 120 PRINT A$
  4. 130 PRINT TAB(27)"REGRESSION ANALYSIS #3001"
  5. 140 COLOR 23
  6. 150 PRINT :PRINT :PRINT TAB(39)"IPCO"
  7. 160 COLOR 7
  8. 170 PRINT :PRINT :PRINT TAB(25)"THE PERSONAL COMPUTER OWNERS GROUP"
  9. 180 PRINT :PRINT :PRINT TAB(19)"p.o. box 10426, pittsburgh,pennsylvania 15234"
  10. 190 PRINT A$
  11. 200 PRINT :PRINT :PRINT :PRINT :PRINT
  12. 210 PRINT TAB(28)"PRESS ANY KEY TO CONTINUE"
  13. 220 A$=INKEY$:IF A$="" THEN 220
  14. 230 CLS
  15. 240 CLS
  16. 250 PRINT TAB(30)"LINEAR REGRESSION":PRINT :PRINT
  17. 260 INPUT "NUMBER OF KNOWN POINTS";N
  18. 270 J=0
  19. 280 K=0
  20. 290 L=0
  21. 300 M=0
  22. 310 R2=0
  23. 320 REM LOOP TO ENTER COORDINATES OF POINTS
  24. 330 FOR I=1 TO N
  25. 340 PRINT "X,Y OF POINT ";I;
  26. 350 INPUT X,Y
  27. 360 REM ACCUMULAT INTERMEDIATE SUMS
  28. 370 J=J+X
  29. 380 K=K+Y
  30. 390 L=L+X^2
  31. 400 M=M+Y^2
  32. 410 R2=R2+X*Y
  33. 420 NEXT I
  34. 430 B=(N*R2-K*J)/(N*L-J^2)
  35. 440 A=(K-B*J)/N
  36. 450 PRINT
  37. 460 PRINT "F(X)=";A;"+(";B;"*X)"
  38. 470 REM COMPUTE REGRESSION ANALYSIS
  39. 480 J=B*(R2-J*K/N)
  40. 490 M=M-K^2/N
  41. 500 K=M-J
  42. 510 PRINT
  43. 520 R2=J/M
  44. 530 PRINT "COEFFICIENT OF DETERMINATION (R^2)=";R2
  45. 540 PRINT "COEFFICIENT OF CORRELATION= ";SQR(R2)
  46. 550 PRINT "STANDARD ERROR OF ESTIMATE=";SQR(K/(N-2))
  47. 560 PRINT
  48. 570 REM ESTIMATE Y-COORDINATES OF POINTS WITH ENTERED X-COORDINATES
  49. 580 PRINT "INTERPOLATION:(ENTER X=0 TO END PROGRAM)"
  50. 590 PRINT "X=";
  51. 600 INPUT X
  52. 610 REM RESTART OR END PROGRAM
  53. 620 IF X=0 THEN 660
  54. 630 PRINT "Y=";A+B*X
  55. 640 PRINT
  56. 650 GOTO 590
  57. 660 END
  58. =";
  59. 600 INPUT X
  60. 610 REM RESTART OR END PROGRAM
  61. 620 IF X=0 THEN 660
  62. 630 PRINT "Y=";A+B*X
  63. 640