Skip to contents

model_selection Function

model_selection_h.Rd

Perform model selection among the models fit with varying number of mixture components (number of clocks).

Usage

model_selection_h(
  data,
  draws_and_summary,
  log_lik_matrix_list,
  elbo_iterations,
  n_components = 0
)

Arguments

data

list of k_max of lists: $input_data:list: List of 7: $S: int, $N: int, $karyotype: num (0 or 1), $seg_assignment: num, $peaks:List of N of num (1:2), $NV: num, $DP: num

draws_and_summary

list of lenght k_max of draws form the variational method with the summary statistics of the draws from the approximate posterior

log_lik_matrix_list

list of lenght k_max

elbo_iterations

list of lenght k_max

n_components

number of components specified from user

Value

result_model_selection: list $best_fit, $best_K, $model_selection_tibble, $entropy_list)