[論文レビュー] Improved AGN light curve analysis with the z-transformed discrete correlation function

Tal Alexander|arXiv (Cornell University)|Feb 6, 2013

Remote Sensing and Land Use参考文献 5被引用数 52

ひとこと要約

本論文では、AGN研究で一般的なスパarsely、不規則にサンプリングされた光曲線における相互相関関数の推定に向け、等頻度バインニングとフィッシャーのz変換を用いることでバイアスを低減し、誤差推定を可能にする、ロバストな手法としてz変換を施した離散相互相関関数（ZDCF）を導入する。ZDCFは補間法や従来のDCF手法を凌駆し、特にデータポイントが少ない場合に優れた性能を示し、実際のAGN光曲線における信頼性の高い遅延推定と変動性解析を可能にする。

ABSTRACT

The cross-correlation function (CCF) is commonly employed in the study of AGN, where it is used to probe the structure of the broad line region by line reverberation, to study the continuum emission mechanism by correlating multi-waveband light curves and to seek correlations between the variability and other AGN properties. The z -transformed discrete correlation function (ZDCF) is a new method for estimating the CCF of sparse, unevenly sampled light curves. Unlike the commonly used interpolation method, it does not assume that the light curves are smooth and it does provide errors on its estimates. The ZDCF corrects several biases of the discrete correlation function method of Edelson & Krolik (1988) by using equal population binning and Fisher's z -transform. These lead to a more robust and powerful method of estimating the CCF of sparse light curves of as few as 12 points. Two examples of light curve analysis with the ZDCF are presented. 1) The ZDCF estimate of the auto-correlation function is used to uncover a correlation between AGN magnitude and variability time scale in a small simulated sample of very sparse and irregularly sampled light curves. 2) A maximum likelihood function for the ZDCF peak location is used to estimate the time-lag between two light curves. Fortran 77 and 95 code implementations of the ZDCF and the maximum likelihood peak location (PLIKE) algorithms are freely available (see http://www.weizmann.ac.il/weizsites/tal/research/software/).

研究の動機と目的

• スパarselyで不規則にサンプリングされたAGN光曲線の解析において、補間法や離散相互相関関数（DCF）手法の限界を解決すること。
• 特に不規則なサンプリングと誤差推定の欠如に起因する、元のDCF手法におけるバイアスを克服すること。
• 12ポイント程度の少数データでも、信頼性のある相関推定と誤差区間を提供する統計的にロバストな手法を開発すること。
• ZDCFのピーク位置を最尤法で特定することで、AGN光曲線間の正確な遅延推定を可能にすること。
• 実世界のスパarselyなデータを用いた、AGNの変動性、連続スペクトル発生機構、リバーブレーションマッピングのより信頼性の高い解析を促進すること。

提案手法

• 時間遅れ分布に対して等頻度バインニングを適用し、各バインに同じ数のデータポイントが含まれるようにする（等時間幅ではなく）。
• フィッシャーのz変換を用いて相関係数推定の分散を安定化させ、正規性を向上させ、バイアスを低減する。
• 光曲線の重複するセグメントにおける共分散と分散を用いて、各時間遅れにおける相関係数を推定する。
• 各バイン内での時間遅れの分布に基づき、ZDCF推定値の誤差区間を計算し、非対称信頼区間を提供する。
• ZDCFピーク位置を最尤関数を用いて定義することで、時間遅れ推定の精度を向上させる。
• FORTRAN 77およびFORTRAN 95で実装し、コミュニティ利用可能な公開コードを提供する。

実験結果

リサーチクエスチョン

• RQ1 スパarselyで不規則にサンプリングされたAGN光曲線において、相互相関関数を信頼性高く推定する方法は何か？
• RQ2 従来のDCF手法に内在する統計的バイアスは何か？ そして、z変換と等頻度バインニングを用いることで、それらのバイアスはどのように是正できるか？
• RQ3 ZDCF手法は、補間法やDCFが失敗する状況でも、相関係数の正確な誤差推定を提供できるか？
• RQ4 ZDCFは、既存の手法と比較して、AGN光曲線間の時間遅れ推定をどの程度改善するか？
• RQ5 非常にスパarselyなシミュレーテッド光曲線において、ZDCFはAGNの明るさと変動時標の相関を検出できるか？

主な発見

• 等頻度バインニングを用いることで、等時間バインニングに比べて相関推定のバイアスが低減され、特にスパarselyな光曲線において性能が向上する。
• フィッシャーのz変換を用いることで、相関係数推定の分散が安定化し、より信頼性が高く正規分布に近い標本分布が得られる。
• ZDCFは、各バインに使用された時間遅れの実際の分布を反映した非対称誤差区間を提供するため、対称誤差区間よりも不確実性をより正確に表現できる。
• 12ポイント程度のシミュレーションでも、ZDCFは真の相関を的確に回復し、AGNの明るさと変動時標の間の相関を明らかにする。
• ZDCFのピーク位置を最尤法で特定することで、単純なバインニングに比べ、より正確で頑健な時間遅れ推定が可能になる。
• 特にS/Nが低く、サンプリングが疎な状況下でも、補間法や標準DCFに比べ、統計的信頼性と誤差推定の両面でZDCFが優れている。