Learn R Programming

PhylogeneticEM (version 1.8.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 positions and the mean values are compatible.

Usage

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

Value

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

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.