Trying to understand distributions
Keith Sloan <[email protected]> Tue, 13 Jul 2021 05:05:13 +0100
| Newsgroups | gmane.comp.python.scientific.user |
|---|---|
| Message-ID | <[email protected]> |
If I try to follow Robert Kern's advice to us use all of the fitted parameters but without success ax3 = fig.add_subplot(4, 1, 3) ag1, bg1, cg1 =stats.gamma.fit(counts) print(ag1, bg1, cg1) x = np.linspace(stats.gamma.ppf(0.1, ag1),stats.gamma.ppf(0.99, ag1), int(ag1)) ax3.plot(x, stats.gamma.pdf(x, ag1, scale = cg1),'r-', lw=5, alpha=0.6, label='gamma pdf') ax4 = fig.add_subplot(4, 1, 4) ag2, bg2, cg2 =stats.gamma.fit(RErange1[xfield].data) print(ag2, bg2, cg2) x = np.linspace(stats.gamma.ppf(0.1, ag2),stats.gamma.ppf(0.99, ag2), int(ag2)) ax4.plot(x, stats.gamma.pdf(x, ag2, scale = cg2),'r-', lw=5, alpha=0.6, label='gamma pdf') RK -- ========== Art & Ceramics =========== https://www.instagram.com/ksloan1952/ _______________________________________________ SciPy-User mailing list [email protected] https://mail.python.org/mailman/listinfo/scipy-user