find_rotation {PhylogeneticEM}

Test for rotation invariant datasets

Description

find_rotation takes two fits from from PhyloEM, and test if their datasets are equal up to a rotation.

Usage

find_rotation(res1, res2, tol = NULL)

Arguments

res1

an object of class PhyloEM.

res2

an object of class PhyloEM.

tol

relative numerical tolerance. Default to .Machine$double.eps^(0.5).

Value

If appropriate, the rotation matrix rot such that dat1 = rot