r - 绘制 `density` 对象时出现问题

Question

我在尝试绘制 density 对象时遇到问题。考虑例如

require(grDevices)

set.seed(43)
d0 = density(rexp(1e5,rate=1))
d1 = density(rexp(1e5,rate=1.8))


plot(d1, col="white", xlab = "x", ylab="Density", main = "")
polygon(d1, col=adjustcolor("blue", alpha.f=0.2))
lines(d0, col="white")
polygon(d0, col=adjustcolor("red", alpha.f=0.2))

目前，它看起来符合我的预期。缩放 Y 轴的低值时会出现问题。考虑例如

plot(d1, col="white", xlim=c(2.5,3), xlab = "x", ylab="Density", main = "", 
ylim=c(0,0.02))
polygon(d1, col=adjustcolor("blue", alpha.f=0.2))
lines(d0, col="white", xlim=c(2.5,3), ylim=c(0,0.02))
polygon(d0, col=adjustcolor("red", alpha.f=0.2))

奇怪的是，多边形的下部没有达到密度 = 0。而且，一个多边形的末端低于另一个。设置 yaxs="i" 和 xaxs="i" 时问题仍然存在。

这是怎么回事，如何解决这个问题？

---

有了一些个人数据，我得到了类似的东西

Answer 1

另一种选择是将密度估计的 x 范围扩展到数据范围之外一点，这样密度估计实际上在两端都为零。这避免了人为更改密度估计中的任何值的需要。例如:

set.seed(43)
d0 = density(rexp(1e5,rate=1))
d1 = density(rexp(1e5,rate=1.8))

d1$y[c(1, length(d1$y))]

[1] 2.987316e-03 1.235864e-06

set.seed(43)
d1 = rexp(1e5,rate=1.8)
d1 = density(d1, from=min(d1) - 0.05*diff(range(d1)), to=max(d1) + 0.05*diff(range(d1)))

d1$y[c(1, length(d1$y))]

[1] 6.334144e-17 3.797333e-17