Added dump data files for car-iar-poisson models

knitr/car-iar-poisson/bym2.data.R

@@ -0,0 +1,48 @@
+N <-
+56L
+N_edges <-
+132
+node1 <-
+c(1, 1, 1, 1, 2, 2, 3, 3, 4, 4, 4, 5, 5, 5, 5, 6, 7, 7, 7, 7, 
+9, 9, 9, 9, 9, 10, 10, 11, 13, 13, 14, 14, 14, 15, 15, 15, 16, 
+16, 16, 16, 17, 17, 18, 18, 18, 18, 18, 20, 21, 21, 23, 23, 23, 
+23, 23, 24, 24, 24, 24, 24, 24, 24, 24, 25, 25, 26, 26, 26, 27, 
+27, 27, 28, 28, 29, 29, 29, 30, 30, 30, 30, 30, 31, 31, 31, 31, 
+32, 33, 33, 34, 34, 34, 34, 34, 34, 34, 35, 35, 36, 36, 36, 37, 
+37, 38, 38, 38, 38, 38, 39, 39, 40, 40, 40, 41, 41, 41, 42, 42, 
+44, 44, 45, 46, 46, 47, 47, 47, 48, 49, 49, 49, 51, 52, 55)
+node2 <-
+c(5, 9, 11, 19, 7, 10, 6, 12, 18, 20, 28, 11, 12, 13, 19, 8, 
+10, 13, 16, 17, 11, 17, 19, 23, 29, 16, 22, 12, 17, 19, 31, 32, 
+35, 25, 29, 50, 17, 21, 22, 29, 19, 29, 20, 28, 33, 55, 56, 55, 
+29, 50, 29, 34, 36, 37, 39, 27, 30, 31, 44, 47, 48, 55, 56, 26, 
+29, 29, 42, 43, 31, 32, 55, 33, 45, 34, 43, 50, 38, 42, 44, 45, 
+56, 32, 35, 46, 47, 35, 45, 56, 39, 40, 42, 43, 51, 52, 54, 37, 
+46, 37, 39, 41, 41, 46, 42, 44, 49, 51, 54, 40, 41, 41, 49, 52, 
+46, 49, 53, 43, 51, 48, 49, 56, 47, 53, 48, 49, 53, 49, 52, 53, 
+54, 54, 54, 56)
+y <-
+c(9, 39, 11, 9, 15, 8, 26, 7, 6, 20, 13, 5, 3, 8, 17, 9, 2, 7, 
+9, 7, 16, 31, 11, 7, 19, 15, 7, 10, 16, 11, 5, 3, 7, 8, 11, 9, 
+11, 8, 6, 4, 10, 8, 2, 6, 19, 3, 2, 3, 28, 6, 1, 1, 1, 1, 0, 
+0)
+K <-
+2
+x <-
+structure(c(1.6, 1.6, 1, 2.4, 1, 2.4, 1, 0.7, 0.7, 1.6, 0.7, 
+1.6, 1, 2.4, 0.7, 1.6, 1, 0.7, 0.7, 1, 0.7, 1.6, 1, 0.7, 0.1, 
+0.1, 0.7, 0.7, 1, 1, 0.7, 2.4, 1, 0.7, 0.7, 0, 1, 0.1, 1.6, 0, 
+0.1, 1.6, 1.6, 0, 0.1, 0.7, 0.1, 0.1, 0, 0.1, 0.1, 0, 0.1, 0.1, 
+1.6, 1, 3.2, 3.2, 2, 4.8, 2, 4.8, 2, 1.4, 1.4, 3.2, 1.4, 3.2, 
+2, 4.8, 1.4, 3.2, 2, 1.4, 1.4, 2, 1.4, 3.2, 2, 1.4, 0.2, 0.2, 
+1.4, 1.4, 2, 2, 1.4, 4.8, 2, 1.4, 1.4, 0, 2, 0.2, 3.2, 0, 0.2, 
+3.2, 3.2, 0, 0.2, 1.4, 0.2, 0.2, 0, 0.2, 0.2, 0, 0.2, 0.2, 3.2, 
+2), .Dim = c(56L, 2L))
+E <-
+c(1.4, 8.7, 3, 2.5, 4.3, 2.4, 8.1, 2.3, 2, 6.6, 4.4, 1.8, 1.1, 
+3.3, 7.8, 4.6, 1.1, 4.2, 5.5, 4.4, 10.5, 22.7, 8.8, 5.6, 15.5, 
+12.5, 6, 9, 14.4, 10.2, 4.8, 2.9, 7, 8.5, 12.3, 10.1, 12.7, 9.4, 
+7.2, 5.3, 18.8, 15.8, 4.3, 14.6, 50.7, 8.2, 5.6, 9.3, 88.7, 19.6, 
+3.4, 3.6, 5.7, 7, 4.2, 1.8)
+scaling_factor <-
+0.485317494985636

knitr/car-iar-poisson/bym2_islands.stan

@@ -28,10 +28,10 @@ transformed data {
   for (i in 1:N_components) {
     component_starts[i] = c_offset;
     c_offset = c_offset + nodes_per_component[i];
-    component_ends[i] = c_offset - 1;
+    component_ends[i] = c_offset - 2;
   }
   for (i in 1:N_con_comp) {
-    for (j in component_starts[i]:component_ends[i]) {
+    for (j in component_starts[i]:component_ends[i]+1) {
       scaling_factor[j] = scales[i];
     }
   }

knitr/car-iar-poisson/simple_iar.data.R

@@ -0,0 +1,23 @@
+N <-
+56
+N_edges <-
+132
+node1 <-
+c(1, 1, 1, 1, 2, 2, 3, 3, 4, 4, 4, 5, 5, 5, 5, 6, 7, 7, 7, 7, 
+9, 9, 9, 9, 9, 10, 10, 11, 13, 13, 14, 14, 14, 15, 15, 15, 16, 
+16, 16, 16, 17, 17, 18, 18, 18, 18, 18, 20, 21, 21, 23, 23, 23, 
+23, 23, 24, 24, 24, 24, 24, 24, 24, 24, 25, 25, 26, 26, 26, 27, 
+27, 27, 28, 28, 29, 29, 29, 30, 30, 30, 30, 30, 31, 31, 31, 31, 
+32, 33, 33, 34, 34, 34, 34, 34, 34, 34, 35, 35, 36, 36, 36, 37, 
+37, 38, 38, 38, 38, 38, 39, 39, 40, 40, 40, 41, 41, 41, 42, 42, 
+44, 44, 45, 46, 46, 47, 47, 47, 48, 49, 49, 49, 51, 52, 55)
+node2 <-
+c(5, 9, 11, 19, 7, 10, 6, 12, 18, 20, 28, 11, 12, 13, 19, 8, 
+10, 13, 16, 17, 11, 17, 19, 23, 29, 16, 22, 12, 17, 19, 31, 32, 
+35, 25, 29, 50, 17, 21, 22, 29, 19, 29, 20, 28, 33, 55, 56, 55, 
+29, 50, 29, 34, 36, 37, 39, 27, 30, 31, 44, 47, 48, 55, 56, 26, 
+29, 29, 42, 43, 31, 32, 55, 33, 45, 34, 43, 50, 38, 42, 44, 45, 
+56, 32, 35, 46, 47, 35, 45, 56, 39, 40, 42, 43, 51, 52, 54, 37, 
+46, 37, 39, 41, 41, 46, 42, 44, 49, 51, 54, 40, 41, 41, 49, 52, 
+46, 49, 53, 43, 51, 48, 49, 56, 47, 53, 48, 49, 53, 49, 52, 53, 
+54, 54, 54, 56)