4

我正在寻找在数据框中创建一个现有列的副本，该列偏移了许多行。

例如，如果 column2 是 column1 偏移 1 的副本，则

> dataframe
$column1
[1] 1 2 3 4 5

$column2
[1] 0 1 2 3 4

我使用以下代码取得了一些成功：

offset7 <- rep(0, 7)
dataframe$column1.prev7 = c(offset7, dataframe$column1[1:(length(dataframe$column1)-7)])

但是，如果我偏移 30 或更多，它就会开始出错。我的数据足够长，这不会成为偏移量大于行数的问题。错误是：

Error in dataframe$column1[1:(length(dataframe$column1) - 30)] : 
  only 0's may be mixed with negative subscripts

提前致谢！与 plyr 配合得很好的快速无循环版本将是首选。这里的目的是将时间序列数据分解为长达一年的各种滞后，然后以各种方式分析结果。

4

1 回答 1

9

请使用适当的时间序列类进行时间序列操作。最受欢迎的是zoo和xts，它们都有大量的文档。

作为一个简单的例子，考虑

> library(xts)
> foo <- xts(100:109, order.by=Sys.Date()+0:9)
> merge(foo,  l1=lag(foo,1), lm1=lag(foo,-1))
           foo  l1 lm1
2010-11-18 100  NA 101
2010-11-19 101 100 102
2010-11-20 102 101 103
2010-11-21 103 102 104
2010-11-22 104 103 105
2010-11-23 105 104 106
2010-11-24 106 105 107
2010-11-25 107 106 108
2010-11-26 108 107 109
2010-11-27 109 108  NA
>

但不要用手去做。并在此处搜索 '[r] xts' 或 [r] zoo' 以在 R 标签中进行搜索。

于 2010-11-18T21:27:18.700 回答