标准化核密度图的宽度＆amp;传奇规模

Question

我有以下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的一种方式 正如评论中所提到的，诀窍是将每年可能存在的几个数据帧合并到一个数据框中，理想情况下是一个长格式，其中一列包含＆＃39;年＆＃39;信息。

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