Learn R Programming

PhylogeneticEM (version 1.0.0)

find_shift_values: 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 postions 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.
m_Y
the vector of values of the means at the tips.

Value

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

Details

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