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

View source: R/clone_id.R

get_logLikR Documentation

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.