[论文解读] Simulation Study Comparing Two Tests of Second-order Stationarity and Confidence Intervals for Localized Autocovariance

Guy P. Nason|arXiv (Cornell University)|Mar 21, 2016

Image and Signal Denoising Methods参考文献 2被引用 5

一句话总结

本文通过模拟比较了两种二阶平稳性检验方法——Priestley-Subba Rao (PSR) 检验与基于哈尔小波的演化谱检验 (HWTOS)，结果表明 HWTOS 在重尾分布和近单位根过程中的表现更优。研究进一步验证了局部自协方差估计的近似正态性，并提供了使用 locits 包实现可复现性的 R 命令。

ABSTRACT

This report compares two tests of second-order stationarity through simulation. It also provides several examples of localised autocovariances and their approximate confidence intervals on different real and simulated data sets. An empirical verification of an asymptotic Gaussianity result is also demonstrated. The commands use to produce figures in a companion paper are also described.

研究动机与目标

• 比较两种平稳性检验方法（Priestley-Subba Rao (PSR) 与基于演化谱的哈尔小波检验 (HWTOS)）的实证性能。
• 评估非平稳时间序列中局部自协方差置信区间的准确性。
• 通过实证方法验证在局部平稳条件下，局部自协方差估计量的渐近高斯性。
• 提供可复现的 R 代码，使用 locits 包重现配套论文中的图表。
• 评估波动聚集效应及参数取值对检验势与大小的影响，特别是在高绝对值系数的 AR(1) 过程中。

提出的方法

• 从七个平稳时间序列模型（S1–S7）中各模拟 1,000 次实现，每次样本量 T=512，用于估计经验检验大小。
• 通过 R 语言 fractal 包中的 `stationarity` 函数实施 PSR 检验，通过 locits 包中的 `hwtos2` 函数实施 HWTOS 检验。
• 在 HWTOS 检验中应用多重比较校正（FDR 与 Bonferroni 方法），以控制多重检验下的第一类错误率。
• 使用 `Rvarlacv` 函数在特定时间点估计局部自协方差（LACV），采用不同窗口大小与滞后阶数。
• 基于模拟中估计值的样本均值与样本方差，构建 LACV 的近似置信区间。
• 可视化 LACV 估计的抽样分布，并与理论正态密度进行比较，以评估其渐近正态性。

实验结果

研究问题

• RQ1 在不同平稳模型下，PSR 与 HWTOS 检验的经验大小如何比较，特别是在重尾或近单位根过程中？
• RQ2 当应用于高绝对值系数的 AR(1) 过程（如 α = ±0.9）时，HWTOS 检验是否能保持适当的 I 类错误率？
• RQ3 局部自协方差估计在多大程度上服从正态分布？不同样本大小下置信区间的性能如何？
• RQ4 AR(1) 过程中 α = -0.9 时的波动聚集效应如何影响平稳性检验的解释？
• RQ5 哪些 R 命令与工作流能够使用 locits 包复现配套论文中的关键图表？

主要发现

• HWTOS 检验比 PSR 更为保守，经验大小通常低于名义显著性水平 5%，尤其在轻尾噪声下更为明显。
• 对于 α = 0.9 的 AR(1) 过程，PSR 检验经验大小显著膨胀（12.4%），而 HWTOS 在多重比较校正下保持接近名义水平（4.0–4.7%）。
• 对于 α = -0.9 的情况，HWTOS 的经验大小偏高（Bonferroni 校正下为 20.3%），表明对波动聚集敏感，但仅有 17% 的拒绝是由单一显著小波系数导致。
• PSR 在重尾数据中表现较差，约 60% 的情况下拒绝原假设，而 HWTOS 保持更好控制（经验大小约 10–20%）。
• 对于近单位根 AR(2) 过程（S7），PSR 的经验大小高达 23.9%，而 HWTOS 表现更优（7.3–7.4%），但仍高于名义水平。
• 随着样本量增加（T=512 至 4096），局部自协方差估计的抽样分布与正态分布高度吻合，密度估计与理论正态曲线一致。