python - Python 中 R 的 qqplot 的等价物

Question

我正在尝试使用 qqplots 比较 python 中的 2 个分布。在 R 中，使用 qqplot() 非常简单，如下所示:

qqplot(beaver2$temp[beaver2$activ==1],beaver2$temp[beaver2$activ==0])

这将比较 2 个发行版。在Python中是否有类似的方法？

Y=np.random.normal(0.5,0.1,1000)
Y_hat=np.random.gamma(0.5,0.2,1000)

一种选择是获取 2 组的残差，并针对正态分布绘制 qqplot

Answer 1

不确定您到底想要什么，但请检查这是否对您有帮助。您可以使用 statsmodel.api 中的 sm 来绘制 qq 绘图。

import statsmodels.api as sm
pp_x = sm.ProbPlot(Y, fit=True)
pp_y = sm.ProbPlot(Y_hat, fit=True)
fig = pp_x.qqplot(line='45', other=pp_y)
plt.show()

结果: