Calculate the celda_G log likelihood for user-provided cluster assignments

This function calculates the log likelihood of each clustering of genes generated over multiple iterations of Gibbs sampling.

calculateLoglikFromVariables.celda_G(counts, y, L, beta, delta, gamma)

Arguments

counts	Integer matrix. Rows represent features and columns represent cells.
y	Numeric vector. Denotes feature module labels.
L	Integer. Number of feature modules.
beta	Numeric. Concentration parameter for Phi. Adds a pseudocount to each feature module in each cell. Default 1.
delta	Numeric. Concentration parameter for Psi. Adds a pseudocount to each feature in each module. Default 1.
gamma	Numeric. Concentration parameter for Eta. Adds a pseudocount to the number of features in each module. Default 1.
...	Additional parameters.

Value

The log likelihood of the provided cluster assignment, as calculated by the celda_G likelihood function

Examples

celda.sim = simulateCells(model="celda_G")
loglik = calculateLoglikFromVariables(celda.sim$counts, model="celda_G",
                                      y=celda.sim$y, L=celda.sim$L,
                                      beta=celda.sim$beta, delta=celda.sim$delta,
                                      gamma=celda.sim$gamma)