Test the SciPy built-in function for One-Way ANOVA

oneway_anova_test.py

"""Test the SciPy built-in function for One-Way ANOVA"""
import matplotlib.pyplot as plt
import numpy as np
import scipy.stats as sstats

numargs = sstats.f.numargs
[ dfn, dfd ] = [0.9,] * numargs
rv = sstats.f(dfn, dfd)

# Display frozen pdf

plt.figure()
x = np.linspace(0, np.minimum(rv.dist.b, 3))
h = plt.plot(x, rv.pdf(x))
plt.show()

# Check accuracy of cdf and ppf

plt.figure()
prb = sstats.f.cdf(x, dfn, dfd)
h = plt.semilogy(np.abs(x - sstats.f.ppf(prb, dfn, dfd)) + 1e-20)
plt.show()