---
title: Variable Bandpass Periodic Block Bootstrap
url: https://www.emergentmind.com/topics/variable-bandpass-periodic-block-bootstrap-vbpbb
type: topic
---

# Variable Bandpass Periodic Block Bootstrap

Variable Bandpass Periodic Block Bootstrap (VBPBB) is a nonparametric resampling method for periodically correlated (PC) time series that combines narrowband frequency isolation with period-aware bootstrap resampling. Introduced by Valachovic to address the failure of ordinary bootstrap and conventional block-bootstrap schemes to preserve periodic dependence, VBPBB first separates a target periodic component from interference such as noise or other harmonics and then resamples the isolated component in a way that preserves its periodic correlation structure. Across the literature, the method is used for two closely related purposes: inference for periodic means and harmonics, and structure-preserving imputation for incomplete time series [2402.03491; 2505.02008].

## 1. Emergence and problem formulation

VBPBB was developed for settings in which a time series contains one or more periodic or cyclostationary components, such as annual, weekly, daily, or harmonic structure. In this setting, the central statistical difficulty is not merely temporal dependence, but *periodic* dependence: observations separated by one full cycle or by harmonic multiples can remain strongly correlated, and that correlation depends on phase within the cycle. The literature therefore treats VBPBB as a method for PC series rather than for generic weakly dependent processes [2402.03491].

This development is motivated by limitations of earlier bootstrap procedures. Ordinary bootstrap destroys temporal dependence. The classical Moving-Block Bootstrap preserves short-range dependence up to the block length but destroys correlations at lags beyond that block length, so it cannot reproduce long-lag periodic correlations. Seasonal or Periodic Block Bootstrap methods improve on this by aligning resampling with the cycle, but when they operate on the unfiltered series they also resample other frequencies and noise, which inflates bootstrap variance and widens confidence bands [2509.03710].

The method quickly became a framework for empirical work on environmental, epidemiological, and energy time series. Applications in the supplied literature include US COVID-19 mortality, New York State COVID-19 hospitalizations, Manhattan PM$_{2.5}$, Turkey electricity consumption, ischemic heart disease hospitalizations in New York, and time-series imputation in daily temperature and simulation settings [2403.06343; 2404.11006; 2404.08738; 2404.03786; 2509.03710; 2505.02008]. A related extension, the Variable Multiple Bandpass Periodic Block Bootstrap (VMBPBB), was introduced for time series with multiple periodically correlated components, with separate filtering and resampling for each component before recombination [2502.07462].

## 2. Periodic correlation and bandpass separation

The mathematical setting is a PC process whose mean and covariance repeat with period $p$. One common formulation is
$$
\mu_{t+p}=\mu_t,\qquad \gamma(t+p,\tau)=\gamma(t,\tau),
$$
while another writes
$$
\mathrm{Cov}\bigl(X(t),X(t+h)\bigr)=\gamma_h\bigl(t \bmod p\bigr).
$$
Both formulations express the same idea: second-order structure depends on phase within the cycle rather than only on lag [2402.03491; 2509.03710].

VBPBB uses frequency separation to isolate the PC component of interest before resampling. In the original exposition, the filtered series is obtained with the Kolmogorov–Zurbenko Fourier Transform (KZFT) filter,
$$
Y_t=\mathrm{KZFT}_{m,k,\nu}(X)_t
=\sum_{u=-k(m-1)/2}^{k(m-1)/2} a_u^{(m,k)} e^{-2\pi i \nu u} X_{t+u},
$$
where $m$ is an odd window length, $k$ is the number of iterations, and $\nu=1/p$ is the target frequency. Its energy transfer function is
$$
\bigl|B(\lambda-\nu)\bigr|^2
=
\Bigl|
\frac{1}{m}\frac{\sin(\pi m(\lambda-\nu))}{\sin(\pi(\lambda-\nu))}
\Bigr|^{2k}.
$$
The filter is designed so that energy near the target frequency passes while nearby interference is attenuated [2402.03491].

A closely related exposition writes the bandpass component as the difference of two low-pass Kolmogorov–Zurbenko filters,
$$
B_p(\omega)=H_{m_1,k_1}(\omega)-H_{m_2,k_2}(\omega),
$$
with
$$
H_{m,k}(\omega)=
\Biggl(\frac{1}{m}\sum_{j=-(m-1)/2}^{(m-1)/2} e^{-i\omega j}\Biggr)^k
=
\Biggl(\frac{\sin(\tfrac{m\omega}{2})}{m\,\sin(\tfrac{\omega}{2})}\Biggr)^k.
$$
In time domain form,
$$
P_t=(B_p * X)_t=\sum_j b_{p,j}X_{t-j},
$$
and the residual is $R_t=X_t-P_t$ [2505.02008].

