Using R and gdalUtils is there a way to apply a 3x3 low-pass filter to a virtual raster while writing it out as a .tif?

Question

I am mosaicking many small multi-band raster chunks in R using gdalUtils after running a complicated processing scheme. I am essentially looking for a way to apply a 3x3 low-pas