get_logLik: Log likelihood of clone_id model It returns P(A, D | C, I,...

In PMBio/cardelino: Clone Identification from Single Cell Data

Log likelihood of clone_id model It returns P(A, D | C, I, theta0, theta1)

Description

Log likelihood of clone_id model It returns P(A, D | C, I, theta0, theta1)

Usage

get_logLik(A1, B1, Config, Assign, theta0, theta1)

Arguments

A1	variant x cell matrix of integers; number of alternative allele reads in variant i cell j
B1	variant x cell matrix of integers; number of reference allele reads in variant i cell j
Config	variant x clone matrix of float values. The clone-variant configuration probability, averaged by posterior samples
Assign	cells x clone matrix of float values. The cell-clone assignment probability, averaged by posterior samples
theta0	the binomial rate for alternative allele from config = 0
theta1	the binomial rate for alternative allele from config = 1

Value

logLik, a float of log likelihood

Author(s)

Yuanhua Huang

PMBio/cardelino documentation built on Nov. 21, 2022, 4:52 a.m.