update.partitionsNumber.gen: Update formula in the general case

In pbastide/PhylogeneticEM: Automatic Shift Detection using a Phylogenetic EM

Update formula in the general case

Description

update.partitionsNumber.gen apply the actualization formula to get Nk and Ak of a node given its daughters.

Usage

update.partitionsNumber.gen(daughtersParams, ...)

Arguments

daughtersParams

: matrix with 2*npart columns. Each row contains the result of partitionsNumber for the children of a given node

Details

Uses functions sum_partitions and sum_simplex to compute the needed sums.

Value

vector of size 2*npart. For k=1:npart it contains the number of partitions with k groups compatible with the sub tree starting at the current node ; and for k=(npart+1):2*npart, it contains the number of "marked" partitions with (k-npart) groups compatible with the sub tree starting at the current node.

pbastide/PhylogeneticEM documentation built on Feb. 12, 2024, 1:27 a.m.