擅长:python、mysql、java
<p>使用scipy怎么样?您可以从<a href="http://docs.scipy.org/doc/scipy/reference/stats.html#continuous-distributions" rel="noreferrer">continuous distributions in the scipy.stats library</a>中选择所需的分发。</p>
<p>广义gamma函数有非零的歪斜和峰度,但是您需要做一点工作来确定要使用哪些参数来指定分布以获得特定的均值、方差、歪斜和峰度。这里有一些代码让你开始。</p>
<pre><code>import scipy.stats
import matplotlib.pyplot as plt
distribution = scipy.stats.norm(loc=100,scale=5)
sample = distribution.rvs(size=10000)
plt.hist(sample)
plt.show()
print distribution.stats('mvsk')
</code></pre>
<p>这将显示正态分布中10000个元素样本的直方图,均值为100,方差为25,并打印分布的统计信息:</p>
<p><code>(array(100.0), array(25.0), array(0.0), array(0.0))</code></p>
<p>用广义伽马分布代替正态分布</p>
<pre><code>distribution = scipy.stats.gengamma(100, 70, loc=50, scale=10)
</code></pre>
<p>你得到统计数据[均值,方差,偏斜,峰度]
<code>(array(60.67925117494595), array(0.00023388203873597746), array(-0.09588807605341435), array(-0.028177799805207737))</code>。</p>