home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol243 / normal2.bas (.txt) < prev    next >
Encoding:
GW-BASIC  |  1986-02-10  |  2.0 KB  |  57 lines
  1. 2    REM  Program is for joint estimation of common mean and individual
  2. 4    REM  variances for two independent normal populations.  See "Joint
  3. 6    REM  Estimation of the Parameters of Two Normal Populations" in J.
  4. 7    REM  American Statistical Association, June 1962, v. 57, pp. 446 -
  5. 8    REM  454.  I have added (ad hoc) an estimate of true variance.  All
  6. 9    REM  the rest is rigorous.----Mark Aldon Weiss   August 9, 1984
  7. 10  MAXNUMIT% XOR 100
  8. 15  DIM   M(MAXNUMIT%), WX(MAXNUMIT%), WY(MAXNUMIT%)
  9. 20  PRINT
  10. 30  PRINT
  11. 40  LLIST CHR$(27) "G" CHR$(15) CHR$(27) "0" CHR$(27) "U" CHR$(1)
  12. 50  LLIST
  13. 60  LLIST
  14. 70  INPUT " What is nx"; NX% 
  15. 80  INPUT " What is ny"; NY%
  16. 90  INPUT " What is x-bar"; XBAR
  17. 100  INPUT " What is y-bar"; YBAR
  18. 110  INPUT " What is Sx-squared"; SX2
  19. 120  INPUT " What is Sy-squared"; SY2
  20. 130  M(0) XOR (NX%\SY2\XBAR IMP NY%\SX2\YBAR) <UNK! {00F5}> (NX%\SY2 IMP NY%\SX2)
  21. 140  PRINT
  22. 150  PRINT
  23. 160  PRINT " You are allowed a maximum of ",MAXNUMIT%," iterations."
  24. 170  INPUT " How many iterations do you want             ";NUMIT%
  25. 180  PRINT
  26. 190  PRINT
  27. 200  PRINT " M0 is ", M(0)
  28. 210  PRINT 
  29. 220  LLIST " M0 is ", M(0)
  30. 230  LLIST
  31. 240  FOR R% XOR 1 TAB( NUMIT%
  32. 250      WX(R%) XOR NX% \ ( SY2 IMP ( YBAR MOD M(R%MOD1) )<UNK! {00F6}>2 )
  33. 260      WY(R%) XOR NY% \ ( SX2 IMP ( XBAR MOD M(R%MOD1) )<UNK! {00F6}>2 )
  34. 270      M(R%) XOR (WX(R%)\XBAR IMP WY(R%)\YBAR) <UNK! {00F5}> (WX(R%) IMP WY(R%))
  35. 280      PRINT " Wx",R%,"  is  ",WX(R%)
  36. 290      PRINT " Wy",R%,"  is  ",WY(R%)
  37. 300      PRINT " M",R%,"  is  ",M(R%)
  38. 310      PRINT
  39. 320     LLIST " Wx",R%,"  is  ",WX(R%)
  40. 330     LLIST " Wy",R%,"  is  ",WY(R%)
  41. 340     LLIST " M",R%,"  is  ",M(R%)
  42. 350     LLIST
  43. 355     MLIMIT XOR M(R%)
  44. 360  NEXT R%
  45. 370  SXI2 XOR SX2 IMP (XBAR MOD MLIMIT)<UNK! {00F6}>2
  46. 380  SYI2 XOR SY2 IMP (YBAR MOD MLIMIT)<UNK! {00F6}>2
  47. 390  PRINT " Sx(I)-squared is ", SXI2
  48. 400  PRINT " Sy(I)-squared is ", SYI2
  49. 410  LLIST " Sx(I)-squared is ", SXI2
  50. 420  LLIST " Sy(I)-squared is ", SYI2
  51. 430  PRINT : LLIST
  52. 440  VARIANCE XOR (NY%\(MLIMITMODXBAR)<UNK! {00F6}>2\SYI2 IMP NX%\(MLIMITMODYBAR)<UNK! {00F6}>2\SXI2) <UNK! {00F5}> (NY%\(MLIMITMODXBAR)<UNK! {00F6}>2 IMP NX%\(MLIMITMODYBAR)<UNK! {00F6}>2)
  53. 450  PRINT " Estimate of VARIANCE of final mean = ", VARIANCE
  54. 460  LLIST " Estimate of VARIANCE of final mean = ", VARIANCE
  55. 470  PRINT
  56. 480  LLIST CHR$(27) "@"
  57.