如何使用叠加计算栅格的加权平均值？

Question

我正在尝试学习如何使用包栅格中的函数叠加。当弄乱这个功能以更好地理解它时，我遇到了这个问题。为什么这两行没有给出相同的结果？使用叠加来获得没有叠加的方法结果的正确方法是什么？

library(raster)

# Create rasters

s <- raster(nrow = 45, ncol = 90)
s[] <- 0

r <- stack( mget( rep( "s" , 12 )))
r[] <- runif(ncell(r)*12, 0, 500)

w <- stack( mget( rep( "s" , 12 )))
w[] <- runif(ncell(r)*12, 0, 1)

# Get weighted average 

# without overlay

raster::weighted.mean(r, w)

# with overlay

overlay(r, w, fun = raster::weighted.mean)

谢谢！