問題タブ [biwavelet]

この質問はかなり漠然としていますが、R で biwavelet パッケージを使用して成功した人はいますか? 次のコードがあります。

そしてそれはエラーを返します：

image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)), のエラー: 'x' と 'y' の値は有限であり、欠落していない必要があります

私のデータには欠損値が含まれておらず、すべて有限であるため、このエラーを理解できません。

例を再現しようとすると：

面倒くさいので、このパッケージでデータが失敗する理由がわかりません。上記の例のデータは、私のデータ、つまり「マトリックス」と同じクラスです。ヘルプやアドバイスをいただければ幸いです。

変更されました：

次のコード例が機能しない理由を誰かが提案できますか:

次のようなエラーが生成されます。

この問題を解決するにはどうすればよいですか?

追加：

Gavin Simpsons の回答への回答: 1 つの周波数のみを含むようにデータを保持し、時間ベクトルを変更すると、コードは正常に動作します。

次の例では:

特定の時系列の日周サイクルを見ている場合、使用可能な最低の解像度は何ですか? この例では現在、1 時間間隔で測定されたデータを使用しているため、1 日に 24 のデータ ポイントがあります。これは、2 時間、3 時間、さらには 4 時間の時間ステップを使用した場合と同じくらい重要でしょうか?

また、データの解像度を下げたときに結果の有意性が変化するかどうかを計算する方法はありますか

base グラフィックと ggplot グラフィックを組み合わせた図を生成したいと思います。次のコードは、R のベース プロット関数を使用した図を示しています。

生成するもの

これらのパネルのほとんどは、レポートに含めるのに十分に見えます。ただし、自己相関を示すプロットは改善する必要があります。これは、ggplot を使用すると、はるかによく見えます。

ただし、ggplot はベース グラフィックではないため、ggplot をレイアウトまたは par(mfrow) と組み合わせることはできません。ベースグラフィックスから生成された自己相関プロットをggplotによって生成されたものに置き換えるにはどうすればよいですか? すべての図が ggplot で作成されている場合、 grid.arrange を使用できることはわかっていますが、ggplot で生成されたプロットが 1 つだけの場合、どうすればよいですか?

多くのステーションの 30 年間のデータを含む毎日の時系列があります。biwavelet次のコードを使用して、データの周期性をテストするためにRのパッケージを使用しています

コードの実行時にこのエラーが発生しました

私が使用しているデータセットの中には、毎年の毎日の値が欠けているものがあります。私はかなり新しいRとプログラミングです。biwaveletパッケージを使用すると、値が欠落しているとこのようなエラーが発生するようです? 私の質問は、R でこの問題を解決する方法はありますか?

I cannot seem to find the wt() function in R used on this page: Combine base and ggplot graphics in R figure window

I could not find it in ??wt and on searching on google. It must be in some package. There is no library statement on the above page. Thanks for your help.

R のクロス ウェーブレット関数を理解しようとしていますが、biwavelet パッケージを使用して位相遅れ矢印を時間遅れに変換する方法がわかりません。例えば：

これらは私の 2 つの時系列です。どちらも同じですが、一方が他方より遅れています。矢印は 45 度の位相角を示しています。下向きまたは上向きは 90 度 (同相または異相) を意味しているように見えるので、私の解釈では 45 度の遅れを見ています。

これをタイムラグに変換するにはどうすればよいですか。つまり、これらの信号間のタイムラグをどのように計算しますか?

これは特定の波長に対してのみ実行できることをオンラインで読みました（これは、特定の期間を意味すると思いますか？）。では、365 の期間に関心があり、シグナル間のタイム ステップが 1 日であるとすると、タイム ラグはどのように計算されるでしょうか?

ウェーブレット コヒーレンス解析に R の biwavelet パッケージを使用しています。すべての時間ステップで 2 つの時系列間の相関 (rsq) を取得することに興味があります。私は使用しています

このコマンドは、列が時系列のタイム ステップ数と等しいが 35 行の行列を返します。各行が何を示しているのか、またはこのマトリックスをどのように解釈するのかがわかりません。どんな助けでも大歓迎です。

編集: dput 時系列 1 のデータ (最初の 10 行) - wtc 関数の実行時に変換をログに記録します

時系列 2 (最初の 10 行)

そのサブセットからの rsq 行列

biwavelet パッケージのプロット関数を使用していますが、このメッセージが表示され続けます

輪郭(): 円形/長いセグリスト -- オプション("max.contour.segments") > 25000 を設定しますか?

Rで次のコードを試しましたが、機能しません

R で max.contour.segments を増やすにはどうすればよいですか。