slow_matrix {nodeSub}

this function calculates the p matrix within R this is slower than the C++ implementation in get_p_matrix but provides a way to debug and verify

Description

this function calculates the p matrix within R this is slower than the C++ implementation in get_p_matrix but provides a way to debug and verify

Usage

slow_matrix(eig, branch_length, rate)

Arguments

eig

eigen object

branch_length

branch length

rate

substitution rate

Value

p matrix