我想使用rollapply或rollapplyr将modwt函数应用于我的时间序列数据。
我熟悉rollapply / r的工作方式,但我需要一些帮助来设置输出,以便在使用rollapply时可以正确存储我的结果。
wavelim包中的modwt函数采用时间序列并将其分解为J级,对于我的特定问题J = 4,这意味着我将从我的单个时间序列中存储4组系数,存储在5的列表中。这个列表我只关心d1,d2,d3& D4。
modwt函数的输出如下所示
> str(ar1.modwt)
List of 5
$ d1: num [1:200] -0.223 -0.12 0.438 -0.275 0.21 ...
$ d2: num [1:200] 0.1848 -0.4699 -1.183 -0.9698 -0.0937 ...
$ d3: num [1:200] 0.5912 0.6997 0.5416 0.0742 -0.4989 ...
$ d4: num [1:200] 1.78 1.86 1.85 1.78 1.65 ...
$ s4: num [1:200] 4.64 4.42 4.19 3.94 3.71 ...
- attr(*, "class")= chr "modwt"
- attr(*, "wavelet")= chr "la8"
- attr(*, "boundary")= chr "periodic"
在上面的示例中,我已将modwt函数应用于长度为200的全长时间序列,但我希望使用rollapply将其应用于30的小滚动窗口。
我已经尝试了以下但输出是一个大矩阵,我不能轻易识别哪些值属于d1,d2,d3或d4
roller <- rollapplyr(ar1, 30,FUN=modwt,wf="la8",n.levels=4,boundary="periodic")
这是一个大矩阵,具有以下结构:
> str(roller)
List of 855
$ : num [1:30] 0.117 -0.138 0.199 -1.267 1.872 ...
$ : num [1:30] -0.171 0.453 -0.504 -0.189 0.849 ...
$ : num [1:30] 0.438 -0.3868 0.1618 -0.0973 -0.0247 ...
$ : num [1:30] -0.418 0.407 0.639 -2.013 1.349 ...
...省略了很多行...
$ : num [1:30] 0.307 -0.658 -0.105 1.128 -0.978 ...
[list output truncated]
- attr(*, "dim")= int [1:2] 171 5
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:5] "d1" "d2" "d3" "d4" ...
如何设置一个变量,以便为每个刻度d1,d2,d3和d4存储(200-30)+1列表及其中的列表?
有关可重复的示例,请使用以下内容:
library(waveslim)
data(ar1)
ar1.modwt <- modwt(ar1, "la8", 4)
答案 0 :(得分:3)
定义调用modwt的modwt2,获取前4个组件并将它们串联成数字向量。然后使用rollapplyr
与rollr
一起使用,其中rollr
的每一行都是对modwt2
的一次调用的结果。最后,将rollr
的每一行重新整形为一个单独的矩阵,并创建这些矩阵的列表L
:
modwt2 <- function(...) unlist(head(modwt(...), 4))
rollr <- rollapplyr(ar1, 30, FUN = modwt2, wf = "la8", n.levels = 4, boundary = "periodic")
L <- lapply(1:nrow(rollr), function(i) matrix(rollr[i,], , 4))
如果需要30 x 4 x 171阵列,则以下内容将其简化为3d数组:
simplify2array(L)
或作为清单列表:
lapply(L, function(x) as.list(as.data.frame(x)))
2)这是一种替代解决方案,它直接使用lapply
并返回一个列表,每个列表的组件都是由d1,d2,d3和d4组成的列表。
lapply(1:(200-30+1), function(i, ...) head(modwt(ar1[seq(i, length = 30)], ...), 4),
wf = "la8", n.levels = 4, boundary = "periodic")
更新:代码改进,展开(1)和添加(2)。