find_shift_values: Find values given edges. OU stationary case. Ultrametric...

View source: R/parsimonyNumber.R

find_shift_valuesR Documentation

Find values given edges. OU stationary case. Ultrametric tree.

Description

find_actualized_shift_values computes the values of the shifts their positions, and the mean values at the tips. Warning : this function does not check for consistency. Please make sure that the shifts positions and the mean values are compatible.

Usage

find_shift_values(shifts_edges, T_tree_ac, vect_Y, p, ntaxa)

Arguments

shifts_edges

a vector of positions of shifts on the tree.

T_tree_ac

matrix of incidence of the tree, result of function incidence.matrix.

Details

This function uses qr.solve for rectangular linear system solving.

Value

vector, with first entry the values at the root, and other entries the values of the shifts.


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