次のコード例が機能しない理由を誰かが提案できますか:
require(biwavelet)
t <- seq(1/24, 365, 1/24)
A <- 2
fs <- 1/24
y <- A + sin(2*pi*fs*t)
d = cbind(t,y)
wt.t1 <- wt(d)
plot(wt.t1)
次のようなエラーが生成されます。
Error in image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)), :
invalid z limits
この問題を解決するにはどうすればよいですか?
追加:
Gavin Simpsons の回答への回答: 1 つの周波数のみを含むようにデータを保持し、時間ベクトルを変更すると、コードは正常に動作します。
require(biwavelet)
A <- 2
fs <- 1/24
y <- A + sin(2*pi*fs*t)
d <- cbind(seq(1,8760), y)
wt.t1 <- wt(d)
plot(wt.t1)