Build index for applicable dates to year, month parameters (and consider designing this as a date_beg, date_end range)
If multiple dates, then handle raster brick with some default averaging (method = mean), with option to use a different method like sd() (standard deviation); This can be handled with simple raster algebra, that should allow for na.rm = T as an argument; eg sum(raster_brick, na.rm=T).
A follow-up enhancement to #57 would be to:
year,monthparameters (and consider designing this as adate_beg,date_endrange)method = mean), with option to use a different method likesd()(standard deviation); This can be handled with simple raster algebra, that should allow forna.rm = Tas an argument; egsum(raster_brick, na.rm=T).