获取向量的最后n个元素.有没有比使用length()函数更好的方法？

Question

如果为了论证,我想要Python中10长度向量的最后五个元素,我可以在范围索引中使用" - "运算符,这样:

>>> x = range(10)
>>> x
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9]
>>> x[-5:]
[5, 6, 7, 8, 9]
>>>

R中最好的方法是什么？有没有比我现在的技术更简洁的方法,即使用length()函数？

> x <- 0:9
> x
 [1] 0 1 2 3 4 5 6 7 8 9
> x[(length(x) - 4):length(x)]
[1] 5 6 7 8 9
> 

这个问题与btw的时间序列分析有关,通常只对最近的数据有用.

Answer 1

看到?tail和?head一些方便的功能:

> x <- 1:10
> tail(x,5)
[1]  6  7  8  9 10

为了论证的缘故:除了最后五个元素之外的所有元素都是:

> head(x,n=-5)
[1] 1 2 3 4 5

正如@Martin Morgan在评论中所说,还有另外两种可能性比尾部解决方案更快,以防你必须在1亿个值的向量上进行一百万次.为了便于阅读,我会选择尾巴.

test                                        elapsed    relative 
tail(x, 5)                                    38.70     5.724852     
x[length(x) - (4:0)]                           6.76     1.000000     
x[seq.int(to = length(x), length.out = 5)]     7.53     1.113905     

基准代码:

require(rbenchmark)
x <- 1:1e8
do.call(
  benchmark,
  c(list(
    expression(tail(x,5)),
    expression(x[seq.int(to=length(x), length.out=5)]),
    expression(x[length(x)-(4:0)])
  ),  replications=1e6)
)

Answer 2

你可以在 R 中用另外两个字符做同样的事情：

x <- 0:9
x[-5:-1]
[1] 5 6 7 8 9

或者

x[-(1:5)]

Answer 3

tail这里基于速度的不赞同似乎并没有真正强调速度较慢的部分来自于尾部更安全的事实,如果你不确定x的长度会超过n这个数字要分组的元素:

x <- 1:10
tail(x, 20)
# [1]  1  2  3  4  5  6  7  8  9 10
x[length(x) - (0:19)]
#Error in x[length(x) - (0:19)] : 
#  only 0's may be mixed with negative subscripts

Tail将简单地返回最大元素数而不是生成错误,因此您不需要自己进行任何错误检查.使用它的一个很好的理由.更安全的代码,如果额外的微秒/毫秒对您的使用无关紧要.

Answer 4

怎么样rev(x)[1:5]？

x<-1:10
system.time(replicate(10e6,tail(x,5)))
 user  system elapsed 
 138.85    0.26  139.28 

system.time(replicate(10e6,rev(x)[1:5]))
 user  system elapsed 
 61.97    0.25   62.23