Function used in simulate
for OU initialization.
init.simulate.OU(phy, p, root.state, optimal.value, simulate_random, ...)
paramSimu: array p x Nnode x 3, filled with NAs. Slice paramSimu[, ntaxa + 1, ] (array p x 3) is initialized with simulated states, root expectations, and optimal values for all the traits.
Input tree.
dimension of the trait simulated
(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)