extract_limits

zdm.analyze_cube.extract_limits(x, y, p, method=1, dumb=False, interp=False)

Parameters:

• x (np.ndarray) – xvalues of data (independent variable on which we place limits)

• y (np.ndarray) – yvalues of data (probability distribution)

• p – probability at each tail of the distribution (e.g. 0.16 for ~1 sigma, 0.025 for 2 sigma) Use “dumb = True” to interpet this as the C.I. (e.g. 0.68 for 1 sigma, 0.95 for 2 sigma)

• p – probability of confidence interval (e.g. (1-0.68)/2 for ~1 sigma)

• method – 1: includes most-likely points in C.I. first 2: sets equal probability to either side

• interp – if True, interpolates between points

Returns:

lower and upper bounds of range

Return type:

v0,v1