This function takes in a model object, observed data, and the configuration file used to produce the real-data stan data object. It then runs the model under the prior and returns the observable (generated quantities) for some number of draws.
Arguments
- model
a CmdStanModel object defining the model
- real_data
observed data, formated for stan, needed to match, for example, the number of wastewater sites
- config
the output of get_config_vals()
- ndraws
integer, number of draws from the prior predictive distribution to return
- ...
further arguments passed to model$sample(), e.g. # chains
- mclapply_cores
integer specifying number of threads to use for posterior-predictive-analysis-level parallelization, separate from stan parallelization