Generates the JAGS code from the Hyde network and uses it to
create an object representing a Bayesian graphical model.
Usage
compileJagsModel(network, data = NULL, ...)
Arguments
network
An object of class HydeNetwork
data
A list of data values to be observed in the nodes. It is
passed to the data argument of rjags::jags. Alternatively,
a data frame representing a policy matrix may be provided to compile
multiple JAGS models.
...
Additional arguments to be passed to jags.model
Details
compileJagsModel is a partial wrapper for
jags.model. Running compileJagsModel(network) is
equivalent to running jags.model(textConnection(writeNetworkModel(network))).
Value
Returns a compiledHydeNetwork object. The jags element
of this object is suitable to pass to coda.samples. Otherwise,
the primary function of the object is plotting the network with
observed data shown.
Author(s)
Benjamin Nutter
See Also
jags.model
Examples
data(PE, package="HydeNet")
Net <- HydeNetwork(~ wells +
pe | wells +
d.dimer | pregnant*pe +
angio | pe +
treat | d.dimer*angio +
death | pe*treat,
data = PE)
compiledNet <- compileJagsModel(Net, n.chains=5)
#* Generate the posterior distribution
Posterior <- HydePosterior(compiledNet,
variable.names = c("d.dimer", "death"),
n.iter = 1000)
Posterior
#* For a single model (ie, not a decision model), the user may choose to
#* use the code{rjags} function code{coda.samples}.
#* However, this does not have a succinct print method
s <- coda.samples(compiledNet$jags,
variable.names = c("d.dimer", "death"),
n.iter=1000)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(HydeNet)
Loading required package: nnet
Loading required package: rjags
Loading required package: coda
Linked to JAGS 4.1.0
Loaded modules: basemod,bugs
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/HydeNet/compileJagsModel.Rd_%03d_medium.png", width=480, height=480)
> ### Name: compileJagsModel
> ### Title: Compile Jags Model from a Hyde Network
> ### Aliases: compileJagsModel
>
> ### ** Examples
>
> data(PE, package="HydeNet")
> Net <- HydeNetwork(~ wells +
+ pe | wells +
+ d.dimer | pregnant*pe +
+ angio | pe +
+ treat | d.dimer*angio +
+ death | pe*treat,
+ data = PE)
>
>
> compiledNet <- compileJagsModel(Net, n.chains=5)
Compiling model graph
Resolving undeclared variables
Allocating nodes
Graph information:
Observed stochastic nodes: 0
Unobserved stochastic nodes: 7
Total graph size: 90
Initializing model
Warning messages:
1: In bind_rows_(x, .id) :
binding factor and character vector, coercing into character vector
2: In bind_rows_(x, .id) :
binding factor and character vector, coercing into character vector
3: In bind_rows_(x, .id) :
binding factor and character vector, coercing into character vector
4: In bind_rows_(x, .id) :
binding factor and character vector, coercing into character vector
>
> #* Generate the posterior distribution
> Posterior <- HydePosterior(compiledNet,
+ variable.names = c("d.dimer", "death"),
+ n.iter = 1000)
> Posterior
d.dimer death chain_index obs_index
1 240.0554 No 1 1
2 228.9688 No 1 2
3 236.8301 No 1 3
4 235.9967 No 1 4
5 236.8104 No 1 5
6 243.2212 No 1 6
7 241.8320 No 1 7
8 234.8408 No 1 8
9 241.5265 No 1 9
10 237.7027 No 1 10
11 236.6772 No 1 11
12 236.6804 No 1 12
13 237.6084 No 1 13
14 243.5160 No 1 14
15 236.8758 No 1 15
16 233.9627 No 1 16
17 242.0520 No 1 17
18 240.3758 No 1 18
19 242.9079 No 1 19
20 233.9463 No 1 20
21 237.7165 No 1 21
22 240.2489 No 1 22
23 244.0848 No 1 23
24 248.4026 No 1 24
25 229.5003 No 1 25
26 233.0004 No 1 26
27 251.4121 No 1 27
28 241.6711 No 1 28
29 263.1135 No 1 29
30 238.9652 No 1 30
31 236.6187 Yes 1 31
32 238.3798 No 1 32
33 247.0881 Yes 1 33
34 356.6785 No 1 34
35 239.9796 No 1 35
36 235.6587 No 1 36
37 240.4534 No 1 37
38 235.3943 No 1 38
39 264.6220 No 1 39
40 321.6461 No 1 40
41 240.2583 No 1 41
42 240.3618 No 1 42
43 316.9179 No 1 43
44 240.5303 No 1 44
45 238.2582 No 1 45
46 239.0557 No 1 46
47 233.9876 No 1 47
48 236.1256 No 1 48
49 237.2944 No 1 49
50 239.6672 No 1 50
51 233.0045 No 1 51
52 236.8800 No 1 52
53 243.7963 No 1 53
54 233.2028 No 1 54
55 240.3231 No 1 55
56 246.8852 No 1 56
57 264.9559 No 1 57
58 300.2878 No 1 58
59 244.0375 No 1 59
60 236.5748 No 1 60
61 235.8641 No 1 61
62 242.9417 No 1 62
63 242.2320 No 1 63
64 230.0906 No 1 64
65 234.1900 No 1 65
66 235.6543 Yes 1 66
67 240.8381 No 1 67
68 258.0946 No 1 68
69 240.6204 No 1 69
70 233.0573 No 1 70
71 243.2229 No 1 71
72 237.2404 No 1 72
73 235.0120 No 1 73
74 233.0425 No 1 74
75 243.8169 No 1 75
76 279.6995 No 1 76
77 263.5627 No 1 77
78 304.4741 No 1 78
79 305.0717 No 1 79
80 251.1467 No 1 80
81 237.7738 No 1 81
82 303.5272 No 1 82
83 266.3612 No 1 83
84 236.7266 No 1 84
85 236.6102 No 1 85
86 313.4194 No 1 86
87 240.2179 No 1 87
88 238.6489 No 1 88
89 239.2500 No 1 89
90 245.2279 No 1 90
91 240.8836 No 1 91
92 239.9402 No 1 92
93 240.7574 Yes 1 93
94 237.1998 No 1 94
95 303.2442 No 1 95
96 228.0362 No 1 96
97 313.8207 No 1 97
98 244.0377 No 1 98
99 245.4594 No 1 99
100 333.2056 No 1 100
101 305.3218 No 1 101
102 233.2238 No 1 102
103 292.1871 No 1 103
104 234.4522 Yes 1 104
105 240.5637 No 1 105
106 238.0909 No 1 106
107 307.2861 No 1 107
108 247.8391 No 1 108
109 313.9613 No 1 109
110 239.4559 No 1 110
111 240.7908 No 1 111
112 236.2179 No 1 112
113 232.7067 No 1 113
114 308.6162 No 1 114
115 246.1673 No 1 115
116 304.4818 No 1 116
117 229.8164 No 1 117
118 237.5734 No 1 118
119 246.6249 Yes 1 119
120 234.8245 No 1 120
121 231.9740 No 1 121
122 239.9252 No 1 122
123 250.1992 No 1 123
124 235.4092 No 1 124
125 309.1576 Yes 1 125
126 240.0302 No 1 126
127 236.6945 No 1 127
128 308.0390 No 1 128
129 306.1418 No 1 129
130 249.5918 No 1 130
131 308.6025 Yes 1 131
132 236.1138 No 1 132
133 236.4675 Yes 1 133
134 337.8868 No 1 134
135 238.1912 No 1 135
136 233.2491 No 1 136
137 240.3126 No 1 137
138 301.6226 Yes 1 138
139 244.8686 No 1 139
140 230.3803 No 1 140
141 244.4653 No 1 141
142 267.7222 No 1 142
143 268.1041 No 1 143
144 239.2603 No 1 144
145 239.1956 No 1 145
146 307.1732 Yes 1 146
147 241.4490 No 1 147
148 306.9009 No 1 148
149 232.0773 No 1 149
150 237.4446 No 1 150
151 321.4513 No 1 151
152 243.7489 No 1 152
153 240.5240 No 1 153
154 247.2318 No 1 154
155 302.0321 No 1 155
156 235.5921 No 1 156
157 261.1949 No 1 157
158 234.9331 No 1 158
159 246.2207 No 1 159
160 302.0851 Yes 1 160
161 235.7651 No 1 161
162 240.5852 No 1 162
163 274.5652 No 1 163
164 248.4437 No 1 164
165 246.0420 No 1 165
166 239.5404 No 1 166
167 246.8227 No 1 167
168 295.8808 No 1 168
169 241.2402 No 1 169
170 252.6831 No 1 170
171 291.6038 Yes 1 171
172 231.4190 No 1 172
173 236.9986 No 1 173
174 333.4928 Yes 1 174
175 238.2346 No 1 175
176 240.4645 No 1 176
177 347.7333 Yes 1 177
178 234.0712 No 1 178
179 238.9691 No 1 179
180 232.2598 No 1 180
181 314.0694 No 1 181
182 234.8495 No 1 182
183 233.5854 No 1 183
184 267.7642 No 1 184
185 314.5363 No 1 185
186 228.9505 No 1 186
187 236.4884 No 1 187
188 234.2897 No 1 188
189 243.2411 No 1 189
190 237.8269 No 1 190
191 237.6370 No 1 191
192 306.4623 No 1 192
193 236.9449 No 1 193
194 242.6451 No 1 194
195 311.6332 No 1 195
196 241.0308 No 1 196
197 273.5399 No 1 197
198 251.8375 Yes 1 198
199 235.0187 No 1 199
200 234.0214 No 1 200
201 232.4204 No 1 201
202 244.7440 No 1 202
203 243.6298 No 1 203
204 245.1712 No 1 204
205 240.3253 No 1 205
206 241.7494 No 1 206
207 239.3493 Yes 1 207
208 236.2473 No 1 208
209 232.0085 No 1 209
210 233.2606 No 1 210
211 248.8653 No 1 211
212 253.9008 No 1 212
213 232.8948 No 1 213
214 311.0262 No 1 214
215 242.9728 No 1 215
216 267.6113 Yes 1 216
217 241.9458 Yes 1 217
218 239.6723 No 1 218
219 341.9919 No 1 219
220 272.2647 No 1 220
221 245.4182 No 1 221
222 241.8452 No 1 222
223 236.1154 No 1 223
224 334.4767 Yes 1 224
225 238.6235 No 1 225
226 309.7394 Yes 1 226
227 239.8097 No 1 227
228 247.2449 No 1 228
229 234.0084 No 1 229
230 238.0224 No 1 230
231 235.0480 No 1 231
232 236.8967 No 1 232
233 245.2227 No 1 233
234 232.2395 No 1 234
235 312.5140 No 1 235
236 308.4960 Yes 1 236
237 236.7503 No 1 237
238 234.6722 No 1 238
239 237.6541 No 1 239
240 240.8878 No 1 240
241 312.0095 Yes 1 241
242 231.5671 No 1 242
243 260.7936 No 1 243
244 236.7478 No 1 244
245 240.7593 No 1 245
246 231.0260 No 1 246
247 239.5327 No 1 247
248 242.8355 No 1 248
249 239.9659 No 1 249
250 243.6496 No 1 250
251 308.6733 No 1 251
252 247.3427 No 1 252
253 243.3965 No 1 253
254 236.4139 No 1 254
255 237.4108 No 1 255
256 251.4074 No 1 256
257 274.2904 No 1 257
258 240.5974 No 1 258
259 307.1792 No 1 259
260 252.9318 No 1 260
261 247.6187 No 1 261
262 299.3641 No 1 262
263 241.7724 No 1 263
264 233.2476 No 1 264
265 245.8048 No 1 265
266 310.9647 No 1 266
267 246.5929 No 1 267
268 234.1665 No 1 268
269 254.0425 No 1 269
270 239.0180 No 1 270
271 305.0225 No 1 271
272 243.4413 No 1 272
273 242.2036 No 1 273
274 317.3357 No 1 274
275 270.8245 No 1 275
276 226.8104 No 1 276
277 230.9723 No 1 277
278 275.9748 No 1 278
279 250.4573 No 1 279
280 240.6471 No 1 280
281 250.6662 No 1 281
282 235.5732 No 1 282
283 244.4798 No 1 283
284 234.0997 No 1 284
285 306.4812 No 1 285
286 248.5174 No 1 286
287 244.9128 No 1 287
288 241.8046 No 1 288
289 302.8245 No 1 289
290 240.0165 No 1 290
291 273.7086 No 1 291
292 231.7180 No 1 292
293 243.3711 No 1 293
294 235.5345 No 1 294
295 236.4194 No 1 295
296 303.9472 No 1 296
297 236.0634 No 1 297
298 312.8428 No 1 298
299 243.0947 No 1 299
300 238.8568 No 1 300
301 235.8516 No 1 301
302 303.1704 No 1 302
303 243.2326 No 1 303
304 231.1800 No 1 304
305 341.7638 No 1 305
306 244.9423 No 1 306
307 246.5463 No 1 307
308 306.6596 Yes 1 308
309 242.7666 No 1 309
310 235.1357 No 1 310
311 232.2319 No 1 311
312 232.5253 Yes 1 312
313 240.7247 No 1 313
314 235.2562 No 1 314
315 242.0085 No 1 315
316 302.7058 Yes 1 316
317 240.5972 No 1 317
318 243.0508 No 1 318
319 310.6307 Yes 1 319
320 239.3277 No 1 320
321 244.1874 No 1 321
322 242.7881 No 1 322
323 241.0625 No 1 323
324 237.4599 No 1 324
325 233.0354 No 1 325
326 239.7789 No 1 326
327 303.6891 No 1 327
328 245.5611 No 1 328
329 240.4593 No 1 329
330 307.9586 Yes 1 330
331 234.6394 No 1 331
332 267.8695 Yes 1 332
333 240.3092 No 1 333
334 232.0038 No 1 334
335 239.3490 No 1 335
336 233.8993 No 1 336
337 235.8179 No 1 337
338 340.6429 No 1 338
339 236.3423 No 1 339
340 311.3741 Yes 1 340
341 228.1563 No 1 341
342 235.6080 No 1 342
343 237.1641 No 1 343
344 278.8061 No 1 344
345 301.4899 No 1 345
346 240.3732 No 1 346
347 339.7007 No 1 347
348 238.1614 No 1 348
349 309.4602 No 1 349
350 240.8313 No 1 350
351 233.9137 No 1 351
352 304.6886 No 1 352
353 233.7468 No 1 353
354 240.4547 No 1 354
355 304.6932 Yes 1 355
356 245.1033 No 1 356
357 242.6141 No 1 357
358 242.6034 No 1 358
359 308.4275 No 1 359
360 239.0835 No 1 360
361 240.9989 No 1 361
362 236.5532 No 1 362
363 265.3860 No 1 363
364 242.5568 No 1 364
365 241.8813 No 1 365
366 247.4449 No 1 366
367 308.5176 No 1 367
368 307.1464 No 1 368
369 269.2151 No 1 369
370 236.7683 No 1 370
371 232.7996 No 1 371
372 243.8774 Yes 1 372
373 263.9413 No 1 373
374 314.2993 No 1 374
375 241.9901 No 1 375
376 243.3460 No 1 376
377 240.7871 No 1 377
378 236.4621 No 1 378
379 245.2165 No 1 379
380 281.9419 No 1 380
381 234.2493 No 1 381
382 242.0849 No 1 382
383 240.9983 No 1 383
384 236.0189 No 1 384
385 240.9716 No 1 385
386 269.3173 No 1 386
387 237.4571 No 1 387
388 239.0852 No 1 388
389 236.7594 No 1 389
390 307.2145 No 1 390
391 243.5391 No 1 391
392 242.0075 No 1 392
393 249.9576 No 1 393
394 240.9574 No 1 394
395 291.8259 No 1 395
396 261.2630 No 1 396
397 255.0137 No 1 397
398 332.8382 No 1 398
399 243.5118 No 1 399
400 242.6284 No 1 400
401 306.5232 Yes 1 401
402 240.4316 No 1 402
403 272.2747 Yes 1 403
404 316.4411 No 1 404
405 270.4623 No 1 405
406 241.5676 No 1 406
407 307.1282 No 1 407
408 272.3460 No 1 408
409 232.9105 No 1 409
410 302.1516 No 1 410
411 238.9420 No 1 411
412 240.0428 No 1 412
413 243.9045 No 1 413
414 311.1270 Yes 1 414
415 262.2377 No 1 415
416 226.1526 No 1 416
417 239.4075 No 1 417
418 262.7464 No 1 418
419 237.1992 No 1 419
420 239.1140 No 1 420
421 269.8893 No 1 421
422 242.4586 No 1 422
423 246.5022 No 1 423
424 224.3018 No 1 424
425 233.1978 No 1 425
426 244.9693 No 1 426
427 307.2683 Yes 1 427
428 229.1313 No 1 428
429 236.8495 No 1 429
430 231.2805 No 1 430
431 274.9995 No 1 431
432 236.2897 No 1 432
433 244.8840 No 1 433
434 304.4227 Yes 1 434
435 239.1514 No 1 435
436 239.1729 No 1 436
437 247.3503 No 1 437
438 224.9750 No 1 438
439 232.2488 No 1 439
440 270.8350 No 1 440
441 239.1676 No 1 441
442 240.3677 No 1 442
443 232.2874 No 1 443
444 241.9341 No 1 444
445 262.8751 No 1 445
446 243.1060 No 1 446
447 240.7895 No 1 447
448 308.4634 Yes 1 448
449 248.9335 No 1 449
450 249.0457 No 1 450
451 247.0744 No 1 451
452 239.7668 No 1 452
453 247.6795 No 1 453
454 306.6280 No 1 454
455 235.1047 No 1 455
456 303.9795 No 1 456
457 239.7722 No 1 457
458 262.6511 No 1 458
459 239.2911 No 1 459
460 245.8000 No 1 460
461 242.3339 No 1 461
462 231.1512 No 1 462
463 240.8658 No 1 463
464 241.8280 No 1 464
465 236.3678 No 1 465
466 279.4997 No 1 466
467 243.0329 No 1 467
468 275.9841 No 1 468
469 241.5224 No 1 469
470 239.2879 No 1 470
471 239.5034 No 1 471
472 245.1415 No 1 472
473 304.5580 No 1 473
474 237.1774 No 1 474
475 240.8737 No 1 475
476 244.9624 No 1 476
477 243.7947 No 1 477
478 251.8780 No 1 478
479 232.4435 No 1 479
480 234.5323 No 1 480
481 247.2508 No 1 481
482 244.2422 No 1 482
483 241.5419 No 1 483
484 318.8758 No 1 484
485 244.3784 No 1 485
486 299.2707 Yes 1 486
487 239.2751 No 1 487
488 240.9240 No 1 488
489 242.9651 No 1 489
490 243.6111 No 1 490
491 234.5892 No 1 491
492 244.0379 No 1 492
493 269.0530 No 1 493
494 265.8738 No 1 494
495 242.0566 No 1 495
496 240.9841 No 1 496
497 246.2239 No 1 497
498 302.3615 No 1 498
499 231.9864 No 1 499
500 248.4218 No 1 500
501 244.4063 No 1 501
502 244.8031 No 1 502
503 235.4623 No 1 503
504 309.4585 Yes 1 504
505 239.3600 No 1 505
506 247.1160 No 1 506
507 270.1305 No 1 507
508 236.3439 Yes 1 508
509 242.2468 No 1 509
510 239.0807 No 1 510
511 270.9809 No 1 511
512 240.0600 Yes 1 512
513 269.4741 No 1 513
514 308.1178 Yes 1 514
515 240.2096 No 1 515
516 241.3405 No 1 516
517 233.5454 No 1 517
518 310.0776 Yes 1 518
519 234.1974 No 1 519
520 242.0505 No 1 520
521 241.9222 No 1 521
522 238.2666 No 1 522
523 234.0851 No 1 523
524 308.0969 Yes 1 524
525 307.6915 No 1 525
526 237.6406 No 1 526
527 304.6651 Yes 1 527
528 240.5769 No 1 528
529 260.8309 No 1 529
530 243.9493 No 1 530
531 238.7649 No 1 531
532 236.0359 No 1 532
533 268.3616 No 1 533
534 241.8040 Yes 1 534
535 234.4815 No 1 535
536 243.7400 No 1 536
537 238.6310 No 1 537
538 243.9569 No 1 538
539 247.6912 No 1 539
540 306.5151 No 1 540
541 239.2594 No 1 541
542 307.4474 No 1 542
543 234.8883 No 1 543
544 248.0528 No 1 544
545 243.1140 No 1 545
546 241.0461 No 1 546
547 231.3711 No 1 547
548 245.2213 No 1 548
549 236.4754 No 1 549
550 271.8141 No 1 550
551 235.5341 No 1 551
552 235.6722 No 1 552
553 234.0819 No 1 553
554 244.5249 No 1 554
555 236.4663 No 1 555
556 303.7446 No 1 556
557 240.7544 No 1 557
558 243.7189 No 1 558
559 231.9391 No 1 559
560 309.7779 No 1 560
561 337.0425 No 1 561
562 237.5868 No 1 562
563 242.7798 No 1 563
564 236.2231 No 1 564
565 306.7160 No 1 565
566 232.3664 No 1 566
567 248.0793 No 1 567
568 237.3606 Yes 1 568
569 309.4453 No 1 569
570 234.8863 No 1 570
571 277.6786 No 1 571
572 234.5649 No 1 572
573 232.0807 No 1 573
574 238.3002 No 1 574
575 237.7041 No 1 575
576 237.0918 No 1 576
577 242.9315 No 1 577
578 262.2703 No 1 578
579 316.4800 No 1 579
580 265.8577 No 1 580
581 251.3631 No 1 581
582 229.6346 No 1 582
583 250.2157 No 1 583
584 338.5449 No 1 584
585 241.5222 No 1 585
586 241.7962 No 1 586
587 271.1830 No 1 587
588 238.7672 No 1 588
589 233.0706 No 1 589
590 236.0410 No 1 590
591 233.9203 No 1 591
592 298.2036 Yes 1 592
593 243.1053 No 1 593
594 315.8495 No 1 594
595 230.8304 No 1 595
596 244.5070 No 1 596
597 235.1002 No 1 597
598 234.0205 No 1 598
599 237.0283 No 1 599
600 249.0888 No 1 600
601 237.2412 No 1 601
602 233.9538 No 1 602
603 234.5847 No 1 603
604 241.5171 No 1 604
605 238.2759 No 1 605
606 309.1580 Yes 1 606
607 233.3175 No 1 607
608 239.5378 No 1 608
609 229.4378 No 1 609
610 247.8263 No 1 610
611 249.2420 No 1 611
612 307.9031 No 1 612
613 267.9032 No 1 613
614 238.9962 No 1 614
615 245.8479 No 1 615
616 235.8448 No 1 616
617 228.5330 No 1 617
618 238.7552 No 1 618
619 279.2717 No 1 619
620 248.2473 No 1 620
621 245.0854 No 1 621
622 271.2500 No 1 622
623 234.0937 No 1 623
624 238.3286 No 1 624
625 309.4446 Yes 1 625
626 239.8316 No 1 626
627 241.6313 No 1 627
628 300.0900 No 1 628
629 242.6895 No 1 629
630 234.5550 No 1 630
631 272.7999 No 1 631
632 308.5198 No 1 632
633 231.7385 No 1 633
634 269.5142 No 1 634
635 244.6430 No 1 635
636 238.4095 No 1 636
637 253.1682 No 1 637
638 241.4307 No 1 638
639 267.8274 No 1 639
640 312.6152 No 1 640
641 240.8667 No 1 641
642 244.1488 No 1 642
643 314.4889 No 1 643
644 245.4820 No 1 644
645 242.9923 No 1 645
646 253.8542 No 1 646
647 306.8319 No 1 647
648 242.8578 No 1 648
649 233.9700 No 1 649
650 236.9383 No 1 650
651 238.8203 No 1 651
652 312.3609 Yes 1 652
653 308.1697 No 1 653
654 239.8117 No 1 654
655 247.9824 No 1 655
656 242.8285 No 1 656
657 235.7331 No 1 657
658 239.4835 No 1 658
659 241.0515 No 1 659
660 306.7133 Yes 1 660
661