compute_actualization_matrix_ultrametric: Compute Matrix W of actualization (Ultrametric case)

In PhylogeneticEM: Automatic Shift Detection using a Phylogenetic EM

Compute Matrix W of actualization (Ultrametric case)

Description

compute_actualization_matrix_ultrametric computes a squares p*Nedge bloc diagonal matrix of the (I_p - exp(-A * (h - t_pa(j))))_j node.

Usage

compute_actualization_matrix_ultrametric(
  tree,
  selection.strength,
  times_shared = compute_times_ca(tree)
)

Arguments

tree	a phylogenetic tree.
selection.strength	the selection strength of the process.
times_shared	a matrix, result of function compute_times_ca.

Details

Careful: the root is not taken into account in this function.

Value

Matrix of size p*Nedge

PhylogeneticEM documentation built on Aug. 31, 2022, 9:16 a.m.