只有0可以与负下标 - W2CWM2C小波混合使用

时间:2017-07-11 21:54:42

标签: r time-series wavelet

我正在使用W2CWM2C软件包作为我的项目的一部分,用于调查应用于时间序列数据的小波。

我使用W2CWM2C软件包对我的数据执行双变量小波相关,然后将它们绘制到屏幕上。

pg 6上的W2CWM2C文档类似,我使用的是2个时间序列变量

library(data.table) # load the necessary packages
test = data.table(a=seq(1,10), b=seq(10,19)) #create the data.table
test[,c:=a+b]#create a new column named c and fill it with the sum of a and b

对于那些感兴趣的人,我的原始数据存储在一个数据框中,代表以15分钟为间隔记录的2周网络数据。数据具有24小时的季节性,即峰值大约每24小时发生一次,因为每个变量的测量值为1:1343,因此频率为4 x 24 = 96,我希望这是正确的!

我已将以下代码用于我的数据样本

> class(wideRawTS)
[1] "mts"    "ts"     "matrix"

> head(wideRawTS[,c(1:2)],10)
Time Series:
Start = c(1, 1) 
End = c(1, 10) 
Frequency = 96 
         DO0182U09A3 DO0182U09B3
1.000000     -101.50     -103.37
1.010417     -101.32     -102.75
1.020833     -101.45     -103.30
1.031250     -100.91      -95.92
1.041667     -100.91     -103.04
1.052083     -100.97     -103.67
1.062500     -101.48     -103.26
1.072917     -100.85     -103.86
1.083333     -100.87     -103.96
1.093750     -100.77     -103.33

当我运行此代码时,我收到以下错误:

wideRawTS <- ts(wideRawDF[,2:13], frequency = 96)
U09A3_B3 <- wideRawTS[,c(1,2)]
Wname <- "la8"
J <- 8
lmax <- 10
compWCC <- WCC(U09A3_B3, Wname, J, lmax, device = "screen", NULL,NULL,NULL,NULL,NULL)

这与我的数据是负面的吗?我已经用Google搜索了这个错误,但是我发现了这个帖子,但是我的头脑已经过了。通过阅读一些评论,我认为这个问题与正面和负面的子集有关,我不确定它是否适用于此。

如果有任何人遇到类似问题的经验,我会很高兴听到你如何解决它。

提供可重现的示例here是托管在Google云端硬盘上的wideRawDF数据框的副本。

0 个答案:

没有答案