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