getReadMatrix.NB {HeritSeq} | R Documentation |
Simulate a count matrix from negative binomial mixed effect models (NBMM).
Description
Simulate a (possibly unbalanced) count matrix from NBMM.
Under NBMM, an observed number of reads aligned to feature/gene ,
, follows a negative binomial (NB) distribution with mean
and variance
, where
is the dispersion parameter, shared across strains. The
generalized linear model uses a
-link:
Usage
getReadMatrix.NB(vec.num.rep, alphas, sigma2s, phis)
Arguments
vec.num.rep |
A vector of replicate numbers for each strain. |
alphas |
Intercept vector |
sigma2s |
Random effect variance vector |
phis |
Dispersion parameter in NB models, |
Value
A matrix with NB reads.
is the
total number of samples;
is the number of features. Column names
are sample names of the form "Ss_r", where S stands for sample, s is the
strain number, r is the replicate number within the strain. Row names
are the feature names of the form "Gene g", where g is the feature index.
Examples
## Generate a sequencing dataset with 5 features and 6 strains.
## Assign parameter values.
rep.num <- c(3, 5, 2, 3, 4, 2)
a0s <- c(-1, 1, 2, 5, 10)
sig2s <- c(10, 0.2, 0.1, 0.03, 0.01)
phis <- c(0.5, 1, 0.05, 0.01, 0.1)
set.seed(1234)
## Generate reads:
nbData <- getReadMatrix.NB(rep.num, a0s, sig2s, phis)