Skip to content
Prev 2727 / 12125 Next

[R-pkg-devel] mvrnorm, eigen, tests, and R CMD check

> Your explanation needs to be a bit more general in the
    > case of identical eigenvalues - each distinct eigenvalue
    > has an associated subspace, whose dimension is the number
    > repeats of that eigenvalue and the eigenvectors for that
    > eigenvalue are an orthonormal basis for that subspace.
    > (With no repeated eigenvalues this gives your 'unique up
    > to sign'.)

Thank you, Bill, notably for the concrete example of non-trivial
eigenspaces (per eigenvector). 
Note I did say

  "... such that at least for the good cases where all eigenspaces
   are 1-dimensional, ..."

knowing well that only in that case it "is easy".
I have a gut feeling but may be wrong that such simplistic post
processing may also help (to get cross-platform reproducibility)
in the case of MASS::mvrnorm() where repeated eigenvalues will
be common in practice.

Martin