stats_computation.py 文件源码

def median_buffer_range(mag, magerr):
    """This function returns the ratio of points that are between plus or minus 10% of the
    amplitude value over the mean

    :param mag: the time-varying intensity of the lightcurve. Must be an array.
    :param magerr: photometric error for the intensity. Must be an array.

    :rtype: float
    """

    mag, magerr = remove_bad(mag, magerr)
    n = np.float(len(mag))
    amp = amplitude(mag, magerr) 
    #mean = meanMag(mag, magerr)
    mean = np.median(mag)
    a = mean - amp/10. 
    b = mean + amp/10. 

    median_buffer_range = len(np.argwhere((mag > a) & (mag < b))) / n

    return median_buffer_range