Pith. sign in

REVIEW 2 major objections 6 minor 66 references

AdaptSPEC-X: Covariate Dependent Spectral Modeling of Multiple Nonstationary Time Series

T0 review · 2 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read AdaptSPEC-X jointly estimates time-varying means and spectra for a panel of nonstationary time series, handles missing values, and makes predictions at unobserved covariate values.

desk verdict A competent synthesis that extends AdaptSPEC to panels with covariates and missingness; the missing-data imputation is the load-bearing assumption that needs testing. read the letter →

arxiv 1908.06622 v2 pith:6CT5VANI submitted 2019-08-19 stat.ME

Michael Bertolacci , Ori Rosen , Edward Cripps , Sally Cripps This is my paper · ORCID
classification stat.ME MSC 62M1562F1562G08
keywords nonstationarytimeseriesspectralestimationcovariate-dependentmixturelogisticstickbreakingprocessmissingdatasmoothingsplinesBayesiannonparametricslocallystationaryprocesses
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

AdaptSPEC-X is a Bayesian method for analyzing a panel of possibly nonstationary time series as a single object. The paper claims that modeling each series with a covariate-dependent infinite mixture of AdaptSPEC components, piecewise stationary processes whose segment log spectra are smoothing splines, allows time-varying means and spectra to be estimated jointly across series, with missing values handled inside the model and predictive inference available at unobserved covariate values. These features matter because single-series spectral estimates cannot be improved by finer sampling once the record is fixed, whereas a panel can borrow strength across similar series. In a simulation with four regimes and 10% missingness, the method recovers the true mean and log spectrum, with median mean-squared error below 0.02 for the mean and below 0.08 for the spectrum at almost all observed and unobserved locations. Applications to Australian rainfall and US measles incidence reproduce established drought and post-vaccine signals and add a new one: estimated rainfall variability has declined since 1950.

What carries the argument

The central object is the AdaptSPEC-X model: a covariate-dependent infinite mixture (Equation 6) whose components are AdaptSPEC models, piecewise stationary processes in which the unknown number of segments and their boundaries are sampled by reversible-jump MCMC and each segment's log spectrum is a smoothing spline. Covariate dependence enters through logistic stick-breaking weights (Equation 8), whose log odds are modeled by a thin-plate Gaussian process prior on the covariates (Equation 9). Missing values are handled by treating the Whittle likelihood as an exact multivariate normal with a symmetric circulant precision matrix built from the estimated spectrum, so the conditional distribution of missing observations is multivariate normal (Equations 11-13). The MCMC scheme combines data augmentation for missing values, Polya-Gamma sampling for the stick-breaking coefficients, Riemann manifold Hamiltonian Monte Carlo for the spline coefficients, and a label-swapping move to improve mixing.

What would settle it

Simulate panels from a known non-Gaussian, zero-inflated process, such as monthly aggregated rainfall with a known time-varying spectrum and 10-30% missingness, then compare AdaptSPEC-X's imputed values and estimated log spectra against the truth; systematic bias in imputations or median log-spectrum MSE well above the Gaussian-case values reported in Section 5 would show that the multivariate normal missing-data assumption is the failing component.

Watch

Extended reading notes

Core claim

The central claim is that a panel of nonstationary time series can be modeled jointly by an infinite mixture of AdaptSPEC components, with mixture weights driven by time-independent covariates through a logistic stick-breaking process. Each component segments a series into an unknown number of locally stationary pieces, estimates a time-varying mean and a smoothing-spline log spectrum per piece, and imputes missing values from the multivariate normal conditional implied by the Whittle likelihood. The paper argues that this structure simultaneously solves four problems that had been treated separately: multiple series, nonstationarity in both mean and spectrum, multiple covariates, and missing data. In the simulation study, the model recovers the true time-varying mean and log spectrum at observed and unobserved covariate values; in the applications, it detects the World War II and Millennium droughts and the post-1963 measles decline, and it estimates that rainfall variability across Australia has declined since 1950.

Load-bearing premise

The paper's load-bearing premise is that the Whittle likelihood can be treated as an exact multivariate normal distribution for the entire series, so that missing values are drawn from a Gaussian whose covariance is built from the estimated spectrum; if that approximation is poor, for example with strongly non-Gaussian or zero-inflated rainfall, or with short segments where the Whittle error is large, the imputed values and every downstream estimate inherit the bias.

Editorial extensions

