我有一些数据在左侧明显被截断。我希望用密度估计来拟合它,以某种方式处理它,而不是试图让它平滑。
什么已知的方法(像往常一样,在 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)

谢谢 :)
我有一些数据在左侧明显被截断。我希望用密度估计来拟合它,以某种方式处理它,而不是试图让它平滑。
什么已知的方法(像往常一样,在 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)

谢谢 :)
R 的 logspline 包具有 oldlogspline 函数,它将使用观察数据和审查数据的混合来估计密度。
密度函数还有一个参数来指示“from要估计密度的网格的最左侧”。继续上面的例子:
lines(density(x, from = 0), col = 4, lwd = 3)
from但是,如您所见,这与没有上述参数的分布完全相同。它只是从0开始,仅此而已。