¿Cómo puedo crear un mosaico ráster utilizando la lista de rásteres?

Necesito crear varios mosaicos rasterizados. Estoy usando la versión raster 2.0-31 del paquete en una computadora con Windows de 64 bits. Creo que hice mi tarea revisando todos los blogs posibles y haciendo esta pregunta a algunos colegas, pero aún no puedo encontrar una solución.

El problema que tengo es que no puedo crear un mosaico si mis cuadrículas se enumeran en un objeto ráster. encontréeste ejemplo aunque pensé que podía aplicar, pero no, recibo un mensaje de error extraño. El siguiente ejemplo representa mi problema:

r <- raster()
r1 <- crop(r, extent(-10, 10, -10, 10))
r2 <- crop(r, extent(0, 20, 0, 20))
r3 <- crop(r, extent(10, 30, 10, 30))

r1[] <- 1:ncell(r1)
r2[] <- 1:ncell(r2)
r3[] <- 1:ncell(r3)
rasters1 <- list(r1, r2, r3)

mos <- mosaic(rasters1,fun=mean)

Este es el error que recibo:

Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘mosaic’ for signature ‘"list", "missing"’

También probé la función sugerida enaquí, pero tampoco funcionó.

fmerge <- function(rasters1, fun, ...){
  ex <- raster(union(rasters1))
  res(ex) <- res(rasters1[[1]])
  for( i in 1:length(rasters1) )
    rasters[[i]] <- merge(rasters1[[i]], ex)
  rasters <- stack(rasters1)
  fun(rasters1, ...)
}

rfm <- fmerge(rasters1, mean, na.rm=T)

Este es el mensaje de error:

Error in raster(union(rasters1)) : 
  error in evaluating the argument 'x' in selecting a method for function 'raster': Error in as.vector(y) : argument "y" is missing, with no default

Respuestas a la pregunta(1)

Su respuesta a la pregunta