Correlation Analysis between environmental channels

Projects

• Hot correlation map
• Investigation into correlators
  • finding linear correlation
  • finding non-linear correlation

Project description

KAGRAには主チャンネル以外に環境チャンネルが1万程度ある。チャンネルデータの相関を取ることで効率的にノイズ源の特定を行うことを目的とする。 KAGRAには重力波チャンネル以外に環境チャンネルが1万程度の補助チャンネルの設置が予定されている。 本研究ではそれらの補助チャンネル同士の相関を取ることによって、補助チャンネルや重力波チャンネルに混入したノイズ源の特定ならびに混入したノイズがチャンネル間を伝搬メカニズムするメカニズムを明らかにすることを目的とする。

ノイズが伝搬するメカニズムは一般には線形または非線形メカニズムが考えられるが、過去の重力波検出器診断(detector characterization)では主に線形メカニズムの特定に重きを置いていた。本研究ではさらに踏み込んで、非線形なメカニズムをも明らかにできるような相関解析を目指している。 それらのメカニズムを解き明かすための手段として、ピアソンの積率相関係数とMaximum Information Coefficientを用いる。前者は線形的な相関を明らかにし、後者は線形以外の相関をも明らかにすることができる。

Project schedule and Task List

1タスクの定義は1週間程度で終わらせることができる研究

他研究グループのDetector Characterizationによる先行研究の把握

• 以下の論文を読む(1タスク)

  1 http://arxiv.org/abs/1111.4516

  2 http://arxiv.org/abs/1107.2948

  3 http://arxiv.org/abs/1403.1431

Hot correlation mapの作成

• テストデータの作成, Haskell上でのデータのIO(1タスク)
• ピアソンの積率相関係数の計算を行う関数の作成(1タスク)
• 2つのチャンネルデータを読み込み、あるデータ長Nについてデータの相関係数を求める。 チャンネル間に混入したノイズが同時刻だけ限らないので、相関係数を計算する開始時刻をずらして計算を行い、2つのデータの中で相関係数の最大値を求める(2タスク)
• HTMLを用いて表形式で相関マップ(Hot correlation map)の作成する。横軸と縦軸がchannel名(データの名前)、各channel名がクロスするセルに相関係数の値を書き、わかりやすさのためにその値の大きさを色合いで表現する (1タスク)
• 上項目をHaskellの関数として作成する、入力 相関値のテーブル、出力はhtmlソース(1タスク, ~7/14)
• それぞれのデータがガウシアンの場合は、ピアソンの積率相関係数の分布はデータ数をnとして自由度(n-2)のt分布に従うことが知られている。しかしデータがガウシアンでない場合の分布は計算されたピアソンの積率相関係数の分布は簡単にはモデル化することができないと予想される。そこで、過去の計算した相関値からcumlativeには何%の位置なるのかを計算して、評価したい。このような手法を用いることで、例えばchannel1とchannel2の相関値R_12とchannel1とchannel3の相関値R_13は直接は比較することはできないが、全体の何%に所属しているかで相関があるorないを評価することができる。

  • -> 間野さんより、permutation testを用いてはどうかというとともに、C言語で書かれたサンプルコードをいただいた。

• permutation testについては理解できたのでこれをHaskellで実装する(1タスク)

混入したノイズの同定

• upconverted noise に関する先行研究の理解 http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-20-8-8329

* 混入したノイズをトイモデルとして数式化する(1タスク) * たたみこみを用いて上項目のノイズデータを時系列で生成する(1タスク) * 生成したノイズに対して、相関マップを計算する。その統計量の振る舞いから、ピアソンの積率相関係数以外にどのような解析手法を用いると、効果的に混入したノイズをデータの中から特定できるかを考察する。既存の統計量では評価できない場合は、新しい統計量を定義する。(3タスク)*/ /*非ガウスノイズが混入した場合についても同様の評価を行う(時間と相談)*/

• ピアソンの積率相関係数以外に、相関を調べる統計量の導入を検討する(1タスク)

シミュレーションノイズの生成

* (2タスク) Virgoではエンドミラーを透過した散乱光による非線形ノイズの論文が報告されているので、そのモデルの理解、モデルに沿って非線形ノイズの例となるシミュレーションデータの生成 http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-20-8-8329

