Get the complete log likelihood for a given celda model.

Get the complete log likelihood for a given celda model.

completeLogLikelihood(celda.mod)

Arguments

celda.mod	Celda object of class "celda_C", "celda_G", or "celda_CG".

Value

Numeric Vector. The log-likelihood of the model's cluster assignments during each iteration.

Examples

celda.sim = simulateCells("celda_CG")
celda.mod = celda_CG(celda.sim$counts, K=celda.sim$K, L=celda.sim$L,
                     nchains=1, max.iter=1)
#> --------------------------------------------------------------------
#> Starting Celda_CG: Clustering cells and genes.
#> --------------------------------------------------------------------
#> Thu Sep 06 12:56:45 2018 .. Initializing chain 1 with 'random' (seed=12345)
#> Thu Sep 06 12:56:45 2018 .... Determining if any cell clusters should be split.
#> Thu Sep 06 12:56:45 2018 .... Cluster 4 was reassigned and cluster 3 was split in two.
#> Thu Sep 06 12:56:45 2018 .... Determining if any gene clusters should be split.
#> Thu Sep 06 12:56:46 2018 .... Cluster 2 was reassigned and cluster 9 was split in two.
#> Thu Sep 06 12:56:46 2018 .... Completed iteration: 1 | logLik: -1272652.93151523
#> Thu Sep 06 12:56:46 2018 .. Finished chain 1 with seed 12345
#> --------------------------------------------------------------------
#> Completed Celda_CG. Total time: 0.6726151 secs
#> --------------------------------------------------------------------
complete.loglik = completeLogLikelihood(celda.mod)