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)
Run the code above in your browser using DataLab