Prepare parameters setup for PRECAST model fitting.

model_set_PRECAST(
  Sigma_equal = FALSE,
  Sigma_diag = TRUE,
  mix_prop_heter = TRUE,
  error_heter = TRUE,
  Sp2 = TRUE,
  wpca_int = FALSE,
  int.model = "EEE",
  coreNum = 5,
  coreNum_int = coreNum,
  beta_grid = seq(0.2, 4, by = 0.2),
  maxIter_ICM = 6,
  maxIter = 20,
  epsLogLik = 1e-05,
  verbose = TRUE,
  seed = 1
)

Arguments

Sigma_equal

a logical value, whether set the mixture covariance matrices equal, default as FALSE.

Sigma_diag

a logical value, whether set the mixture covariance matrices diagonal, default as TRUE.

mix_prop_heter

a logical value, whether set the smoothing parameter of each data batch to be different, default as TRUE.

error_heter

a logical value, whether use the heterogenous error for PRECAST model, default as TRUE.

Sp2

a logical value, whether add the intrisical CAR component in the PRECAST model, default as TRUE.

wpca_int

an optional logical value, means whether use the weighted PCA to obtain the initial values of loadings and other paramters, default as FALSE which means the ordinary PCA is used.

int.model

an optional string, specify which Gaussian mixture model is used in evaluting the initial values for DR-SC, default as "EEE"; and see Mclust for more models' names.

coreNum

an integer value to decide how many cores are used to run in parallel, default as 5.

coreNum_int

an integer value to decide how many cores are used in parallel computing initial values.

beta_grid

a numeric vector, specify the smoothness parameter of Random Markov Field. The default is seq(0.2,4,0.2).

maxIter_ICM

the maximum iteration of ICM step, default as 6.

maxIter

the maximum iteration of ICM-EM algorithm. The default is 25

epsLogLik

an optional positive vlaue, tolerance vlaue of relative variation rate of the observed pseudo log-loglikelihood value, defualt as '1e-5'.

verbose

a logical value, whether output the information in iteration.

seed

a postive integer, the random seed to be set.

Examples

model_set_PRECAST(Sigma_equal=TRUE, coreNum=5)
#> $Sigma_equal
#> [1] TRUE
#> 
#> $Sigma_diag
#> [1] TRUE
#> 
#> $mix_prop_heter
#> [1] TRUE
#> 
#> $error_heter
#> [1] TRUE
#> 
#> $Sp2
#> [1] TRUE
#> 
#> $wpca_int
#> [1] FALSE
#> 
#> $int.model
#> [1] "EEE"
#> 
#> $coreNum
#> [1] 5
#> 
#> $coreNum_int
#> [1] 5
#> 
#> $beta_grid
#>  [1] 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8
#> [20] 4.0
#> 
#> $maxIter_ICM
#> [1] 6
#> 
#> $maxIter
#> [1] 20
#> 
#> $epsLogLik
#> [1] 1e-05
#> 
#> $verbose
#> [1] TRUE
#> 
#> $seed
#> [1] 1
#>