num_obs <- 18
r <- 
structure(c(0, 1, 2, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 2, 1),
.Dim = c(18))
nr <- 
structure(c(4, 7, 3, 5, 10, 10, 4, 5, 11, 4, 3, 3, 6, 3, 4, 5, 3, 0),
.Dim = c(18))
num_comp <- 4
num_inter <- 7
X_comp <- 
structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
-0.22314355131421, 0.470003629245736, -2.52572864430826, -0.693147180559945,
-0.22314355131421, 0.470003629245736, -1.83258146374831, 0, -0.22314355131421,
0.470003629245736, -1.42711635564015, 0.405465108108164, -Inf, -0.7339691750802,
-Inf, -0.693147180559945, -Inf, -0.22314355131421, -Inf, 0, -Inf,
0.470003629245736, -Inf, 0.405465108108164, -0.7339691750802, -1.83258146374831,
-Inf, -0.693147180559945, -0.22314355131421, -0.22314355131421, -Inf, 0,
-0.22314355131421, 0.470003629245736, -Inf, 0.405465108108164,
-0.22314355131421, 0.693147180559945, -Inf, -0.693147180559945, -Inf, -Inf,
-2.52572864430826, 0, -Inf, -Inf, -1.83258146374831, 0.405465108108164, -Inf,
-Inf, -1.13943428318836, -0.693147180559945, -Inf, -Inf, -0.7339691750802, 0,
-Inf, -Inf, -0.446287102628419, 0.405465108108164, -0.22314355131421, -Inf,
-1.83258146374831, -0.693147180559945, -0.22314355131421, -Inf,
-1.13943428318836, 0, -0.22314355131421, -Inf, -1.42711635564015,
0.405465108108164),
.Dim = c(4, 18, 2))
X_inter <- 
structure(c(1.28, 1.28, 1.28, 0, 0, 0, 0.0768, 0.64, 1.28, 1.6, 0, 0, 0, 0, 0,
0, 0, 0, 0.064, 0.128, 0.192, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.128, 0.256,
0.192, 0.4, 0.8, 1.2, 0, 0, 0, 0.24, 0.8, 1.2, 0.4, 0, 0, 0, 0, 0, 0.4, 0.8,
1.2, 0.128, 0.256, 0.384, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.8, 1.6,
2.4, 0.24, 0.8, 2.4, 0.08, 0.8, 2.4, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0.04, 0.16,
0.36, 0, 0, 0, 0, 0, 0, 0, 0.08, 0.24, 0.16, 0.48, 0.96, 0.08, 0.32, 0.36,
0.0512, 0.2048, 0.4608, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
.Dim = c(18, 7))
num_groups <- 3
num_strata <- 2
group <- 
structure(c(1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3),
.Dim = c(18))
stratum <- 
structure(c(1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2),
.Dim = c(18))
group_stratum_cid <- 
structure(c(1, 2, 2),
.Dim = c(3))
prior_tau_dist <- 1
prior_EX_prob_comp <- 
structure(c(0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9),
.Dim = c(3, 4))
prior_EX_prob_inter <- 
structure(c(0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9,
0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9),
.Dim = c(3, 7))
prior_EX_mu_mean_comp <- 
structure(c(-0.693147180559945, -0.693147180559945, -0.693147180559945,
-0.693147180559945, 0, 0, 0, 0),
.Dim = c(4, 2))
prior_EX_mu_sd_comp <- 
structure(c(1, 1, 1, 1, 0.5, 0.5, 0.5, 0.5),
.Dim = c(4, 2))
prior_EX_tau_mean_comp <- 
structure(c(-1.38629436111989, -0.693147180559945, -1.38629436111989,
-0.693147180559945, -1.38629436111989, -0.693147180559945, -1.38629436111989,
-0.693147180559945, -2.07944154167984, -1.38629436111989, -2.07944154167984,
-1.38629436111989, -2.07944154167984, -1.38629436111989, -2.07944154167984,
-1.38629436111989),
.Dim = c(2, 4, 2))
prior_EX_tau_sd_comp <- 
structure(c(0.707293041387699, 0.707293041387699, 0.707293041387699,
0.707293041387699, 0.707293041387699, 0.707293041387699, 0.707293041387699,
0.707293041387699, 0.707293041387699, 0.707293041387699, 0.707293041387699,
0.707293041387699, 0.707293041387699, 0.707293041387699, 0.707293041387699,
0.707293041387699),
.Dim = c(2, 4, 2))
prior_EX_corr_eta_comp <- 
structure(c(1, 1, 1, 1),
.Dim = c(4))
prior_EX_mu_mean_inter <- 
structure(c(0, 0, 0, 0, 0, 0, 0),
.Dim = c(7))
prior_EX_mu_sd_inter <- 
structure(c(0.707106781186548, 0.707106781186548, 0.707106781186548,
0.707106781186548, 0.707106781186548, 0.707106781186548, 0.707106781186548),
.Dim = c(7))
prior_EX_tau_mean_inter <- 
structure(c(-1.38629436111989, -1.38629436111989, -1.38629436111989,
-1.38629436111989, -1.38629436111989, -1.38629436111989, -1.38629436111989,
-1.38629436111989, -1.38629436111989, -1.38629436111989, -1.38629436111989,
-1.38629436111989, -1.38629436111989, -1.38629436111989),
.Dim = c(2, 7))
prior_EX_tau_sd_inter <- 
structure(c(0.35364652069385, 0.35364652069385, 0.35364652069385,
0.35364652069385, 0.35364652069385, 0.35364652069385, 0.35364652069385,
0.35364652069385, 0.35364652069385, 0.35364652069385, 0.35364652069385,
0.35364652069385, 0.35364652069385, 0.35364652069385),
.Dim = c(2, 7))
prior_EX_corr_eta_inter <- 1
prior_NEX_mu_mean_comp <- 
structure(c(-0.693147180559945, -0.693147180559945, -0.693147180559945,
-0.693147180559945, 0, 0, 0, 0),
.Dim = c(4, 2))
prior_NEX_mu_sd_comp <- 
structure(c(1, 1, 1, 1, 0.5, 0.5, 0.5, 0.5),
.Dim = c(4, 2))
prior_NEX_mu_mean_inter <- 
structure(c(0, 0, 0, 0, 0, 0, 0),
.Dim = c(7))
prior_NEX_mu_sd_inter <- 
structure(c(0.707106781186548, 0.707106781186548, 0.707106781186548,
0.707106781186548, 0.707106781186548, 0.707106781186548, 0.707106781186548),
.Dim = c(7))
prior_is_EXNEX_comp <- 
structure(c(1, 1, 1, 1),
.Dim = c(4))
prior_is_EXNEX_inter <- 
structure(c(1, 1, 1, 1, 1, 1, 1),
.Dim = c(7))
prior_PD <- 0