Как добавить легенду к 2-мерной гистограмме, созданной с помощью hist2d
?
Я строю график плотности двух наложенных параметров (температуры и осадков) климатологических данных (rdata).Файл доступен здесь .
Вот мой код:
# load libraries
library(raster)
library(gplots)
library(fields)
# load data
load("clim_data.rdata")
# assign colors
col <- (colorRampPalette(c("blue", "green", "yellow", "red")))(8)
# plot
par(mar=c(5,5,2,15), xpd=T)
hist2d(clim_data$Prec, clim_data$Temp,
nbins=250,
col=col,
FUN=function(x) log(length(x)),
main="Title",
xlab="Precipitation [mm]",
ylab="Temperature [°C]",
xlim=c(0, 8000), ylim=c(-40, 40), legend=T)
Это сделает сюжет (без легенды) и приведет к этой ошибке:
"legend" is not a graphical parameter
Я не нашелчто-нибудь о легенде в hist2d
RDocumentation .Но все равно попробовал, потому что согласно 2-му ответу этого поста генерация легенды для сюжета исторических2 должна быть возможной. представленный результат - это именно то, что я ищу.
Как я могу добавить подобную цветовую полосу, указывающую плотность на моем графике?Мне известны варианты, использующие ggplot2
или hexbinplot
, но я бы предпочел базовый сюжет.Я также знаю о colorbar.plot
, add.color.bar
, color.bar
и legend
, но изо всех сил пытаюсь найти значения плотности, сгенерированные hist2d
.Большое спасибо!