From af72b63ae519c3cffc84cacadb7300ffceaa39df Mon Sep 17 00:00:00 2001 From: shivesh Date: Sat, 13 Apr 2019 12:37:18 -0500 Subject: add discovery check --- utils/plot.py | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'utils') diff --git a/utils/plot.py b/utils/plot.py index 4f9d961..954c24b 100644 --- a/utils/plot.py +++ b/utils/plot.py @@ -1141,6 +1141,10 @@ def plot_x(data, outfile, outformat, args): scales, statistic = ma.compress_rows(r_data[idim][itex][isrc]).T print 'scales', scales print 'statistic', statistic + max_st = np.max(statistic) + if args.stat_method is StatCateg.BAYESIAN: + if (statistic[0] - max_st) > np.log(10**(BAYES_K)): + raise AssertionError('Discovered LV!') try: tck, u = splprep([scales, statistic], s=0) -- cgit v1.2.3