Извлечение одного графика из функции графика в GeoR
Я использую geoR
пакет, и я хотел бы отобразить только один из графиков, полученных при использовании plot
функция. Используя набор данных, поставляемый с пакетом:
library(geoR); data(elevation)
plot(elevation)
Это дает 4 графика на сетке 2 х 2, как показано ниже. Я хотел бы использовать только правый нижний график, но я не уверен, как получить этот график в одиночку.
Поэтому я попытался построить это с нуля:
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. Любые предложения будут приветствоваться.