KLDvglm function

Kullback-Leibler Divergence