Calculate the celda_C log likelihood for user-provided cluster assignments

Calculate the celda_C log likelihood for user-provided cluster assignments

calculateLoglikFromVariables.celda_C(counts, sample.label, z, K, alpha, beta)

Arguments

counts

Integer matrix. Rows represent features and columns represent cells.
sample.label

Vector or factor. Denotes the sample label for each cell (column) in the count matrix.
z

Numeric vector. Denotes cell population labels.
K

Integer. Number of cell populations.
alpha

Numeric. Concentration parameter for Theta. Adds a pseudocount to each cell population in each sample. Default 1.
beta

Numeric. Concentration parameter for Phi. Adds a pseudocount to each feature in each cell population. Default 1.
...

Additional parameters.

Examples

celda.sim = simulateCells(model="celda_C")
loglik = calculateLoglikFromVariables(celda.sim$counts, model="celda_C",
                                      sample.label=celda.sim$sample.label,
                                      z=celda.sim$z, K=celda.sim$K,
                                      alpha=celda.sim$alpha, beta=celda.sim$beta)