rnegbinmix_syn_nonsyn: Simulating synonymous and non-synonymous mutation counts

In johannabertl/SelectionMix: Negative binomial mixture model

Description

rnegbinmix_syn_nonsyn simulates counts of synonymous and non-synonymous mutations, using a negative binomial and a negative binomial mixture, respectively, with the same parameters alpha and beta.

Usage

rnegbinmix_syn_nonsyn(n, alpha, beta, c, p)

Arguments

n	number of genes
alpha	scalar, positive
beta	scalar, positive
c	numeric, positive, length k
p	numeric, probability vector, length k

Details

The function MASS::rnegbin from package MASS is used for the simulation of negative binomial data. The function rnegbinmix is used for the negative binomial mixture.

Author(s)

Johanna Bertl

Examples

# Mixture with 3 components

c1 = 0.01; c2 = 1; c3=100
p1 = 0.33; p2 = 0.33; p3 = 0.34
alpha = 10
beta = 5

simulated.mixture = rnegbinmix_syn_nonsyn(1000, alpha = alpha, beta = beta, c = c(c1, c2, c3), p = c(p1, p2, p3))

johannabertl/SelectionMix documentation built on May 3, 2019, 4:03 p.m.