Arithmetic operations on multi-dimensional arrays with Missings

[fabiangans]

Or using mapslices:

mapslices(xi->mean(skipmissing(xi)),x,1)

although this is probably slower than a hand-coded function.