|
9 | 9 | },
|
10 | 10 | {
|
11 | 11 | "cell_type": "code",
|
12 |
| - "execution_count": 36, |
| 12 | + "execution_count": 1, |
13 | 13 | "metadata": {},
|
14 | 14 | "outputs": [],
|
15 | 15 | "source": [
|
|
18 | 18 | },
|
19 | 19 | {
|
20 | 20 | "cell_type": "code",
|
21 |
| - "execution_count": 1, |
| 21 | + "execution_count": 2, |
22 | 22 | "metadata": {},
|
23 | 23 | "outputs": [
|
24 | 24 | {
|
|
46 | 46 | },
|
47 | 47 | {
|
48 | 48 | "cell_type": "code",
|
49 |
| - "execution_count": 2, |
| 49 | + "execution_count": 3, |
50 | 50 | "metadata": {},
|
51 | 51 | "outputs": [],
|
52 | 52 | "source": [
|
|
58 | 58 | },
|
59 | 59 | {
|
60 | 60 | "cell_type": "code",
|
61 |
| - "execution_count": 3, |
| 61 | + "execution_count": 4, |
62 | 62 | "metadata": {},
|
63 | 63 | "outputs": [
|
64 | 64 | {
|
|
99 | 99 | " <tbody>\n",
|
100 | 100 | " <tr>\n",
|
101 | 101 | " <th>0</th>\n",
|
102 |
| - " <td>3542</td>\n", |
103 |
| - " <td>61</td>\n", |
104 |
| - " <td>M</td>\n", |
105 |
| - " <td>1997</td>\n", |
106 |
| - " <td>1.690</td>\n", |
107 |
| - " <td>1.69</td>\n", |
108 |
| - " <td>8</td>\n", |
109 |
| - " <td>1.1</td>\n", |
| 102 | + " <td>5401</td>\n", |
| 103 | + " <td>59</td>\n", |
| 104 | + " <td>F</td>\n", |
| 105 | + " <td>1996</td>\n", |
| 106 | + " <td>1.40</td>\n", |
| 107 | + " <td>1.34</td>\n", |
| 108 | + " <td>5</td>\n", |
| 109 | + " <td>1.0</td>\n", |
110 | 110 | " <td>0</td>\n",
|
111 |
| - " <td>4526</td>\n", |
| 111 | + " <td>4655</td>\n", |
112 | 112 | " <td>0</td>\n",
|
113 | 113 | " <td>NaN</td>\n",
|
114 | 114 | " </tr>\n",
|
115 | 115 | " <tr>\n",
|
116 | 116 | " <th>1</th>\n",
|
117 |
| - " <td>4770</td>\n", |
118 |
| - " <td>57</td>\n", |
119 |
| - " <td>F</td>\n", |
120 |
| - " <td>1998</td>\n", |
121 |
| - " <td>1.260</td>\n", |
122 |
| - " <td>1.17</td>\n", |
123 |
| - " <td>4</td>\n", |
124 |
| - " <td>0.9</td>\n", |
| 117 | + " <td>771</td>\n", |
| 118 | + " <td>75</td>\n", |
| 119 | + " <td>M</td>\n", |
| 120 | + " <td>1996</td>\n", |
| 121 | + " <td>3.23</td>\n", |
| 122 | + " <td>3.61</td>\n", |
| 123 | + " <td>10</td>\n", |
| 124 | + " <td>1.3</td>\n", |
125 | 125 | " <td>0</td>\n",
|
126 |
| - " <td>4089</td>\n", |
| 126 | + " <td>4698</td>\n", |
127 | 127 | " <td>0</td>\n",
|
128 | 128 | " <td>NaN</td>\n",
|
129 | 129 | " </tr>\n",
|
130 | 130 | " <tr>\n",
|
131 | 131 | " <th>2</th>\n",
|
132 |
| - " <td>3374</td>\n", |
133 |
| - " <td>66</td>\n", |
| 132 | + " <td>1502</td>\n", |
| 133 | + " <td>71</td>\n", |
134 | 134 | " <td>F</td>\n",
|
135 | 135 | " <td>1996</td>\n",
|
136 |
| - " <td>1.010</td>\n", |
137 |
| - " <td>1.34</td>\n", |
| 136 | + " <td>1.16</td>\n", |
| 137 | + " <td>1.29</td>\n", |
138 | 138 | " <td>4</td>\n",
|
139 |
| - " <td>0.8</td>\n", |
| 139 | + " <td>0.9</td>\n", |
| 140 | + " <td>0</td>\n", |
| 141 | + " <td>4837</td>\n", |
140 | 142 | " <td>0</td>\n",
|
141 |
| - " <td>2615</td>\n", |
142 |
| - " <td>1</td>\n", |
143 |
| - " <td>Digestive</td>\n", |
| 143 | + " <td>NaN</td>\n", |
144 | 144 | " </tr>\n",
|
145 | 145 | " <tr>\n",
|
146 | 146 | " <th>3</th>\n",
|
147 |
| - " <td>7867</td>\n", |
148 |
| - " <td>53</td>\n", |
| 147 | + " <td>4107</td>\n", |
| 148 | + " <td>61</td>\n", |
149 | 149 | " <td>F</td>\n",
|
150 | 150 | " <td>1997</td>\n",
|
151 |
| - " <td>0.705</td>\n", |
152 |
| - " <td>1.25</td>\n", |
153 |
| - " <td>2</td>\n", |
154 |
| - " <td>0.8</td>\n", |
| 151 | + " <td>2.72</td>\n", |
| 152 | + " <td>2.51</td>\n", |
| 153 | + " <td>10</td>\n", |
| 154 | + " <td>0.9</td>\n", |
155 | 155 | " <td>0</td>\n",
|
156 |
| - " <td>4547</td>\n", |
| 156 | + " <td>4458</td>\n", |
157 | 157 | " <td>0</td>\n",
|
158 | 158 | " <td>NaN</td>\n",
|
159 | 159 | " </tr>\n",
|
160 | 160 | " <tr>\n",
|
161 | 161 | " <th>4</th>\n",
|
162 |
| - " <td>7674</td>\n", |
163 |
| - " <td>50</td>\n", |
| 162 | + " <td>5355</td>\n", |
| 163 | + " <td>55</td>\n", |
164 | 164 | " <td>F</td>\n",
|
165 |
| - " <td>1996</td>\n", |
166 |
| - " <td>0.793</td>\n", |
| 165 | + " <td>1997</td>\n", |
167 | 166 | " <td>1.20</td>\n",
|
168 |
| - " <td>2</td>\n", |
| 167 | + " <td>1.19</td>\n", |
| 168 | + " <td>4</td>\n", |
169 | 169 | " <td>0.8</td>\n",
|
170 | 170 | " <td>0</td>\n",
|
171 |
| - " <td>4761</td>\n", |
| 171 | + " <td>4514</td>\n", |
172 | 172 | " <td>0</td>\n",
|
173 | 173 | " <td>NaN</td>\n",
|
174 | 174 | " </tr>\n",
|
|
178 | 178 | ],
|
179 | 179 | "text/plain": [
|
180 | 180 | " index age sex sample.yr kappa lambda flc.grp creatinine mgus t \\\n",
|
181 |
| - "0 3542 61 M 1997 1.690 1.69 8 1.1 0 4526 \n", |
182 |
| - "1 4770 57 F 1998 1.260 1.17 4 0.9 0 4089 \n", |
183 |
| - "2 3374 66 F 1996 1.010 1.34 4 0.8 0 2615 \n", |
184 |
| - "3 7867 53 F 1997 0.705 1.25 2 0.8 0 4547 \n", |
185 |
| - "4 7674 50 F 1996 0.793 1.20 2 0.8 0 4761 \n", |
| 181 | + "0 5401 59 F 1996 1.40 1.34 5 1.0 0 4655 \n", |
| 182 | + "1 771 75 M 1996 3.23 3.61 10 1.3 0 4698 \n", |
| 183 | + "2 1502 71 F 1996 1.16 1.29 4 0.9 0 4837 \n", |
| 184 | + "3 4107 61 F 1997 2.72 2.51 10 0.9 0 4458 \n", |
| 185 | + "4 5355 55 F 1997 1.20 1.19 4 0.8 0 4514 \n", |
186 | 186 | "\n",
|
187 |
| - " event chapter \n", |
188 |
| - "0 0 NaN \n", |
189 |
| - "1 0 NaN \n", |
190 |
| - "2 1 Digestive \n", |
191 |
| - "3 0 NaN \n", |
192 |
| - "4 0 NaN " |
| 187 | + " event chapter \n", |
| 188 | + "0 0 NaN \n", |
| 189 | + "1 0 NaN \n", |
| 190 | + "2 0 NaN \n", |
| 191 | + "3 0 NaN \n", |
| 192 | + "4 0 NaN " |
193 | 193 | ]
|
194 | 194 | },
|
195 |
| - "execution_count": 3, |
| 195 | + "execution_count": 4, |
196 | 196 | "metadata": {},
|
197 | 197 | "output_type": "execute_result"
|
198 | 198 | }
|
|
251 | 251 | "Auto-assigning NUTS sampler...\n",
|
252 | 252 | "INFO:pymc3:Auto-assigning NUTS sampler...\n",
|
253 | 253 | "Initializing NUTS using adapt_diag...\n",
|
254 |
| - "INFO:pymc3:Initializing NUTS using adapt_diag...\n", |
255 |
| - "Multiprocess sampling (4 chains in 4 jobs)\n", |
256 |
| - "INFO:pymc3:Multiprocess sampling (4 chains in 4 jobs)\n", |
257 |
| - "NUTS: [tau_raw_log__, tau_s_raw_log__, mu, a0, b0]\n", |
258 |
| - "INFO:pymc3:NUTS: [tau_raw_log__, tau_s_raw_log__, mu, a0, b0]\n", |
259 |
| - " 0%| | 0/2000 [00:00<?, ?it/s]" |
| 254 | + "INFO:pymc3:Initializing NUTS using adapt_diag...\n" |
260 | 255 | ]
|
261 | 256 | }
|
262 | 257 | ],
|
263 | 258 | "source": [
|
264 | 259 | "with weibull:\n",
|
265 |
| - " trace = pm.sample(1000, tune=1000, init='adapt_diag',threads=1)" |
266 |
| - ] |
267 |
| - }, |
268 |
| - { |
269 |
| - "cell_type": "code", |
270 |
| - "execution_count": null, |
271 |
| - "metadata": {}, |
272 |
| - "outputs": [], |
273 |
| - "source": [ |
274 |
| - "pm.traceplot(trace, \n", |
275 |
| - " varnames=['mu','alpha', 'beta'], \n", |
276 |
| - " lines=dict(mu=mu_stan, alpha=alpha_stan, beta=beta_stan));" |
| 260 | + " trace = pm.sample(1000, tune=1000, init='adapt_diag',cores=1)" |
277 | 261 | ]
|
278 | 262 | },
|
279 | 263 | {
|
|
0 commit comments