[論文レビュー] MIARMA: An information preserving method for filling gaps in time series. Application to CoRoT light curves

J. Pascual-Granado, R. Garrido+1 more|arXiv (Cornell University)|Oct 3, 2014

Stellar, planetary, and galactic studies参考文献 21被引用数 8

ひとこと要約

MIARMA は、星の光曲線の周波数スペクトルを保持するための新しいギャップ補完手法であり、前向き・後向きの自己回帰移動平均（ARMA）モデルを用いて、特に CoRoT データにおける時間系列の周波数成分を正確に再現する。線形補間よりも優れており、スペクトルに偽のアリエスを生じさせず、特に低周波数パルサーである Be 星や δ Scuti 星において、従来検出できなかった微細構造を明らかにする。

ABSTRACT

The method here presented intends to minimize the effect of the gaps in the power spectra by gap-filling preserving the original information, that is, in the case of asteroseismology, the stellar oscillation frequency content. We make use of a forward-backward predictor based on autoregressive moving average modelling (ARMA) in the time domain. The method MIARMA is particularly suitable for replacing invalid data such as those present in the light curves of the CoRoT satellite due to the pass through the South Atlantic Anomaly, and eventually for the data gathered by the NASA planet hunter Kepler. We select a sample of stars from the ultra-precise photometry collected by the asteroseismic camera on board the CoRoT satellite: the δ Scuti star HD 174966, showing periodic variations of the same order as the CoRoT observational window, the Be star HD 51193, showing longer time variations, and the solar-like HD 49933, with rapid time variations. We showed that in some cases linear interpolations are less reliable to what was believed. In particular: the power spectrum of HD 174966 is clearly aliased when this interpolation is used for filling the gaps; the light curve of HD 51193 presents a much more aliased spectrum than expected for a low frequency harmonic signal; and finally, although the linear interpolation does not affect noticeably the power spectrum of the CoRoT light curve of the solar-like star HD 49933, the ARMA interpolation showed rapid variations previously unidentified that ARMA interprets as a signal. In any case, the ARMA interpolation method provides a cleaner power spectrum, that is, less contaminated by spurious frequencies. In conclusion, MIARMA appears to be a suitable method for filling gaps in the light curves of pulsating stars observed by CoRoT since the method preserves their frequency content, which is a necessary condition for asteroseismic studies.

研究の動機と目的

• ギャップによって生じる偽の周波数アリエスの問題、特に CoRoT 光曲線における問題を解決すること。
• 低周波数パルサーにおいて、元の周波数成分を保持できない線形補間の限界を克服すること。
• 信頼できる星震動解析が可能な真の信号成分を保持するギャップ補完手法を開発すること。
• δ Scuti 星、Be 星、太陽型振動星を含む多様な星の種類に適用可能な手法を確保すること。
• 既存のインpaint技術が誤ってノイズを信号として再構築する可能性があるのを防ぐ、堅牢な代替手法を提供すること。

提案手法

• 自己回帰移動平均（ARMA）モデルに基づく前向き・後向き予測子を用いて、時間系列の欠損値を推定する。
• 有効なデータセグメントを ARMA プロセスとしてモデル化し、解析的関数形を仮定せずに時間的ダイナミクスを捉える。
• 元の信号のスペクトル特性を保持しながら、繰り返し適用してギャップを埋める。
• モデルの適合度を評価し、ノイズと真の信号を区別するためにアカイケ情報量基準（AIC）を用いる。
• 人工的な白色ノイズ時系列データを用いて検証し、人工的なパワーを導入しないことを確認する。
• 信号の忠実度とスペクトル忠実度を評価するために、k-inpainting および線形補間と比較する。

実験結果

リサーチクエスチョン

• RQ1 線形補間は、低周波数パルサーのギャップのある光曲線において、真の周波数成分をどれほど不正確に保っているのか？
• RQ2 ARMA を用いたギャップ補完は、CoRoT の観測窓関数によって生じる偽のアリエスを排除できるか？
• RQ3 MIARMA は、従来の補間法で隠れていた、光曲線における以前に検出されなかった微細構造を明らかにできるか？
• RQ4 合成データにおける白色ノイズの統計的性質を保つ観点で、MIARMA は k-inpainting よりも優れているか？
• RQ5 ARMA 補間は、ギャップのある時系列において、内在的信号とノイズを信頼性高く区別できるか？

主な発見

• 線形補間は、δ Scuti 星や Be 星（例：HD 174966 や HD 51193）のパワー・スペクトルに顕著なスペクトルアリエスを導入し、真の振動周波数を歪めている。
• 太陽型星 HD 49933 では、線形補間はパワー・スペクトルに顕著な影響を及ぼさないが、MIARMA は以前に識別されなかった急速な変動を明らかにし、それが真の信号であると解釈できる。
• MIARMA は、星の種類に関係なく、パワー・スペクトルからのすべてのアリエスを効果的に除去する。長周期変光星に対しても、線形補間が失敗する状況でも同様に有効である。
• 完全にランダムな白色ノイズ時系列に適用した場合、MIARMA は標準偏差 σ = 0.0753 の補間値を生成し、ゼロに近い値であるため、人工的な信号再構築は行われていない。
• これに対して、k-inpainting アルゴリズムは、同一の白色ノイズデータに対して標準偏差 σ = 0.2730 を生成し、ノイズを信号として再構築しているため、信号保持の原則に反する。
• アカイケ情報量基準（AIC）は、MIARMA で補間された HD 49933 の光曲線における微細構造が白色ノイズではないことを確認しており、物理的信号としての解釈を支持する。