Test for rotation invariant datasets

Source: R/estimateEM.R

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

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