The common decomposition in later work is
$$
X(t)=\sum_{v\in\mathcal V} X^{(v)}(t)+\varepsilon(t),
$$
where each $X^{(v)}$ is a narrowband PC component and $\varepsilon(t)$ is remainder noise. Reported practical guidance ties the filter width to frequency separation. Examples in the literature include choosing $m$ as the smallest odd integer exceeding a rule based on the gap between adjacent target frequencies, or choosing $m$ so that the passband extends no farther than halfway to the nearest unwanted frequency. Reported iteration counts vary by study, including $k=1$, $k=2$, and settings in which $k_j$ is often $3$ [2509.03710; 2404.03786; 2505.02008].

## 3. Bootstrap construction

Once a periodic component has been isolated, VBPBB resamples it in a way that respects periodic phase. The core idea is that resampling must occur on the filtered component, not on the raw series. The original formulation partitions the filtered series into seasonal subsequences
$$
S_i=\{Y_i,Y_{i+p},Y_{i+2p},\dots\},\qquad i=1,\dots,p,
$$
then cycles through these subsequences while sampling with replacement to form a bootstrap replicate. Other expositions describe nonoverlapping blocks of length exactly equal to the period, or phase strata
$$
I_r=\{t:t\equiv r \pmod p\},\qquad r=1,\dots,p,
$$
within which observations are resampled and then written back to their phase-aligned positions. These formulations differ operationally, but they share the same structural principle: preserve the phase-specific correlation pattern of the periodic component [2402.03491; 2509.03710].

In the imputation literature, the algorithm is often presented componentwise. For each significant period $p_j$, the method filters the series to obtain $P_t^{(j)}$, partitions the filtered component into blocks of length $p_j$, samples those blocks with replacement, concatenates the sampled blocks, and reconstructs
$$
X_t^{*(b)}=R_t+\sum_j P_t^{(j),*(b)}
$$
for bootstrap replicate $b$. Amplitude summaries such as $\hat A_j^{*(b)}$ may then be computed from the resampled blocks [2505.02008].

A related bias study allows block length $b_j$ to be an integer multiple of the period, such as $T$ or $2T$, and resamples overlapping blocks from each filtered component. This formulation emphasizes the bias-variance trade-off in block choice: larger $b_j$ reduces truncation bias but can increase bootstrap variance [2509.08647].

Implementation guidance in the supplied papers is consistent on several points. The block length is ordinarily set equal to the period of the component under study. Reported numbers of bootstrap replications include $500$–$1{,}000$, and one practical recommendation is $B\ge 1{,}000$ for stable median estimates. The original computational summary gives total cost
$$
O(nkm+Bn),
$$
combining one KZFT pass and $B$ linear-time bootstrap replicates [2402.03491; 2505.02008; 2508.19535].

## 4. Inference for periodic means and significant components

The main inferential target in early VBPBB work is the periodic mean of a filtered component or of a sum of significant components. For a component with period $p$, one computes phase-wise or time-wise bootstrap summaries and obtains confidence bands from empirical quantiles. One representative expression is
$$
\bigl[\widehat{L}(t),\widehat{U}(t)\bigr]
=
\bigl[
Q_{0.025}\{X^{\mathrm{(tot)},b}(t)\},
Q_{0.975}\{X^{\mathrm{(tot)},b}(t)\}
\bigr],
$$
where
$$
X^{\mathrm{(tot)},b}(t)=\sum_{k=1}^K X^{(v_k),b}(t)
$$
aggregates the bootstrap replicates of the significant components [2509.03710].

The literature reports several closely related significance rules. In one formulation, a component is declared significant if a bootstrap confidence interval for its amplitude excludes zero,
$$
\min \mathrm{CI}_j>0
\quad\Longleftrightarrow\quad
0\notin \mathrm{CI}_j.
$$
An equivalent test statistic is
$$
T_j=\frac{\hat A_j}{\widehat{\sigma}_j},
$$
compared with $z_{1-\alpha/2}$ [2505.02008].

In another formulation, significance is based on the periodic-mean confidence band not admitting a horizontal line, equivalently
$$
\min_t \widehat{U}(t)<\max_t \widehat{L}(t).
$$
Application papers often state the rule more simply: a component is significant when the 95% pointwise confidence band for the bootstrapped periodic mean excludes zero [2509.03710; 2404.08738].

This inferential architecture generalizes naturally to multiple components. VMBPBB formalizes the separate filtering, separate componentwise resampling, and final recombination needed to preserve several periodic correlations simultaneously in a multiple-periodic setting [2502.07462].

## 5. Use in missing-data imputation

A major later use of VBPBB is missing-data imputation for univariate time series with strong seasonal or cyclic structure. In this setting, significant periodic components extracted by VBPBB are supplied to the Amelia II multiple-imputation framework as auxiliary variables. Ahmad et al. describe using median vectors of periodic blocks and/or the bootstrap replicates of the significant components, after which Amelia II’s bootstrap-EM algorithm imputes missing values while preserving the overall sample mean and the sample variance-covariance structure. Their formulation states that, in practice, the resulting imputations match the original periodic means and variances up to Monte Carlo error [2505.02008].

