set.seed(123)
y <- c(rnorm(40, 0, 0.3), rnorm(20, 5, 0.3))
g <- c(rep(1:6, each = 10))
out <- fit_fSAN(y = y, group = g, est_method = "MCMC",
mcmc_param = list(nrep = 500, burn = 200))
get_model(out)
get_time(out)
hp <- get_params(out)
sims <- get_sim(out)
Run the code above in your browser using DataLab