@inherit_docstring_from(rv_discrete_meco)<EOL><INDENT>def fit_mle(self, data, b=None):<DEDENT>
<EOL>mu = np.mean(data)<EOL>if not b:<EOL><INDENT>b = np.sum(data)<EOL><DEDENT>p = _geom_solve_p_from_mu_vect(mu, b)<EOL>if len(np.atleast_1d(p)) == <NUM_LIT:1>:<EOL><INDENT>return float(p), b<EOL><DEDENT>else:<EOL><INDENT>return p, b<EOL><DEDENT>
%(super)s
        In addition to data, requires ``b``, the upper limit of the
        distribution.