这个问题很模糊,但是有人在R中使用过biwavelet包并成功吗?我有以下代码:
require(biwavelet)
t1 <- cbind(DecTime,Temp)
## continuous wavelet transform
wt1 <- wt(t1)
plot(wt1)
它返回错误:
image.default(x $ t,yvals,t(zvals),zlim = zlims,ylim = rev(range(yvals)),错误:
“x”和“y”值必须是有限且不丢失的
我不明白此错误,因为我的数据不包含任何缺失值,而且它们都是有限的。
尝试重现示例时:
require(biwavelet)
Date = seq(from=as.POSIXct("2011-01-01 00:00"),
to=as.POSIXct("2011-12-31 23:00"), length=8760)
DecTime = julian(Date, Date[1])
data=cbind(as.numeric(DecTime), rnorm(8760))
## Continuous wavelet transform
wt.t1=wt(data)
plot(wt.t1)
它很烦人,因此我不知道为什么我的数据在此程序包中失败。上面示例中的数据与我的数据属于同一类,即“矩阵”。任何帮助或建议,将不胜感激。
已更改:
require(biwavelet)
Date = seq(from=as.POSIXct("2011-01-01 00:00"),
to=as.POSIXct("2011-12-31 23:00"), length=8760)
DecTime = julian(Date, Date[1])
D <- c(4.0267, 4.0211, 4.0005,4.0042,4.0042,4.0191)
data=cbind(as.numeric(DecTime[1:6]),as.numeric(D))
## Continuous wavelet transform
wt.t1=wt(data)
plot(wt.t1)
> data
[,1] [,2]
[1,] 0.00000000 4.0267
[2,] 0.04166667 4.0211
[3,] 0.08333333 4.0005
[4,] 0.12500000 4.0042
[5,] 0.16666667 4.0042
[6,] 0.20833333 4.0191
> class(data)
[1] "matrix"
> class(data[,1])
[1] "numeric"
> class(data[,2])
[1] "numeric"
> wt.t1=wt(data)
> plot(wt.t1)
Error in image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)), :
invalid z limits
最佳答案
万一您错过了我对其他帖子的答复,您发现wt.R函数中有一个错误(错误的括号)。该错误已在biwavelet软件包的0.12版中得到了解决,因此您的代码现在可以使用了。
感谢您发现错误!
关于r - 连续小波变换,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/11192291/