规范核密度图的宽度和规模传奇(Standardising Kernel Density Plot&#

我有以下的4核密度图，但想传说规模以及情节宽/高是在所有4相同的比较。

我的代码是：

kde_pipit_2014_bw <- density(Pipit_ppp_2016, sigma=4.18, edge=TRUE, kernel="gaussian") 
kde_pipit_2015_bw <- density(Pipit_ppp_2016, sigma=4.18, edge=TRUE, kernel="gaussian") 
kde_pipit_2016_bw <- density(Pipit_ppp_2016, sigma=4.18, edge=TRUE, kernel="gaussian") 
kde_pipit_2016_bw <- density(Pipit_ppp_2016, sigma=4.18, edge=TRUE, kernel="gaussian")

par(mfrow=c(2,2),cex=0.7, mai=c(0.1,0.1,0.2,0.2))
plot(kde_pipit_2014_bw) 
plot(kde_pipit_2015_bw) 
plot(kde_pipit_2016_bw) 
plot(kde_pipit_2017_bw)

下面是该地块的图像。 有没有什么办法可以将它调整到相同大小的比较？

Answer 1:

您还没有提供的样本数据，所以我用mtcars 。 这是一种方式GGPLOT2的诀窍是，在注释中，从每年与包含“年份”信息的列合并的可能存在的几个dataframes，到一个数据帧，进入理想的长格式。

library(ggplot2)

ggplot(mtcars, aes(x = mpg, y = wt)) + 

  stat_density2d(aes(fill = ..density..), geom = "raster", contour = FALSE) +
  scale_fill_gradient(low="green",high="red")+
  facet_grid(~cyl)  ## you can replace this with your year

#the colors are just examples, I did not try to reproduce exactly your color scale 

文章来源: Standardising Kernel Density Plot's Width & Legend Scale