Nothing
DataSimu <-
function(){
cat("Generating covariance matrix......")
set1_matrix<-matrix(0,500,500)
for(i in 1:500){
for(j in 1:500){
if(i==j)set1_matrix[i,j]=1
}
}
set2_matrix<-matrix(0,500,500)
for(i in 1:500){
for(j in 1:500){
if(i==j)set2_matrix[i,j]=1
if(i<=50 & j<=50 & i!=j)set2_matrix[i,j]=0.6
}
}
set3_matrix<-set1_matrix
set4_matrix<-set2_matrix
set5_matrix<-matrix(0,500,500)
for(i in 1:500){
for(j in 1:500){
if(i==j)set5_matrix[i,j]=1
if(i<=50 & j<=50 & i!=j){
if((i<=25 & j<=25)|(i>25 & j>25) )set5_matrix[i,j]=0.6
else set5_matrix[i,j]=-0.6
}
}
}
cat("Simulating data......")
set1_data<-mvrnorm(n = 20, rep(0,500), Sigma=set1_matrix)
set2_data<-mvrnorm(n = 20, c(rep(0.75,50),rep(0,450)), Sigma=set2_matrix)
set3_data<-mvrnorm(n = 20, c(rep(0.75,50),rep(0,450)), Sigma=set3_matrix)
set4_data<-mvrnorm(n = 20, rep(0,500),Sigma=set4_matrix)
set5_data<-mvrnorm(n = 20, c(rep(0.75,25),rep(-0.75,25),rep(0,450)),Sigma=set5_matrix)
set6_data<-cbind(set2_data[,1:10],set3_data[,1:10],set4_data[,1:10],set5_data[,1:10],set1_data[,1:460])
control_data<-mvrnorm(n = 20, rep(0,500), Sigma=set1_matrix)
cat("Done!")
return(list(set1_data,set2_data,set3_data,set4_data,set5_data,set6_data,control_data))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.