Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations
read the original abstract
This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it - an improved and generalized version of Bayesian Blocks (Scargle 1998) - that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piecewise linear and piecewise exponential representations, multi-variate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by (Arias-Castro, Donoho and Huo 2003). In the spirit of Reproducible Research (Donoho et al. 2008) all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.
This paper has not been read by Pith yet.
Forward citations
Cited by 3 Pith papers
-
EP260321a/SN 2026gzf: The Faintest Shock Breakout Associated with a Broad-Lined Supernova
EP260321a is the faintest observed shock breakout tied to a broad-lined Type Ic supernova, interpreted as a choked weak outflow from a stripped star.
-
GRB 250706B/C: Insight-HXMT Discovery of a High-Luminosity Burst as a Candidate for Fallback-Regulated Accretion in the Prompt Emission
GRB 250706B/C exhibits temporal features consistent with fallback-regulated accretion operating on a high-luminosity branch in a collapsar.
-
Temporal evolution of the periodic GeV signal from 4FGL J1913.2+0512 and analysis of the SS 433 / W50 lobes
Fermi LAT data analysis detects a time-evolving ~162-day periodic GeV signal from 4FGL J1913.2+0512 linked to SS 433, prominent in 2008-2018 but diminished afterward, indicating multi-year evolution in gamma-ray produ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.