[論文レビュー] Quantile-Based Spectral Analysis in an Object-Oriented Framework and a Reference Implementation in R: The quantspec Package

Tobias Kley|arXiv (Cornell University)|Aug 28, 2014

Control Systems and Identification参考文献 10被引用数 10

ひとこと要約

本稿では、時系列の分位数ベーススペクトル解析のためのオブジェクト指向フレームワークである quantspec Rパッケージを紹介する。このフレームワークにより、ロバストで、正規分布でなく、2次モーメントに依存しない解析が可能となる。コプシラとランクを用いた分位数パワースペクトル密度推定器の体系的実装が行われ、スムージングにより一貫性が保たれるとともに性能が向上し、従来のL2手法を超える包括的でオープンソースのツールボックスを提供する。

ABSTRACT

Quantile-based approaches to the spectral analysis of time series have recently attracted a lot of attention. Despite a growing literature that contains various estimation proposals, no systematic methods for computing the new estimators are available to date. This paper contains two main contributions. First, an extensible framework for quantile-based spectral analysis of time series is developed and documented using object-oriented models. A comprehensive, open source, reference implementation of this framework, the R package quantspec, was recently contributed to CRAN by the author of this paper. The second contribution of the present paper is to provide a detailed tutorial, with worked examples, to this R package. A reader who is already familiar with quantile-based spectral analysis and whose primary interest is not the design of the quantspec package, but how to use it, can read the tutorial and worked examples (Sections 3 and 4) independently.

研究の動機と目的

• 時系列の分位数ベーススペクトル解析のための体系的で拡張可能かつオブジェクト指向のフレームワークの開発。
• 実用的用途を想定した、分位数スペクトル推定器の包括的でオープンソースのR実装（quantspecパッケージ）の提供。
• 最近提唱された分位数ベーススペクトル推定器のための体系的計算手法の不足を解消すること。
• 重い尾を持つ時系列、外れ値、非正規的依存構造、さらには2次モーメントを持たないプロセスを含む、ロバストな解析を可能にすること。
• 将来の拡張を支援するための高度な機能、例えばバンド幅選択、多スペクトル推定、グラフィカル可視化のサポート。

提案手法

• 分位数および確率レベルの依存性を定義するため、ラプラスおよびコプシラ交差自己共分散カーネルに基づくフレームワークが構築される。それぞれ、Cov(I{X_t ≤ q1}, I{X_{t-k} ≤ q2}) および Cov(I{F(X_t) ≤ τ1}, I{F(X_{t-k}) ≤ τ2}) で表される。
• フーリエ変換を用いてスペクトル表現を導出する：ラプラススペクトル密度カーネル 𝔽_{q1,q2}(ω) = (1/2π) Σ_k γ_k(q1,q2) e^{-i k ω} およびコプシラスペクトル密度カーネル 𝔽_{τ1,τ2}(ω) = (1/2π) Σ_k γ_k^U(τ1,τ2) e^{-i k ω}。
• quantspecパッケージは、ユーザーが指定可能なバンド幅を用いたEpanechnikovカーネルスムージングを用いて、分位数パワースペクトル（quantilePG）およびそのスムージング版（smoothedPG）を実装する。
• シミュレーション研究や性能評価のための周波数ドメイン値を抽出可能な柔軟なアクセサメソッド（例：getValues）をサポートする。
• モンテカルロ法を用いたシミュレーション研究では、繰り返し実現（例：R=50000）により、バイアスと分散の評価のための平均二乗誤差のルート積分（RIMSE）を計算する。
• 拡張性を重視して設計されており、GitHubのイシュー追跡およびトピックブランチを活用して、将来的な拡張（例：局所定常プロセス、高次多スペクトル）をサポートする。

実験結果

リサーチクエスチョン

• RQ1 時系列の分位数ベーススペクトル解析のための、オブジェクト指向的かつ拡張可能なフレームワークを体系的に設計するにはどうすればよいか？
• RQ2 バイアスと分散の観点から、スムージングされた分位数パワースペクトルと未スムージング版の性能はどのように比較されるか？
• RQ3 quantspecパッケージは、非正規的、重い尾、または2次モーメントを持たないプロセスのコプシラスペクトル密度を信頼性高く推定できるか？
• RQ4 有限標本特性を評価するためのシミュレーション研究をどのように構築すればよいか？
• RQ5 将来的な拡張（例：データ駆動型バンド幅選択、多スペクトル推定）を可能にするための主要な設計原則は何か？

主な発見

• quantspecパッケージは、Rにおける分位数ベーススペクトル解析の包括的でオープンソースかつ拡張可能なフレームワークを成功裏に実装した。
• スムージングされた分位数パワースペクトルは、未スムージング版に比べて顕著に低い平均二乗誤差のルート積分（RIMSE）を示し、さまざまな分位数ペアにおいてRIMSEが約0.03から0.005に低下した。
• コプシラスペクトル密度カーネルは単調変換に対して不変であるため、系列的依存性と周辺分布の特徴を分離可能である。
• フレームワークは2次モーメントを要件としないため、重い尾や非対称なプロセスに対してもロバストな解析を可能にする。
• シミュレーション結果により推定器の一貫性が確認され、標本サイズ（例：N=2^9 から N=2^12）および反復回数（例：R=100 から R=50000）が増加するにつれてRIMSEが減少した。
• パッケージはプロダクション環境で使用可能であり、拡張性も高く、GitHubで開発中の計画された機能（例：データ駆動型バンド幅選択、高次多スペクトル推定）がすでに開発段階にある。