get_D_KL_continuous_highD {singleCellHaystack}

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

Description

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

Usage

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.

Value

A numerical value, the Kullback-Leibler divergence