mcem_algorithm {DGP4LCF}

Monte Carlo Expectation Maximization (MCEM) algorithm to return the Maximum Likelihood Estimate (MLE) of DGP Parameters.

Description

This function is used to return the MLE of DGP parameters.

Usage

mcem_algorithm(
  ind_x,
  ig_parameter = 10^-2,
  increasing_rate = 0.5,
  prob_conf_interval = 0.9,
  iter_count_num = 5,
  x,
  mcem_parameter_setup_result,
  ipt_x = FALSE,
  missing_list = NULL,
  missing_num = NULL
)

Arguments

ind_x

A logical value. ind_x = TRUE uses the model including the intercept term for subject-gene mean in within-MCEM-Gibbs sampler; otherwise uses the model without the intercept term.

ig_parameter

A numeric scalar. Hyper-parameters for the prior Inverse-Gamma distribution.

increasing_rate

A numeric scalar. Rate of increasing the sample size.

prob_conf_interval

A numeric scalar. The probability that the true change in the Q-function is larger than the lower bound.

iter_count_num

A numeric scalar. Maximum number of increasing the sample size; a larger number than this would end the algorithm.

x

A list of n elements. Each element is a matrix of dimension (p, q_i), storing the gene expression observed at q_i time points for the ith subject.

mcem_parameter_setup_result

A list of objects returned from the function 'mcem_parameter_setup'.

ipt_x

A logical value. ind_x = TRUE denotes the need to impute for NAs of gene expression. The default value is ind_x = FALSE.

missing_list

A list of n elements. Each element is a matrix of dimension (missing_num, 2): each row corresponds to the position of one NA that needs imputation; first and second columns denote the row and column indexes, respectively, of the NA in the corresponding person's matrix of gene expression.

missing_num

A vector of n elements. Each element corresponds to a single person's number of NAs that needs imputation.

Value

The MLE of DGP parameters.

Examples

# See examples in vignette
vignette("bsfadgp_regular_data_example", package = "DGP4LCF")
vignette("bsfadgp_irregular_data_example", package = "DGP4LCF")