data(icecream)
#run MCMC sampler (use way more than 20 draws for actual use
icecream_scr_est <- icecream %>% dplyr::filter(id<20) %>% vd_est_vdm_screen(R=20)
ec_trace_screen(icecream_scr_est, burnin=1)
Run the code above in your browser using DataLab