Question

我有一个向量

v = [..., -10, -10, -10, ..., 1, 2, 5, 6, 7, 9, ...]

geom_density像密度函数一样，以平滑的方式绘制此向量的直方图！

如何使用auc（库MESS的函数，曲线下面积）计算给定时间间隔内此类矢量的密度图的曲线下面积，例如{ {1}}？

Answer 1

“ geom_density以平滑的方式绘制此向量的直方图，就像密度函数一样！” 好吧，这是因为geom_density执行了kernel density estimation！因此它不是“像密度函数” ，而是密度函数。

在geom_density的掩盖下，实际上是stats::density执行密度估计。给出内核密度估计值，以便它们定义曲线下单位面积合适的概率密度函数。

我们可以通过

确认

x <- rnorm(100)
dens <- density(x)
df <- data.frame(x = dens$x, y = dens$y)
sum(df$y) * diff(df$x)[1]
#[1] 1.000952

足够近。

通过将df中的相应值相加，可以直接在特定范围内集成密度函数；由于您不提供样本数据，因此请您自行决定。