如何在时间序列中自动创建n个滞后？

maximusdooku 发表于 Dev

极大主义

我有一个带有列t的数据框。我想创建n个滞后列，其名称类似t-1，t-2等。

  year      t  t-1 t-2
19620101    1   NA  NA
19630102    2   1   NA
19640103    3   2   1
19650104    4   3   2
19650104    5   4   3
19650104    6   5   4

我的想法是，我将分四个步骤进行操作：

使用“粘贴”的列名循环
使用“粘贴”的滞后列的临时数据帧的循环
创建滞后列的循环
捆绑他们。

但是我无法继续执行代码。有点粗糙：

df_final<-lagged(df="odd",n=3)

lagged<-function(df,n){
   df<-zoo(df)
   lags<-paste("A", 1:n, sep ="_")
   for (i in 1:5) {
     odd<-as.data.frame(lag(odd$OBS_Q,-1*i,na.pad =  TRUE))

   #Cbind here
   }

我一直在写这个函数。你能显示一些吗？或更简单的方法

参考：R向量/数据帧中的基本滞后

附录：

真实数据：

x<-structure(list(DATE = 19630101:19630104, PRECIP = c(0, 0, 0,0), 
               OBS_Q = c(1.61, 1.48, 1.4, 1.33), swb = c(1.75, 1.73, 1.7,1.67), 
               gr4j = c(1.9, 1.77, 1.67, 1.58), isba = c(0.83, 0.83,0.83, 0.83), 
               noah = c(1.31, 1.19, 1.24, 1.31), sac = c(1.99,1.8, 1.66, 1.57), 
               swap = c(1.1, 1.05, 1.08, 0.99), vic.mm.day. = c(2.1,1.75, 1.55, 1.43)), 
          .Names = c("DATE", "PRECIP", "OBS_Q", "swb","gr4j", "isba", "noah", "sac", "swap", "vic.mm.day."), 
          class = c("data.table","data.frame"), row.names = c(NA, -4L))

要滞后的列是OBS_Q。

大卫·阿伦堡

如果您正在寻找效率，请尝试data.table新shift功能

library(data.table) # V >= 1.9.5
n <- 2
setDT(df)[, paste("t", 1:n) := shift(t, 1:n)][]
#    t t 1 t 2
# 1: 1  NA  NA
# 2: 2   1  NA
# 3: 3   2   1
# 4: 4   3   2
# 5: 5   4   3
# 6: 6   5   4

在这里，您可以为新列设置任何名称（在内paste），也不需要将其绑定回原始名称，因为这会使用:=运算符通过引用更新您的数据集。

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。