ts.intersect不适用于xts对象

2023年2月13日 297次阅读

以下产生错误

a1 = as.xts(ts(rnorm(20), start=c(1980,1), freq=4))
a2 = as.xts(ts(rnorm(30), start=c(1983,1), freq=4))
a = ts.intersect(a1,a2)

Error in .cbind.ts(list(...), .makeNamesTs(...), dframe = dframe, union = FALSE) : 
  no time series supplied

文件说应该是

两个或多个单变量或多变量时间序列,或者可以强制转换为时间序列的对象.

最佳答案 ts.intersect通过查找tsp属性来确定对象是否是ts对象. as.xts.ts删除了tsp属性,这就是它不被强制回ts对象的原因.

这看起来像xts-> ts-> xts转换中的错误,但我需要仔细看看.

作为解决方法,您可以手动将tsp属性添加到xts对象(请注意,这可能会导致其他xts方法出现问题,例如str.xts)并添加.tsp属性.

set.seed(21)
A1 <- ts(rnorm(20), start=c(1980,1), freq=4)
A2 <- ts(rnorm(30), start=c(1983,1), freq=4)
# convert to xts
a1 <- as.xts(A1)
a2 <- as.xts(A2)
# add tsp attribute
# (so stats:::.cbind.ts will think these are coercible to ts objects)
tsp(a1) <- tsp(A1)
tsp(a2) <- tsp(A2)
# add .tsp attribute
# (needed for as.ts.xts to work)
attr(a1,'.tsp') <- tsp(A1)
attr(a2,'.tsp') <- tsp(A2)

a <- ts.intersect(a1,a2)
点赞