I want to calculate the Kullback-Leibler divergence
It’s implemented in Distances.jl (not the other cool package, just used there): this package (one of the metrics): https://github.com/turingtest37/SequencerJ.jl/blob/master/docs/src/index.md
and in the original Python version. See: sequencer.org