FacPad (version 2.0)

gibbs2: A Collapsed Gibbs Sampling Algorithm for the Inference of Sparse Bayesian Factor Models_version2

Description

In each iteration, the algorithm iteratively updates each entry in the binary matrix Z, loading matrix W and factor activity matrix X, as well as other model parameters.

Usage

gibbs2(matrixY, matrixL,eta0,eta1,alpha_tau = 1, 
beta_tau = 0.01, tau_sig = 1, max_iter = 10000, 
thin = 10, alpha_sigma = 0.7, beta_sigma = 0.3, file_name)

Value

The algorithm will store the inferred binary indicator matrix Z, loading matrix W and factor activity matrix X (as well as tau_g if not pre-difined) in each thinned iteration and write them into .RData file with name defined by the user.

Examples

data(matrixY)
data(matrixL)
result<-gibbs2(matrixY,matrixL,eta0=0.2,eta1=0.2,
max_iter=50,thin=10,file_name="test_v2_50iter.RData")