R/rmix.R

rmix <-
function (pars, n) 
{
    n1 = round(n * pars[1])
    n2 = n - n1
    x1 = rgamma(n1, shape = pars[2], rate = pars[3])
    x2 = rnorm(n2, mean = pars[4], sd = pars[5])
    return(c(x1, x2))
}
Vivianstats/scImpute documentation built on Aug. 21, 2019, 1:49 p.m.