Get Posterior Predictive Samples
Draw predicted frequencies based on posterior distribution of (a) individual estimates (default) or (b) for a new participant (if numItems
is provided; does not consider continuous or discrete predictors in traitMPT).
posteriorPredictive( fittedModel, M = 100, numItems = NULL, expected = FALSE, nCPU = 4 )
fittedModel |
|
M |
number of posterior predictive samples. As a maximum, the number of posterior samples in |
numItems |
optional: a vector with the number of items per MPT tree to sample predicted data for a new participant (first, a participant vector θ is sampled from the hierarchical posterior; second, frequencies are generated). |
expected |
if |
nCPU |
number of CPUs used for parallel sampling. For large models and many participants, this requires considerable computer-memory resources (as a remedy, use |
by default, a list of M
posterior-predictive samples (i.e., matrices) with individual frequencies (rows=participants, columns=MPT categories). For M=1
, a single matrix is returned. If numItems
is provided, a matrix with samples for a new participant is returned (rows=samples)
## Not run: # add posterior predictive samples to fitted model # (facilitates plotting using ?plotFit) fittedModel$postpred$freq.pred <- posteriorPredictive(fittedModel, M=1000) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.