l_gpois: calculate joint likelihood of a gamma-poisson model
Description
calculate joint likelihood of a gamma-poisson model
Usage
l_gpois(Y_obs, lambda_ref, d, alpha, beta, phi = 1)
Value
numeric Joint log likelihood
Arguments
- Y_obs
numeric vector Gene expression counts
- lambda_ref
numeric vector Reference expression levels
- d
numeric Total library size
- alpha
numeric Alpha parameter of Gamma-Poisson distribution
- beta
numeric Beta parameter of Gamma-Poisson distribution