`plot.density`延伸＆＃34; xlim＆＃34;超出我的数据范围。为什么以及如何解决它？

Question

使用下面的代码，我试图得到不同分布的密度图。

dens <- apply(df[,c(7,9,12,14,16,18)], 2, density)

plot(NA, xlim=range(sapply(dens, "[", "x")), ylim=range(sapply(dens, "[", "y")))
mapply(lines, dens, col=1:length(dens))

legend("topright", legend=names(dens), fill=1:length(dens),bty = "n",lwd=1, cex=0.7)

所有变量的最大上限为5.但是我的行超过了5.我需要在代码中更改以修复绘图？

Answer 1

默认情况下，from会扩展范围，使密度曲线在极端值接近0。是否要将曲线限制在数据范围内？如果是这样，您需要在to内使用density()和x_range <- range(df[,c(7,9,12,14,16,18)]) dens <- apply(df[,c(7,9,12,14,16,18)], 2, density, from = x_range[1], to = x_range[2]) 个参数。

set.seed(0); X <- matrix(rnorm(300), 100, 3)

## range of your data
x_range <- range(X)
# [1] -2.904899  2.658658

## default handling of `density`
dens <- apply(X, 2, density)
range(sapply(dens, "[[", "x"))
# [1] -3.922346  3.696451

或许最好提供一个可重复的例子。

## customized `from` and `to`
dens <- apply(X, 2, density, from = x_range[1], to = x_range[2])
range(sapply(dens, "[[", "x"))
# [1] -2.904899  2.658658

try