Estoy trabajando con la observación diaria de datos climáticos organizados en archivos .nc. Los leí usando el comando stack del paquete ráster. Cada archivo (correspondiente a un año) es un elemento RasterStack con las siguientes características:

class       : RasterStack 
dimensions  : 360, 720, 259200, 365  (nrow, ncol, ncell, nlayers)
resolution  : 0.5, 0.5  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax) 

cada capa es la trama de valores de un día.
Me gustaría sumar las capas para calcular los valores mensuales. Creo que la solución debería usar calc o stackApply {raster}, pero no pude encontrar una manera de sumar de la capa x a la capa y o una forma de subconjugar el RasterStack antes de sumar.

Preparé unarchivo de ejemplo con solo 12 capas (para reducir el tamaño).

No sé exactamente cómo proponer un código, lo siento, pero debería ser algo como:

library(raster)
setwd("myfolder")
data<-stack(mydata.nc)

datasum<- stackApply(data, ??? ,fun=sum)

Gracias