def mean(a, axis=None):
    if a.mask is numpy.ma.nomask:
        return super(numpy.ma.MaskedArray, a).mean(axis=axis)
counts = numpy.logical_not(a.mask).sum(axis=axis)
if counts.shape:
    pass
`numpy.mean` используется с кортежем в качестве аргумента оси: не работает с замаскированным массивом