Learn R Programming

PhylogeneticEM (version 1.0.0)

compute_expectations.scOU: Compute the expected states of a scOU

Description

compute_expectations.scOU use the matrix formulation to compute the expected values at all the nodes. Assumes a stationary root.

Usage

compute_expectations.scOU(phylo, root.state, shifts, alpha, U_tree = NULL, times_shared = NULL)

Arguments

root.state
(list): state of the root, with: random : random state (TRUE) or deterministic state (FALSE) value.root : if deterministic, value of the character at the root exp.root : if random, expectation of the character at the root var.root : if random, variance of the character at the root (pxp matrix)
shifts
(list) position and values of the shifts : edges : vector of the K id of edges where the shifts are values : matrix p x K of values of the shifts on the edges (one column = one shift) relativeTimes : vector of dimension K of relative time of the shift from the parentnode of edges
phylo:
Input tree.

Value

paramSimu: array p x nNodes x 2 (BM). For each trait t, 1 <= t="" <="p," paramsimu[t,="" ,="" ]="" has="" two="" columns,="" both="" containing="" the="" expected="" values="" for="" all="" nodes.="" dl="">