How do I prevent writeRaster() from writing bad metadata statistics (mean and stddev) in R?

Question

I am trying to output a raster which contains the combined bands from two multi-band rasters with a couple simple statements in R using the raster package. eg.: