我正在使用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数据框的副本。