r - 使用 colwise 的密度图 - 每条线的颜色不同？

Question

我需要一张不同密度线的图，每条线都是另一种颜色。这是一个示例代码（但要小得多），使用内置的 data.fame USArrests。我希望它可以使用吗？

colors <- heat.colors(3)  
plot(density(USArrests[,2], bw=1, kernel="epanechnikov", na.rm=TRUE),col=colors[1])     
lines1E <- function(x)lines(density(x,bw=1,kernel="epanechnikov",na.rm=TRUE))    
lines1EUSA <- colwise(lines1E)(USArrests[,3:4])`  

目前，该代码colwise()仅使用一种颜色生成。我怎样才能让每一行都有另一种颜色？或者有没有更好的方法来绘制几条不同颜色的密度线？

Answer 1

我不太遵循您的示例，因此我创建了自己的示例数据集。首先，创建一个包含三列的矩阵：

m = matrix(rnorm(60), ncol=3)

然后绘制第一列的密度：

plot(density(m[,1]), col=2)

使用您的lines1E函数作为模板：

lines1E = function(x) {lines(density(x))}

我们可以在图中添加多条曲线：

colwise(lines1E)(as.data.frame(m[ ,2:3]))

---

就个人而言，我只会使用：

##Added in NA for illustration
m = matrix(rnorm(60), ncol=3)
m[1,] = NA
plot(density(m[,1], na.rm=T))
sapply(2:ncol(m), function(i) lines(density(m[,i], na.rm=T), col=i))

要得到：