ctr <- rsctrl(burn_in = 10, n_eff = 3, step=10, seed = 0, tfixed = FALSE)
mat <- matrix(sample(c(0,1), 50, replace = TRUE), nr = 10)
all_m <- rsampler(mat, ctr)
summary(all_m)
some_m <- rsextrobj(all_m, 1, 2)
summary(some_m)
if (FALSE) {
save(some_m, file = "some.RSobj")
some_new <- load("some.RSobj")
summary(some_new)
}
Run the code above in your browser using DataLab