get_D_KL_highD {singleCellHaystack}

Calculates the Kullback-Leibler divergence between distributions for the high-dimensional version of haystack().

Description

Calculates the Kullback-Leibler divergence between distributions for the high-dimensional version of haystack().

Usage

get_D_KL_highD(classes, density.contributions, reference.prob, pseudo = 0)

Arguments

classes

A logical vector. Values are T is the gene is expressed in a cell, F is not.

density.contributions

A matrix of density contributions of each cell (rows) to each center point (columns).

reference.prob

A reference distribution to calculate the divergence against.

pseudo

A pseudocount, used to avoid log(0) problems.

Value

A numerical value, the Kullback-Leibler divergence