1

我正在使用该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轴。任何建议都会受到欢迎。

4

0 回答 0