Displaying 1 result from an estimated 1 matches for "dorum".
Did you mean:
forum
2010 Dec 27
1
R-code to generate random rotation matrix for rotation testing
...(i in 1:10) {
# R is random matrix of independent standard-normal entries
R <- matrix(rnorm(16),ncol=4)
R <- qr.Q(qr(R, LAPACK = TRUE))
# Z shall be a random matrix with the same mean and covariance structure as Y
Z <- crossprod(R,Y)
}
#####
A suggestion for the procedure exists (in Dorum et al. http://www.bepress.com/sagmb/vol8/iss1/art34/ , end of chapter 2.1), but a hint to a (fast) implementation would be greatly appreciated.
Best regards and a happy new year,
Martin Krautschke
-----------------------
Martin Krautschke
Student at University of Vienna
--
Sicherer, schnelle...