Добавить легенду к сюжету Hist2D в R - PullRequest
0 голосов
/ 09 мая 2019

Как добавить легенду к 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.Большое спасибо!

...