get_D_KL_continuous_highD function

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

get_D_KL_continuous_highD(
  weights,
  density.contributions,
  reference.prob,
  pseudo = 0
)

Arguments

• weights: A numerical vector with expression values of a gene.
• 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