@inherit_docstring_from(rv_discrete_meco)<EOL><INDENT>def fit_mle(self, data, b=None):<DEDENT>
data = np.array(data)<EOL>length = len(data)<EOL>if not b:<EOL><INDENT>b = np.sum(data)<EOL><DEDENT>return _trunc_logser_solver(length, b), b<EOL>
%(super)s
b : float
    The upper bound of the distribution. If None, fixed at sum(data)