Python中R的qqplot等价

Question

我尝试使用qplots比较python中的两个发行版。在R中，使用qqplot()很简单，如下所示：

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

它将比较两个分布。在Python中有类似的方法吗？在

^{pr2}$

一种选择是得到2个集合的残差，并根据Normal Distribution绘制qqplot

Answer 1

我发现statsmodels.api提供了一个函数qqplot_2samples

Y=np.random.normal(0.5,0.1,1000)
Y_hat=np.random.normal(0.5,0.2,1000)
import statsmodels.api as sm
sm.qqplot_2samples(Y,Y_hat,line='45')
plt.show()

Answer 2

不确定你到底想要什么，但看看这是否对你有帮助。你可以使用sm中的statsmodel.api来绘制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()

结果是：