如何解决R:在向量中找到在一定数量的迭代中不会减少的位置
numbers <- c(1,0.9,0.8,0.7,0.71,0.72,0.69,0.696,0.697,0.55,0.6,0.66,0.56,0.58)
鉴于numbers
,我想找到第一个索引,该索引的下一个numbers
值在n = 5
中的值不减小。在上面的示例中,我正在寻找的索引是8,因为当numbers[8] <= numbers[9:(8 + 5)]
时。这是我的尝试:
myfun <- function(numbers,n){
for(i in 1:length(numbers)){
if(all(numbers[i] <= numbers[(i + 1):(i + n)])){
return(i)
}
}
}
> myfun(numbers,5)
[1] 8
是否有一种无需编写for循环即可更快地获得答案的方法?
解决方法
编辑
我想我之前误解了这个问题(感谢@thelatemail引起注意)。您想找出小于所有下一个n
值的值。
您可以通过滚动操作来做到这一点。
n <- 5
which(zoo::rollapply(numbers,n,function(x) all(x >= x[1])))[1]
#[1] 8
更早的答案
这将返回numbers
中连续增加的序列的索引。
您可以使用rle
:
n <- 5
with(rle(diff(numbers) > 0),sum(lengths[seq_len(which(lengths >= n & values)[1] - 1)])) + 1
#[1] 8
您可以将其分解以更好地理解:
diff
给出连续数字之间的差异。
diff(numbers)
# [1] -0.100 -0.100 -0.100 0.010 -0.010 0.020 -0.030 0.006 0.001
#[10] 0.003 0.010 0.010 -0.170 0.050 0.060 -0.110 0.010 0.020
我们将其与> 0
进行比较,以得出TRUE
用于增加值,FALSE
用于减小值。
diff(numbers) > 0
# [1] FALSE FALSE FALSE TRUE FALSE TRUE FALSE TRUE TRUE TRUE TRUE
#[12] TRUE FALSE TRUE TRUE FALSE TRUE TRUE
我们在其上应用rle
:
tmp <- rle(diff(numbers) > 0)
tmp
#Run Length Encoding
# lengths: int [1:10] 3 1 1 1 1 5 1 2 1 2
# values : logi [1:10] FALSE TRUE FALSE TRUE FALSE TRUE ...
我们找到一个位置,其中递增序列的长度大于等于n
tmp$lengths >= n & tmp$values
#[1] FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE
使用which
获取索引,如果有多个,则使用[1]
选择第一个:
which(tmp$lengths >= n & tmp$values)[1]
[1] 6
sum
之前此索引的所有长度,因此比数字大-1
sum(tmp$lengths[seq_len(which(tmp$lengths >= n & tmp$values)) - 1])
#[1] 7
现在在上面的数字上加+1,以获得下一个索引。
如果使用这种分步方法,则可以比顶部的单一衬套更轻松地处理不同的边缘情况。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。