[論文レビュー] Fitting and Comparison of Models of Radio Spectra

Bojan Nikolic|ArXiv.org|Dec 11, 2009

Gaussian Processes and Bayesian Inference参考文献 6被引用数 3

ひとこと要約

本論文は、ネストド・サブセット法を用いたベイジアンスペクトルフィッティングフレームワークを提示し、放射線スペクトルモデルを客観的に比較可能にし、完全な事後分布推定と証拠に基づくモデル選択を可能にする。近傍銀河に適用したところ、NGC 7331に対しては連続的注入に自己吸収を含むモデルが、高い証拠により強く支持された一方、他の銀河に対しては単純なモデルで十分であり、従来の最尤推定法に比べてより高いロバスト性を示した。

ABSTRACT

I describe an approach to fitting and comparison of radio spectra based on Bayesian analysis and realised using a new implementation of the nested sampling algorithm. Such an approach improves on the commonly used maximum-likelihood fitting of radio spectra by allowing objective model selection, calculation of the full probability distributions of the model parameters and provides a natural mechanism for including information other than the measured spectra through priors. In this paper I cover the theoretical background, the algorithms used and the implementation details of the computer code. I also briefly illustrate the method with some previously published data for three near-by galaxies. In forthcoming papers we will present the results of applying this analysis larger data sets, including some new observations, and the physical conclusions that can be made. The computer code as well as the overall approach described here may also be useful for analysis of other multi-chromatic broad-band observations and possibly also photometric redshift estimation. All of the code is publicly available, licensed under the GNU General Public License, at http://www.mrao.cam.ac.uk/~bn204/galevol/speca/index.html

研究の動機と目的

• 希な、ノイズの多い、非ガウス的なデータが存在する状況下で、統計的に厳密かつ客観的な方法で電波スペクトルモデルのフィッティングと比較を実現すること。
• 従来の最尤推定法の限界、例えば客観的なモデル選択ができないこと、不完全な誤差評価があることへの対処。
• 非ガウス的かつ相関のある不確実性を含む、モデルパラメータの完全な事後分布の提供。
• 物理的事前分布を分析に自然に統合し、外部制約を用いたパラメータ推定を改善すること。
• 一般用途で利用可能で、多波長および光度赤方偏移解析に応用可能な、公開済みの汎用ツールの開発。

提案手法

• ベイズの定理を用いて、データと仮説が与えられたもとでのモデルパラメータの事後分布を計算するベイズ推論を採用。
• ベイズ的証拠 Z と完全な事後分布を効率的に計算するために、ネストド・サブセット法（Skilling, 2006）を採用。
• シンクロトロン、熱的、自己吸収放射の解析的表現を用いて電波スペクトルをモデル化し、カスタマイズ可能な事前分布を導入。
• 観測された周波数におけるフラックス密度と、モデルとパラメータから予測されたフラックス密度を比較して尤度を計算。
• ファン図を用いて、事後分布全体にわたる予測スペクトルの性能を可視化。
• 公開済みのオープンソースコードとして実装し、GNU GPLライセンスのもとで、電波および多バンド天文学への広範な応用を目的としている。

実験結果

リサーチクエスチョン

• RQ1 近傍銀河の観測フラックス密度を最もよく説明するのは、べき乗則、連続的注入、または自己吸収を含むどのスペクトルモデルか？
• RQ2 データが希であり、誤差が非ガウス的である状況下で、どのようにしてモデル選択を客観的かつ統計的に厳密に行えるか？
• RQ3 物理的事前分布は、スペクトルフィッティングにおけるパラメータ推定の精度を向上させ、パラメータのデグレネシーをどれほど軽減できるか？
• RQ4 連続的注入にシンクロトロン自己吸収を含む複雑なモデルが、単純な代替モデルよりもベイズ的証拠によって定量的に優位に評価されるか？
• RQ5 ファン図は、モデルフィットの質とパラメータ不確実性を直感的かつ情報豊かに可視化するのにどのように役立つか？

主な発見

• NGC 7331では、連続的注入にシンクロトロン自己吸収を含む（CI+SSA）モデルのベイズ的証拠値が、他のモデルと比べて数個のオーダーも高いことから、強く支持された。
• NGC 7331におけるCI+SSAモデルのパラメータの周辺事後分布は、特にスペクトル指数と折り返し周波数において顕著な非ガウス性を示した。
• NGC 3627では、自己吸収を含まない連続的注入モデルが好まれており、データには低周波数での自己吸収の兆候は認められなかった。
• 最初に分析した銀河ではべき乗則モデルが最良のフィットを示し、4つのテストモデルの中で最も高いベイズ的証拠を示した。
• 自己吸収を含むモデル（例：30 MHzで低周波数の折り返しを含む）は、事前分布による制限により、可能な折り返し周波数の範囲が制限され、結果として性能が劣った。
• ファン図は、事後分布全体にわたる予測スペクトルの範囲を効果的に可視化し、モデルがデータの特徴を再現できない領域を明らかにした。