If this is right

  • Joint modeling across a panel improves local spectral estimation where the asymptotics of single-series locally stationary processes cannot be improved by further observation.
  • Panels with substantial missingness, 10% in simulation and 26% in the measles application, can be analyzed by nonparametric spectral methods without discarding or ad-hoc gap-filling the missing periods.
  • Spectral quantities can be predicted at unobserved covariate values, so features like the seasonal rainfall cycle can be interpolated to locations with no gauge record.
  • Time-varying means are estimated within the model, so series need not be pre-centered, and low-frequency information is not removed by a preliminary detrending step.
  • The method scales to large panels, about 192,000 observations in the measles application, through a combination of basis truncation and efficient sampling.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The missing-data assumption is the part most likely to limit the method: for zero-inflated rainfall-like series, a transformation or copula layer over the Gaussian imputation is a natural extension that could be tested by masking known values and comparing imputations.
  • The thin-plate GP prior on the log odds is smooth in covariate space, so regimes that change abruptly across space will be smoothed over; a piecewise or nonstationary GP would sharpen those boundaries.
  • Extending the LSBP to time-varying covariates, such as climate indices or policy changes, would let the mixture weights react to external drivers; the paper lists this as future work, but the model machinery suggests the extension is direct.
  • A head-to-head comparison with single-series AdaptSPEC on the same simulated panels would isolate how much accuracy is gained from borrowing strength across series, which the paper motivates but does not separately quantify.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. This paper introduces AdaptSPEC-X, a Bayesian covariate-dependent infinite mixture model for a panel of nonstationary time series. Each mixture component is an AdaptSPEC model; the mixture weights are driven by time-independent covariates through a logistic stick-breaking process with a thin-plate Gaussian process prior on the log odds; and AdaptSPEC is extended with (i) an unknown time-varying segment mean and (ii) missing-data imputation based on treating the Whittle likelihood as a multivariate normal density (Section 3.3, Eqs. 11-13). Estimation is by MCMC combining data augmentation, reversible-jump moves, a label-swapping step, and a Riemann manifold Hamiltonian Monte Carlo update of the spline coefficients, for which the metric tensor is shown to be constant in the coefficients (Appendix A.1.2). The method is evaluated on a simulation with Gaussian AR(2) data from four known covariate regions, with 10% missingness, and applied to Australian monthly rainfall and US measles incidence. The simulation shows recovery of the true time-varying mean and log spectrum at most locations, and the applications produce the reported drought and vaccination signals.

Significance. If the claims hold, the paper makes a substantial methodological contribution: it packages the per-series AdaptSPEC model into a panel model with covariate-dependent weights, which is a natural and useful extension, and it demonstrates predictive inference at unobserved covariate values. The technical execution is careful: the truncated-normal conditional for the segment mean (Appendix A.1.1) and the constant metric tensor of the RMHMC step (Eq. A.4) are derived correctly, including the observation that only the zero-frequency periodogram depends on the mean. The paper also ships reproducible code and data on GitHub, and the Discussion honestly lists the main limitations (smooth stationary GP for the log odds, Whittle inefficiency for short segments, measurement error, discrete counts, time-varying covariates). The main weakness is that the headline claim of handling missing values is validated only for Gaussian data, while the applications involve nonnegative, zero-inflated series for which the Gaussian imputation assumption is most questionable.

major comments (2)
  1. [Section 3.3, Eqs. (11)-(13)] The imputation step treats the Whittle likelihood as an exact multivariate normal density for the whole segment and draws missing values from the implied conditional Gaussian. This is a genuine modeling assumption, not an asymptotic convenience, for finite n. In the two applications the data are nonnegative and zero-inflated (9,933 zero months among the rainfall series and 30,439 zero weekly counts among the measles series), so the unbounded Gaussian imputation can produce negative rainfall or incidence values with positive probability and cannot represent point masses at zero. Because the simulation study of Section 5 uses Gaussian AR(2) data with only 10% missingness, it does not exercise the assumption for the data types used in the applications, and the claim that AdaptSPEC-X 'handles missing values' is therefore unverified for exactly the cases where it matters most. I would like to see either a simulation with non-Gaussian and/or zero-inflated data under artificially induced missingness, or a diagnostic in the applications (for example, the proportion of imputed draws that are negative or the distribution of imputed values near zero), together with a discussion of the consequences if that proportion is non-negligible. Section 7 acknowledges Whittle small-sample inefficiency but does not acknowledge this specific assumption or the validation gap.
  2. [Section 5, Figure 3] The simulation study reports absolute recovery errors but contains no comparison with any existing method, even though the paper itself cites close competitors (Bruce et al. 2018, Krafty et al. 2017) and the natural baseline of fitting AdaptSPEC separately to each series. Without such a comparison it is difficult to judge what the covariate-dependent mixture and the LSBP weights add over simpler alternatives, and the notably larger MSE at D2/T2 (median 0.09 for the mean and 0.34 for the spectrum versus below 0.02 and 0.08 elsewhere) is the only quantitative evidence about where the modeling assumptions bite. Adding at least one baseline to the replicated simulation would make the empirical evaluation proportionate to the model's complexity and would substantially strengthen the paper's claims.
