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.

Value

A numerical value, the Kullback-Leibler divergence