r - 按列名从 xts 中选择

Question

我试图在函数中按名称对 xts 对象中的特定列进行操作，但我不断收到错误消息:

if (length(c(year, month, day, hour, min, sec)) == 6 && all(c(year, : 在需要 TRUE/FALSE 的地方缺失值 此外: 警告信息: 1:在 as_numeric(YYYY) 中:强制引入的 NA 2:在 as_numeric(YYYY) 中:强制引入的 NA

如果我有一个 xts 对象:

xts1 <- xts(x=1:10, order.by=Sys.Date()-1:10)
xts2 <- xts(x=1:10, order.by=Sys.Date()+1:10)
xts3 <- merge(xts1, xts2)

然后我可以选择一个特定的列:

xts3$xts1

使用数据框，我可以将 xts3 传递给另一个函数，然后选择一个特定的列:

xts3['xts1']

但是如果我尝试对 xts 对象做同样的事情，我会得到上面的错误。例如

testfun <- function(xts_data){
  print(xts_data['xts1'])
}

调用方式:

testfun(xts3)

这个有效:

testfun <- function(xts_data){
  print(xts_data[,1])
}

但我真的很想按名称选择，因为我无法确定列顺序。

谁能建议如何解决这个问题？

谢谢!

Answer 1

xts - 对象具有类 c("xts", "zoo")，这意味着它们是具有由其创建函数分配的特殊属性的矩阵。虽然 $ 不会成功地处理矩阵，但由于 $.zoo，它可以与 xts 和 zoo 对象一起使用方法。 (也不建议在函数内部使用 $，因为名称评估混淆和部分名称匹配的可能性。)参见:?xts 并检查 sample.xts 对象用第一个例子用 str 创建:

> ?xts
starting httpd help server ... done
> data(sample_matrix)
> sample.xts <- as.xts(sample_matrix, descr='my new xts object')
> 
> str(sample.xts)
An ‘xts’ object on 2007-01-02/2007-06-30 containing:
  Data: num [1:180, 1:4] 50 50.2 50.4 50.4 50.2 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:4] "Open" "High" "Low" "Close"
  Indexed by objects of class: [POSIXct,POSIXt] TZ: 
  xts Attributes:  
List of 1
 $ descr: chr "my new xts object"

 class(sample.xts)
# [1] "xts" "zoo"

这解释了为什么早期的答案建议使用 xts3[ , "x"] 或等效的 xts3[ , 1] 应该成功。 [.xts 函数首先提取“数据”元素，然后返回由 j 参数指定的命名列或编号列。

 str(xts3)
An ‘xts’ object on 2018-05-24/2018-06-13 containing:
  Data: int [1:20, 1:2] 10 9 8 7 6 5 4 3 2 1 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:2] "xts1" "xts2"
  Indexed by objects of class: [Date] TZ: UTC
  xts Attributes:  
 NULL
> xts3[ , "xts1"]
           xts1
2018-05-24   10
2018-05-25    9
2018-05-26    8
2018-05-27    7
2018-05-28    6
2018-05-29    5
2018-05-30    4
2018-05-31    3
2018-06-01    2
2018-06-02    1
2018-06-04   NA
2018-06-05   NA
2018-06-06   NA
2018-06-07   NA
2018-06-08   NA
2018-06-09   NA
2018-06-10   NA
2018-06-11   NA
2018-06-12   NA
2018-06-13   NA

merge.xts 操作可能没有达到您的预期，因为日期范围没有重叠。您似乎可能想要:

> xts4 <- rbind(xts1, xts2)
> str(xts4)
An ‘xts’ object on 2018-05-24/2018-06-13 containing:
  Data: int [1:20, 1] 10 9 8 7 6 5 4 3 2 1 ...
  Indexed by objects of class: [Date] TZ: UTC
  xts Attributes:  
 NULL

请注意，rbind.xts 操作未能交付具有共享列名的对象，因此需要进行数字访问。 (我希望有一个命名的“数据”元素，但您/我们还需要阅读 ?rbind.xts。)