minor comments (6)
  1. [Section 3.2, Eq. (9)] The text says '1N is an n×1 vectors of ones', but the vector should be N×1; this is a typographical error that should be corrected.
  2. [Section 3.2] The notation wh is used both for the log-odds function wh(·) and for the stacked vector (wh(u1),...,wh(uN))'; a bold or subscript convention would avoid ambiguity in Eqs. (9)-(10).
  3. [Footnotes 1-2] The abstract states that software is available in the R package BayesSpec, while the footnotes say the latest CRAN version does not contain AdaptSPEC-X and the code is 'available from the authors'; since the Reproducibility section provides a GitHub link, the footnotes should direct readers there and clarify the CRAN status.
  4. [Section 6.2, Figures 12-14] The authors note that the full-period spectra in Figure 12 look almost identical across states because of the wide power range; since Figure 14 is needed to see geographic heterogeneity, a per-panel normalization or a shared scale restricted to the pre-vaccine period would make Figure 12 more informative.
  5. [Section 5, Process (14)] The four regions have very different membership counts (41, 8, 18, and 33 time series), so the larger MSE at D2/T2 conflates cluster size with covariate structure; a balanced design, or reporting results conditional on cluster size, would separate these two explanations for the degraded performance.
  6. [Section 7, Discussion] The Discussion candidly lists Whittle inefficiency and measurement error as limitations but does not mention the Gaussian missing-data imputation assumption of Section 3.3; a sentence acknowledging this assumption and pointing to the suggested validation would help readers calibrate the method's scope.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: the central simulation is external, and the cited building blocks are not used as self-supporting evidence.

full rationale

We walked the derivation chain from the Whittle likelihood (Eq. 1) through the AdaptSPEC-X mixture (Eqs. 5-10), the missing-data imputation (Eqs. 11-13), and the MCMC scheme (Section 4). The simulation study (Section 5) generates 100 replicates from the known AR(2) process (Eq. 14) with specified time-varying means and spectra and 10% missingness, then compares estimates to that known truth via MSE (Eqs. 15-16). This is an external benchmark; the reported recovery does not reduce to fitted constants or to predictions that are equal to model inputs. The missing-value step treats the Whittle likelihood as an exact multivariate normal and draws missing values from the conditional distribution (11)-(13); this is a statistical modeling assumption and a potential correctness risk for non-Gaussian applications, but it is not circular, since no fitted parameter is renamed as a prediction and the imputation distribution is derived openly from the model's own likelihood. Citations to Rosen et al. (2012) and Rigon and Durante (2020) supply the AdaptSPEC and LSBP building blocks, and some are authored by members of the present team, but the paper does not invoke a uniqueness theorem or forbid alternatives by self-citation; the central claims are tested against synthetic data with known truth. We therefore find no step in which a claimed prediction is identical, by construction, to its inputs.

Assumptions & free parameters 8 free parameters · 6 assumptions · 0 invented entities

The central claim rests on several hand-chosen tuning parameters and on the Whittle-based multivariate normal assumption for missing data. The method introduces no new physical or mathematical entities.

free parameters (8)
  • H (mixture truncation) = 25 (simulation, rainfall), 10 (measles)
    Truncation of the infinite stick-breaking representation; chosen by hand. Authors state higher values made no difference, but no sensitivity analysis is shown.
  • M (maximum number of segments) = 4 (simulation), 17 (rainfall), 18 (measles)
    Prior for the number of segments m is uniform between 1 and M; M is determined by tmin and series length and chosen by hand.
  • tmin (minimum segment length) = 40 (simulation), 60 months (rainfall), 208 weeks (measles)
    Set to ensure enough observations per segment for the Whittle approximation; this choice directly controls segmentation flexibility.
  • J (spline basis functions for log spectrum) = 25 (simulation), 60 (applications)
    Balances prior flexibility and computational cost; chosen by hand without sensitivity analysis.
  • B (number of LSBP basis functions) = 10 (simulation), 20 (applications)
    Truncated basis expansion for the log-odds GP; chosen to capture more than 95% of prior variation, but no sensitivity analysis is reported.
  • Prior range (µ-, µ+) for segment means = (-10, 10) simulation; (0, 30) rainfall; (0, 20) measles
    Support of the uniform prior on segment means; in applications the lower bound enforces positivity of the mean.
  • Prior hyperparameters µβ, Σβ, ντ, Aτ = 0, 100I, 3, 10
    Settings for the LSBP regression coefficients and the half-t prior on τ; fixed without sensitivity analysis.
  • RMHMC leapfrog tuning = step size uniform in [0.1, 1], steps uniform in [1, 10]
    Algorithmic tuning chosen uniformly at random per update; affects mixing of the spline coefficients.
assumptions (6)
  • domain assumption The Whittle likelihood (Eq. 1) is a valid approximation to the true Gaussian likelihood for large n and is used as the likelihood for each segment.
    Invoked throughout Section 2 and for missing-data imputation in Section 3.3. Known to be inaccurate for small samples; the paper cites Contreras-Cristán et al. (2006) and Sykulski et al. (2019) as related limitations.
  • ad hoc to paper Missing values follow the multivariate normal distribution implied by the Whittle likelihood (Eq. 11).
    Section 3.3. This is a modeling assumption justified by reference to Guinness (2019) and parsimony; it is not a theorem and could bias imputations if the process is non-Gaussian or Whittle error is large.
  • domain assumption Each time series is partitioned into independent stationary segments, with independence across segments and across series given component indicators.
    Equations (4)-(5) and Section 3. Piecewise stationarity is standard for AdaptSPEC but ignores dependence across segments and spatial dependence between series.
  • domain assumption Covariates are time-independent and the mixture weights depend on them only through a smooth stationary Gaussian process.
    Section 3.2 and Discussion: the thin-plate GP is smooth and stationary, which may fail for abrupt changes in covariate effects.
  • standard math Standard Bayesian and MCMC results used without proof (Polya-Gamma augmentation, reversible jump, RMHMC).
    Appendices A.2-A.4 rely on cited results from Polson et al. (2013), Rigon and Durante (2020), and Girolami and Calderhead (2011). These are accepted background.
  • domain assumption The model is identifiable up to label switching, managed by the label swapping step of Section 4.
    The label swapping move from Hastie et al. (2015) is assumed to improve mixing; no formal identifiability proof is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of AdaptSPEC-X: Covariate Dependent Spectral Modeling of Multiple Nonstationary Time Series." pith.science (2026). https://pith.science/paper/6CT5VANI

@misc{pith2026190806622,
  author       = {Pith},
  title        = {Pith review of: AdaptSPEC-X: Covariate Dependent Spectral Modeling of Multiple Nonstationary Time Series},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6CT5VANI}},
  note         = {Machine review of arXiv:1908.06622}
}
read the original abstract

We present a method for the joint analysis of a panel of possibly nonstationary time series. The approach is Bayesian and uses a covariate-dependent infinite mixture model to incorporate multiple time series, with mixture components parameterized by a time varying mean and log spectrum. The mixture components are based on AdaptSPEC, a nonparametric model which adaptively divides the time series into an unknown number of segments and estimates the local log spectra by smoothing splines. We extend AdaptSPEC to handle missing values, a common feature of time series which can cause difficulties for nonparametric spectral methods. A second extension is to allow for a time varying mean. Covariates, assumed to be time-independent, are incorporated via the mixture weights using the logistic stick breaking process. The model can estimate time varying means and spectra at observed and unobserved covariate values, allowing for predictive inference. Estimation is performed by Markov chain Monte Carlo (MCMC) methods, combining data augmentation, reversible jump, and Riemann manifold Hamiltonian Monte Carlo techniques. We evaluate the methodology using simulated data, and describe applications to Australian rainfall data and measles incidence in the US. Software implementing the method proposed in this paper is available in the R package BayesSpec.

Figures

Figures reproduced from arXiv: 1908.06622 by the authors.

