home *** CD-ROM | disk | FTP | other *** search
/ Piper's Pit BBS/FTP: ibm 0020 - 0029 / ibm0020-0029 / ibm0028.tar / ibm0028 / MCAD254.ZIP / STAT / GEN_POIS.MCD < prev    next >
Encoding:
Text File  |  1988-12-05  |  2.1 KB  |  101 lines
  1. .MCD 20000 0
  2. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  3. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  4. .CMD SET ORIGIN 0
  5. .CMD SET TOL 0.001000
  6. .CMD MARGIN 0
  7. .CMD LINELENGTH 78
  8. .CMD SET PRNCOLWIDTH 8
  9. .CMD SET PRNPRECISION 4
  10. .TXT 0 40 1 39 
  11. a1,38,49,37
  12. Copyright (c) 1988 by MathSoft, Inc.
  13. .TXT 1 40 1 45 
  14. a1,44,78,43
  15. /EQUATIONS FOR GENERATING POISSON DEVIATES
  16. .TXT 1 -62 1 44 
  17. a1,43,78,42
  18. RANDOM DEVIATES WITH POISSON DISTRIBUTION
  19. .TXT 2 -18 1 75 
  20. a1,74,72,73
  21. This document generates a vector of Poisson distributed random deviates.
  22. .EQN 2 1 1 13 
  23. NRAND~400
  24. .TXT 0 16 1 48 
  25. a1,47,60,46
  26. <-- number of random deviates to be generated
  27. .TXT 1 63 2 75 
  28. a2,74,78,101
  29. Define interval vector using the cumulative distribution and define the 
  30. corresponding count vector:
  31. .EQN 1 -79 1 7 
  32. ▐~1
  33. .TXT 0 8 1 26 
  34. a1,25,69,24
  35. <-- Poisson parameter ▐
  36. .EQN 0 32 1 7 
  37. M~6
  38. .TXT 0 8 1 22 
  39. a1,21,29,20
  40. <-- largest deviate
  41. .TXT 2 -48 2 77 
  42. a2,76,77,147
  43. The tables give the fraction of deviates equal to K and the corresponding 
  44. Poisson probabilities P(K).  The plot shows the comparison graphically.
  45. .EQN 0 79 1 15 
  46. I~0;M+1
  47. .EQN 0 16 1 11 
  48. K~0;M
  49. .EQN 2 -16 5 28 
  50. INT[I~K$(K<I)*e^-▐*▐^K/K!
  51. .EQN 1 -68 9 9 
  52. D[K=
  53. .EQN 1 -10 8 5 
  54. K=
  55. .EQN 0 21 8 9 
  56. P(K){19027}=
  57. .EQN 0 18 10 38 
  58. .8&0&D[K,P(K){1,1,9,28,s}@M&0&K
  59. .EQN 0 74 2 12 
  60. INT[(M+1)~1
  61. .EQN 5 -35 2 10 
  62. OUT[K~K
  63. .TXT 4 -79 1 15 
  64. a1,14,77,13
  65. [Ctrl][PgDn]
  66. .TXT 0 80 1 15 
  67. a1,14,78,13
  68. [Ctrl][PgDn]
  69. .TXT 2 -80 1 44 
  70. a1,43,77,42
  71. The Poisson deviates are in the vector V:
  72. .TXT 0 81 1 45 
  73. a1,44,78,43
  74. Define a vector of uniform random numbers:
  75. .EQN 2 -81 201 11 
  76. V=?
  77. .EQN 0 81 1 19 
  78. J~0;NRAND-1
  79. .EQN 3 0 2 13 
  80. R[J~rnd(1)
  81. .TXT 4 0 1 34 
  82. a1,33,78,32
  83. V holds NRAND Poisson deviates:
  84. .EQN 2 0 3 30 
  85. V[J~(hist(INT,(R{51}){52}J))*OUT
  86. .TXT 5 -1 1 15 
  87. a1,14,78,13
  88. [Ctrl][PgDn]
  89. .TXT 2 0 1 30 
  90. a1,29,78,28
  91. Prepare the histogram plot:
  92. .EQN 2 0 2 11 
  93. BINS[I~I
  94. .EQN 3 1 3 18 
  95. D~hist(BINS,V)/NRAND
  96. .TXT 4 -2 1 39 
  97. a1,38,78,37
  98. Define the Poisson density function.
  99. .EQN 2 25 4 15 
  100. P(k)~e^-▐*▐^k/k!
  101.