擅长:python、mysql、java
<p>看起来您已经正确地设置了问题;<code>rv_continuous</code>的<a href="https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.rv_continuous.html" rel="nofollow noreferrer">documentation</a>是<code>levy_stable</code>的超类,它有所有函数的链接(例如,<a href="https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.rv_continuous.fit.html" rel="nofollow noreferrer">^{<cd3>}</a>)。我的预感是,运行时非常慢是一个不稳定的bug。在</p>
<p>使用<a href="https://github.com/josemiotto/pylevy" rel="nofollow noreferrer">pylevy</a>的<a href="https://pylevy.readthedocs.io/en/latest/levy.html#levy.fit_levy" rel="nofollow noreferrer">^{<cd4>}</a>似乎可以:</p>
<pre><code>import scipy.stats as st, levy
points = 1000
jennys_constant = 8675309
alpha, beta = 1.8, -0.5
draw = st.levy_stable.rvs(alpha, beta, size=points, random_state=jennys_constant)
print(levy.fit_levy(draw))
</code></pre>
<p>结果似乎相当好(而且<code>fit_levy()</code>相当快):</p>
^{pr2}$