R/simCNVdataNorm.R

In CNVassoc: Association Analysis of CNV Data and Imputed SNPs

simCNVdataNorm<-function(n,mu.surrog,sd.surrog,w,mu.y,sd.y,cnv.random=FALSE){
  w<-w/sum(w)
  k<-length(mu.surrog)
  if (cnv.random){
    nj <- rmultinom(n = 1, size = n, prob = w)
  }else{
    nj <- round(n*w)
    nj[1] <- nj[1]+n-sum(nj)
  }
  surrog<-unlist(sapply(seq(along = nj), function(j) rnorm(nj[j], mean = mu.surrog[j], sd = sd.surrog[j]),simplify=FALSE))
  cnv<-rep(1:k,nj)
  resp<-unlist(sapply(1:k,function(j) rnorm(nj[j],mu.y[j],sd.y[j]),simplify=FALSE))
  out<-data.frame(resp,cnv,surrog)
  out<-out[sample(1:n),]
  out
}