if (FALSE) {
data(simdiss, package = "dmbc")
G <- 3
p <- 2
prm.prop <- list(z = 1.5, alpha = .75)
burnin <- 2000
nsim <- 1000
seed <- 2301
set.seed(seed)
control <- list(burnin = burnin, nsim = nsim, z.prop = prm.prop[["z"]],
alpha.prop = prm.prop[["alpha"]], random.start = TRUE, verbose = TRUE,
nchains = 2, thin = 10, store.burnin = TRUE, threads = 2,
parallel = "snow")
sim.dmbc <- dmbc(simdiss, p, G, control)
z <- dmbc_get_configuration(sim.dmbc, chain = 1, est = "mean")
summary(z)
library(bayesplot)
library(ggplot2)
color_scheme_set("mix-pink-blue")
graph <- plot(z, size = 2, size_lbl = 3)
graph + panel_bg(fill = "gray90", color = NA)
}
Run the code above in your browser using DataLab