Displaying 1 result from an estimated 1 matches for "rotationtest".
2010 Dec 27
1
R-code to generate random rotation matrix for rotation testing
...ge is an option, but as far as I can tell I need to draw rotation matrices with determinant -1 as well. Roast and Romer in the limma-bioconductor package appear to have implemented something similar, which seems not to be general enough for my purposes, however. Inspired by the code in the ffmanova-rotationtest function I thought of the following, but it appears to me that there only the covariance, not the mean, is preserved:
#####
# a given Y has independent, multivariate normal rows
library(mvtnorm)
Y <- rmvnorm(4,mean=1:10,sigma=diag(1:10)+3)
# Generation of a set of random matrices Z
for (i in 1...