MBASIC.full: Simultaneously fitting MBASIC model for different numbers of...

In chandlerzuo/mbasic: Matrix Based Analysis for State-space Inference and Clustering.

Description

Simultaneously fitting MBASIC model for different numbers of clusters.

Usage

MBASIC.full(J = NULL, struct = NULL, out = NULL, ncores = 1, ...)

Arguments

J	A vector of the numbers of clusters to be identified.
struct	A list of K by J matrix indicating the structures of each cluster. Default: NULL.
out	The prefix of the file names to write model output information. Default: NULL, when the information is written to the R output.
ncores	The number of CPUs to be used in parallele. Default: 1.
...	Other parameters that are passed to MBASIC.

Details

This function fits MBASIC models with different numbers of clusters and/or different structures in parallel.

Value

A list object including the following fields:

allFits	A list of MBASICFit objects with different numbers of clusters and/or structures.
BestFit	The best model with minimum BIC.
Time	The time used for the model fitting.

Author(s)

Chandler Zuo [e-mail]

See Also

MBASIC

Examples

## Simulate a dataset
dat.sim <- MBASIC.sim(xi = 2, family = "lognormal", I = 1000, fac = rep(1:10, each = 2), J = 3, S = 3, zeta = 0.1)
## Fit the model
dat.sim.fit <- MBASIC.full(Y = dat.sim$Y, S = 3, fac = rep(1:10, each = 2), J = 3:5, maxitr = 3, para = NULL, family = "lognormal", method = "MBASIC", zeta = 0.1, tol = 1e-6, tol.par = 0.001)

chandlerzuo/mbasic documentation built on May 13, 2019, 3:24 p.m.