#############################################################
# --- generate simulation data 
# --- sample size: 50; number of covariates 500 
# --- first 5% covariates differentially expressed (with different mean)
#############################################################

samplesize<-50; num.cov<-500;

z<-matrix(rnorm(samplesize*num.cov, mean=0, sd=1), nrow=num.cov);
z<-round(z, digits=3);

nonzero<-as.integer(0.05*num.cov);

indicator<-rep(0, samplesize);

for(i in 1:as.integer(samplesize/2)){
   indicator[i]<-1;
}

for(i in 1:nonzero){
   for(j in 1:as.integer(samplesize/2)){
      z[i, j]<-z[i, j]+1.5;
   }
}

write.table(t(z), file="covariate.txt", quote=F, append=F,
            row.name=F, col.name=F);

write.table(indicator, file="indicator.txt",
            quote=F, append=F, row.name=F, col.name=F);