r - R 中的不对称 fiddle 图或背靠背密度图

Question

我正在尝试绘制一个不对称的 fiddle 图来比较两个 fiddle 图(我想这基本上只是两个背靠背的核密度估计图)。有没有办法使用 R 的 Vanilla fiddle 绘图工具来做到这一点？

我知道我可以像这样创建一个密度图:

x <- rnorm(1000)
d <- density(x)
plot(d)

而且我知道 density 返回密度估计的 x 和 y 分量，但我似乎无法将这些部分拼凑起来在一起。

Answer 1

我会为此使用density，它应该可以正常工作。我不认为你通过做“ fiddle ”风格的事情得到太多，反射(reflect)在 0 线上。我个人认为相互比较要容易得多。

a <- rnorm(20)
b <- rnorm(50)


ad <- density(a)
bd <- density(b)

abd <- list(x = c(ad$x, bd$x),
            y = c(ad$y, bd$y))

# "violin" style comparison (psuedo-mirrored), switch x and y to make vertical.
plot(range(abd$x), c(-max(abd$y), max(abd$y)), type = "n")
lines(ad$x, ad$y, type = "l")
lines(bd$x, -bd$y, type = "l")
abline(h = 0)

# on top of each other comparison, would nicely generalize for more distributions
plot(range(abd$x), range(abd$y), type = "n")
lines(ad)
lines(bd)