从 xts 对象的子集中减去 xts 对象答案

【问题标题】：Subtracting a xts object from a subset of xts object从 xts 对象的子集中减去 xts 对象
【发布时间】：2013-05-27 10:30:12
【问题描述】：

con = gzcon(url('http://www.systematicportfolio.com/sit.gz', 'rb'))
source(con)
close(con)
load.packages("TTR,PerformanceAnalytics,quantmod,lattice")

#######################################################
#Get and Prep Data
#######################################################
data <- new.env()
tickers<-spl("VTI,IEF,TLT,DBC,VNQ,GLD")

getSymbols(tickers, src = 'yahoo', from = '1980-01-01', env = data)
for(i in ls(data)) data[[i]] = adjustOHLC(data[[i]], use.Adjusted=TRUE)

bt.prep(data, align='remove.na', dates='1990::2013')

我遇到了从 xts 对象中减去特定列的问题。

prices = data$prices
ret = prices / mlag(prices) - 1
ret - ret[,3]  #subtract column three from every other column don't seem to work

有没有快速的解决方案？

我试过了：

apply(ret,2,function(x) x - x[,3]) #doesn't seem to work

有什么想法吗？

【问题讨论】：

SIT 使用 XTS 对象来包装时间序列，所以我认为这不是问题。如果你做 str(ret)，它是一个 xts 对象。
SIT 使用 XTS 对象包装时间序列。上面的示例只是下载数据并让 SIT 将其包装在 XTS 对象周围。但我现在将其编辑为“最小”版本..
这里还有一点点 apply-like: sweep(x, 1, x[, 3])

标签： r xts

【解决方案1】：

下一次，请提供一个最小可重现的例子。例如：

> library(xts)
> data(sample_matrix)
> x <- as.xts(sample_matrix)
> x-x[,1]
Error in `-.default`(x, x[, 1]) : non-conformable arrays
> apply(x, 2, function(y) y-x[,1])
Error in array(ans, c(len.a%/%d2, d.ans), if (!all(vapply(dn.ans, is.null,  : 
  length of 'dimnames' [1] not equal to array extent

问题是 xts 对象默认具有dim 属性，并且在子集时不会像矩阵和动物园类对象那样删除它。您可以通过在子集调用中设置 drop=TRUE 来强制删除它。

> head(x-x[,1,drop=TRUE])
           Open       High         Low       Close
2007-01-02    0 0.07799532 -0.08936727  0.07799532
2007-01-03    0 0.19137980  0.00000000  0.16717014
2007-01-04    0 0.00000000 -0.15681864 -0.08859811
2007-01-05    0 0.00000000 -0.15243423 -0.03887316
2007-01-06    0 0.00000000 -0.13311797 -0.06320448
2007-01-07    0 0.08349916 -0.14025780 -0.14025780

这是可行的，因为x[,1,drop=TRUE] 返回一个“向量 xts”（即无量纲 xts 对象），并且在 - 调用期间，该向量沿 x 被回收。

【讨论】：