Q-Q plot with probabilities

docs/source/plots/graphics_gofplots_qqplot.py

@@ -3,7 +3,7 @@
 Created on Sun May 06 05:32:15 2012
 
 Author: Josef Perktold
-
+editted by: Paul Hobson (2012-08-19)
 """
 from scipy import stats
 from matplotlib import pyplot as plt
@@ -49,3 +49,22 @@
 fig.tight_layout()
 
 plt.gcf()
+
+
+# example with the new ProbPlot class
+import numpy as np
+x = np.random.normal(loc=8.25, scale=3.5, size=37)
+y = np.random.normal(loc=8.00, scale=3.25, size=37)
+pp_x = sm.ProbPlot(x, fit=True)
+pp_y = sm.ProbPlot(y, fit=True)
+
+# probability of exceedance
+fig2 = pp_x.probplot(exceed=True)
+
+# compare x quantiles to y quantiles
+fig3 = pp_x.qqplot(other=pp_y, line='45')
+
+# same as above with probabilities/percentiles
+fig4 = pp_x.ppplot(other=pp_y, line='45')
+
+

statsmodels/api.py

@@ -13,7 +13,7 @@
 from __init__ import test
 from . import version
 from info import __doc__
-from graphics.gofplots import qqplot
+from graphics.gofplots import qqplot, qqplot_2samples, qqline, ProbPlot
 from .graphics import api as graphics
 from .stats import api as stats
 
@@ -25,9 +25,5 @@ def open_help(chmpath=chmpath):
         from subprocess import Popen
         p = Popen(chmpath, shell=True)
 
-
 del os
 del chmpath
-
-
-