R/simCNVdataPois.R

In CNVassoc: Association Analysis of CNV Data and Imputed SNPs

simCNVdataPois<-function (n, mu.surrog, sd.surrog, w, lambda, 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) rpois(nj[j], lambda[j]), simplify = FALSE))
    out <- data.frame(resp, cnv, surrog)
    out <- out[sample(1:n), ]
    out
}