使用时间窗口在R中运行测试

我试图从运行的lepage测试中获得lepage测试统计数据的时间序列

时间窗为9年

以下是我使用的数据（1961年至2009年）

我正在使用nonpar软件包进行lepage测试

我不知道如何在R中实现这一点。 以下是我目前的脚本（这是错误的）：

库（非AR）
当你这样做的时候

test<-rollapplyr(x, 9, function(u) lepage.test(u,u+1)$obs.stat)

令人惊叹的！非常感谢你的帮助！。我不能编辑文章，但应该是：pairs2
 library(nonpar)
 dat<-read.csv("test.csv",header=T)
 x <- as.numeric(dat$Winter)

 for (i in 1:length(x)){
 #window
 n = 9
     stat<-lepage.test(x[i],x[i:i+n])$obs.stat
  }
 return(stat) 
 } 

library(zoo)
test<-rollapplyr(x, 9, function(u) lepage.test(u,u+1)$obs.stat)

test<-rollapplyr(x, 9, function(u) lepage.test(u,u+1)$obs.stat)

library(zoo)
library(nonpar)
x <- zoo(dat$Winter)
y <- rollapply(x, 9, c)
pairs <- lapply(1:(nrow(y)-9), function(i) list(x1 = y[i,], x2 = y[i+9,]))
pairs2 <- lapply(pairs, function(x) lapply(x, as.numeric))
test <- sapply(1:32, function(i) lepage.test(pairs2[[i]]$x1, pairs2[[i]]$x2)$obs.stat)
test
#  [1] 1.8182 2.5782 0.8046 0.2214 0.9449 0.6128 2.4022 7.0955 6.8560 7.8285 5.5361
# [12] 3.2331 2.6603 2.1228 0.9856 0.6616 1.3692 2.1103 1.1509 0.7352 0.5057 1.5349
# [23] 3.8019 4.7344 2.2081 2.3770 0.1579 1.1430 3.8876 3.1820 1.8182 0.1435