Question

我想在xts矩阵的所有周期上运行一个函数。 apply（）非常快，但返回的矩阵与原始对象相比具有转置尺寸：

> dim(myxts)
[1] 7429   48
> myxts.2 = apply(myxts, 1 , function(x) { return(x) })
> dim(myxts.2)
[1]   48 7429
> str(myxts)
An 'xts' object from 2012-01-03 09:30:00 to 2012-01-30 16:00:00 containing:
  Data: num [1:7429, 1:48] 4092500 4098500 4091500 4090300 4095200 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:48] "Open" "High" "Low" "Close" ...
  Indexed by objects of class: [POSIXlt,POSIXt] TZ: 
  xts Attributes:  
 NULL
> str(myxts.2)
 num [1:48, 1:7429] 4092500 4098500 4091100 4098500 0 ...
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:48] "Open" "High" "Low" "Close" ...
  ..$ : chr [1:7429] "2012-01-03 09:30:00" "2012-01-03 09:31:00" "2012-01-03 09:32:00" "2012-01-03 09:33:00" ...
> nrow(myxts)
[1] 7429
> head(myxts)
                       Open    High     Low   Close
2012-01-03 09:30:00 4092500 4098500 4091100 4098500
2012-01-03 09:31:00 4098500 4099500 4092000 4092000
2012-01-03 09:32:00 4091500 4095000 4090000 4090200 
2012-01-03 09:33:00 4090300 4096400 4090300 4094900
2012-01-03 09:34:00 4095200 4100000 4095200 4099900
2012-01-03 09:35:00 4100000 4100000 4096500 4097500

如何保留myxts尺寸？

Answer 1

这是apply记录的内容。来自?apply：

值：

 If each call to ‘FUN’ returns a vector of length ‘n’, then ‘apply’
 returns an array of dimension ‘c(n, dim(X)[MARGIN])’ if ‘n > 1’.

在您的情况下，'n'=48（因为您正在循环行），因此apply将返回维度c(48, 7429)的数组。

另请注意，myxts.2 不是 xts对象。这是一个常规阵列。你有几个选择：

在重新创建xts对象之前转置apply的结果：

data(sample_matrix)
myxts <- as.xts(sample_matrix)
dim(myxts)    # [1] 180   4
myxts.2 <- apply(myxts, 1 , identity)
dim(myxts.2)  # [1]   4 180
myxts.2 <- xts(t(apply(myxts, 1 , identity)), index(myxts))
dim(myxts.2)  # [1] 180   4

向量化您的函数，使其在xts的所有行上运行 object并返回一个xts对象。那你就不用担心了约apply。

最后，请开始提供可重复的示例。它并不那么难，它让人们更容易提供帮助。我在上面提供了一个例子，希望你能在以下问题中使用它。

为什么apply（）返回一个转置的xts矩阵？

1 个答案: