@inherit_docstring_from(rv_discrete_meco)<EOL><INDENT>def fit_mle(self, data, k_array=np.arange(<NUM_LIT:0.1>, <NUM_LIT:100>, <NUM_LIT:0.1>)):<DEDENT>
<EOL>data = np.array(data)<EOL>mu = np.mean(data)<EOL>return mu, _solve_k_from_mu(data, k_array, nbinom_nll, mu)<EOL>
%(super)s

        In addition to data, gives an optional keyword argument k_array
        containing the values to search for k_agg. A brute force search is then
        used to find the parameter k_agg.