Studies in Astronomical Time Series Analysis: V. Bayesian Blocks, A New Method to Analyze Structure in

I describe a new time-domain algorithm for detecting localized structures (bursts), revealing pulse shapes, and generally characterizing intensity variations. The input is raw counting data, in any of three forms: time-tagged photon events (TTE), binned counts, or time-to-spill (TTS) data. The output is the most probable segmentation of the observation into time intervals during which the photon arrival rate is perceptibly constant, i.e., has no statistically significant variations. The idea is not that the source is deemed to have this discontinuous, piecewise constant form, rather that such an approximate and generic model is often useful. Since the analysis is based on Bayesian statistics, I call the resulting structures Bayesian blocks. Unlike most, this method does not stipulate time bins—instead the data determine a piecewise constant representation. Therefore the analysis procedure itself does not impose a lower limit to the timescale on which variability can be detected. Locations, amplitudes, and rise and decay times of pulses within a time series can be estimated independent of any pulse-shape model—but only if they do not overlap too much, as deconvolution is not incorporated. The Bayesian blocks method is demonstrated by analyzing pulse structure in BATSE γ-ray data.

[1]  D. Donoho,et al.  The quasi-periodic oscillations and very low frequency noise of Scorpius X-1 as transient chaos - A dripping handrail? , 1993 .

[2]  Emanuel Parzen,et al.  Change-point approach to data analytic wavelet thresholding , 1996, Stat. Comput..

[3]  Alan Jeffrey,et al.  Handbook of mathematical formulas and integrals , 1995 .

[4]  E. Kolaczyk.  Nonparametric Estimation of Gamma-Ray Burst Intensities Using Haar Wavelets , 1997 .

[5]  Jeffrey D. Scargle,et al.  Wavelet methods in astronomical time series analysis , 1997 .

[6]  John Skilling,et al.  Data analysis : a Bayesian tutorial , 1996 .

[7]  Walter R. Gilks,et al.  Bayesian model comparison via jump diffusions , 1995 .

[8]  Webster West,et al.  Continuous-time estimation of A change-point in a poisson process , 1997 .

[9]  Alan Jeffrey.  CHAPTER 14 – Fresnel Integrals , 1995 .

[10]  K. Worsley.  Confidence regions and tests for a change-point in a sequence of exponential family random variables , 1986 .

[11]  T. Mattfeldt.  Stochastic Geometry and Its Applications , 1996 .

[12]  E. Parzen,et al.  Data dependent wavelet thresholding in nonparametric regression with change-point applications , 1996 .

[13]  S. Chib.  Estimation and comparison of multiple change-point models , 1998 .

[14]  Gerald J. Fishman,et al.  Attributes of Pulses in Long Bright Gamma-Ray Bursts , 1996 .

[15]  A. W. Kemp,et al.  Kendall's Advanced Theory of Statistics. , 1994 .

[16]  P. Bernaola-Galván,et al.  Compositional segmentation and long-range fractal correlations in DNA sequences. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  P. Abry,et al.  Wavelets, spectrum analysis and 1/ f processes , 1995 .

[18]  R. T. Ogden,et al.  Testing change-points with linear trend , 1994 .

[19]  J. Scargle,et al.  The Dripping Handrail Model: Transient Chaos in Accretion Systems , 1996 .

[20]  D. Brillinger,et al.  Some wavelet analyses of point process data , 1997, Conference Record of the Thirty-First Asilomar Conference on Signals, Systems and Computers (Cat. No.97CB36136).

[21]  P. R. Krishnaiah,et al.  Developments in Statistics , 1979 .

[22]  C. Chui,et al.  Wavelets : theory, algorithms, and applications , 1994 .

[23]  F. Haight.  Handbook of the Poisson Distribution , 1967 .

[24]  J. Scargle,et al.  Simulated observations of gamma-ray bursts with GLAST , 1998 .

[25]  Fredrik Gustafsson,et al.  A Change Detection and Segmentation Toolbox for Matlab , 1994 .

[26]  David L. Donoho,et al.  WaveLab and Reproducible Research , 1995 .

[27]  Patrick Billingsley,et al.  Probability and Measure. , 1986 .

[28]  R. Todd Ogden.  Wavelets in the Bayesian change-point analysis , 1999, Optics & Photonics.

[29]  K OrJ.  Numerical Bayesian methods applied to signal processing , 1996 .

[30]  D. Donoho.  On Minimum Entropy Segmentation , 1994 .

[31]  D. Brillinger.  Comparative Aspects of the Study of Ordinary Time Series and of Point Processes , 1978 .

[32]  F. Gustafsson.  Segmentation of Signals Using Piecewise Constant Linear Regression Models , 1994 .

[33]  T. Loredo.  Promise of Bayesian Inference for Astrophysics , 1992 .

[34]  K. Pearson.  Biometrika , 1902, The American Naturalist.

[35]  A. Aldroubi,et al.  Wavelets in Medicine and Biology , 1997 .

[36]  William J. Fitzgerald,et al.  Multiple-order Markov chain Monte Carlo sampling methods with application to a changepoint model , 1997 .

[37]  L. Joseph,et al.  Bayesian Statistics: An Introduction , 1989 .

[38]  D. Donoho,et al.  Translation-Invariant De-Noising , 1995 .

[39]  J. K. Ord,et al.  Handbook of the Poisson Distribution , 1967 .

[40]  P. Billingsley,et al.  Probability and Measure , 1980 .

[41]  K. Riedel.  Numerical Bayesian Methods Applied to Signal Processing , 1996 .

[42]  A. F. Smith.  A Bayesian approach to inference about a change-point in a sequence of random variables , 1975 .

[43]  Gutti Jogesh Babu,et al.  Statistical Challenges in Modern Astronomy , 1992 .

[44]  A. Raftery,et al.  Bayesian analysis of a Poisson process with a change-point , 1986 .

[45]  Michel Loève,et al.  Probability Theory I , 1977 .