Dear BUGS friends,

I received some replies to my "index out of range". Thanks to the replies (really thank you all!), I was able to determine that the problem was that the length of my adj (adjacency) vector did not agree with the num (counts of adjacencies) vector.

However, I now receive the following error on the bottom of the window when I load the 'inits': "adjacency matrix for normal CAR b[1] must be symmetric". I originally used Arcview sw (ISTAT official shapefile) and then imported the geographic file into Geobugs by Maps2winbugs 0.11. I used Geobugs to calculate the adj matrix (please, let me know if there is some better sw to calculate adj matrix) and then I also used a crude macro in SAS to check the symmetry of the matrix (some errors were found). The matrix should be symmetric and I canít understand where is the lack.

Here my model and adj data:

model

{

for (i in 1 : N) {

 O[i]  ~ dpois(mu[i])

 log(mu[i]) <- log(E[i]) + alpha + beta * depriv[i]  + b[i] + h[i]

 RR[i] <- exp(alpha + beta * depriv[i]  + b[i] + h[i])   # Area-specific relative risk (for maps)

 h[i] ~ dnorm(0, tau.h)        # Unstructured random effects

 }

 # CAR prior distribution for spatial random effects:

 b[1:N] ~ car.normal(adj[], weights[], num[], tau.b)

 for(k in 1:sumNumNeigh) {weights[k] <- 1}

 # Other priors:

 alpha  ~ dflat() 

 beta ~ dnorm(0.0, 1.0E-5)

 tau.b  ~ dgamma(0.5, 0.0005)

 sigma.b <- sqrt(1 / tau.b)

 tau.h  ~ dgamma(0.5, 0.0005)

 sigma.h <- sqrt(1 / tau.h)}

 

#..........

 

