1

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

require(gamair)
data(cairo)
data_1 <- within(cairo, Date <- as.Date(paste(year, month, day.of.month, sep = "-")))
data_1 <- data_1[,c('Date','temp')]
data_2 <- data_1

# add a lag
n <- nrow(data_1)
nn <- n - 49
data_1 <- data_1[1:nn,]
data_2 <- data_2[50:nrow(data_2),]
data_2[,1] <- data_1[,1]

require(biwavelet)
d1 <- data_1[,c('Date','temp')]
d2 <- data_2[,c('Date','temp')]
xt1 <- xwt(d1,d2)
plot(xt1, plot.phase = TRUE)

ここに画像の説明を入力 ここに画像の説明を入力

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

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

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

4

2 に答える 2