From 3623d761023be5e41d522fb73f46b1bf6613bb53 Mon Sep 17 00:00:00 2001
From: "Per.Andreas.Brodtkorb"
 <Per.Andreas.Brodtkorb@96a810ac-173e-11df-b9d6-e102a254d9bb>
Date: Thu, 15 Jul 2010 11:28:47 +0000
Subject: [PATCH] Removed reference to ppimport

---
 pywafo/src/wafo/stats/estimation.py | 91 +++++++++++++++--------------
 1 file changed, 46 insertions(+), 45 deletions(-)

diff --git a/pywafo/src/wafo/stats/estimation.py b/pywafo/src/wafo/stats/estimation.py
index de0553e..81a000b 100644
--- a/pywafo/src/wafo/stats/estimation.py
+++ b/pywafo/src/wafo/stats/estimation.py
@@ -865,51 +865,52 @@ class FitDistribution(rv_frozen):
  
 
 def main():
-    _WAFODIST = ppimport('wafo.stats.distributions')
-     #nbinom(10, 0.75).rvs(3)
-    import matplotlib
-    matplotlib.interactive(True)
-    t = _WAFODIST.bernoulli(0.75).rvs(3)
-    x = np.r_[5, 10]
-    npr = np.r_[9, 9]
-    t2 = _WAFODIST.bd0(x, npr)
-    #Examples   MLE and better CI for phat.par[0]
-    R = _WAFODIST.weibull_min.rvs(1, size=100);
-    phat = _WAFODIST.weibull_min.fit(R, 1, 1, par_fix=[nan, 0, nan])
-    Lp = phat.profile(i=0)
-    Lp.plot()
-    Lp.get_CI(alpha=0.1)
-    R = 1. / 990
-    x = phat.isf(R)
-
-    # CI for x
-    Lx = phat.profile(i=0, x=x)
-    Lx.plot()
-    Lx.get_CI(alpha=0.2)
-
-    # CI for logSF=log(SF)
-    Lpr = phat.profile(i=0, logSF=log(R), link=phat.dist.link)
-    Lpr.plot()
-    Lpr.get_CI(alpha=0.075)
-
-    _WAFODIST.dlaplace.stats(0.8, loc=0)
-#    pass
-    t = _WAFODIST.planck(0.51000000000000001)
-    t.ppf(0.5)
-    t = _WAFODIST.zipf(2)
-    t.ppf(0.5)
-    import pylab as plb
-    _WAFODIST.rice.rvs(1)
-    x = plb.linspace(-5, 5)
-    y = _WAFODIST.genpareto.cdf(x, 0)
-    #plb.plot(x,y)
-    #plb.show()
-
-
-    on = ones((2, 3))
-    r = _WAFODIST.genpareto.rvs(0, size=100)
-    pht = _WAFODIST.genpareto.fit(r, 1, par_fix=[0, 0, nan])
-    lp = pht.profile()
+    pass
+#    _WAFODIST = ppimport('wafo.stats.distributions')
+#     #nbinom(10, 0.75).rvs(3)
+#    import matplotlib
+#    matplotlib.interactive(True)
+#    t = _WAFODIST.bernoulli(0.75).rvs(3)
+#    x = np.r_[5, 10]
+#    npr = np.r_[9, 9]
+#    t2 = _WAFODIST.bd0(x, npr)
+#    #Examples   MLE and better CI for phat.par[0]
+#    R = _WAFODIST.weibull_min.rvs(1, size=100);
+#    phat = _WAFODIST.weibull_min.fit(R, 1, 1, par_fix=[nan, 0, nan])
+#    Lp = phat.profile(i=0)
+#    Lp.plot()
+#    Lp.get_CI(alpha=0.1)
+#    R = 1. / 990
+#    x = phat.isf(R)
+#
+#    # CI for x
+#    Lx = phat.profile(i=0, x=x)
+#    Lx.plot()
+#    Lx.get_CI(alpha=0.2)
+#
+#    # CI for logSF=log(SF)
+#    Lpr = phat.profile(i=0, logSF=log(R), link=phat.dist.link)
+#    Lpr.plot()
+#    Lpr.get_CI(alpha=0.075)
+#
+#    _WAFODIST.dlaplace.stats(0.8, loc=0)
+##    pass
+#    t = _WAFODIST.planck(0.51000000000000001)
+#    t.ppf(0.5)
+#    t = _WAFODIST.zipf(2)
+#    t.ppf(0.5)
+#    import pylab as plb
+#    _WAFODIST.rice.rvs(1)
+#    x = plb.linspace(-5, 5)
+#    y = _WAFODIST.genpareto.cdf(x, 0)
+#    #plb.plot(x,y)
+#    #plb.show()
+#
+#
+#    on = ones((2, 3))
+#    r = _WAFODIST.genpareto.rvs(0, size=100)
+#    pht = _WAFODIST.genpareto.fit(r, 1, par_fix=[0, 0, nan])
+#    lp = pht.profile()
                 
 if __name__ == '__main__':
    main()