data(icecream)
#run MCMC sampler (use way more than 50 draws for actual use)
icecream_est <- icecream %>% dplyr::filter(id<20) %>% vd_est_vdm(R=10, keep=1)
#Generate demand predictions
icecream_predicted_demand=
icecream %>% dplyr::filter(id<20) %>%
vd_dem_vdm(icecream_est)
#aggregate
brand_lvl_pred_demand <-
icecream_predicted_demand %>% ec_dem_aggregate("Brand")
Run the code above in your browser using DataLab