get_D_KL function

Calculates the Kullback-Leibler divergence between distributions.

get_D_KL(classes, parameters, reference.prob, pseudo)

Arguments

• classes: A logical vector. Values are T is the gene is expressed in a cell, F is not.
• parameters: Parameters of the analysis, as set by function 'get_parameters_haystack'
• 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