Figure 1
Figure 1. A graphical representation of AdaptSPEC-X. The bottom row lists the main hyperparameters. [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. (a) Underlying surface mapping uj to zj for Process (14), where the 100 sampled locations are shown as crosses. The color of a region indicates the corresponding cluster. Four of the crosses are used as examples in the paper and are colored green and labeled D1 to D4. Four test points labeled T1 to T4 are shown with red diamonds. (b) Example realizations from Process (14) corresponding to zj = 1 at the top through z… view at source ↗
Figure 3
Figure 3. MSE across 100 replications for estimates of mean (top) and spectrum (bottom) for Pro [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figures from the paper (11 more)
Figure 4
Figure 4. Figure 4: Estimated mean ˆµ(t,u) corresponding to the median MSEmean(u) (red) and true mean µ(t,u) (blue) for Process (14). The first row shows the estimates for u = D1–D4 from left to right, respectively, while the second row shows the estimates for test points u = T1–T4 [PITH…
Figure 5
Figure 5. Figure 5: Estimated time varying log spectra log ˆf(t, ω,u) corresponding to the median MSEspec(u) and true time varying log spectra log f(t, ω,u) for Process (14). The first row shows log f(t, ω,u) for u = D1–D4 from left to right, respectively, while the second row shows the e…
Figure 6
Figure 6. Figure 6: (a) Locations of the 151 rainfall sites. Four example sites are marked with green circles, and [PITH_FULL_IMAGE:figures/full_fig_p015_6.png]
Figure 7
Figure 7. Figure 7: Estimated time varying means ˆµ(t,u) for four monthly rainfall sites (first two rows), and four locations without observations (last two rows). Site (u) Event ˆp(· | x) 10525 14042 47053 69018 WW2 drought µ(1940-01,u) < µ(1930-01,u) 0.539 0.442 0.700 0.802 σ 2 (1940-01…
Figure 8
Figure 8. Figure 8: Estimated time varying spectra log ˆf(t, ω,u) for four monthly rainfall sites (first two rows) and four locations without observations (last two rows). The color indicates the log power at the corresponding time and frequency. The ω-axis is on a square-root scale. The …
Figure 9
Figure 9. Figure 9: The same estimated time varying spectra log [PITH_FULL_IMAGE:figures/full_fig_p019_9.png]
Figure 10
Figure 10. Figure 10: Measles incidence rate per 100,000 population for the continental US (that is, excluding [PITH_FULL_IMAGE:figures/full_fig_p021_10.png]
Figure 11
Figure 11. Figure 11: Estimated time varying mean ˆµ(t,u) for measles incidence, where each panel shows the estimate for one state. 23 [PITH_FULL_IMAGE:figures/full_fig_p023_11.png]
Figure 12
Figure 12. Figure 12: Estimated time varying log spectra, log ˆf(t, ω,u) for measles incidence, where each panel shows the estimate for one state. Colors indicate the log power at the corresponding date and ω. The ω-axis is on a square-root scale. The top axis displays the period (1/ω). 24…
Figure 13
Figure 13. Figure 13: Estimated time varying log spectra, log ˆf(t, ω,u), for four states. Estimates for the full study period (as in [PITH_FULL_IMAGE:figures/full_fig_p025_13.png]
Figure 14
Figure 14. Figure 14: Estimated time varying log spectra, log ˆf(t, ω,u) for measles incidence for the pre-vaccine period (<1963) [PITH_FULL_IMAGE:figures/full_fig_p026_14.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

66 extracted references · 64 canonical work pages

  1. [1]

    (1998), Time-Dependent Spectral Analysis of Nonstationary Time Series , Journal of the American Statistical Association, 93, 1488--1501

    Adak, S. (1998), Time-Dependent Spectral Analysis of Nonstationary Time Series , Journal of the American Statistical Association, 93, 1488--1501

    work page 1998

  2. [2]

    (2018), An Introduction to Categorical Data Analysis, 3rd Edition , Wiley

    Agresti, A. (2018), An Introduction to Categorical Data Analysis, 3rd Edition , Wiley

    work page 2018

  3. [3]

    (1924), Application of Schuster's periodogram to long rainfall records, beginning 1748 , Monthly Weather Review, 52, 479--487

    Alter, D. (1924), Application of Schuster's periodogram to long rainfall records, beginning 1748 , Monthly Weather Review, 52, 479--487

    work page 1924

  4. [4]

    J., Reason, C

    Ansell, T. J., Reason, C. J. C., Smith, I. N., and Keay, K. (2000), Evidence for decadal variability in southern Australian rainfall and relationships with regional pressure and sea surface temperature , International Journal of Climatology, 20, 1113--1129

    work page 2000

  5. [5]

    L., Orenstein, W

    Atkinson, W. L., Orenstein, W. A., and Krugman, S. (1992), The Resurgence of Measles in the United States, 1989-1990 , Annual Review of Medicine, 43, 451--463

    work page 1992

  6. [6]

    W., and Cripps, S

    Bertolacci, M., Cripps, E., Rosen, O., Lau, J. W., and Cripps, S. (2019), Climate inference on daily rainfall across the Australian continent, 1876--2015 , Annals of Applied Statistics, 13, 683--712

    work page 2019

  7. [7]

    A., Hall, M

    Bruce, S. A., Hall, M. H., Buysse, D. J., and Krafty, R. T. (2018), Conditional adaptive Bayesian spectral analysis of nonstationary biomedical time series , Biometrics, 74, 260--269

    work page 2018

  8. [8]

    (2018), Bayesian spectral modeling for multiple time series, Journal of the American Statistical Association, 0, 1--38

    Cadonna, A., Kottas, A., and Prado, R. (2018), Bayesian spectral modeling for multiple time series, Journal of the American Statistical Association, 0, 1--38

    work page 2018

Show all 66 references
  1. [9]

    J., Jones, D

    Cai, W. J., Jones, D. A., Harle, K., Cowan, T. D., Power, S. B., Smith, I. N., Arblaster, J. M., and Abbs, D. J. (2007), Past climate change, in Climate Change in Australia, ed. Pearce, K. B. e. a., Aspendale, Victoria, Australia: CSIRO Marine and Atmospheric Research, pp. 17--28

    2007

  2. [10]

    CDC (2019), Centers for Disease Control and Prevention: Measles Cases and Outbreaks , https://www.cdc.gov/measles/cases-outbreaks.html, accessed: 2019-03-21

    2019

  3. [11]

    (2004), Bayesian Estimation of the Spectral Density of a Time Series , Journal of the American Statistical Association, 99, 1050--1059

    Choudhuri, N., Ghosal, S., and Roy, A. (2004), Bayesian Estimation of the Spectral Density of a Time Series , Journal of the American Statistical Association, 99, 1050--1059

    2004

  4. [12]

    and Dunson, D

    Chung, Y. and Dunson, D. B. (2009), Nonparametric B ayes conditional distribution modeling with variable selection, Journal of the American Statistical Association, 104, 1646--1660

    2009

  5. [13]

    Contreras-Crist\' a n, A., Guti\' e rrez-Pe\ n a, E., and Walker, S. G. (2006), A Note on Whittle's Likelihood , Communications in Statistics - Simulation and Computation, 35, 857--875

    2006

  6. [14]

    (1997), Fitting time series models to nonstationary processes, The Annals of Statistics, 25, 1--37

    Dahlhaus, R. (1997), Fitting time series models to nonstationary processes, The Annals of Statistics, 25, 1--37

    1997

  7. [15]

    A., Lee, T

    Davis, R. A., Lee, T. C. M., and Rodriguez-Yam, G. A. (2006), Structural Break Estimation for Nonstationary Time Series Models , Journal of the American Statistical Association, 101, 223--239

    2006

  8. [16]

    D., Pendleton, B

    Duane, S., Kennedy, A. D., Pendleton, B. J., and Roweth, D. (1987), Hybrid Monte Carlo , Physics Letters B, 195, 216--222

    1987

  9. [17]

    Gallant, A. J. E., Reeder, M. J., Risbey, J. S., and Hennessy, K. J. (2013), The characteristics of seasonal-scale droughts in Australia, 1911–2009 , International Journal of Climatology, 33, 1658--1672

    2013

  10. [18]

    (2006), Prior distributions for variance parameters (Comment on article by B rowne and D raper), Bayesian Analysis, 1, 515--534

    Gelman, A. (2006), Prior distributions for variance parameters (Comment on article by B rowne and D raper), Bayesian Analysis, 1, 515--534

    2006

  11. [19]

    and Calderhead, B

    Girolami, M. and Calderhead, B. (2011), Riemann manifold Langevin and Hamiltonian Monte Carlo methods , Journal of the Royal Statistical Society: Series B (Statistical Methodology), 73, 123--214

    2011

  12. [20]

    Green, P. J. (1995), Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , Biometrika, 82, 711--732

    1995

  13. [21]

    T., Bj rnstad, O

    Grenfell, B. T., Bj rnstad, O. N., and Kappey, J. (2001), Travelling waves and spatial hierarchies in measles epidemics, Nature, 414, 716--723

    2001

  14. [22]

    (2019), Spectral density estimation for random fields via periodic embeddings, Biometrika, 106, 267--286

    Guinness, J. (2019), Spectral density estimation for random fields via periodic embeddings, Biometrika, 106, 267--286

    2019

  15. [23]

    C., and von Sachs, R

    Guo, W., Dai, M., Ombao, H. C., and von Sachs, R. (2003), Smoothing Spline ANOVA for Time-Dependent Spectral Analysis , Journal of the American Statistical Association, 98, 643--652

    2003

  16. [24]

    Hannan, E. J. (1973), The asymptotic theory of linear time-series models, Journal of Applied Probability, 10, 130–--145

    1973

  17. [25]

    I., Liverani, S., and Richardson, S

    Hastie, D. I., Liverani, S., and Richardson, S. (2015), Sampling from Dirichlet process mixture models with unknown concentration parameter: mixing issues in large data implementations , Statistics and Computing, 25, 1023--1037

    2015

  18. [26]

    R., Brandling-Bennett, A

    Hinman, A. R., Brandling-Bennett, A. D., and Nieburg, P. I. (1979), The Opportunity and Obligation to Eliminate Measles From the United States , J. Am. Med. Assoc, 242, 1157--1162

    1979

  19. [27]

    K., Drosdowsky, W., and Nicholls, N

    Hope, P. K., Drosdowsky, W., and Nicholls, N. (2006), Shifts in the synoptic systems influencing southwest Western Australia , Climate Dynamics, 26, 751--764

    2006

  20. [28]

    and James, L

    Ishwaran, H. and James, L. F. (2001), Gibbs Sampling Methods for Stick-Breaking Priors , Journal of the American Statistical Association, 96, 161--173

    2001

  21. [29]

    A., Jordan, M

    Jacobs, R. A., Jordan, M. I., Nowlan, S. J., and Hinton, G. E. (1991), Adaptive Mixtures of Local Experts , Neural Comput., 3, 79--87

    1991

  22. [30]

    Joshi, M. K. and Pandey, A. C. (2011), Trend and spectral analysis of rainfall over India during 1901–2000 , Journal of Geophysical Research: Atmospheres, 116

    2011

  23. [31]

    Katz, S. L. and Hinman, A. R. (2004), Summary and Conclusions: Measles Elimination Meeting, 16–17 March 2000 , The Journal of Infectious Diseases, 189, S43--S47

    2004

  24. [32]

    and Gersch, W

    Kitagawa, G. and Gersch, W. (1996), Smoothness Priors Analysis of Time Series , Springer

    1996

  25. [33]

    T., Hall, M., and Guo, W

    Krafty, R. T., Hall, M., and Guo, W. (2011), Functional mixed effects spectral analysis, Biometrika, 98, 583--598

    2011

  26. [34]

    T., Rosen, O., Stoffer, D

    Krafty, R. T., Rosen, O., Stoffer, D. S., Buysse, D. J., and Hall, M. H. (2017), Conditional Spectral Analysis of Replicated Multiple Time Series With Application to Nocturnal Physiology , Journal of the American Statistical Association, 112, 1405--1416

    2017

  27. [35]

    (1992), A historical rainfall data set for Australia , Aust

    Lavery, B., Kariko, A., and Nicholls, N. (1992), A historical rainfall data set for Australia , Aust. Met. Mag., 40, 33--39

    1992

  28. [36]

    and Krafty, R

    Li, Z. and Krafty, R. T. (2018), Adaptive Bayesian Time–Frequency Analysis of Multivariate Time Series , Journal of the American Statistical Association, 0, 1--13

    2018

  29. [37]

    and Daum, R

    Lynfield, R. and Daum, R. S. (2014), The Complexity of the Resurgence of Childhood Vaccine-Preventable Diseases in the United States , Current Pediatrics Reports, 2, 195--203

    2014

  30. [38]

    and Prado, R

    Macaro, C. and Prado, R. (2014), Spectral Decompositions of Multiple Time Series: A Bayesian Non-parametric Approach , Psychometrika, 79, 105--129

    2014

  31. [39]

    Moreno, M. A. (2018), Measles, JAMA Pediatrics, 172, 896--896

    2018

  32. [40]

    Neal, R. M. (1993), Probabilistic inference using Markov Chain Monte Carlo methods , Tech. Rep. CRG-TR-93-1, Department of Computer Science, University of Toronto

    1993

  33. [41]

    Brooks, S., Gelman, A., Jones, G

    --- (2011), MCMC using Hamiltonian dynamics , in Handbook of Markov chain Monte Carlo , eds. Brooks, S., Gelman, A., Jones, G. L., and L., M. X

    2011

  34. [42]

    C., Raz, J

    Ombao, H. C., Raz, J. A., von Sachs, R., and Malow, B. A. (2001), Automatic Statistical Analysis of Bivariate Nonstationary Time Series , Journal of the American Statistical Association, 96, 543--560

    2001

  35. [43]

    A., Staggs, W., Dayan, G

    Parker, A. A., Staggs, W., Dayan, G. H., Ortega-S\` a nchez, I. R., Rota, P. A., Lowe, L., Boardman, P., Teclaw, R., Graves, C., and LeBaron, C. W. (2006), Implications of a 2005 Measles Outbreak in Indiana for Sustained Elimination of Measles in the United States , New Englan...

    2006

  36. [44]

    K., Bednarczyk, R

    Phadke, V. K., Bednarczyk, R. A., Salmon, D. A., and Omer, S. B. (2016), Association Between Vaccine Refusal and Vaccine-Preventable Diseases in the United States: A Review of Measles and Pertussis , JAMA, 315, 1149--1158

    2016

  37. [45]

    G., Scott, J

    Polson, N. G., Scott, J. G., and Windle, J. (2013), Bayesian Inference for Logistic Models Using P\' o lya-Gamma Latent Variables , Journal of the American Statistical Association, 108, 1339--1349

    2013

  38. [46]

    J., Risbey, J

    Pook, M. J., Risbey, J. S., and McIntosh, P. C. (2012), The Synoptic Climatology of Cool-Season Rainfall in the Central Wheatbelt of Western Australia , Monthly Weather Review, 140, 28--43

    2012

  39. [47]

    Priestley, M. B. (1965), Evolutionary Spectra and Non-Stationary Processes, Journal of the Royal Statistical Society. Series B (Methodological), 27, 204--237

    1965

  40. [48]

    and Lall, U

    Rajagopalan, B. and Lall, U. (1998), Interannual variability in western US precipitation , Journal of Hydrology, 210, 51 -- 67

    1998

  41. [49]

    and Durante, D

    Rigon, T. and Durante, D. (2020), Tractable Bayesian Density Regression via Logit Stick-Breaking Priors , Journal of Statistical Planning and Inference, in press

    2020

  42. [50]

    S., McIntosh, P

    Risbey, J. S., McIntosh, P. C., and Pook, M. J. (2013), Synoptic components of rainfall variability and trends in southeast Australia , International Journal of Climatology, 33, 2459--2472

    2013

  43. [51]

    and Stoffer, D

    Rosen, O. and Stoffer, D. S. (2007), Automatic estimation of multivariate spectra via smoothing splines, Biometrika, 94, 335--345

    2007

  44. [52]

    S., and Wood, S

    Rosen, O., Stoffer, D. S., and Wood, S. (2009), Local Spectral Analysis via a Bayesian Mixture of Smoothing Splines , Journal of the American Statistical Association, 104, 249--262

    2009

  45. [53]

    Rosen, O., Wood, S., and Stoffer, D. S. (2012), AdaptSPEC: Adaptive spectral estimation for nonstationary time series , Journal of the American Statistical Association, 107, 1575--1589

    2012

  46. [54]

    M., Olhede, S

    Sykulski, A. M., Olhede, S. C., Guillaumin, A. P., Lilly, J. M., and Early, J. J. (2019), The debiased Whittle likelihood , Biometrika, 106, 251--266

    2019

  47. [55]

    C., England, M

    Ummenhofer, C. C., England, M. H., McIntosh, P. C., Meyers, G. A., Pook, M. J., Risbey, J. S., Gupta, A., and Taschetto, A. S. (2009), What causes southeast Australia's worst droughts? Geophysical Research Letters, 36, L04707,, doi:10.1029/2008GL036801

    2009 doi

  48. [56]

    G., Grefenstette, J., Jung, S

    van Panhuis, W. G., Grefenstette, J., Jung, S. Y., Chok, N. S., Cross, A., Eng, H., Lee, B. Y., Zadorozhny, V., Brown, S., Cummings, D., and Burke, D. S. (2013), Contagious Diseases in the United States from 1888 to the Present , New England Journal of Medicine, 369, 2152--2158

    2013

  49. [57]

    (1990), Spline Models for Observational Data , Society for Industrial and Applied Mathematics

    Wahba, G. (1990), Spline Models for Observational Data , Society for Industrial and Applied Mathematics

    1990

  50. [58]

    P., Ormerod, J

    Wand, M. P., Ormerod, J. T., Padoan, S. A., and Fr\" u hworth, R. (2012), Mean field variational Bayes for elaborate distributions , Bayesian Analysis, 7, 847--900

    2012

  51. [59]

    West, M., Prado, R., and Krystal, A. D. (1999), Evaluation and Comparison of EEG Traces: Latent Structure in Nonstationary Time Series , Journal of the American Statistical Association, 94, 1083--1095

    1999

  52. [60]

    (1953), Estimation and information in stationary time series, Arkiv för Matematik, 2, 423--434

    Whittle, P. (1953), Estimation and information in stationary time series, Arkiv för Matematik, 2, 423--434

    1953

  53. [61]

    Series B (Methodological), 19, 38--63

    --- (1957), Curve and Periodogram Smoothing , Journal of the Royal Statistical Society. Series B (Methodological), 19, 38--63

    1957

  54. [62]

    (2013), Applications of Bayesian Smoothing Splines , in Bayesian Theory and Applications , eds

    Wood, S. (2013), Applications of Bayesian Smoothing Splines , in Bayesian Theory and Applications , eds. Damien, P., Dellaportas, P., Polson, N. G., and Stephens, D. A., Oxford University Press, pp. 309--336

    2013

  55. [63]

    and Zhou, Z

    Yang, J. and Zhou, Z. (2020), Spectral Inference under Complex Temporal Dynamics , Journal of the American Statistical Association, 0, 1--50

    2020

  56. [64]

    H., and Wikle, C

    Yang, W.-H., Holan, S. H., and Wikle, C. K. (2016), Bayesian Lattice Filters for Time-Varying Autoregression and Time–Frequency Analysis , Bayesian Anal., 11, 977--1003

    2016

  57. [65]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence a...

  58. [66]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.