|
23 | 23 | },
|
24 | 24 | {
|
25 | 25 | "cell_type": "code",
|
26 |
| - "execution_count": 87, |
| 26 | + "execution_count": 156, |
27 | 27 | "outputs": [],
|
28 | 28 | "source": [
|
29 | 29 | "import numpy as np\n",
|
|
52 | 52 | },
|
53 | 53 | {
|
54 | 54 | "cell_type": "code",
|
55 |
| - "execution_count": 88, |
| 55 | + "execution_count": 157, |
56 | 56 | "outputs": [],
|
57 | 57 | "source": [
|
58 | 58 | "def get_pb1():\n",
|
|
84 | 84 | },
|
85 | 85 | {
|
86 | 86 | "cell_type": "code",
|
87 |
| - "execution_count": 89, |
| 87 | + "execution_count": 158, |
88 | 88 | "outputs": [],
|
89 | 89 | "source": [
|
90 | 90 | "def fair_die_emission():\n",
|
|
119 | 119 | },
|
120 | 120 | {
|
121 | 121 | "cell_type": "code",
|
122 |
| - "execution_count": 90, |
| 122 | + "execution_count": 159, |
123 | 123 | "outputs": [
|
124 | 124 | {
|
125 | 125 | "name": "stdout",
|
|
185 | 185 | },
|
186 | 186 | {
|
187 | 187 | "cell_type": "code",
|
188 |
| - "execution_count": 91, |
| 188 | + "execution_count": 160, |
189 | 189 | "outputs": [
|
190 | 190 | {
|
191 | 191 | "name": "stdout",
|
|
218 | 218 | },
|
219 | 219 | {
|
220 | 220 | "cell_type": "code",
|
221 |
| - "execution_count": 137, |
| 221 | + "execution_count": 161, |
222 | 222 | "outputs": [
|
223 | 223 | {
|
224 | 224 | "name": "stdout",
|
|
273 | 273 | },
|
274 | 274 | {
|
275 | 275 | "cell_type": "code",
|
276 |
| - "execution_count": 151, |
| 276 | + "execution_count": 162, |
277 | 277 | "outputs": [
|
278 | 278 | {
|
279 | 279 | "name": "stderr",
|
|
288 | 288 | "output_type": "stream",
|
289 | 289 | "text": [
|
290 | 290 | "pi [0.71477891 0.28522109]\n",
|
291 |
| - "a [[2.72155818e-179 4.85324561e-181]\n", |
292 |
| - " [3.87594275e-180 2.70099412e-179]]\n", |
| 291 | + "a [[nan nan]\n", |
| 292 | + " [nan nan]]\n", |
293 | 293 | "b [[nan nan]\n",
|
294 | 294 | " [nan nan]\n",
|
295 | 295 | " [nan nan]\n",
|
|
320 | 320 | " pi = gamma[0]\n",
|
321 | 321 | " for i in range(2):\n",
|
322 | 322 | " for j in range(2):\n",
|
323 |
| - " a[i,j] = np.sum(eta[:,i,j])\n", |
| 323 | + " a[i,j] = np.sum(eta[:,i,j])/np.sum(gamma[:,i])\n", |
324 | 324 | " for i in range(2):\n",
|
325 | 325 | " for k in range(1,7):\n",
|
326 | 326 | " emissions[k-1,i] = np.sum(gamma[sequence==k, i])/np.sum(gamma)\n",
|
|
0 commit comments