截断分布的密度估计?

机器算法验证 r 密度函数 直方图 内核平滑
2022-03-11 15:37:22

我有一些数据在左侧明显被截断。我希望用密度估计来拟合它,以某种方式处理它,而不是试图让它平滑。

什么已知的方法(像往常一样,在 R 中)可以解决这个问题?

示例代码:

set.seed(1341)
x <- c(runif(30, 0, 0.01), rnorm(100,3))
hist(x, br = 10, freq = F)
lines(density(x), col = 3, lwd = 3)

在此处输入图像描述

谢谢 :)

2个回答

R 的 logspline 包具有 oldlogspline 函数,它将使用观察数据和审查数据的混合来估计密度。

密度函数还有一个参数来指示from要估计密度的网格的最左侧”。继续上面的例子:

lines(density(x, from = 0), col = 4, lwd = 3)

from但是,如您所见,这与没有上述参数的分布完全相同。它只是从0开始,仅此而已。