def median_buffer_range2(mag, magerr):
"""This function returns the ratio of points that are more than 20% 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/5.
median_buffer_range = len(np.argwhere((mag < a))) / n
return median_buffer_range
评论列表
文章目录