R/rmpois.R

In XMRF: Markov Random Fields for High-Throughput Genetics Data

rmpois <-
function(n,lambda.v){
  tmp=c();
  # for(i in lambda.v){
  # tmp=rbind(tmp,rpois(n,i))
  # }
  ## require(multicore)
  tmp =  do.call(rbind, parallel::mclapply(lambda.v,function(i) { rpois(n,i)}))
  return(tmp)
}