(http://www.ocean.caos.tohoku.ac.jp/~merge/sstbinary/actvalbm.cgi)(分解能1/20度)を読むモジュールである。ファイル名としてはgzipされたファイルを指定でき、解凍する必要はない。(ftpのURLも直接指定できるがネットワークの負荷を考えること。)メインプログラムとして走らせた時は領域一部をmatplotlibでプロットする。data2010年5月26日水曜日Rでwavelet前回までのようにRにNINO3のデータを渡し、今回は連続wavelet解析を行ってみる。Wavelet解析にはwmtsaパッケージを用いる。以下がプログラムである。プロットする時、series=TRUEでもとの時系列もプロットする。60.3333333333333x軸は年。Y軸は0が1年(2^0)、2が4年(2^2)に対応する。前回に続いてRにNINO3のデータを渡し、今回は自己相関を求めてみる。type="partial"も求めることができる。%locals())図の横の点線は95%の信頼区間をしめす(これより絶対値の小さい相関は95%の信頼区間で無相関ではないとは言えない)。これを90%などにしたいときはci=0.9などとする。以前につくったpythonで取り込むNINO監視領域のデータを、統計ソフトRに持って行き、時系列データとする。Rに取り込むとRを使って統計解析ができるから便利だ。Pythonの世界からRの世界に持っていくにはrpy2を用いる。NINO3データを渡して、Rではtsを使って時系列オブジェクトを作り、plot.tsでプロットする。