def test_erfinv():
from scipy import stats as sstats
chart = Line(show_dots=False)
chart.add('scipy', [sstats.norm.ppf(x / 1000) for x in range(1, 999)])
chart.add('approx', [stats.ppf(x / 1000) for x in range(1, 999)])
# chart.add('approx', [
# special.erfinv(x/1000) - erfinv(x/1000)
# for x in range(-999, 1000)])
return chart.render_response()
def test_erfinv():
from scipy import stats as sstats
chart = Line(show_dots=False)
chart.add('scipy', [sstats.norm.ppf(x / 1000) for x in range(1, 999)])
chart.add('approx', [stats.ppf(x / 1000) for x in range(1, 999)])
# chart.add('approx', [
# special.erfinv(x/1000) - erfinv(x/1000)
# for x in range(-999, 1000)])
return chart.render_response()