e_cpp: function to perform one step of the E-M algorithm

In franciscorichter/emphasis: Evolutionary Modeling on Phylogenetic Applications with Simulations and Importance Sampling

function to perform one step of the E-M algorithm

Description

function to perform one step of the E-M algorithm

function to perform one step of the E-M algorithm

Usage

e_cpp(
  brts,
  init_pars,
  sample_size,
  maxN,
  soc,
  max_missing,
  max_lambda,
  lower_bound,
  upper_bound,
  xtol_rel,
  num_threads = 1
)

e_cpp(
  brts,
  init_pars,
  sample_size,
  maxN,
  soc,
  max_missing,
  max_lambda,
  lower_bound,
  upper_bound,
  xtol_rel,
  num_threads = 1
)

Arguments

brts	vector of branching times
init_pars	vector of initial parameter files
sample_size	number of samples
maxN	maximum number of failed trees
soc	number of lineages at the root/crown (1/2)
max_missing	maximum number of species missing
max_lambda	maximum speciation rate
lower_bound	vector of lower bound values for optimization, should be equal in length to the vector of init_pars
upper_bound	vector of upper bound values for optimization, should be equal in length to the vector of init_pars
xtol_rel	relative tolerance for optimization
num_threads	number of threads used.

Value

a list with the following components:

• treeslist of trees

• rejectednumber of rejected trees

• rejected_overrunsnumber of trees rejected due to too large size

• rejected_lambdanumber of trees rejected due to lambda errors

• rejected_zero_weightsnumber of trees rejected to zero weight

• time_elapsedtime used

• weightsvector of weights

• fhatvector of fhat values

• logfvector of logf values

• loggvector of logg values

a list with the following components:

• treeslist of trees

• rejectednumber of rejected trees

• rejected_overrunsnumber of trees rejected due to too large size

• rejected_lambdanumber of trees rejected due to lambda errors

• rejected_zero_weightsnumber of trees rejected to zero weight

• time_elapsedtime used

• weightsvector of weights

• fhatvector of fhat values

• logfvector of logf values

• loggvector of logg values

franciscorichter/emphasis documentation built on Feb. 19, 2024, 7:36 p.m.