summerizePgmmRJMCMC | R Documentation |
summerizePgmmRJMCMC
summerizePgmmRJMCMC(pgmmResList, trueCluster = NULL)
pgmmResList |
result list from pgmmRJMCMC |
trueCluster |
true cluster allocation |
library("fabMix") library("mclust") library("pgmm") library("mvtnorm") library("mcmcse") library("MASS") library("gtools") n <- 50 p <- 10 q <- 4 K <- 10 syntheticDataset <- simData( sameLambda = TRUE, sameSigma = TRUE, K.true = K, n = n, q = q, p = p, sINV_values = 1 / ((1:p)) ) nsim <- 5 burn <- 0 X <- t(syntheticDataset$data) qnew <- 4 Mstep <- 1 Vstep <- 1 constraint <- c(0, 0, 0) mInit <- 20 mVec <- c(1, 20) res <- pgmmRJMCMC(X, mInit, mVec, qnew, niter = nsim, burn = burn, constraint = constraint, Mstep = Mstep, Vstep = Vstep ) summerizePgmmRJMCMC(res, syntheticDataset$class)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.