在 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 版本中得到修复，因此您上面的两个版本的代码现在都应该可以工作了。

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

Answer 2

我怀疑这是因为您在这里只有一个频率并且没有为此设置功能。我可以通过将白噪声添加到y 来获得绘图：

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 + rnorm(length(y))) ## add some white noise to y
wt.t1 <- wt(d)
plot(wt.t1)

您可能希望联系维护人员报告该问题。当我调试它并反转zlim 时，我得到了做某事的情节，因此diff(zlim) 是肯定的，所以plot() 方法的作者可能做出了一个并非在所有情况下都成立的假设。