Предположим, у меня есть следующие x,y,z данные, которые создают трехмерный массив kde3d_ar:
x <- c(1,4,7)
y <- c(2,10,11)
z <- c(9,15,20)
kde3d_ar <- misc3d::kde3d(x, y, z, n = 3)
Формула для одномерного расхождения Кульбака-Лейблера определяется как 
, где p(t) и q(t) - одномерные распределения вероятностей. Чтобы вычислить дивергенцию Кульбака-Лейблера в трех измерениях, можно использовать тройной интеграл, где t=(x,y,z) и dt=dxdydz. Как я могу рассчитать KL-дивергенцию в трех измерениях в R?