getNBGaussianLikelihood: Calculate the log likelihood of observed read counts

In klarman-cell-observatory/MAUDE: Mean Alterations Using Discrete Expression

Description

Uses a normal distribution (N(mu,sigma)) to estimate how many reads are expected per bin under nullModel, and calculates the log likelihood under a negative binomial model.

Usage

getNBGaussianLikelihood(x, mu, k, sigma = 1, nullModel, libFract)

Arguments

x	a vector of guide counts per bin
mu	the mean for the normal expression distribution
k	the vector of total counts per bin
sigma	for the normal expression distribution (defaults to 1)
nullModel	the bin bounds for the null model (for no change in expression)
libFract	the fraction of the unsorted library this guide comprises (e.g. from unsorted cells, or sequencing the vector)

Value

the log likelihood

Examples

#usually not used directly

klarman-cell-observatory/MAUDE documentation built on Nov. 4, 2019, 3:53 p.m.