Question

我正在尝试使用lag函数在xts对象中创建滞后向量。它在使用$表示法（例如x.ts$r1_lag）在xts对象中定义新向量时有效，但在使用方括号定义新变量时也是如此，即xts[,"r1_lag"]。请参阅以下代码：

library(xts)
x <- data.frame(date=seq(as.Date('2015-01-01'), by='days', length=100),
                runif(1e2), runif(1e2), runif(1e2))
colnames(x) <- c("date", "r1", "r2", "r3")

#the following command works
x.ts <- xts(x, order.by=x$date)
x.ts$r1_lag <- lag(x.ts$r1)
# but the following does not (says subscript is out of bounds)
x.ts <- xts(x, order.by=x$date)
x.ts[,"r1_lag"] <- lag(x.ts[,"r1"])

我需要使用[]符号而不是$符号来引用向量，因为如果我想对多个xts对象中的向量运行滞后转换（多个xts列表中的向量）对象），我无法使用$表示法在对象中定义新向量，即我无法使用下面的程式化循环中的符号来定义新向量：

for (i in letters) {
  for (j in variables) {
    macro.set.ts$i$paste(j,"_L1",sep="") <- lag(macro.set.ts[[i]][,j])
    macro.set.ts$i$paste(j,"_L2",sep="") <- lag(macro.set.ts[[i]][,j], 2)
    macro.set.ts$i$paste(j,"_L4",sep="") <- lag(macro.set.ts[[i]][,j], 4)
  }
}

谢谢！

Answer 1

您没有需要来使用cmake . || exit 1 make || exit 1 ./pcl_visualizer_demo。您可以改为使用[<-.xts：

merge

Answer 2

该错误与lag功能无关。您收到错误，因为您尝试使用另一个xts对象分配xts对象。此示例重现错误：

x.date= seq(as.Date('2015-01-01'), 
                        by = 'days' , length = 5)
x1 <- xts(data.frame(c1=runif(5)), order.by=x.date)
x2 <- xts(data.frame(c2=runif(5)), order.by=x.date)
x1[,'r2'] <- x2

## Error in `[<-.default`(`*tmp*`, , "r2", 
##   subscript out of bounds

我发现这在xts逻辑中是连贯的，因为xts是索引对象。因此，最好合并对象或加入并保留时间序列的索引性质。

merge(x1,x2)

这将cbind 2次系列并修复任何索引问题。事实上，cbind只是merge：

identical(cbind(x1,x2),merge(x1,x2)

我认为这是一种错误，这适用于$<-运算符，而不适用于[<-运算符。

Answer 3

我得到了相同的输出：

x.ts <- cbind(x.ts,lag(x.ts[,"r1"]))

和

x.ts <- transform(x.ts, r1_lag = lag(x.ts[,'r1']))

但是，要小心输出。它看起来可能相同，但结构也有所改变。

Answer 4

这应该有效：

x.ts <- merge(x.ts,lag(x.ts[,"r1"]))

然后，您可能希望重命名添加的最后一列：

dimnames(x.ts)[[2]][5] <- "r1_lag"

结果如下：

> head(x.ts)
           date         r1           r2           r3             r1_lag      
2015-01-01 "2015-01-01" "0.23171030" "0.44174424" "0.3396816640" NA          
2015-01-02 "2015-01-02" "0.97292220" "0.74909452" "0.2793033421" "0.23171030"
2015-01-03 "2015-01-03" "0.52320743" "0.49288463" "0.0193637393" "0.97292220"
2015-01-04 "2015-01-04" "0.36574297" "0.69571803" "0.6411834760" "0.52320743"
2015-01-05 "2015-01-05" "0.37563137" "0.13841216" "0.3087215754" "0.36574297"
2015-01-06 "2015-01-06" "0.48089356" "0.32702759" "0.3967609401" "0.37563137"

> class(x.ts)
[1] "xts" "zoo"

希望这有帮助。

使用$ vs. []表示法创建xts对象的滞后变量

4 个答案: