blatentPPMC {blatent}

blatentPPMC

Description

Simulates data using parameters from posterior distribution of blatent Markov chain.

Usage

blatentPPMC(
  model,
  nSamples,
  seed = model$options$seed,
  parallel = TRUE,
  nCores = 4,
  type = c("mean", "covariance", "univariate", "bivariate", "tetrachoric", "pearson"),
  lowPPMCpercentile = c(0.025, 0.025, 0, 0, 0.025, 0.025),
  highPPMCpercentile = c(0.975, 0.975, 1, 1, 0.975, 0.975)
)

Arguments

model

A blatent MCMC model object.

nSamples

The number of PPMC samples to be simulated.

seed

The random number seed. Defaults to the seed set in the blatent model object.

parallel

If parallelization should be used in PPMC. Defaults to "TRUE".

nCores

If "parallel == TRUE", then specifies the number of cores to use. Defaults to four.

type

The type of statistic to generate, submitted as a character vector. Options include:

• "mean" computes and tabulates the mean for the posterior simulated data for each observed variable.

• "covariance" computes and tabulates the covariance for the posterior simulated data for all pairs of observed variables.

• "univariate" computes and tabulates a Pearson Chi-Square comparing the counts for an observed variable with the counts for a variable from the posterior simulated data, for each observed variable.

• "bivariate"computes and tabulates a Pearson Chi-Square comparing the counts for a pair of observed variables with the counts for a pair of variables from the posterior simulated data, for each pair of observed variables.

• "tetrachoric" computes and tabulates the tetrachoric correlation for the posterior simulated data for all pairs of observed variables.

• "pearson" computes and tabulates the Pearson correlation for the posterior simulated data for all pairs of observed variables.

lowPPMCpercentile

A vector of the lower bound percentiles used for flagging statistics against PPMC predictive distributions. Results are flagged if the observed statistics percentile is lower than the number in the vector. Provided in order of each term in "type". Defaults to "c(.025, .025, 0, 0, .025, .025)".

highPPMCpercentile

A vector of the upper bound percentiles used for flagging statistics against PPMC predictive distributions. Results are flagged if the observed statistics percentile is higher than the number in the vector. Provided in order of each term in "type". Defaults to "c(.975, .975, 1, 1, .975, .975)".