Thin 'echoice2'-vd draw objects
vd_thin_draw(est, burnin_perc = 0.5, total_draws = NULL)
is an 'echoice2' draw object (list)
how much burn-in to remove
how many draws to keep after thinning
thinned 'echoice2' draw object (list)
data(icecream)
#run MCMC sampler (use way more than 50 draws for actual use)
icecream_est <- icecream %>% dplyr::filter(id<100) %>% vd_est_vdm(R=50, keep = 1)
#> Using 16 cores
#> MCMC in progress
#> MCMC complete
#> Total Time Elapsed: 0.00 minutes
#without thinning, yields R=50 draWs
dim(icecream_est$MUDraw)
#> [1] 50 21
icecream_est_thinned <- vd_thin_draw(icecream_est,.5)
#> Draws post burnin: 26
#26 draws left after thinning about half
dim(icecream_est_thinned$MUDraw)
#> [1] 26 21