REVIEW 4 major objections 5 minor 53 references
A Bayesian spike-and-slab prior over Fourier coefficients can select true periodicities from dense frequency grids while remaining sparse enough for wearable-device time series.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
A Bayesian spike-and-slab frequency-selection method with a refined grid and hierarchical multivariate extension recovers oscillatory components from wearable time series.
T0 review reviewed 2026-08-01 challenge →
load-bearing objection Useful Bayesian frequency-selection method with a solid simulation study and a real multivariate twist, but the MCMC acceptance probability as written is incomplete and the main comparison is unfairly stacked. the 4 major comments →
Frequency Selection in Bayesian Spectral Modeling of Time Series Data with Applications to Wearable Device Measurements
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The central discovery is that frequency selection can be formulated as a sparse regression problem over a dense grid of candidate frequencies and solved with a structured spike-and-slab prior. Each candidate frequency has two coefficients (sine and cosine) that are included or excluded jointly through a latent indicator; a minimum-separation constraint embedded in the prior forbids nearby active frequencies, which keeps the selected set parsimonious and identifiable. A Metropolis-within-Gibbs sampler with add/delete/swap moves explores the posterior, and posterior inclusion probabilities quantify each frequency's relevance. In simulations the method identifies the true frequencies and their
What carries the argument
The load-bearing object is the bivariate spike-and-slab prior on each frequency's sine-cosine coefficient pair, together with the minimum-separation constraint that gives zero prior mass to configurations with two active frequencies closer than d grid bins. The refined frequency grid (spacing smaller than the standard Fourier spacing) is what buys resolution beyond Fourier bins. The sampler couples these pieces: add/delete/swap proposals are restricted to valid configurations, coefficients are drawn from the conjugate Gaussian conditional, and the multivariate version places a Dirichlet prior over joint inclusion patterns so that shared rhythms are encouraged only when the data support them.
Load-bearing premise
The stochastic-search MCMC is assumed to converge to the correct constrained posterior, but the paper asserts irreducibility without proof and the acceptance ratio in Eq. (5) is written without explicit proposal or Jacobian terms; if the chain does not mix over the dense candidate grid, every posterior inclusion probability and power estimate is unreliable.
What would settle it
Take a short time series with a small number of candidate frequencies, enumerate the exact posterior over inclusion configurations (small enough to compute directly), then compare it with the empirical inclusion frequencies from the MCMC sampler; any systematic discrepancy would indicate the sampler is not targeting the stated posterior. Alternatively, simulate two frequencies separated by less than d times the grid spacing and check whether the procedure ever places posterior mass on both, revealing how the separation constraint rather than the data governs resolution.
If this is right
- If the central claim holds, researchers can replace manual periodogram peak-picking with automatically selected frequencies that come with posterior probabilities of inclusion.
- A dense grid rather than the standard Fourier grid becomes a practical default, since the paper shows the extra computational cost buys better frequency localization and power estimates.
- Joint modeling of multiple wearable signals can separate rhythms that are genuinely shared across modalities from those specific to one signal, useful for connecting thermoregulation, activity, and other physiological rhythms.
- The method offers a parsimony check on penalized alternatives: in the real actigraphy data, LASSO and Elastic Net selected hundreds of frequencies while spectralSS selected a handful of interpretable components.
- Frequency-specific inclusion probabilities and the posterior distribution over model size give a principled, interpretable measure of uncertainty about which periodicities matter.
Where Pith is reading between the lines
- The minimum-separation distance d acts as a resolution parameter; treating it as random would require handling a combinatorial normalizing constant, so a systematic sensitivity analysis over d, as the paper's simulations partially provide, is the practical way to learn how many distinct rhythms the data can actually support.
- Because the multivariate sampler evaluates likelihood ratios only for affected components, the machinery could scale beyond two channels if the joint inclusion-pattern prior were compressed; as written, the 2^D configuration space is the main bottleneck.
- The line-spectrum model is a deliberate simplification, so on broad-band or non-stationary signals the method will interpret spectral mass as discrete peaks; it may be best used as a screening tool before smooth-spectrum or time-varying spectral analysis.
- The real-data results generate testable clinical hypotheses, such as whether ultradian components near 8 or 12 hours differ systematically across epilepsy patients or track seizure frequency; these could be checked in larger cohorts.
- The paper's mechanism for selecting frequencies, if valid, could be adapted to other settings where a sparse set of sinusoidal components is sought, such as detecting oscillatory structure in neural, financial, or environmental time series.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes spectralSS, a Bayesian spike-and-slab method for selecting oscillatory components in time series. The model places a spike-and-slab prior on paired sine/cosine Fourier coefficients, enforces a minimum-separation constraint in the prior, and uses a refined grid of candidate frequencies to improve resolution. Inference is carried out with a Metropolis-within-Gibbs sampler using add/delete/swap moves, and posterior inclusion probabilities are used to select frequencies. The framework is extended to multivariate series through a Dirichlet prior over joint inclusion patterns. The paper reports extensive simulations, comparisons with Lasso/EN/SCAD/MCP, and two wearable-data applications. The central claim is that spectralSS recovers the true set of frequencies more accurately and parsimoniously than penalized alternatives and that the multivariate extension identifies shared and component-specific rhythms.
Significance. The modeling framework is attractive and likely useful for chronobiology and wearable-device data. It provides a transparent generative model, explicit uncertainty quantification for selected frequencies, and a flexible multivariate extension. The simulation suite is broad — including closely spaced frequencies, model misspecification, heavy-tailed errors, and correlated bivariate errors — and the authors provide code. However, the manuscript's central empirical claims are currently not fully supported: the only explicit Metropolis–Hastings acceptance formula appears incomplete, the headline comparison confounds the selection mechanism with grid resolution, and the frequency-error metrics are aggregate sums that can mask large per-frequency errors. If these issues are addressed, the paper would be a useful contribution to Bayesian spectral analysis and digital-health applications.
major comments (4)
- [Section 3, Eq. (5)] The acceptance probability for the add move is incomplete. β_j^prop is said to be sampled from the conjugate Normal posterior in Eq. (6), but Eq. (5) contains no proposal density q(β_j^prop | z_prop, …) and no ratio of add/delete move probabilities. A valid Metropolis–Hastings ratio for this transdimensional move must include these terms; the denominator should also contain p(β_curr | z_curr) if the joint posterior in Eq. (4) is the target. The assertion that irreducibility implies correctness is not sufficient: irreducibility alone does not imply that the chain has the intended stationary distribution unless detailed balance (or an equivalent stationarity condition) holds. Since every posterior inclusion probability, power estimate, and model-size posterior in Sections 4 and 5 is averaged over this chain, this is a load-bearing issue. The same omission appears in the multivariate accept
- [Sections 4.1.1 and 4.1.2] The comparison between spectralSS and the penalized baselines is not on an equal footing. The paper states that spectralLASSO/EN/SCAD/MCP 'operate on the same set of Fourier frequencies,' while spectralSS is evaluated on a refined grid with Δω = 0.0001 and L = 5000 candidates. In the simulation of Section 4.1, the true frequencies 1/67, 1/21, and 1/13 are not exactly on the 512-point Fourier grid, so a large portion of the reported AE_F advantage for spectralSS is simply the benefit of being able to place a candidate at (or very near) the true frequency. To support the 'superior performance' claim, either give all methods the same refined grid, or additionally report spectralSS on the Fourier grid, or compare using a nearest-grid baseline. Otherwise the comparison conflates grid resolution with the statistical selection mechanism.
- [Section 4.1.1, AE_F and AE_P definitions] The two error metrics used in the headline comparisons are aggregate sums: AE_F = |∑_j ω_j − ∑_j ω̂_j| and AE_P = |∑_j P(ω_j) − ∑_j P̂(ω_j)|. These are not per-frequency matching errors. A model that selects a completely wrong set of frequencies can have a small AE_F if the selected frequencies straddle the true ones so that the sums nearly cancel. The metrics also do not define how to handle unequal model sizes. This affects Figures 2 and S1 and the corresponding claims about 'accuracy in frequency estimation.' Please replace these with matching-based metrics (e.g., nearest-neighbor assignment, Hungarian matching, or a comparable F-measure) and report per-frequency power errors under the same matching.
- [Sections 4.1.1 and 5.1.3] The penalized baselines are implemented with glmnet and ncvreg, but the paper does not state how the tuning parameter λ is selected. The number of frequencies selected by Lasso, EN, SCAD, and MCP depends critically on λ; without a reproducible selection rule (cross-validation, BIC, oracle, or a path summary), the claims that EN/LASSO select hundreds of frequencies and are less parsimonious than spectralSS are not reproducible. The real-data claim that spectralLASSO selects on average 322 frequencies and spectralEN 353 frequencies across six subjects needs the same clarification. Please specify the tuning procedure and provide sensitivity of the sparsity comparison to it.
minor comments (5)
- [Section 2.2, Eqs. (2)–(3)] The notation M_max is imported from the Fourier-frequency setup of Section 2.1 but is then used for the dense grid of Section 2.1.1; L is not defined in Eq. (3). Please define the candidate set once and use consistent notation.
- [Section 4.1] Typo 'PPI ¿ 0.5' should be 'PPI > 0.5'.
- [Figure 3] The labels '(c) PPI (dim 2)' and '(d) PPI (dim 2)' appear to be duplicated; one panel should refer to dimension 1 and the other to dimension 2.
- [Section 5.1.2] The text says 'we selected a grid of Fourier frequencies' but the method in Section 2.1.1 uses a refined grid; please clarify whether the real-data grid is the standard Fourier grid or a finer grid, and report its spacing.
- [Section 3.2] Minor wording: 'a more stricter rule' should be 'a more stringent rule.'
Circularity Check
No significant circularity: posterior frequency selection is derived from an explicit likelihood and priors; the minor self-citation for a separation heuristic is not load-bearing.
full rationale
The paper's central outputs—posterior inclusion probabilities, selected frequencies, and power estimates—are derived from the stated saturated Fourier regression likelihood (Eq. 1), the spike-and-slab prior (Eq. 2), the constrained inclusion prior (Eq. 3), and stated hyperparameters (Beta, inverse-gamma, slab variance). None of these inputs contains the target frequencies or the true model size; in simulations, the true frequencies are used only to score outputs, not to set priors. The one self-citation, "A similar principle is discussed by Hadj-Amar et al. (2020)" (Section 2.2), merely notes a precedent for the minimum-separation heuristic; the paper independently specifies d and studies its effect (Section 2.2.1 and Supplementary S3), so the citation is not load-bearing. The refined-grid claim is tested directly against the Fourier-grid variant and the method is benchmarked against external penalized competitors (spectralLASSO, spectralEN, spectralSCAD, spectralMCP) implemented with glmnet and ncvreg, so the empirical claims are self-contained relative to external benchmarks. The skeptical concern about Eq. (5)—that the written acceptance ratio omits proposal-density and reverse-move terms, and that the text asserts "the Markov chain is irreducible over the support of the posterior. Hence, the sampler is valid and converges to the correct target distribution" (Section 3)—is a real methodological/completeness issue, but it is not circular: the claimed target distribution is not constructed from the sampler's outputs, and no fitted parameter is renamed as a prediction. The paper also openly discloses its own limitations (fixed d and exponential 2^D scaling in Section 6), further indicating that no load-bearing claim is concealed as an input. No equation reduces to its own output by construction, and no central result is forced by a self-citation chain.
Axiom & Free-Parameter Ledger
free parameters (5)
- minimum separation d =
d = 5 in applications; d = 3 in illustrative simulation
- grid spacing Δω =
Δω = 0.0001 (L = 5000) in simulations; Fourier-resolution grid in applications
- slab variance σβ² =
σβ² = 10
- Beta sparsity hyperparameters a,b =
a=1, b=10 in illustrative simulation; b=10,000 in real-data applications
- Dirichlet hyperparameter α for inclusion patterns =
α = (10,3,3,3) in bivariate simulation and application
axioms (4)
- domain assumption Time series are adequately represented as a finite sum of sinusoids at candidate frequencies plus independent Gaussian noise (Eq. 1 and multivariate analogue).
- ad hoc to paper The stochastic search sampler converges to the posterior; proposal distributions/ratios are correct and the chain is irreducible over the constrained support.
- domain assumption Refined grid basis functions plus minimum separation d yield an identifiable and well-conditioned design.
- standard math Fourier decomposition of a finite series is exact at Fourier frequencies.
Cite this review
Pith. "Pith review of Frequency Selection in Bayesian Spectral Modeling of Time Series Data with Applications to Wearable Device Measurements." pith.science (2026). https://pith.science/paper/JKODXOKH
@misc{pith2026260715157,
author = {Pith},
title = {Pith review of: Frequency Selection in Bayesian Spectral Modeling of Time Series Data with Applications to Wearable Device Measurements},
year = {2026},
howpublished = {\url{https://pith.science/paper/JKODXOKH}},
note = {Machine review of arXiv:2607.15157}
}
read the original abstract
This paper introduces a Bayesian spike-and-slab framework for spectral analysis of time series data. The proposed method combines frequency selection and dimensionality reduction with a refined grid of candidate frequencies, enabling high-resolution recovery of oscillatory components while promoting sparsity through a structured spike-and-slab prior. A stochastic search algorithm efficiently explores the posterior space, yielding posterior inclusion probabilities that quantify the relevance of each frequency. We extend the framework to multivariate signals via a hierarchical prior on frequency inclusion patterns, allowing the model to capture both shared and component-specific rhythms across multiple time series. Extensive simulation studies demonstrate the method's robustness and superior performance in frequency estimation and spectral power reconstruction compared to existing approaches. Applied to actigraphy data from individuals with partial-onset seizures, the univariate model identifies clinically relevant circadian and ultradian rhythms. In a second application, for the joint analysis of physical activity and skin temperature from a healthy individual, the multivariate model reveals partially overlapping rhythmic components consistent with known physiological coupling. This work establishes a powerful and interpretable approach to spectral analysis, with broad applicability to wearable data, chronobiology, and personalized health monitoring.
Figures
Reference graph
Works this paper leans on
-
[1]
Handbook of Bayesian Variable Selection , pages=
Discrete Spike-and-Slab Priors: Models and Computational Aspects , author=. Handbook of Bayesian Variable Selection , pages=. 2021 , publisher=
2021
-
[2]
and Vannucci, M
Brown, P.J. and Vannucci, M. and Fearn, T. , journal=. Multivariate
-
[3]
2000 , publisher=
Time series analysis and its applications , author=. 2000 , publisher=
2000
-
[4]
1991 , publisher=
Time series: theory and methods , author=. 1991 , publisher=
1991
-
[5]
Variable selection for nonparametric
Savitsky, Terrance and Vannucci, Marina and Sha, Naijun , journal=. Variable selection for nonparametric. 2011 , publisher=
2011
-
[6]
Annals of Statistics
Optimal predictive model selection , author=. Annals of Statistics
-
[7]
Scientific reports , volume=
Ultradian rhythms in heart rate variability and distal body temperature anticipate onset of the luteinizing hormone surge , author=. Scientific reports , volume=. 2020 , publisher=
2020
-
[8]
Computational statistics & data analysis , volume=
Brain waves analysis via a non-parametric Bayesian mixture of autoregressive kernels , author=. Computational statistics & data analysis , volume=. 2022 , publisher=
2022
-
[9]
Journal of the American Statistical Association , volume=
Frequency detection and change point estimation for time series of complex oscillation , author=. Journal of the American Statistical Association , volume=. 2024 , publisher=
2024
-
[10]
International Statistical Review/Revue Internationale de Statistique , pages=
Bayesian graphical models for discrete data , author=. International Statistical Review/Revue Internationale de Statistique , pages=. 1995 , publisher=
1995
-
[11]
2007 , publisher=
Data analysis using regression and multilevel/hierarchical models , author=. 2007 , publisher=
2007
-
[12]
Journal of The Royal Society Interface , volume=
Hidden Markov models for monitoring circadian rhythmicity in telemetric activity data , author=. Journal of The Royal Society Interface , volume=. 2018 , publisher=
2018
-
[13]
Journal of the American Statistical Association , volume=
Bayesian model search for nonstationary periodic time series , author=. Journal of the American Statistical Association , volume=. 2020 , publisher=
2020
-
[14]
1993 , publisher=
Spectral analysis for physical applications , author=. 1993 , publisher=
1993
-
[15]
Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=
Regularization and variable selection via the elastic net , author=. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=. 2005 , publisher=
2005
-
[16]
Excerpts from
Bretthorst, G Larry , booktitle=. Excerpts from. 1988 , publisher=
1988
-
[17]
Journal of the American Statistical Association , volume=
Bayesian estimation of the spectral density of a time series , author=. Journal of the American Statistical Association , volume=. 2004 , publisher=
2004
-
[18]
Psychophysiology , volume=
Recommendations and publication guidelines for studies using frequency domain and time-frequency domain analyses of neural time series , author=. Psychophysiology , volume=. 2022 , publisher=
2022
-
[19]
Handbook of
Tadesse, Mahlet G and Vannucci, Marina , year=. Handbook of. Chapman and Hall/CRC. , publisher=
-
[20]
Annals of neurology , volume=
Seizure Cycles under Pharmacotherapy , author=. Annals of neurology , volume=. 2024 , publisher=
2024
-
[21]
Brain Communications , volume=
Ultradian rhythms in accelerometric and autonomic data vary based on seizure occurrence in paediatric epilepsy patients , author=. Brain Communications , volume=. 2024 , publisher=
2024
-
[22]
Sleep , volume=
Sleep--wake behavioral characteristics associated with depression symptoms: findings from the Multi-Ethnic Study of Atherosclerosis , author=. Sleep , volume=. 2024 , publisher=
2024
-
[23]
Epilepsia , volume=
Actigraphic correlates of neuropsychiatric symptoms in adults with focal epilepsy , author=. Epilepsia , volume=. 2023 , publisher=
2023
-
[24]
Journal of Machine Learning Research , volume=
From fourier to koopman: Spectral methods for long-term time series prediction , author=. Journal of Machine Learning Research , volume=
-
[25]
Advances in Neural Information Processing Systems , volume=
Fourier-transform-based attribution priors improve the interpretability and stability of deep learning models for genomics , author=. Advances in Neural Information Processing Systems , volume=
-
[26]
1998 , publisher=
Spectral analysis of time-series data , author=. 1998 , publisher=
1998
-
[27]
IEEE Control Systems Letters , volume=
Spectral characterization of the multi-seasonal component of the Italian electric load: a LASSO-FFT approach , author=. IEEE Control Systems Letters , volume=. 2019 , publisher=
2019
-
[28]
2019 IEEE International Conference on Software Maintenance and Evolution (ICSME) , pages=
Performance-influence model for highly configurable software with fourier learning and lasso regression , author=. 2019 IEEE International Conference on Software Maintenance and Evolution (ICSME) , pages=. 2019 , organization=
2019
-
[29]
Variable selection via
George, Edward I and McCulloch, Robert E , journal=. Variable selection via. 1993 , publisher=
1993
-
[30]
Approaches for
George, Edward I and McCulloch, Robert E , journal=. Approaches for. 1997 , publisher=
1997
-
[31]
Spike and slab variable selection: frequentist and Bayesian strategies , author=
-
[32]
Journal of the american statistical association , volume=
Bayesian variable selection in linear regression , author=. Journal of the american statistical association , volume=. 1988 , publisher=
1988
-
[33]
medRxiv , year=
Stability and Volatility of Human Rest-Activity Rhythms: Insights from Very Long Actograms (VLAs) , author=. medRxiv , year=
-
[34]
Alzheimer Disease & Associated Disorders , volume=
Impact of disturbed rest-activity rhythms on activities of daily living in moderate and severe dementia patients , author=. Alzheimer Disease & Associated Disorders , volume=. 2021 , publisher=
2021
-
[35]
Sleep Health , volume =
Sleep-wake behaviors associated with cognitive performance in middle-aged participants of the Hispanic Community Health Study/Study of Latinos , author=. Sleep Health , volume =. 2024 , publisher=
2024
-
[36]
Seizure , volume=
More than sleep and wake disturbances: An actigraphic study showing the sleep-wake pattern dysregulation in epilepsy , author=. Seizure , volume=. 2022 , publisher=
2022
-
[37]
1999 , publisher=
Discrete-time signal processing , author=. 1999 , publisher=
1999
-
[38]
Communications of the Association for Computing Machinery , volume=
A spectral method for confidence interval generation and run length control in simulations , author=. Communications of the Association for Computing Machinery , volume=. 1981 , publisher=
1981
-
[39]
American Journal of Physiology-Regulatory, Integrative and Comparative Physiology , volume=
Circadian rhythm of heat production, heart rate, and skin and core temperature under unmasking conditions in men , author=. American Journal of Physiology-Regulatory, Integrative and Comparative Physiology , volume=. 1994 , publisher=
1994
-
[40]
Progress in brain research , volume=
Mechanisms and functions of coupling between sleep and temperature rhythms , author=. Progress in brain research , volume=. 2006 , publisher=
2006
-
[41]
Bayesian Analysis , volume=
Bayesian approximations to hidden semi-Markov models for telemetric monitoring of physical activity , author=. Bayesian Analysis , volume=. 2023 , publisher=
2023
-
[42]
Annual Review of Pharmacology and Toxicology , volume=
Circadian rhythms: mechanisms and therapeutic implications , author=. Annual Review of Pharmacology and Toxicology , volume=. 2007 , publisher=
2007
-
[43]
NPJ Digital Medicine , volume=
Picado: a mobile e-health platform for real-time monitoring of circadian dynamics of symptoms and physiological functions in cancer patients , author=. NPJ Digital Medicine , volume=. 2018 , publisher=
2018
-
[44]
Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=
Regression shrinkage and selection via the lasso , author=. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=. 1996 , publisher=
1996
-
[45]
Journal of the American statistical Association , volume=
Variable selection via nonconcave penalized likelihood and its oracle properties , author=. Journal of the American statistical Association , volume=. 2001 , publisher=
2001
-
[46]
Annals of Statistics , volume=
Nearly unbiased variable selection under minimax concave penalty , author=. Annals of Statistics , volume=
-
[47]
Sensors , volume=
Smart wearables for cardiac monitoring—real-world use beyond atrial fibrillation , author=. Sensors , volume=. 2021 , publisher=
2021
-
[48]
Sleep , volume=
State of the science and recommendations for using wearable technology in sleep and circadian research , author=. Sleep , volume=. 2024 , publisher=
2024
-
[49]
Nature biotechnology , volume=
Wearable biosensors for healthcare monitoring , author=. Nature biotechnology , volume=. 2019 , publisher=
2019
-
[50]
Biology , volume=
Episodic ultradian events—ultradian rhythms , author=. Biology , volume=. 2019 , publisher=
2019
-
[51]
Journal of biological rhythms , volume=
Evidence for a coupled oscillator model of endocrine ultradian rhythms , author=. Journal of biological rhythms , volume=. 2018 , publisher=
2018
-
[52]
Journal of the American Statistical Association , year=
Bayesian spectral modeling for multiple time series , author=. Journal of the American Statistical Association , year=
-
[53]
Computational Statistics & Data Analysis , volume=
Fast Bayesian inference on spectral analysis of multivariate stationary time series , author=. Computational Statistics & Data Analysis , volume=. 2023 , publisher=
2023
This paper was first reviewed by deepseek-v4-flash on August 1, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.