hi,
I'm experiencing problems in running car models.
For example, this bernoulli car.normal model with 2 dependent variables
works in openbugs with a small dataset (25 dependent variables, attached
below), and not work in winbugs (the error message "undefined real result"
appears in the "mask" windows after starting the updates).
The same model does not work in openbugs with a slightly larger dataset (n =
119, attached below).
model {
alpha ~ dnorm(0,0.01)
beta[1] ~ dnorm(0,0.01)
beta[2] ~ dnorm(0,0.01)
vrho ~ dnorm(0, 0.2) I(0,)
tau <- 1/vrho
for (i in 1:n) {
dep[i] ~ dbern(p[i])
logit(p[i]) <- rho[i]+xbeta[i]+ alpha
xbeta[i] <- beta[1] * ind_1[i]+beta[2]*ind_2[i]
}
rho[1:n] ~ car.normal(adj[], weights[], num_adj[], tau)
for(k in 1:sumNumNeigh) {
weights[k] <- 1
}
}
I had similar issues with other example models. For instance, when trying to
replicate the analysis in appendix 3 of the book Haining 2003 (spatial data
analysis theory and practice), "something went wrong in procedure Node.Set".
I do not understand where is the problem.
Somebody has experience with these models and can provide a script that works
well?
many thanks
Francesco Ficetola
THE DATA USED FOR THE ANALYSES
## 25 observations, work with openbugs but not with winbugs
list( n=25,
dep=c(1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 0,
0),
ind_1=c(1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1,
0, 0),
ind_2=c(1, 2, 4, 5, 1, 2, 3, 6, 5, 4, 5, 6, 4, 5, 6, 7, 8, 4, 8, 5, 8, 9, 5,
6, 9),
num_adj=c(3, 5, 5, 5, 3, 5, 8, 8, 8, 5, 5, 8, 8, 8, 5, 5, 8, 8, 8, 5, 3, 5, 5,
5, 3),
adj=c(2, 6, 7, 1, 3, 6, 7, 8, 2, 4, 7, 8, 9, 3, 5, 8, 9, 10, 4, 9, 10, 1, 2,
7, 11, 12, 1, 2, 3, 6, 8, 11, 12, 13, 2, 3, 4, 7, 9, 12, 13, 14, 3, 4, 5, 8,
10, 13, 14, 15, 4, 5, 9, 14, 15, 6, 7, 12, 16, 17, 6, 7, 8, 11, 13, 16, 17, 18,
7, 8, 9, 12, 14, 17, 18, 19, 8, 9, 10, 13, 15, 18, 19, 20, 9, 10, 14, 19, 20,
11, 12, 17, 21, 22, 11, 12, 13, 16, 18, 21, 22, 23, 12, 13, 14, 17, 19, 22, 23,
24, 13, 14, 15, 18, 20, 23, 24, 25, 14, 15, 19, 24, 25, 16, 17, 22, 16, 17, 18,
21, 23, 17, 18, 19, 22, 24, 18, 19, 20, 23, 25, 19, 20, 24),
sumNumNeigh=144)
## 119 observations, does not work
list( n=119,
dep=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1),
ind_1=c(-3.103189404, -0.074161100, 1.343733759, -0.341856188, 0.738782744,
0.046027834, -0.575924581, -0.966636803, 0.588965921, -0.074161100,
-0.422984776, -0.299967642, -0.109557255, 0.383466677, 0.774178900,
0.562524328, 0.588965921, -0.441044916, -0.065899249, 0.005198028,
-1.387951245, -1.707345744, -1.460833778, 1.123878399, 1.081989853,
0.955849737, 1.226190109, 1.020690866, 1.448873669, -0.738633895, -1.261102800,
-1.718730577, -1.463956760, -1.463956760, -1.986425666, -3.016436502,
1.161348981, -0.053455676, -1.361645050, -1.035773956, -1.718730577,
0.888934183, 1.201822943, 0.955310977, -0.726190322, -0.569259458, 1.044519272,
-0.866399519, -2.625724280, 0.259971845, -2.505352429, -0.262497061,
0.811110720, 1.046593697, 0.070039157, -1.138085666, 0.441642682, -0.789261990,
-1.407377482, 0.342453954, 0.588965921, -0.205917782, 0.607026061, 0.070039157,
-0.050332694, 0.005198028, -0.110265584, 0.588965921, 0.977032860, 0.899270459,
0.899270459, 0.136954711, 0.899270459, -0.197655932, 0.441642682, 0.709337771,
0.588965921, 0.441642682, 0.509606793, 0.441642682, 0.588965921, 0.070039157,
-0.387588620, 1.659989480, 0.263094827, -0.050332694, -0.625838736,
0.362283555, -1.172844376, -0.050332694, -0.655283709, -0.194532950,
-0.625838736, -0.778300842, -0.255831937, 1.046593697, 1.536789429,
0.978263559, 0.955849737, 1.152895576, 0.040594184, 0.506483811, 0.070039157,
0.506483811, 0.832354905, -0.262497061, 0.588965921, 1.424862346, -0.901795675,
-0.770038992, 0.259971845, 0.551495339, -1.756086857, 0.717599621, 0.033934672,
1.370390365, 1.081989853, 0.832354905, 0.441642682),
ind_2=c(0.50278458, 0.72259997, 0.60533525, 0.55529948, 0.56689244,
0.39900492, 0.45488576, 0.48029833, 0.48780745, 0.49493996, 0.43257942,
0.43582152, 0.34562137, 0.32437156, 0.32184286, 0.32021504, 0.15828237,
0.55781773, 0.65898266, 0.64521698, 0.51575239, 0.33574613, 0.26062330,
0.16922060, 0.52505389, 0.53728810, 0.45326413, 0.49706954, 0.52139518,
0.69812148, 0.69039481, 0.72860903, 0.02874139, 0.03449957, 0.05503147,
0.42266842, 0.36127642, 0.49120833, 0.50212476, 0.56195451, 0.50567450,
0.37214024, 0.43027076, 0.35887468, 0.33844319, 0.72613089, 0.68265321,
0.65267208, 0.60723962, 0.62034666, 0.34557961, 0.85212421, 0.69462643,
0.39626379, 0.36684025, 0.50586308, 0.32320080, 0.42030626, 0.27007022,
0.52079290, 0.51144758, 0.63270656, 0.60321469, 0.41683553, 0.34860648,
0.41571990, 0.48461640, 0.61877749, 0.17018495, 0.01224653, 0.00000000,
0.00000000, 0.17590250, 0.33181416, 0.32497835, 0.00000000, 0.36212528,
0.73127063, 0.82381242, 0.85863270, 0.44769934, 0.48591420, 0.52840606,
0.54505171, 0.11027521, 0.11147037, 0.11642350, 0.07191562, 0.18606148,
0.25148321, 0.26821279, 0.26590689, 0.27806283, 0.30118350, 0.31554809,
0.59571527, 0.45458219, 0.20300381, 0.71936958, 0.46815348, 0.39867541,
0.38717332, 0.38198487, 0.35496596, 0.60095786, 0.52663490, 0.39017296,
0.17752609, 0.21907811, 0.06177426, 0.12876361, 0.16205503, 0.16439350,
0.16720015, 0.00000000, 0.00000000, 0.00000000, 0.38493514, 0.00000000),
num_adj=c(17, 21, 19, 19, 19, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17,
17, 1, 1, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 23, 23, 23, 21, 21,
5, 5, 5, 5, 5, 5, 23, 23, 25, 25, 25, 25, 25, 20, 20, 8, 8, 8, 8, 8, 8, 8, 8,
8, 6, 7, 10, 10, 7, 6, 6, 14, 14, 14, 14, 13, 11, 11, 13, 17, 18, 19, 19, 14,
14, 14, 14, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 15, 11, 17, 21, 21, 23,
24, 25, 25, 20, 21, 14, 6, 6, 6, 6, 6, 6, 6, 4, 4, 4, 4, 4),
adj=c(2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 1, 3, 4, 5,
6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 63, 64, 65, 66, 1, 2, 4, 5, 6,
7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 64, 65, 1, 2, 3, 5, 6, 7, 8, 9,
10, 11, 12, 13, 14, 15, 16, 17, 18, 64, 65, 1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12,
13, 14, 15, 16, 17, 18, 64, 65, 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 13, 14, 15,
16, 17, 18, 1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 1, 2,
3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8,
10, 11, 12, 13, 14, 15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14,
15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 1,
2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8,
9, 10, 11, 12, 14, 15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,
15, 16, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 1,
2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 1, 2, 3, 4, 5, 6, 7, 8,
9, 10, 11, 12, 13, 14, 15, 16, 18, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,
14, 15, 16, 17, 20, 19, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 23, 24, 25, 26, 27, 28, 29, 30,
31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 24, 25,
26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51,
52, 21, 22, 23, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46,
47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 27, 28, 29, 30,
31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24,
25, 26, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51,
52, 21, 22, 23, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46,
47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 30, 31, 32, 33, 34, 35,
36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 29,
31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24,
25, 26, 27, 28, 29, 30, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 51,
52, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 33, 34, 35, 36, 37, 44, 45, 46,
47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 34, 35,
36, 37, 44, 45, 46, 47, 48, 49, 50, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
32, 33, 35, 36, 37, 44, 45, 46, 47, 48, 49, 50, 21, 22, 23, 24, 25, 26, 27, 28,
29, 30, 31, 32, 33, 34, 36, 37, 44, 45, 46, 47, 48, 49, 50, 21, 22, 23, 24, 25,
26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 46, 47, 48, 49, 50, 21, 22, 23, 24,
25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 46, 47, 48, 49, 50, 39, 40, 41,
42, 43, 38, 40, 41, 42, 43, 38, 39, 41, 42, 43, 38, 39, 40, 42, 43, 38, 39, 40,
41, 43, 38, 39, 40, 41, 42, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33,
34, 35, 45, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30,
31, 32, 33, 34, 35, 44, 46, 47, 48, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27,
28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 47, 48, 49, 50, 51, 52, 21, 22,
23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 48, 49,
50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37,
44, 45, 46, 47, 49, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32,
33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 50, 51, 52, 21, 22, 23, 24, 25, 26, 27,
28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 44, 45, 46, 47, 48, 49, 51, 52, 21, 22,
23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 44, 45, 46, 47, 48, 49, 50, 52, 21, 22,
23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 44, 45, 46, 47, 48, 49, 50, 51, 54, 55,
56, 57, 58, 59, 60, 61, 53, 55, 56, 57, 58, 59, 60, 61, 53, 54, 56, 57, 58, 59,
60, 61, 53, 54, 55, 57, 58, 59, 60, 61, 53, 54, 55, 56, 58, 59, 60, 61, 53, 54,
55, 56, 57, 59, 60, 61, 53, 54, 55, 56, 57, 58, 60, 61, 53, 54, 55, 56, 57, 58,
59, 61, 53, 54, 55, 56, 57, 58, 59, 60, 63, 64, 65, 66, 67, 68, 2, 62, 64, 65,
66, 67, 68, 2, 3, 4, 5, 62, 63, 65, 66, 67, 68, 2, 3, 4, 5, 62, 63, 64, 66, 67,
68, 2, 62, 63, 64, 65, 67, 68, 62, 63, 64, 65, 66, 68, 62, 63, 64, 65, 66, 67,
70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 105, 106, 107, 69, 71, 72, 73, 74,
75, 76, 77, 78, 79, 80, 105, 106, 107, 69, 70, 72, 73, 74, 75, 76, 77, 78, 79,
80, 105, 106, 107, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80, 105, 106, 107,
69, 70, 71, 72, 74, 75, 76, 77, 78, 79, 80, 106, 107, 69, 70, 71, 72, 73, 75,
76, 77, 78, 79, 80, 69, 70, 71, 72, 73, 74, 76, 77, 78, 79, 80, 69, 70, 71, 72,
73, 74, 75, 77, 78, 79, 80, 105, 106, 69, 70, 71, 72, 73, 74, 75, 76, 78, 79,
80, 96, 103, 104, 105, 106, 107, 69, 70, 71, 72, 73, 74, 75, 76, 77, 79, 80,
96, 102, 103, 104, 105, 106, 107, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 80,
96, 101, 102, 103, 104, 105, 106, 107, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78,
79, 96, 101, 102, 103, 104, 105, 106, 107, 82, 83, 84, 85, 86, 87, 88, 89, 90,
91, 92, 93, 94, 95, 81, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 81,
82, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 81, 82, 83, 85, 86, 87, 88,
89, 90, 91, 92, 93, 94, 95, 81, 82, 83, 84, 86, 87, 88, 89, 90, 91, 92, 93, 94,
95, 99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85, 87, 88, 89, 90, 91, 92,
93, 94, 95, 99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85, 86, 88, 89, 90,
91, 92, 93, 94, 95, 99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85, 86, 87,
89, 90, 91, 92, 93, 94, 95, 99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85,
86, 87, 88, 90, 91, 92, 93, 94, 95, 99, 100, 101, 102, 103, 104, 81, 82, 83,
84, 85, 86, 87, 88, 89, 91, 92, 93, 94, 95, 98, 99, 100, 101, 102, 103, 104,
81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 92, 93, 94, 95, 98, 99, 100, 101, 102,
103, 104, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 93, 94, 95, 98, 99, 100,
101, 102, 103, 104, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 94, 95, 98,
99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92,
93, 95, 98, 99, 100, 101, 102, 103, 104, 81, 82, 83, 84, 85, 86, 87, 88, 89,
90, 91, 92, 93, 94, 98, 99, 100, 101, 102, 103, 104, 77, 78, 79, 80, 97, 98,
99, 100, 101, 102, 103, 104, 105, 106, 107, 96, 98, 99, 100, 101, 102, 103,
104, 105, 106, 107, 90, 91, 92, 93, 94, 95, 96, 97, 99, 100, 101, 102, 103,
104, 105, 106, 107, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98,
100, 101, 102, 103, 104, 105, 106, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
96, 97, 98, 99, 101, 102, 103, 104, 105, 106, 79, 80, 85, 86, 87, 88, 89, 90,
91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 102, 103, 104, 105, 106, 78, 79, 80,
85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 103, 104,
105, 106, 77, 78, 79, 80, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97,
98, 99, 100, 101, 102, 104, 105, 106, 77, 78, 79, 80, 85, 86, 87, 88, 89, 90,
91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 105, 106, 69, 70, 71,
72, 76, 77, 78, 79, 80, 96, 97, 98, 99, 100, 101, 102, 103, 104, 106, 107, 69,
70, 71, 72, 73, 76, 77, 78, 79, 80, 96, 97, 98, 99, 100, 101, 102, 103, 104,
105, 107, 69, 70, 71, 72, 73, 77, 78, 79, 80, 96, 97, 98, 105, 106, 109, 110,
111, 112, 113, 114, 108, 110, 111, 112, 113, 114, 108, 109, 111, 112, 113, 114,
108, 109, 110, 112, 113, 114, 108, 109, 110, 111, 113, 114, 108, 109, 110, 111,
112, 114, 108, 109, 110, 111, 112, 113, 116, 117, 118, 119, 115, 117, 118, 119,
115, 116, 118, 119, 115, 116, 117, 119, 115, 116, 117, 118),
sumNumNeigh = 1852)
-------------------------------------------------------------------
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
|