r - mauchly の検定に違反した共分散行列を作成する

Question

Mauchly の検定が帰無仮説を棄却することを示すように、固有値によって共分散行列を作成する方法。rv を生成する実験計画の共分散行列のパターンはあるので、球形性の仮定が破られていますか? 与えられた固有値で任意の共分散行列を計算し、その行列に対して次のソースコードで実験を作成した場合、mauchlys テストは帰無仮説を棄却することがあります。球面性の仮定に違反する正規分布データをシミュレートする方法は?

require(corpcor)
require(reshape2)
require(mvtnorm)

getCovariance<-function(lambda){
   j<-1
   Sigma<--1
   while(any(Sigma<=0)){    
      set.seed(j)
  X<-matrix(rnorm(50*length(lambda)),50,length(lambda))  
  R<-cor(X)           
  P<-eigen(R)$vector 
  Sigma<-t(P)%*%diag(lambda)%*%P   
      j<-j+1
   }             
   Sigma
 }

 subject<-50
 treatment<-4
 set.seed(1)
 lambda<-sample(1:subject,treatment)
 cov<-getCovariance(lambda)
 exp<-rmvnorm(mean=rep(0,treatment),sigma=cov,n=subject,method="chol")
 fit<-lm(exp~1)
 mauchly.test(fit,X=~1)