def mean(a, rep=<NUM_LIT>, **kwargs):
return rfunc(a, ma.mean, rep, **kwargs)<EOL>
Compute the average along a 1D array like ma.mean,
    but with a representativity coefficient : if ma.count(a)/ma.size(a)>=rep,
    then the result is a masked value