ggplot2でgeom_raster（）の余白プロットをプロットする

Question

私はいくつかのデータとこのデータのマージンプロットを含むラスタマップをプロットしようとしています。 種類がlevelplot()rasterVisのように機能しているか、線ではなく色で（プロット内の色を反映して）より優れています。私がしようとしていることを見て、以下で明らかになります。 （密度ラインと）rasterVis例：

しかし、私はggplot2を使用します。 data used in this example（8KB）をダウンロードしてください。

私はgrid.arrange()を使用しようとしました：

library(ggplot2) 
library(raster) 
library(maps) 
library(grid) 
library(scales) #for muted() 
library(gridExtra) 

load("example.Rsave") #load data 

#prepare data for the margin plot 
var=tp_rlim; min=0; max=1; mid=0.5 
varname <- colnames(var)[3] 
zonal <- aggregate(var[,3], list(var$Lat), mean) 
colnames(zonal) <- c("Lat", varname) 

#plot the margin plot 
ggzonal <- ggplot(zonal, aes(x="", y=Lat, fill=tp))+ 
      geom_raster() + 
      theme(legend.position="none", 
        axis.title.x = element_blank(), 
        axis.title.y = element_blank(), 
        axis.ticks = element_blank(), 
        axis.text.x = element_blank(), 
        axis.text.y = element_blank() 
        ) + 
      scale_fill_gradient2(limits=c(min, max), space="Lab", midpoint=mid) 

#plot main plot 
    gg <- ggplot(var, aes_string(x="Lon", y="Lat", fill=varname)) + 
      geom_raster() + 
      coord_quickmap(xlim = range(var$Lon), ylim = range(var$Lat)) + 
      scale_fill_gradient2(limits=c(min, max), space="Lab", midpoint=mid) + 
      theme(legend.position="bottom", 
        legend.key.width = unit(1.5, "cm"), 
        legend.key.height = unit(0.3, "cm"), 
        legend.title=element_blank() 
       ) + 
      borders(region=unique(na.omit(map.where("world", var$Lon, var$Lat))), colour="black") 

    grid.arrange(gg, ggzonal, ncol=2, widths=c(5, 1)) 

しかし、結果は満足できるものではない：

私もggExtra::ggMarginalを使用しようとしましたが、している：

ggExtra::ggMarginal(gg) 

を

または：

ggExtra::ggMarginal(gg, type="histogram") 

私はこれをどのように行うことができますか？

library(cowplot) 

grobs <- ggplotGrob(gg + theme(legend.position="bottom"))$grobs 
legend_b <- grobs[[which(sapply(grobs, function(x) x$name) == "guide-box")]] 

prow <- plot_grid(gg + theme(legend.position="none"), ggzonal, align = 'h', rel_widths = c(5, 1)) 

plot_grid(prow, legend_b, ncol = 1, rel_heights = c(1, .2)) 
ggsave('stack.png', w = 5, h = 5) 

結果：

Answer 1

これはうまくいくかもしれない、

# fix the ylim of annotation plot 
gg2 = ggzonal + coord_cartesian(ylim = range(var$Lat)) 
# devtools::install_github("baptiste/egg") 
library(egg) 
ggarrange(gg, gg2, widths = c(10,1)) 

Answer 2

素晴らしいcowplotパッケージを使用して、命令がhere見つけ、私たちは、次の操作を行うことができます