问题描述
我想对 ks::kde
生成的核密度估计在所有维度上进行积分(在以下示例中为四重积分)。
library(ks)
data(iris)
virginica <- iris[iris$Species == 'virginica',1:4]
lims <- apply(virginica,2,range)
virg_kde <- kde(x = virginica,gridsize = rep(20,ncol(virginica)),xmin = lims[1,],xmax = lims[2,binned = FALSE)
核密度估计在 virg_kde$estimate
中找到,是一个 4 维数组。我想计算它的体积。未导出的函数 ks:::integral.kde
似乎不适用于这种高维数据。我的微积分技能有点生疏,但我决定在 stackoverflow 上而不是 math.stackexchange 上问这个问题,因为我希望有人可以指点我在 R 中的现有函数。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)