update_SigmaINV_faster: Gibbs sampling for Sigma^{-1}

View source: R/fabMix.R

update_SigmaINV_fasterR Documentation

Gibbs sampling for \Sigma^{-1}

Description

Gibbs sampling for updating \Sigma^{-1} for the xCU model.

Usage

update_SigmaINV_faster(x_data, z, y, Lambda, mu, K, alpha_sigma, beta_sigma)

Arguments

x_data

n\times p matrix containing the observed data

z

Allocation vector

y

n\times q matrix containing the latent factors

Lambda

K\times p\times q array with factor loadings

mu

K\times p array containing the marginal means

K

Number of components

alpha_sigma

Prior parameter alpha

beta_sigma

Prior parameter beta

Value

p\times p matrix with the common variance of errors per component \Sigma^{-1}.

Author(s)

Panagiotis Papastamoulis

Examples

library('fabMix')
# simulate some data
n = 8                # sample size
p = 5                # number of variables
q = 2                # number of factors
K = 2                # true number of clusters
sINV_diag = 1/((1:p))    # diagonal of inverse variance of errors
set.seed(100)
syntheticDataset <- simData(sameLambda=TRUE,K.true = K, n = n, q = q, p = p, 
                        sINV_values = sINV_diag)

# use the real values as input and update SigmaINV
update_SigmaINV_faster(x_data = syntheticDataset$data, 
	z = syntheticDataset$class, 
	y = syntheticDataset$factors, 
	Lambda = syntheticDataset$factorLoadings, 
	mu = syntheticDataset$means, 
	K = K, 
	alpha_sigma = 0.5, beta_sigma = 0.5)


fabMix documentation built on May 29, 2024, 2:53 a.m.