Question

我试图绘制一个不对称的小提琴图来比较两个小提琴图（我猜这基本上只是两个核密度估计图背靠背）。有没有办法用R？

的香草小提琴绘图工具做到这一点

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

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

我知道density会返回密度估算的x和y成分，但我似乎无法将各个部分放在一起。

Answer 1

我只会使用density这样做，它应该可以正常工作。我不认为你通过做事来获得很大的收获和小提琴＆＃34;风格，镜像横跨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)

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

1 个答案: