r - 在 R 中使用 biwavelet 包

Question

谁能指出以下示例代码不起作用的原因:

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 的回答:如果我将数据保留为仅包含一个频率但更改时间向量，则代码可以正常工作。

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)

Answer 1

您在 wt.R 函数中发现了一个错误(错误的括号)。该错误已在 biwavelet 包的 0.12 版本中修复，因此您上面的两个版本的代码现在都应该可以工作了。

感谢您发现错误。请不要犹豫，向软件包的维护者(即我)发送电子邮件告知 future 的错误。