计算并绘制多个密度？

Question

我有一个包含多列的矩阵，我想计算每列的密度，然后在一个R基图中绘制这些密度。如果该图具有自动更正的比例，也会更容易。

m <- matrix(rnorm(10), 5, 10))

Answer 1

创建密度d的列表，计算xlim和ylim的值，并使用这些值创建空图。最后，在该图上绘制每个密度，并可选地绘制图例。根据要求，它仅使用基数R。

set.seed(123)
m <- matrix(rnorm(50), 5, 10) # test data

d <- apply(m, 2, density)

xlim <- range(sapply(d, "[[", "x"))
ylim <- range(sapply(d, "[[", "y"))
plot(NA, xlim = xlim, ylim = ylim, ylab = "density")

nc <- ncol(m)
cols <- rainbow(nc)
for(i in 1:nc) lines(d[[i]], col = cols[i])

legend("topright", legend = 1:nc, lty = 1, col = cols, cex = 0.7)

Answer 2

也可以使用ggplot2：

library(reshape2)
library(ggplot2)
#Data
set.seed(123)
m <- as.data.frame(matrix(rnorm(50), 5, 10))
#Melt
meltdata <- melt(m)
#Plot 1
ggplot(meltdata,aes(value,color=variable))+
  geom_density()+ggtitle('Plot 1')
#Plot 2
ggplot(meltdata,aes(value,fill=variable))+
  geom_density(alpha=0.6)+ggtitle('Plot 2')