x <- matrix(rnorm(1000),ncol=10)
y <- matrix(rnorm(500),ncol=10)
normal.cortest(x) #just test if this matrix is an identity
normal.cortest(x,y) #do these two matrices differ?
mat.cortest(x)
mat.cortest(x,y) #twice the degrees of freedom as the Jennrich
jennrich.cortest(x,y) #
Run the code above in your browser using DataLab