slow_matrix: 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