* (2タスク) 生成したシミュレーションデータをMIC, peasonを用いて解析を行う 検出の可否、統計量の分布について調べる。

• => それぞれの手法の相関の度合いを表す検出統計量は単純には比較することができないため、それぞれの手法のreceiver operating characteristic(ROC) curveの比較を行った。

* (2タスク) 論文執筆 (イントロダクション) と同時に結果についての議論、確認

• => 毎日、端山さんと昼食をとりながらの議論を行い研究を進めた。

* (2タスク) 論文執筆 (シミュレーション、結果、まとめ)

• => 執筆した原稿をPRDに投稿した(2016年2月16日)

* (1タスク) レフリーコメントへの対応方針について議論

• => レフリーからコメントを受け取った(2016年3月1日) 実際の重力波検出器から得られたデータに対してMICを適用する、またはより精密なシミュレーションに基づいてノイズデータを生成した方ががよいとのレフリーコメントであった。 今回はより精密なシミュレーションに基づいてノイズを再生成し、相関解析を行う方針を取る。

* (1タスク) Virgo検出器のnoise budgetについて調査する

• => The Virgo Sensitivity Curve にVirgo検出器のnoise budgetがまとめられていた。 それぞれのノイズソースのスペクトルもtextデータとして公開されているのでこれを用いる。 ラインノイズについては単純なガウシアンノイズではないことが知られているため、今回はViolin noiseといったラインノイズはノイズ生成には用いていない。

* (1タスク) noise budgetに沿ったノイズデータのシミュレーション生成

• => 上記のテキストデータをノイズスペクトルの絶対値とし、ガウス乱数を振り、実数部と虚数成分の生成を行った。 それを逆フーリエ変換することで上記noise bugetに沿った時系列ノイズデータの生成を行った。

* (2タスク) 上記ノイズデータにアップコンバージョンノイズを埋め込み、相関解析(Pearson, MIC)を用いて解析する。

• 時系列グラフ、スペクトル、データのホワイトニング、ホワイトニングしたスペクトル、散布図、receiver operating characteristic(ROC) curveの確認をする。

  =>

* (1タスク) レスポンスレターの準備、修正後の原稿の回覧、再投稿

• => 共著者(端山さん、間野さん、Didierさん、神田さん)に回覧する。また、干渉計の散乱光のプロフェッショナルでもある阿久津さんにも原稿に目を通していただき有益なコメントをいただいた。 (2016年7月29日)PRDに再投稿した。

Analysis method

始めにピアソンの積率相関係数を用いて、1つのデータについて解析を行えるようにする。 その後、他の統計量の導入を考えている。その際、Haskellのソースコードでは解析に用いる手法の引数部分を書き換えるだけで良いように、プログラムの開発を行ってゆく。 他の解析手法として、現在検討しているのは、線形相関以外の相関をも調べられるMaximum Information Coefficient(MIC)などである。

(MICの論文orWEBページ) http://www.slideshare.net/logics-of-blue/mic-31810194

最終的には相関統計量だけではなく散布図(scatter plot)も同時に目で確認する必要がある。

Results

Presentation

• Hirotaka Yuzurihara, Kazuhiro Hayama, Shuhei Mano, Didier Verkindt and Nobuyuki Kand, "Unveiling linearly and non-linearly correlated signals between gravitational wave detectors and environmental monitors" Physical Review D(Under article submission) draft of manuscript

• Hirotaka Yuzurihara, Kazuhiro Hayama, Shuhei Mano, Didier Verkindt and Nobuyuki Kand, "A method for non-linear correlated noise in gravitational waves detectors", Hyannis, USA, 17 June 2016 web page slide

• 譲原浩貴他, "KAGRA Detector Characterization : KAGRAテータの重力波探索に向けた補助信号間に現れる相関した特徴の新しい発見手法の確立", 日本物理学会 2014年秋季大会, 講演番号 20aSB-11, 佐賀大学, 2014年9月20日 スライド

• Hirotaka Yuzurihara, Kazuhiro Hayama, Shuhei Mano, Takahiro Yamamoto, Atsushi Nishizawa and Nobuyuki Kanda, "Developing a novel method for finding correlated features between control, environment channels, KAGRA detector characterization", 25 November 2014 スライド web page

Document