按列名从xts中选择

时间:2018-06-03 19:07:13

标签: r xts

我正在尝试在函数中按名称操作xts对象中的特定列,但我一直收到错误:

Error in if (length(c(year, month, day, hour, min, sec)) == 6 && all(c(year, : missing value where TRUE/FALSE needed In addition: Warning messages: 1: In as_numeric(YYYY) : NAs introduced by coercion 2: In as_numeric(YYYY) : NAs introduced by coercion

如果我有一个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对象做同样的事情,我会得到上面的错误。 e.g。

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

跟:

testfun(xts3)

这有效:

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

但我真的想通过名字选择,因为我无法确定列顺序。

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

谢谢!

2 个答案:

答案 0 :(得分:2)

xts - 对象具有类c("xts", "zoo"),这意味着它们是具有由其创建函数指定的特殊属性的矩阵。虽然$使用矩阵不会成功,但由于xts方法,它可与zoo$.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函数提取&#34;数据&#34;首先返回元素,然后返回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 - 操作无法传递具有共享列名的对象,因此需要进行数字访问。 (我原本期望一个名为&#34; Data&#34;元素,但您/我们还需要阅读?rbind.xts。)

答案 1 :(得分:1)

输入?`[.xts`,您会看到该函数有ij参数(等等)。

  

i - 要提取的行。数字,时基或ISO-8601样式(见详情)

     

j - 要提取的列,数字或按名称

您通过'xts1'作为i参数,而它应该是j。所以你的功能应该是

testfun <- function(xts_data){
  print(xts_data[, 'xts1']) # or xts3[j = 'xts1']
}