adj = c(18, 8, 6, 5, 4, 3, 2,

25, 13, 8, 5, 1,

5, 4, 1,

11, 10, 9, 5, 4, 3,

7, 5, 3, 2, 1,

58, 53, 47, 36, 22, 17, 3,

33, 30, 29, 21, 15, 13, 10,

15, 13, 10, 9, 7, 5,

24, 21, 10,

41, 40, 32, 31, 24, 12,

48, 42, 41, 28, 23, 12,

42, 35, 32, 22, 12, 3, 47,

37, 36, 35, 23, 17, 3,

47, 42, 36, 23, 22,

80, 76, 59, 52, 39, 25, 18, 8,

39, 33, 29, 19, 13,

80, 77, 76, 60, 52, 38,

49, 45, 33, 30, 21, 19, 13,

64, 51, 48, 45, 41, 33, 29, 21,

51, 49, 39, 30, 29, 19, 13,

49, 45, 34, 31, 29, 24, 20, 19, 10,

49, 48, 41, 40, 34, 29, 21,

75, 73, 64, 56, 51, 49, 42, 41, 32, 45,

96, 84, 75, 56, 53, 42, 37, 36, 35, 23,

84, 75, 73, 48, 47, 42,

100, 96, 92, 84, 58, 47, 37, 36,

102, 92, 66, 53, 37, 3,

112, 102, 100, 58, 53,

99, 93, 88, 80, 63, 60, 59,

114, 110, 108, 104, 94, 71,

82, 75, 71, 64, 56, 48,

94, 91, 88, 77, 71, 64, 60, 51,

77, 63, 59, 52, 51, 39,

145, 138, 126, 122, 111, 95, 89,

138, 126, 120, 118, 101, 95, 89,

118, 111, 101, 95, 84, 82, 75,

122, 118, 111, 103, 89, 84, 75,

95, 89, 84, 82, 73, 56, 48, 47,

86, 85, 66, 57, 44, 3,

146, 134, 109, 86, 66, 61, 155,

86, 72, 57,

74, 72, 68, 67, 62, 57, 55, 54, 50, 44, 3,

86, 78, 72, 61, 50, 46, 44,

57, 46, 44,

61, 57, 50, 46, 3,

117, 113, 105, 90, 72, 68, 62, 55,

105, 81, 74, 68, 62,

123, 117, 90, 81, 74,

136, 123, 105, 90, 68, 67, 46, 3,

163, 155, 150, 146, 141, 132, 85, 66, 172,

150, 141, 134, 66,

141, 134, 132, 131, 125, 112, 102, 85, 61, 58, 3,

113, 109, 85, 78, 72, 61, 57,

148, 146, 127, 113, 86, 85, 72,

177, 172, 168, 163, 148, 146, 134, 85, 179,

133, 127, 117, 109, 86, 81, 72,

117, 113, 109, 86, 81, 78, 62, 57, 55, 54, 46,

72, 55, 46,

202, 197, 191, 175, 170, 160, 159,

160, 147, 133, 123, 113, 105, 81, 72, 170,

182, 180, 170, 160, 139, 136, 117, 105, 74,

207, 197, 180, 170, 123,

227, 219, 207, 202, 190, 170, 160, 139,

236, 214, 188, 67,

214, 209, 207, 188, 180, 149, 136, 123,

217, 214, 209, 182, 149, 136, 67,

331, 296, 262, 236, 217, 188, 149, 136, 74, 46, 3,

156, 131, 122, 115, 112, 66,

176, 164, 152, 145, 131, 125, 122, 186,

167, 152, 150, 134, 132, 131, 66,

167, 164, 156, 150, 141, 131,

192, 186, 176, 167, 156, 152, 131,

192, 186, 178, 164, 163, 152, 150, 141,

178, 167, 163, 152, 141, 134, 132,

178, 177, 163, 155, 134,

201, 192, 179, 178, 172, 168, 163, 155,

219, 212, 202, 201, 190, 179, 175, 168,

159, 148, 147, 146, 133, 113, 109,

202, 191, 190, 179, 168, 160, 159, 148,

191, 179, 177, 175, 159, 155, 148, 146,

97, 80, 79, 59, 38, 18,

97, 93, 77, 76, 59, 38, 107,

154, 153, 143, 137, 116, 114, 110, 108, 106,

107, 106, 99, 97, 80, 77, 116,

137, 124, 114, 108, 98, 94,

142, 138, 120, 118, 111,

171, 169, 161, 154, 151, 120,

161, 151, 145, 142, 126, 122, 118, 111,

171, 165, 151, 145, 142, 138,

234, 230, 216, 206, 193, 186, 176,

206, 203, 198, 187, 186, 176, 145,

216, 211, 203, 193, 192, 176, 167, 164, 156,

220, 211, 201, 186, 178, 177, 167, 164,

241, 220, 216, 201, 192, 186,

242, 220, 212, 201,

226, 220, 212, 211, 192, 191, 179, 178, 177,

243, 242, 226, 219, 202, 201, 191, 179,

243, 227, 212, 207, 202, 197, 191, 246,

253, 248, 246, 236, 227, 214, 209, 207,

268, 267, 260, 254, 248, 243, 231, 227, 219, 274,

246, 243, 231, 219, 209, 207, 197, 248,

265, 262, 256, 253, 248, 231, 217, 214, 209, 67,

272, 267, 265, 253, 236,

274, 272, 268, 265, 256, 253, 248, 246,

129, 119, 107, 93, 80, 79, 76,

166, 157, 140, 116, 107,

185, 183, 169, 158, 154, 143, 124,

200, 196, 194, 185, 171, 165, 154, 153,

213, 204, 200, 185, 169,

232, 225, 221, 215, 208,

225, 221, 213, 208, 196, 194, 169,

221, 208, 200, 189, 171, 169,

215, 208, 198, 196, 174, 171,

237, 230, 228, 215, 203, 198, 193, 187,

237, 234, 228, 216, 206, 203,

237, 230, 216, 203,

241, 237, 234, 230, 220, 211, 203, 186,

260, 259, 257, 254, 251, 243, 241, 226, 220, 212,

260, 257, 254, 242,

286, 278, 268, 259, 257, 254, 246, 242,

286, 274, 267, 260, 254, 246,

302, 291, 278, 274, 268, 265, 260,

286, 268, 267, 265, 246,

338, 335, 334, 331, 326, 325, 315, 300, 286, 278, 265, 262,

302, 286, 274, 272, 267, 262, 256, 236,

334, 331, 302, 265, 236, 67,

403, 384, 383, 371, 358, 331, 296,

289, 288, 287, 277, 240, 79,

288, 287, 270, 240, 129, 97, 76, 18, 14,

240, 119, 97, 79,

240, 210, 205, 157, 129, 107, 97,

210, 195, 173, 166, 128, 119, 107,

173, 158, 157, 140, 128, 184,

184, 183, 173, 166, 153, 143, 140,

205, 195, 183, 173, 166, 158,

328, 309, 305, 298, 294, 293, 285, 284, 280, 277, 276, 266, 258, 240, 238, 205, 204, 185, 184, 158, 153,

204, 194, 183, 169, 153,

255, 247, 238, 225, 213, 194, 185, 183,

290, 275, 271, 255, 238, 225, 213, 204,

284, 280, 275, 247, 238,

284, 275, 258, 255, 247, 244,

290, 280, 275, 258, 183, 271,

255, 244, 229, 228, 225, 215,

249, 244, 237, 232, 230, 215, 206, 198,

266, 261, 258, 244, 237, 228,

261, 257, 249, 241, 234, 230, 228, 216, 206,

276, 266, 257, 249, 237,

285, 266, 261, 257, 183,

285, 278, 276, 261, 260, 259, 251, 242, 241, 237,

257, 242, 241,

341, 337, 326, 323, 320, 315, 302, 296, 278,

286, 278,

325, 323, 320, 315, 302, 300, 291, 286, 285, 260, 257,

323, 315, 302, 278,

325, 323, 302, 300, 278,

235, 14,

319, 313, 312, 297, 289, 288, 287, 321,

321, 319, 297, 295,

321, 316, 312, 306, 295, 289, 288, 317,

306, 297, 295, 288, 277, 270, 316, 317,

389, 387, 385, 375, 353, 350, 277,

369, 353, 350, 332, 328, 317, 306, 289, 270, 240, 183,

317, 297, 289, 277, 316,

293, 290, 284, 238, 183,

332, 328, 322, 305, 183,

314, 307, 305, 294, 183,

332, 322, 314, 309, 298, 183,

366, 365, 361, 355, 353, 352,

344, 336, 330, 322, 314, 307,

322, 307, 305, 298, 329,

310, 307, 304, 298, 285, 183,

318, 310, 294, 285,

342, 310, 304, 285,

337, 320, 296, 285,

346, 342, 320, 318, 304, 296, 294, 278, 276, 257, 183,

346, 325, 320, 296,

347, 341, 338, 296,

354, 341, 338, 335, 326,

458, 296,

472, 468, 465, 338, 296,

478, 475, 465, 458, 413, 403, 384, 383, 358, 354, 347, 335, 334, 302, 296,

383, 371, 357, 338, 334, 331,

500, 491, 465, 458, 296,

491, 478, 475, 472, 458, 338, 500,

478, 465, 338, 296,

413, 403, 383, 357, 338, 296,

500, 484, 478, 472, 465,

514, 500, 484, 296, 526,

526, 520, 484, 296,

520, 514, 484,

424, 421, 407, 388, 372, 370, 364, 351,

390, 388, 372, 345, 313, 235, 14,

356, 351, 345, 319, 316, 312, 297, 295,

372, 351, 348, 321, 319, 313,

375, 370, 356, 350, 316, 306, 297, 289, 277,

356, 321, 317, 306, 297, 289,

424, 421, 409, 405, 376, 375, 364, 356, 317,

410, 407, 398, 397, 394, 390, 376, 372, 348, 235,

405, 387, 370, 369, 350, 317,

401, 389, 387, 369,

409, 408, 405, 389, 387, 385,

424, 416, 415, 408, 405, 401, 370,

449, 434, 427, 424, 421, 415,

418, 415, 409, 401, 389, 380,

418, 415, 408, 401, 385, 380, 369, 353,

456, 449, 436, 435, 424, 418, 415,

389, 380, 369, 361, 352, 332, 322, 277, 362,

380, 373, 366, 365, 362, 353, 352,

420, 418, 380, 365, 361,

441, 436, 434, 420, 415, 408, 389, 380, 373,

438, 420, 396, 373, 366, 362, 361,

443, 438, 382, 366, 365, 453,

396, 382, 365, 362, 361, 355, 344,

366, 362, 352, 344, 336, 322,

355, 352, 344, 329, 322,

382, 366, 355, 336, 330, 329,

368, 318, 310, 296, 285,

382, 368, 330,

391, 382, 342, 330, 310, 296,

87, 70, 65, 43, 26,

69, 65, 27,

162, 144, 87, 83, 65, 26,

252, 250, 233, 162, 144, 83,

273, 269, 250, 223, 83,

279, 269, 264, 252, 250,

308, 299, 281, 279, 273, 252, 83,

121, 87, 70, 27,

199, 162, 130, 87, 43,

233, 223, 218, 199, 144, 121, 87, 69,

264, 263, 250, 245, 223, 218, 162,

245, 233, 199, 162, 130,

279, 273, 263, 250, 233,

308, 282, 279, 269, 263, 245,

282, 281, 263, 239, 233, 224, 218, 130,

218, 162, 130, 121,

282, 245, 239, 181, 135, 130, 83,

224, 222, 135, 83,

245, 224, 130,

324, 301, 283, 222, 83,

359, 311, 303, 301, 83, 283,

349, 340, 327, 311,

367, 359, 349, 343, 327, 324, 301, 292,

324, 311, 303, 292, 283,

393, 378, 377, 367, 349, 324, 311, 83,

399, 393, 392, 377, 367, 359,

399, 378, 377, 359, 83,

349, 343, 339, 327,

399, 386, 367, 363, 339,

399, 395, 379, 363,

399, 378, 367,

462, 428, 422, 395, 393, 392, 386, 379, 378, 367, 83,

422, 400, 399, 386, 363, 462,

423, 400, 381, 363, 360, 339,

400, 395, 386, 381, 379, 367, 339, 333,

423, 422, 395, 381, 363, 333,

479, 466, 444, 423, 422, 360,

479, 455, 423, 422,

485, 455, 422, 360,

445, 433, 360, 235,

445, 442, 440, 433, 425, 390, 388, 360, 348, 16, 14,

445, 425, 360, 235,

450, 442, 432, 390, 360, 235,

432, 425, 419, 411, 402, 397, 388, 348, 235,

450, 439, 425, 419, 390,

447, 439, 432, 431, 411, 390,

459, 450, 447, 432, 431, 419,

414, 411, 406, 397, 390,

426, 414, 412, 404, 402, 398, 397,

412, 406, 404, 397, 394, 388,

473, 470, 467, 464, 452, 448, 446, 430,

448, 446, 437, 431, 426, 419, 414, 402, 390,

467, 464, 459, 454, 451, 448, 431, 411,

429, 417, 412, 410, 406, 398, 394,

437, 429, 426, 417, 414, 406, 404, 398,

430, 429, 412, 410, 404, 394,

417, 410, 407, 404, 398, 388,

430, 427, 417, 407, 404, 394, 388,

454, 452, 448, 437, 429, 427, 417, 410,

482, 480, 473, 456, 454, 449, 430, 427,

480, 456, 452, 436, 435, 434, 427, 421,

502, 486, 480, 473, 452,

452, 449, 435, 430, 421, 410, 407,

480, 461, 457, 452, 449, 436, 434, 441, 488,

488, 480, 474, 471, 463, 460, 457, 456,

488, 474, 461,

457, 453, 443, 441, 420, 396, 365,

476, 463, 461, 457,

481, 476, 474, 461, 460,

457, 443, 438, 396,

457, 453, 438, 396, 382, 296,

555, 537, 524, 495, 477, 466, 422, 360,

495, 492, 485, 483, 360,

529, 524, 523, 522, 516, 492, 485, 483, 477,

543, 537, 529, 495, 485,

543, 537, 533, 524, 522, 495,

537, 523, 495, 492, 422,

543, 537, 533, 522, 516, 495,

543, 537, 529, 523, 522,

543, 533, 529, 523, 495,

503, 499, 498, 492, 483, 459, 450, 360,

495, 492, 477, 469, 360,

519, 507, 503, 492, 469,

519, 516, 507, 498, 495, 483, 477, 469, 422,

499, 497, 469, 464, 451, 450, 447, 446, 439, 503,

487, 486, 482, 470, 454, 452,

499, 497, 493, 467, 459, 454, 446,

506, 503, 497, 469, 459, 464,

518, 510, 507, 506, 499, 498, 469, 459,

519, 510, 507,

530, 527, 519, 518, 515, 507, 503,

530, 515, 510, 507, 498, 492, 422,

556, 547, 536, 527, 525, 519, 517, 510, 422,

518, 517, 509, 503, 499, 497,

527, 517, 510, 506, 503,

530, 527, 525, 518, 509, 506, 497,

536, 530, 517, 511, 509, 505,

509, 497, 493, 490, 487,

509, 506, 504, 499, 493, 464, 459, 467, 517,

511, 509, 504, 490, 486, 473, 470,

525, 509, 505, 487, 486,

513, 505, 496, 488, 486, 482, 480,

546, 545, 536, 534, 531, 525, 513, 511, 502, 486,

558, 551, 547, 545, 530, 525, 505,

501, 488, 476, 474, 463,

488, 481, 471, 463, 461,

502, 496, 488, 486, 482, 461, 456, 452, 449,

502, 488, 480, 513,

521, 513, 512, 508, 502, 501, 496, 481, 480, 474, 471, 461, 456,

534, 532, 513, 505,

549, 548, 546, 542, 540, 532, 528, 521,

548, 546, 544, 540, 534, 532,

548, 545, 544, 542, 534, 505,

558, 548, 546, 536, 534, 505,

564, 558, 552, 550, 549, 546, 545, 544, 542, 528,

528, 521, 512, 508, 494, 488, 481, 476,

501, 494, 481, 463, 460, 457, 296,

512, 501, 488,

528, 521, 508, 501, 488,

544, 542, 540, 535, 534, 531, 521, 513,

540, 532, 521,

544, 542, 535, 532, 521,

544, 540, 535, 532, 528, 513, 512, 501, 488,

549, 548, 544, 539, 521, 512, 501, 494,

553, 550, 548, 544, 539, 528,

539, 494,

553, 541, 539, 538, 528, 501, 476, 296, 557,

476, 461, 460, 456, 453, 443, 441, 438, 296, 436,

553, 550, 549, 541, 538, 528, 494, 557,

537, 533, 529, 524, 523, 522,

543, 533, 529, 524, 523, 516, 485, 422,

485, 422,

561, 560, 547, 530, 422,

580, 575, 572, 561, 556, 422,

594, 589, 580, 422,

580, 560, 422,

580, 569, 561, 560, 559,

561, 559, 556, 551, 536, 530,

572, 569, 566, 561, 551, 547,

587, 580, 572, 566, 559,

587, 569, 567, 559, 558, 551,

598, 597, 587, 567, 566, 564, 551, 548, 545, 536,

566, 558, 551,

609, 598, 596, 564, 558, 568,

597, 596, 568, 558, 552, 548,

568, 564, 550, 548,

597, 596, 592, 578, 564, 552, 550,

578, 573, 568, 553, 552, 549, 548, 539,

592, 591, 590, 578, 553, 550,

592, 573, 568, 550, 596,

591, 590, 583, 573, 570, 557, 550, 549, 541, 539, 494,

570, 565, 553, 541, 539, 494,

590, 583, 576, 565, 557, 553,

590, 585, 576, 570, 553,

599, 591, 585, 583, 573, 553, 570,

599, 593, 590, 583, 582, 577, 576,

577, 576, 574, 565, 557, 553, 539, 494, 296,

577, 576, 570, 557, 541,

585, 583, 577, 570, 565, 541,

585, 582, 579, 576, 574, 565, 541,

582, 577, 574, 296,

601, 600, 588, 422,

607, 604, 601, 588, 571, 422,

610, 607, 602, 600, 422,

610, 604, 602, 601, 600,

611, 610, 606, 603, 602, 594,

610, 607, 605, 604, 594, 422,

610, 606, 605, 603, 601,

622, 621, 619, 613, 611, 605, 603, 601, 589, 587,

611, 606, 605, 594, 589,

621, 608, 606, 587,

589, 587, 586, 575, 572, 569, 560, 422,

606, 603, 594, 587, 586, 580,

613, 608, 606, 598, 589, 580, 569, 566, 558,

624, 621, 615, 613, 598, 587,

615, 609, 608, 597, 587, 558,

596, 591, 578, 573, 568,

609, 597, 592, 591, 578, 568, 564,

609, 599, 596, 592, 590, 573, 553,

618, 616, 615, 609, 593, 591, 590, 585,

616, 599, 585, 582,

611, 607, 605, 604, 602, 601,

619, 606, 601,

627, 624, 623, 622, 621, 606, 601,

624, 623, 619, 613, 608, 606,

631, 624, 618, 615,

635, 628, 624, 618, 616,

631, 628, 616, 615, 599,

628, 624, 618, 609, 608, 599, 598,

615, 599, 598, 597, 596, 591,

654, 653, 647, 642, 635, 632,

654, 640, 639, 638, 632,

647, 643, 638, 637, 632, 631, 630, 624, 616,

658, 654, 651, 642, 640, 636, 632,

667, 660, 658, 655, 653, 651, 647, 642, 639, 638,

655, 649, 644, 637, 616,

649, 646, 643, 635, 616,

674, 668, 663, 655, 646, 641, 616, 582,

658, 657, 654, 639, 636,

667, 664, 657, 654, 651, 639,

675, 667, 665, 658, 657, 601,

664, 658, 652, 651, 636, 601,

657, 636, 601,

657, 652, 651, 639, 632, 627, 601,

686, 667, 665, 664,

686, 675, 674, 673, 664, 660, 658, 654,

676, 674, 667, 660,

642, 640, 639, 638, 636, 635, 630, 627,

635, 632, 627, 624,

636, 632, 630, 624, 619, 601,

635, 631, 630, 628, 627, 623, 621, 619, 615, 608,

624, 621, 619,

642, 639, 632,

654, 647, 638,

660, 655, 654, 653, 643, 638, 635,

676, 670, 668, 663, 660,

674, 673, 671, 670, 668, 660,

676, 674, 671, 663, 644,

676, 671, 663, 660,

671, 670, 668, 660, 655, 644,

663, 660, 654, 649, 647, 646, 644, 643,

655, 646, 643, 637,

655, 649, 647, 637, 635,

676, 673, 671, 670, 667, 663, 655, 654, 647,

374, 83,

422,

581, 422, 83,

601,

633, 581, 422,

296,

582, 296,

582, 296,

620, 582,

620, 296,

645, 633, 581,

662, 645,

633, 601, 422,

665, 633, 601,

672, 641, 620,

686, 683, 665,

674, 641,

672, 669, 620,

669,

693, 669,

693, 669,

690, 662,

681, 662,

681, 679, 645,

681, 679,

679,

681, 645,

612,

645, 612, 581,

690, 683,

694, 683,

39, 38, 25, 18, 6, 1,

15, 13, 11, 9, 7, 6, 4, 2, 1,

39, 33, 30, 29, 25, 19, 15, 11, 6, 5,

39, 38, 13, 8, 6,

24, 21, 20, 19, 15, 12, 11, 7, 3,

32, 28, 24, 23, 10, 3,

37, 22, 3,

19, 13, 11, 10, 5,

11, 7, 5,

31, 28, 21, 20, 12, 10,

41, 40, 34, 28, 24, 21,

56, 48, 47, 35, 32, 23,

53, 47, 37, 35, 22,

60, 59, 51, 39, 38,

60, 52, 51, 38, 33, 30, 25, 13, 8,

71, 64, 63, 60, 52, 49, 48, 39, 33,

45, 40, 31, 21,

45, 41, 34, 31, 28,

49, 48, 45, 40, 32, 28, 31,

112, 92, 66, 58,

106, 99, 91, 77, 63,

108, 106, 99, 94, 88, 63,

110, 108, 106, 98, 91, 71, 63,

114, 101, 98, 71, 120,

101, 89, 75, 73, 71,

120, 114, 111, 104, 89, 82, 71,

104, 101, 98, 94, 82, 73, 64, 63, 51,

176, 156, 145, 138, 131, 125, 118, 115, 103, 95,

122, 115, 100, 96, 95, 84,

125, 122, 112, 103, 100,

103, 96, 95, 89, 75, 56, 53, 47,

103, 100, 84, 53, 47,

125, 115, 102, 100, 92, 66,

90, 81, 72, 68, 55, 46,

90, 81, 74, 62, 46,

155, 148, 134, 127, 109, 85, 168,

81, 72, 62, 54, 46,

197, 190, 175, 170, 159, 147, 133, 123, 117,

209, 207, 182, 139, 136, 123,

197, 190, 160, 139, 123, 117,

231, 227, 219, 209, 197, 182, 180, 139,

236, 231, 217, 209, 188, 182,

231, 227, 214, 207, 188, 182, 180, 236,

182, 180, 149, 123, 74, 67, 188,

188, 182, 136, 67,

164, 156, 152, 141, 125, 122, 66,

192, 177, 172, 167, 163, 150, 201,

178, 177, 172, 167, 155, 150, 134,

201, 191, 177, 175, 168, 155, 212,

219, 212, 197, 191, 190, 175,

160, 159, 148, 133, 127, 117,

160, 147, 127, 117, 113,

190, 175, 168, 160, 148, 147, 127,

168, 159, 155, 147, 146, 127, 109, 175,

143, 140, 128, 124, 108, 107, 106, 93,

106, 93, 91, 88, 77,

124, 116, 108, 107, 99, 94, 93, 91, 88,

124, 116, 110, 106, 98, 94, 91,

137, 124, 120, 110, 104, 101, 98,

151, 138, 126, 120, 161,

171, 165, 161, 142, 138, 120,

165, 154, 151, 142, 137, 126, 114, 111, 104, 101,

154, 124, 120, 114, 110,

169, 165, 153, 137, 124, 120,

193, 187, 176, 174, 171, 161, 156, 138, 122, 118,

196, 189, 174, 169, 165, 161, 151, 145,

198, 189, 187, 171, 145,

206, 198, 193, 176, 174, 145,

193, 187, 186, 164, 156, 145, 122, 203,

242, 241, 226, 216, 211, 201, 192,

267, 253, 246, 236, 231, 227,

272, 267, 256, 248, 236, 231,

267, 265, 256, 253,

157, 128, 119, 116, 106, 97, 93, 80,

166, 158, 143, 128, 116,

158, 153, 140, 124, 116,

255, 247, 225, 204, 200, 194,

229, 225, 215, 208, 200, 196,

229, 221, 215, 200, 196, 189,

232, 229, 228, 221, 208, 206, 198, 189,

228, 215, 206, 193, 189, 187, 174,

257, 251, 242, 237, 220, 216, 211,

254, 246, 242, 227, 219, 212,

268, 260, 259, 246, 243, 242,

358, 338, 331, 302, 262,

358, 357, 334, 302, 296, 262, 67,

383, 358, 357,

277, 270, 205, 183, 129, 119, 79,

205, 195, 157, 119,

195, 184, 166, 158, 157,

210, 205, 184, 173, 157,

240, 210, 195, 184, 183, 119,

290, 284, 271, 258, 255, 247, 238,

271, 258, 247, 244, 232, 225, 213, 204, 275,

284, 275, 271, 266, 255, 249, 244, 183,

271, 258, 255, 249, 232, 228, 225,

276, 261, 258, 249, 183,

348, 313, 287, 235, 79, 16,

319, 313, 295, 288, 270, 79, 14,

297, 295, 289, 287, 270, 79,

375, 369, 317, 277,

280, 183,

353, 352, 336, 332, 329, 314, 309, 307, 305, 355,

362, 361, 355, 353, 336, 322,

330, 329, 322, 314, 310, 298, 294,

368, 342, 330, 318, 307, 304, 294,

346, 337, 325, 285, 278,

354, 347, 326, 325, 296,

413, 384, 383, 357, 338,

403, 384, 338,

491, 484, 475, 465, 338, 296,

520, 491, 484, 472, 465, 296,

526, 520, 514, 500, 491, 478, 296,

372, 364, 356, 345, 321, 376,

388, 376, 364, 351, 348, 345,

345, 321, 313, 312, 295, 287,

348, 345, 319, 295, 287, 14,

370, 364, 351, 321, 317, 316,

376, 372, 370, 356, 351,

427, 421, 410, 394, 388, 376,

409, 401, 387, 375, 370,

435, 434, 421, 416, 415, 409, 376, 370,

436, 435, 434, 424, 418, 416, 409, 408, 389, 380,

424, 415, 409,

353, 328, 322, 309, 305, 277,

332, 309, 277, 183,

418, 415, 408, 389, 373, 361, 353,

441, 438, 418, 373, 365,

457, 456, 449, 441, 434, 418, 415,

443, 396, 391, 368, 366, 344, 330, 296,

391, 382, 368, 344, 329, 310, 307,

87, 69, 27, 26,

273, 264, 252, 233, 223,

269, 83,

162, 121, 70, 69, 65, 43, 27,

87, 43, 27,

281, 279, 264, 245, 233,

281, 273, 269, 264, 263,

282, 281, 269, 83,

308, 281, 245, 224, 83,

245, 239, 224, 218, 199, 135, 121,

224, 181, 130,

324, 303, 301, 292, 222,

303, 283, 181, 83,

349, 343, 340, 339, 311,

311, 301, 283,

393, 378, 359,

367, 359, 343, 340, 339, 311, 327,

367, 363, 349, 340, 333, 327, 379,

400, 363, 333,

442, 440, 360, 235,

451, 447, 446, 439, 419, 411,

459, 451, 439, 431, 419,

459, 447, 446, 431,

426, 412, 411, 406, 402,

448, 437, 429, 414, 412, 411, 406,

437, 430, 426, 417, 412, 404,

448, 430, 429, 426, 411, 412,

519, 515, 510, 503, 498, 492,

497, 493, 490, 470, 464, 454, 446,

525, 517, 511, 506, 504, 497, 490, 487,

509, 504, 493, 487, 470, 467,

511, 505, 502, 487, 482, 480, 473, 470,

546, 545, 542, 532, 531, 505,

567, 566, 559, 558, 547, 536,

579, 577, 541, 296,

600, 571, 422,

605, 603, 602, 589, 586, 422,

73, 71, 63, 51, 49, 48,

79, 76, 38, 8, 1,

255, 247, 244, 232, 229, 221, 213, 204, 200,

290, 280, 275, 247, 204, 183,

325, 315, 300, 278,

347, 341, 335, 325, 302,

347, 338, 326, 302,

403, 384, 371, 358, 357, 338,

405, 401, 385, 375, 369,

449, 435, 427, 424, 407, 376, 370,

223, 162, 83, 69,

399, 83,

399, 392, 379, 378, 363, 359, 349, 339, 311,

455, 444, 422, 400, 360, 333,

455, 444, 422,

406, 402, 398, 390, 388,

454, 446, 437, 430, 426, 411,

469, 459, 439, 432, 425, 360,

490, 487, 486, 473, 467, 454,

530, 518, 517, 510,

504, 497, 490, 467, 464,

532, 531, 521, 505, 502, 496, 488,

572, 560, 559, 556, 547,

115, 112, 103, 96, 92, 53,

457, 456, 438, 436, 420, 418,

440, 360, 235,

422, 399, 395, 83,

695, 686, 683, 674,

694, 689, 683, 675, 674, 667, 665,

681, 662,

422,

690, 688, 687, 682, 681, 666, 665, 645, 633,

629, 620, 617, 595, 582, 579, 574, 541, 520, 514, 500, 494, 484, 478, 476, 475, 472, 468, 458, 457, 443, 384, 382, 368, 357, 354, 346, 342, 341, 338, 337, 331, 325, 285, 67, 584,

693, 684, 669, 659, 641, 620,

698, 696, 672, 669,

694, 692, 686, 683, 676, 673, 668, 667, 644, 641,

692, 674, 672, 659, 644, 620, 582,

684, 672, 669, 659, 641, 629, 626, 582, 296,

690, 689, 686, 683, 675, 664, 662, 661, 633, 601,

665, 664, 661, 657, 652, 636, 634, 633, 627, 622, 619, 611, 610, 607, 606, 600, 571, 422, 614,

681, 680, 679, 677, 666, 662, 650, 648, 633, 612, 581,

665, 662, 661, 648, 645, 634, 625, 601, 581, 422,

581, 563, 489, 462, 428, 422, 399, 393, 374, 359, 324, 308, 303, 299, 282, 269, 252, 224, 223, 222, 181, 144, 69,

685, 681, 677, 645, 612, 691,

688, 685, 682, 680, 679, 677, 662, 645,

698, 696, 693, 684, 678, 672, 620,

697, 687, 683, 665, 662,

697, 695, 694, 690, 689, 686, 674, 665,

634, 633, 625, 604, 602, 601, 600, 594, 588, 586, 581, 580, 575, 571, 563, 560, 556, 555, 537, 530, 519, 516, 492, 485, 479, 466, 462, 455, 444, 423, 400, 399, 395, 83, 562, 554,

612, 581, 489, 83,

679, 650, 645, 581, 374, 656,

650, 648, 645, 633, 625, 612, 563, 422, 374, 83,

644, 641, 626, 620, 617, 616, 595, 593, 585, 579, 577, 296,

485, 483, 477, 469, 466, 455, 450, 445, 442, 440, 433, 425, 423, 333, 235,

67, 66, 61, 58, 46, 44, 37, 23, 22, 17, 12, 10, 7, 4, 1, 74,

646, 644, 637, 635, 631, 618, 599, 593, 582), sumNumNeigh = 4370)

Thanks,

Marco   

 

 

ARPAM Dipartimento di Ancona
Servizio di Epidemiologia Ambientale
Via Ruggeri, 3N - 60100 Ancona - Italy
Direct line +39 071 2132784
Mobile +39 347 7943415
Fax +39 071 2132731
 
Mailto: [log in to unmask]
Web: www.arpa.marche.it/doc/htm/frameset_epidemiologia.htm
------------------------------------------------------------------- This list is for discussion of modelling issues and the BUGS software. For help with crashes and error messages, first mail [log in to unmask] To mail the BUGS list, mail to [log in to unmask] Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html Please do not mail attachments to the list. To leave the BUGS list, send LEAVE BUGS to [log in to unmask] If this fails, mail [log in to unmask], NOT the whole list