• ベストアンサー

複数の時系列データのコヒーレンシイ計算

複数の時系列データのコヒーレンシイ計算をExcel、SCILAB、フリーソフト等で簡易的に行う方法を教えてください。位相のずれた複数時系列データを解析して、位相差を求めたいと思います。位相差を求めるのにコヒーレンシイが良いとの話があり、コヒーレンシイの計算方法を質問させて頂きました。只、最終目的は位相差を明確にすることですので、他の良い方法があれば、その方法でも結構です。尚、先般、「複数の時系列データの位相計算」で質問しました所、意味が不明とのことで、再度質問をさせて頂きました。

質問者が選んだベストアンサー

  • ベストアンサー
回答No.1

フリーの統計解析ソフトRはご存知と思います。Rのインストール・一般的使用例は多数解説があるので割愛します。例えば http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html コヒーレンスを計算する関数cohはRの基本的な関数の中にはないので、パッケージseewaveをインストールします。 1.Rのパッケージメニューで「パッケージのインストール」を選択 2.seewaveを選択してOKボタンを押す 3.library("seewave") と入力するとcohが使えるようになる。 二つの時系列をts1とts2とすれば > coh(ts1,ts2) でコヒーレンスが計算されます。seewaveについては http://sueur.jerome.perso.neuf.fr/WebPage_PapersPDF/Sueuretal_Bioacoustics_2008.pdf

fisherman-2011
質問者

お礼

早速のご回答ありがとうございました。統計解析用ソフトRについては、今回初めて知りましたが、webで調査しますと、色々の実例が紹介されている有用なソフトと認識しました。現在、windows/base/R-2.13.0-win.exe をインストールしましたが、途中文字化け等が現れている状態です。インストール並びに起動方法が分かり、ご提案の方法でコヒーレンシが計算できましたら、ご連絡いたします。まずは、ご回答のお礼まで。

fisherman-2011
質問者

補足

先般は有難うございました。その後、検討を重ね、ご指摘のcoh関数を設定できる様になりました。只、以下のエラーメッセージが現れて計算が進みません。 > coh(bd,tor) 以下にエラー inputw(wave = wave1, f = f) : 引数 "f" がありませんし、省略時既定値もありません 再度お手数をおかけしますが、宜しければこの対処法を教えてください。

全文を見る
すると、全ての回答が全文表示されます。

その他の回答 (1)

回答No.2

Rについては http://www.okada.jp.org/RWiki/ のQ&A(初級者コース)、Q&A(中級者コース)で質問されたらいかがでしょうか。、すみません私もRについて詳しくはありません。

fisherman-2011
質問者

お礼

色々と有難うございました。

全文を見る
すると、全ての回答が全文表示されます。

