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