plotSampledParameter {batchmix}

Plot sampled vector parameter

Description

Plot the sampled values for a sampled vector from the output of the “mixtureModel“ function. Not recommended for large B or P.

Usage

plotSampledParameter(samples, parameter, R = NULL, thin = 1, burn_in = 0)

Arguments

samples

The output of the “mixtureModel“ function.

parameter

The name of the parameter to be plotted (a string).

R

The number of iterations run. Defaults to the number of samples for the cluster membership.

thin

The thinning factor of the sampler. Defaults to 1.

burn_in

The samples at the beginning of the chain to drop. Defaults to 0.

Value

A ggplot object of the values in each sampled batch mean per iteration.

Examples

# Data in matrix format
X <- matrix(c(rnorm(100, 0, 1), rnorm(100, 3, 1)), ncol = 2, byrow = TRUE)

# Observed batches represented by integers
batch_vec <- sample(seq(1, 3), size = 100, replace = TRUE)

# MCMC iterations (this is too low for real use)
R <- 50
thin <- 1

# MCMC samples and BIC vector
samples <- runBatchMix(X, R, thin, batch_vec, "MVN", K_max = 8)

# Plot the sampled value of the cluster means against MCMC iteration
parameter <- "means"
plotSampledParameter(samples, parameter, R, thin)