|
723 | 723 | "# sample N Bernoulli random variables from Ber(0.05).\n", |
724 | 724 | "# each random variable has a 0.05 chance of being a 1.\n", |
725 | 725 | "# this is the data-generation step\n", |
726 | | - "occurances = mc.rbernoulli(p_true, N)\n", |
| 726 | + "occurrences = mc.rbernoulli(p_true, N)\n", |
727 | 727 | "\n", |
728 | | - "print occurances # Remember: Python treats True == 1, and False == 0\n", |
729 | | - "print occurances.sum()" |
| 728 | + "print occurrences # Remember: Python treats True == 1, and False == 0\n", |
| 729 | + "print occurrences.sum()" |
730 | 730 | ], |
731 | 731 | "language": "python", |
732 | 732 | "metadata": {}, |
|
753 | 753 | "cell_type": "code", |
754 | 754 | "collapsed": false, |
755 | 755 | "input": [ |
756 | | - "# Occurances.mean is equal to n/N.\n", |
757 | | - "print \"What is the observed frequency in Group A? %.4f\" % occurances.mean()\n", |
758 | | - "print \"Does this equal the true frequency? %s\" % (occurances.mean() == p_true)" |
| 756 | + "# Occurrences.mean is equal to n/N.\n", |
| 757 | + "print \"What is the observed frequency in Group A? %.4f\" % occurrences.mean()\n", |
| 758 | + "print \"Does this equal the true frequency? %s\" % (occurrences.mean() == p_true)" |
759 | 759 | ], |
760 | 760 | "language": "python", |
761 | 761 | "metadata": {}, |
|
783 | 783 | "collapsed": false, |
784 | 784 | "input": [ |
785 | 785 | "#include the observations, which are Bernoulli\n", |
786 | | - "obs = mc.Bernoulli(\"obs\", p, value=occurances, observed=True)\n", |
| 786 | + "obs = mc.Bernoulli(\"obs\", p, value=occurrences, observed=True)\n", |
787 | 787 | "\n", |
788 | 788 | "#To be explained in chapter 3\n", |
789 | 789 | "mcmc = mc.MCMC([p, obs])\n", |
|
2111 | 2111 | "plt.xlim(0.995, 1)\n", |
2112 | 2112 | "plt.hist(prob_31, bins=1000, normed=True, histtype='stepfilled')\n", |
2113 | 2113 | "plt.title(\"Posterior distribution of probability of defect, given $t = 31$\")\n", |
2114 | | - "plt.xlabel(\"probability of defect occuring in O-ring\")\n" |
| 2114 | + "plt.xlabel(\"probability of defect occurring in O-ring\")\n" |
2115 | 2115 | ], |
2116 | 2116 | "language": "python", |
2117 | 2117 | "metadata": {}, |
|
0 commit comments