我正在尝试用Python做一些分布绘图和拟合,使用SciPy进行统计,使用matplotlib进行绘图。我在一些事情上运气不错,比如创建直方图:
seed(2)
alpha=5
loc=100
beta=22
data=ss.gamma.rvs(alpha,loc=loc,scale=beta,size=5000)
myHist = hist(data, 100, normed=True)
太棒了!
我甚至可以采用相同的伽马参数,并绘制概率分布函数的线函数(经过一些谷歌搜索):
rv = ss.gamma(5,100,22)
x = np.linspace(0,600)
h = plt.plot(x, rv.pdf(x))
如何将PDF线h
叠加在直方图的顶部来绘制直方图myHist
?我希望这是微不足道的,但我一直无法弄清楚。
3条答案
按热度按时间qvtsj1bj1#
把这两件事联系起来
要确保在任何特定的绘图示例中获得所需的内容,请先尝试创建
figure
对象suzh9iv82#
人们可能对绘制任何直方图的分布函数感兴趣。这可以通过
seaborn kde
函数来实现来自coursera的Python数据可视化课程
eh57zj3b3#
扩展Malik的答案,并试图坚持使用香草NumPy,SciPy和Matplotlib。我引入了Seaborn,但它只是用来提供更好的默认值和小的视觉调整:
给了我们以下的情节:
我试图坚持使用最小的功能集,同时产生相对较好的输出,特别是使用SciPy来估计KDE非常容易。