R中的不对称小提琴图或背靠背密度图

Question

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

我知道我可以像这样创建一个单一的密度图：

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

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

Answer 1

我会为此使用density，它应该可以正常工作。我不认为你通过做“小提琴”风格的事情获得太多收益，反映在 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)