関連するQ&A

  • 時系列データの位相計算法

    複数の時系列データの位相計算をExcel、SCILAB、フリーソフト等で簡易的に行う方法を教えてください

  • ダウ平均の時系列データ

    普段、ヤフーの時系列データを使用してデータ解析などを行ったりしているのですが、最近ダウ平均の時系列データが必要になりました。 そこで質問なんですが、ダウ平均の時系列データが取れるサイトをご存じないでしょうか? もし、ご存じありましたら、HPご紹介をお願いできたらと思います。 よろしくお願いします。

  • scilabを使って時系列データをフーリエ変換するには?

    scilabを使ったフーリエ変換についての質問です。 あるデータ(メトロノームに合わせた指タッピングの時間間隔データでっす)を、フーリエ変換しなければいけないのですが分かりません。゜゜(´□`。)°゜。 具体的には、その時系列データ(512個あってexcelに収まってマス、A列に時間(5432ミリ秒とか)、B列に強さ(200g重とか))を、scilabを用いてフーリエ級数変換・・・スペクトル密度のグラフをつくる、な感じなことをするらしいのですが、プログラムがぜんぜんうまく書けません(今月末までつくらないといけないのに・・・)゜(゜≧□≦゜)゜ 生データの自己相関関数をつくって、フーリエ級数変換、高速フーリエ変換、グラフをプロット、のような流れなのかなぁ?と思うのですが、合ってますかぁ? もしよかったらサンプルプログラムを書いてくださると助かります(*^o^*) 先週から私なりにチャレンジしているんですが、エラーばっかりなんですぅ(涙) どうか、よろしくお願いします☆

  • 時系列の揃っていないデータのグラフについて

    添付画像のようなデータがあります。 同じ時系列で推移するデータですが、サンプリング周期の違い等により、ログとして残る時間は歯抜けデータとなります。 これをそのままグラフにすると、添付グラフのように、系列1,2,3で、時系列が揃わないグラフになってしまいます。 グラフ側の設定で時系列が揃うことが望ましいです。 良い方法をご存知のかたおられましたら、回答よろしくお願いいたします。

  • 時系列データを使ったパラメータ推定

    以前に、時系列データを使ったパラメータ推定について質問させて頂き、有効なソフトを教えて頂いたのですが、そのソフトが周りにないことがわかったうえに、購入する余裕もありません。 私が作った企業の長期利潤に関する式は線形なので、対数をとり、最悪手計算でもできるはずだと考えました。 しかし、時系列データ(80都市×15年)を見ると、例えば一つのx値に対して80都市×15年のy値があるのにパラメータの値が1つだけ推定されるという理論が理解できません。 まず、時系列データを処理して値を絞り、最小二乗法でパラメータを推定するのか、とりあえず80×15のパラメータを出して何か処理をするのかもわかりません。経済学、数学ともに専攻外の大学院生で、困っています。時系列データの処理のしかた、考え方をご存知の方、教えてください。お願いします。

  • エクセルで時系列グラフ

    エクセルで時系列データ(単純な数値)をグラフにしたいと考えています。(現在は折れ線グラフを使っています) ある時点からデータの計算方法が変わったために、データを非連続にしたいのですが方法が分かりません。(計算方法が変わる前と後で線をつなげない。ただし時間軸は継続させる) 系列を追加する方法では、新たなデータとして認識されてしまうため、上手くいきません。どうすればいいでしょうか? よろしくお願いします。 *項目軸は四半期データのため、例えば99年の第一四半期なら99Q1という形式にしています。

  • 時系列データをベクトルとして平均値ベクトルの差の検定?

    統計の問題です。 2群の時系列データを持っています。 (スタートから値が上昇し、ある時点から減少するような形) この場合、群間に差があるかどうかを検定する場合は、 時系列データから何か指標(例えば曲線下面積)を分散分析にかける方法がよく用いられているようです。 しかし、時系列データが持っているデータ情報量が上手く生かされていません。 そこで、 http://aoki2.si.gunma-u.ac.jp/lecture/Wilks/wilks2.html にある方法を変数の項目を時系列に変えて適用できないものかと思いました。 時系列データをベクトルとして扱うことは大丈夫なのでしょうか? どなたかアドバイスいただけると幸いです。

  • 2つの時系列データのコヒーレンスの計算について

    工学系の学生をしております。時系列データの処理について勉強しているのですが、コヒーレンスの計算が理解できないので、質問させていただきます。 2つの時系列x(t),y(t)のコヒーレンスを計算したいのですが、教科書(例えば、日野幹雄:スペクトル解析)によると、 coh(w)=|Sxy(w)|^2/( Sxx(w)・Syy(w) )  coh:x(t)とy(t)のコヒーレンス Sxy:x(t)とy(t)のクロススペクトル  Sxx:x(t)のパワースペクトル  Syy:y(t)のパワースペクトル となっており、コヒーレンスcohは0から1の間の値ととると書かれています。 なのですが、x(t),y(t)のフーリエスペクトルをX(w),Y(w)とすると、  Sxx(w)=X*(w)・X(w)/T = |X(w)|^2/T  Syy(w)=Y*(w)・Y(w)/T = |Y(w)|^2/T  Sxy(w)=X*(w)・Y(w)/T = |X(w)|・|Y(w)|/T・exp(-theta) (ここで、X*,Y*はX,Yの複素共役、thetaはクロススペクトルSxyの偏角) となるため、cohの式に入れると、  coh(w)=|Sxy(w)|^2/( Sxx(w)・Syy(w) ) = ( |X(w)|・|Y(w)| )^2 / ( |X(w)|^2・|Y(w)|^2 ) = 1 となってしまい、周波数によらず常に1になってしまうのですが、 上記の考え方は間違っていますでしょうか? 他にもいくつか教科書を読んでみたのですが、どうにも一人で解決出来ず、こちらに質問させていただきました。 よろしくお願いいたします。

  • 複数の異なる検査項目の時系列折れ線グラフ表示

    複数の検査項目を時系列の折れ線グラフ表示をしたいと思います。 その際の基準値の表示(帯)がうまくいきません。 複数の異なる検査項目の結果を時系列のグラフに表示する際、 値がバラバラなので、見づらい時系列表示になってしまいます。 これを解消するために、下限値、上限値を設定してグラフを 見やすくしようと思います。 ですがどう計算式をたてればよいかわかりません。 説明不足かと思いますが、アドバイスをお願いいします。

  • インデックスの時系列データが知りたい

    「MSCI KOKUSAI INDEX」と「シティグループ世界国債インデックス(除く日本)」の、 分配金込の時系列のデータを探しています。 出来るだけ過去からのデータを、希望しています。 円ベースでもドルベースでも構いません。 入手方法をご存知でしたら教えていただけないでしょうか?