我正在使用该geoR
包,我想只显示使用该plot
功能时获得的图表之一。使用随包提供的数据集:
library(geoR); data(elevation)
plot(elevation)
这在 2 x 2 网格上给出了 4 个图,如下所示。我想单独使用右下角的情节,但我不知道如何单独获得这个情节。
所以我尝试从头开始绘制它:
axExFact <-1.1 # to set fitting y-axis limits
ymax <- max(c(0, signif(max(elevation[[2]])*axExFact, digits=1)))
ymin <- min(c(0, signif(min(elevation[[2]])*axExFact, digits=1)))
with(elevation, hist(data, main='Plot', xlab = 'Value ',cex.lab=1.25, font.lab=2,ylim=c(ymin, ymax)))
虽然我可以调整 y 轴限制以完全覆盖数据范围,但我无法沿直方图添加密度估计。以为可以用 完成lines(density(as.numeric[[elevation]]))
,但它不起作用。
plot
因此,仅获取使用该函数获得的图形会容易得多。那么唯一的问题就是调整y轴。任何建议都会受到欢迎。