get_D_KL_highD function

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

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.

Returns

A numerical value, the Kullback-Leibler divergence