r - 单个图形上的核密度图

Question

我一直在尝试使用 R 绘制简单的密度图:

plot(density(Data$X1),col="red")
plot(density(Data$X2),col="green")

因为我想比较，所以我想把两者都画在一个图中。但是'matplot'不起作用!我还尝试使用 ggplot2 作为:

library(ggplot2)
qplot(Data$X1, geom="density")
qplot(Data$X2, add = TRUE, geom="density")

同样在这种情况下，图单独出现(尽管我写了 add=TRUE)!!请问有人能想出一个简单的解决方案吗？

Answer 1

在 ggplot2 或 lattice 中，您需要 reshape 数据以重新组合它们。

例如:

    dat <- data.frame(X1= rnorm(100),X2=rbeta(100,1,1))
    library(reshape2)
    dat.m <- melt(dat)

使用``格子`

  densityplot(~value , groups = variable, data=dat.m,auto.key = T)

使用``ggplot2`

 ggplot(data=dat.m)+geom_density(aes(x=value, color=variable))

编辑 添加 X1+X2

使用 lattice 和扩展的公式接口(interface)，做到这一点非常容易:

densityplot(~X1+X2+I(X1+X2) , data=dat)   ## no need to reshape data!!