home *** CD-ROM | disk | FTP | other *** search
/ Play and Learn 2 / 19941.ZIP / 19941 / EDUCMATH / STATS / NORMAL.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1994-02-04  |  3.1 KB  |  59 lines
  1. 1  '                      NORMAL DISTRIBUTION
  2. 2  '               Copyright Tracy L. Gustafson, M.D.
  3. 3  '              Round Rock, Texas. Version 3.2, 1985
  4. 4  ON ERROR GOTO 5000:CHAIN MERGE "EPIMRG",5
  5. 15  DIM D(1,1),CS(1,1),T(1),N$(1),X(1),X2(1),MD(1),SD(1)
  6. 22  DATA "NORMAL DISTRIBUTION",30,21
  7. 30  DQ="Do you want to evaluate another ":PRINT:PRINT LEFT$(DQ,14);":":PRINT
  8. 35  PRINT TAB(13);"1.)  Determine if the mean of your sample is significantly":PRINT TAB(19);"different from the population mean.":PRINT
  9. 40  PRINT TAB(13);"2.)  Determine what percent of test values lie within":PRINT TAB(19);"a given range.":PRINT
  10. 45  PRINT TAB(13);"3.)  Evaluate a known Z value."
  11. 50  LOCATE 16,30:PRINT "Enter choice:":AR=16:AC=44:GOSUB 4800:ASUB=VAL(IP$):IF ABS(ASUB-2)>1 THEN BEEP:GOTO 50
  12. 55  ON ASUB GOTO 60,130,195
  13. 60  CLS:PRINT TAB(11);"COMPARING SAMPLE MEAN TO POPULATION MEAN (two-tailed Z test)":PRINT TAB(11);STRING$(60,205)
  14. 65  LOCATE 5,8:PRINT "(Enter RETURN if you want to enter known MEAN and STD. DEV.)"
  15. 70  LOCATE 4,1:GOSUB 4000
  16. 75  LOCATE 7,15:PRINT "What is the SAMPLE NUMBER you wish to analyze?";:AR=7:AC=63:GOSUB 4200
  17. 80  LOCATE 9,20:PRINT "This sample has a MEAN of";X(NS)/T(NS)
  18. 85  LOCATE 11,15:PRINT "What is the POPULATION MEAN for this variable?":AR=11:GOSUB 4800:M=VAL(IP$)
  19. 90  XZ=(X(NS)/T(NS)-M)/SD(NS):GOSUB 220
  20. 95  LOCATE 25,10:PRINT DQ;"MEAN using this DATAFILE?  ";:INPUT;"",A$:IF A$="y" OR A$="Y" THEN CLS:GOTO 75 ELSE 210
  21. 100  PRINT:PRINT:PRINT TAB(25);
  22. 105  INPUT "What is the SAMPLE MEAN?   ",SM:PRINT:PRINT TAB(21);
  23. 110  INPUT "What is the POPULATION MEAN?   ",M:PRINT:PRINT TAB(11);
  24. 115  INPUT "What is the SAMPLE STANDARD DEVIATION?   ",SV
  25. 120  XZ=(SM-M)/SV:GOSUB 220
  26. 125  LOCATE 25,12:PRINT DQ;:INPUT;"known SAMPLE MEAN?  ",A$:IF A$="y" OR A$="Y" THEN CLS:GOTO 100 ELSE 210
  27. 130  CLS:PRINT TAB(21);"PERCENT OF TEST VALUES IN A GIVEN RANGE":PRINT TAB(21);STRING$(39,205)
  28. 135  LOCATE 5,23:INPUT "What is the MEAN of this test?   ",SM
  29. 140  LOCATE 7,15:INPUT "What is the STANDARD DEVIATION of this test?   ",SV
  30. 145  LOCATE 9,23:PRINT "Enter limits of RANGE in question:":PRINT TAB(23);"LOWER LIMIT";TAB(45);"UPPER LIMIT"
  31. 150  AR=CSRLIN:AC=27:GOSUB 4800:LL=VAL(IP$):AC=49:GOSUB 4800:UL=VAL(IP$)
  32. 155  XZ=(LL-SM)/SV:GOSUB 180:PC=P
  33. 160  XZ=(UL-SM)/SV:GOSUB 180:PC=(PC-P)*100
  34. 165  COLOR CLR2,CLR1:GOSUB 250
  35. 170  PRINT TAB(13);"The PERCENT of test values in this range = ";:IF PC<100 THEN PRINT PC; ELSE PRINT "> 99.99999";
  36. 175  PRINT TAB(80):COLOR CLR1,CLR2:LOCATE 25,18:PRINT DQ;:INPUT;"TEST RANGE?  ",A$:IF A$="y" OR A$="Y" THEN 130 ELSE 210
  37. 180  IF ABS(XZ)>6 THEN P=0 ELSE GOSUB 240:P=PT*R
  38. 185  IF XZ<0 THEN P=1-P
  39. 190  RETURN
  40. 195  CLS:PRINT TAB(15);"EVALUATING Z VALUE - NORMAL DISTRIBUTION (two-tailed)":PRINT TAB(15);STRING$(53,205)
  41. 200  LOCATE 5,30:INPUT "Enter Z value:   ",XZ:GOSUB 220
  42. 205  LOCATE 25,20:PRINT DQ;:INPUT;"Z value?  ",A$:IF A$="y" OR A$="Y" THEN 195
  43. 210  LOCATE 25,4:INPUT;"Do you want to perform another analysis using the NORMAL DISTRIBUTION?  ",A$:IF A$="y" OR A$="Y" THEN 20
  44. 215  GOTO 3000
  45. 220  IF ABS(XZ)>6 THEN P=0:GOTO 225 ELSE GOSUB 240:P=2*PT*R
  46. 225  GOSUB 250:COLOR CLR2,CLR1:PRINT TAB(33);"p = ";
  47. 230  IF P<1E-06 THEN PRINT "< 10 (-6)";TAB(71) ELSE PRINT P;TAB(71)
  48. 235  COLOR CLR1,CLR2:RETURN
  49. 240  R=1/SQR(EXP(XZ*XZ)*6.28319):W=1/((ABS(XZ)*0.231642)+1)
  50. 245  W2=W*W:PT=(W*0.319382)-(W2*0.356564)+(W*W2*1.78148)-(W2*W2*1.82126)+(W2*W2*W*1.33027):RETURN
  51. 250  PLAY "MS O3 L64 G O2 GE L9 E"
  52. 255  PRINT:PRINT:PRINT TAB(10);:RETURN
  53. 4010  IF FILE$="" THEN 100
  54. 4025  ERASE D,CS,T,N$,X,X2,MD,SD
  55. 4030  DIM D(A,C),CS(A,C),T(A),N$(A),X(A),X2(A),MD(A),SD(A)
  56. 5000  BEEP:IF ERR<>53 AND ERR<>71 THEN 5010 ELSE LOCATE 2,10:PRINT "Please place EPISTAT in drive A: (or other default).":PRINT TAB(25);"Press any key to continue:"
  57. 5005  A$=INKEY$:IF A$="" THEN 5005 ELSE RESUME
  58. 5010  ON ERROR GOTO 0:END
  59.