A later formalization introduces frequency-specific covariates $\tilde x_t^{(j)}$ and writes the augmented imputation model as
$$
y_t \sim \mathcal{N}\bigl(\alpha_0+\sum_j \alpha_j \tilde x_t^{(j)},\sigma^2\bigr).
$$
The stated objective is to ensure that the imputed values inherit the true underlying periodic structure of the data, especially annual, semiannual, and monthly cycles [2508.19535].

In the daily-temperature comparison reported by Ahmad et al., missing at random was imposed at 13% overall and 30% on weekends. Against standard Amelia II, Amelia II with VBPBB achieved $\mathrm{MAE}=1.326$ and $\mathrm{RMSE}=4.630$, compared with $\mathrm{MAE}=3.015$ and $\mathrm{RMSE}=10.353$ without VBPBB. The same study reports a reduction of MAE by 57%, a reduction of RMSE by 56%, and an increase in correlation with observed data from 0.862 to 0.972 [2505.02008].

| Method | MAE | RMSE |
|---|---:|---:|
| Amelia II with VBPBB | 1.326 | 4.630 |
| Amelia II without VBPBB | 3.015 | 10.353 |

A broader simulation assessment uses $T=6{,}000$ daily observations, annual, harmonic, and monthly signals, noise variances $\sigma^2\in\{0.1,4,10,25,100\}$, and MCAR missingness from 5% to 70%. Across all three signal stages, the VBPBB-enhanced approach consistently outperforms standard imputation. Reported gains are largest when multiple periodic components are retained, and under very low noise with 5% missingness the reductions in MAE and RMSE are about 93%; even at $\sigma^2=100$ and 70% missingness, improvements remain on the order of about 9%–16% depending on the stage [2508.19535].

## 6. Empirical record, extensions, and limitations

The empirical record supplied for VBPBB is dominated by comparisons of confidence-band width and component detectability against unfiltered periodic bootstrap procedures. Across these studies, the recurring result is that pre-filtering removes extraneous variance before resampling, producing much tighter confidence bands while preserving periodic correlation structure.

| Application | Reported significant components | Reported comparative result |
|---|---|---|
| US COVID-19 mortality [2403.06343] | Seasonal, 2nd–5th annual harmonics, weekly | PBB/VBPBB median CI-width ratios 3.49, 5.98, 6.41, 13.87, 32.93, and 5.47 |
| Manhattan PM$_{2.5}$ [2404.08738] | Semi-annual and weekly | GSBB CIs 11.00× wider for semi-annual and 15.87× wider for weekly |
| Turkey electricity [2404.03786] | Daily, weekly, annual-related components | Median CI-width ratios 1.64×, 4.61×, 3.65×, combined 1.58× |
| NY COVID-19 hospitalizations [2404.11006] | Annual, 3rd–5th harmonics | VBPBB CIs 1.64×, 5.26×, 6.55×, and 12.36× narrower |
| IHD hospitalizations in New York [2509.03710] | Weekly, annual, 2nd, 3rd, 5th, and 6th harmonics | GSBB-to-VBPBB band-width ratio 35.6 for annual and up to 179.1 for the fifth annual harmonic |

The method’s theoretical and practical limits are also explicit in the literature. VBPBB is designed for time series that contain narrowband cyclical components, or that can be represented as a sum of a small number of PC components plus residual noise. It depends on filter specification, can suffer from end effects, and requires sufficiently long series so that filtering does not remove too much data at the edges. Several papers note that poor choices of $m$, $k$, or bandwidth may under-filter, over-filter, or damp true peaks. One application summary states that VBPBB reduces to ordinary PBB in the extreme case $m=1$ or, more generally, when the passband is chosen wide enough to pass all frequencies [2402.03491; 2403.06343].

The bias analysis in later work makes these issues more precise. Under conditions including second-order periodic stationarity of each component, $b_j\to\infty$, $b_j/n\to0$, and weak dependence of the residual noise, the overall-mean bias and the pointwise-mean bias of the bootstrap distribution are both
$$
O\!\Bigl(\max_j \frac{b_j}{n}\Bigr).
$$
That study recommends centering each bandpass filter at $\omega_j=2\pi j/T$, choosing half-widths to exclude nearby harmonics, and balancing the bias-variance trade-off in block length. It also states that VBPBB will not reproduce nonperiodic trends or localized events unless they repeat each cycle [2509.08647].

Current future-direction statements in the supplied papers focus on automated selection of frequencies and bandwidths, extension beyond univariate MCAR settings, structural breaks, multivariate VBPBB exploiting cross-series periodic correlations, and parallelized filtering and blocked bootstrap schemes [2508.19535]. Together, these directions place VBPBB within a broader program of frequency-aware resampling for cyclostationary inference and temporally structured imputation.

Source: https://www.emergentmind.com/topics/variable-bandpass-periodic-block-bootstrap-vbpbb