For/while 循环直到保持顺序答案

【问题标题】：For/ while loops until sequence is maintainedFor/while 循环直到保持顺序
【发布时间】：2017-12-04 04:19:43
【问题描述】：

我有一个如下所示的数据框：

input <- structure(list(rank = c(1L, 2L, 3L, 3L, 4L, 4L, 5L, 5L, 6L, 6L, 
7L, 7L, 8L, 8L, 9L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 13L, 14L, 
14L, 15L, 16L, 17L, 18L, 19L), sequence = c("HRIGRGGRYGRKGVAI", 
"TQIDELPVDFAAYLGE", "AFSIGLLQRLDFRHNL", "QNDKIAPQDQDSFLDD", "SMHAEMPKSDRERVMN", 
"AQSVIFANTRRKVDWI", "PGRVSDVIKRGALRTE", "AEVISRIGEFLSNSSK", "GGDIIAQAQSGTGKTG", 
"TFVGGTRVQDDLRKLQ", "QGLVLSPTRELALQTA", "DWIAEKLNQSNHTVSS", "NIVINFDLPTNKENYL", 
"AGVIVAVGTPGRVSDV", "SDRERVMNTFRSGSSR", "GFEKPSSIQQRAIAPF", "SGTGKTGAFSIGLLQR", 
"LDTLMDLYETVSIAQS", "VRPIPSFDDMPLHQNL", "MPEEVLELTKKFMRDP", "QQRAIAPFTRGGDIIA", 
"LHEIEAHYHTQIDELP", "LVARGIDVHHVNIVIN", "ANTRRKVDWIAEKLNQ", "VLVLDEADEMLSQGFA", 
"RGALRTESLRVLVLDE", "PQDQDSFLDDQPGVRP", "YGRKGVAINFVTEKDV", "SSKFCETFVGGTRVQD", 
"RVLVTTDLVARGIDVH"), start_position = c(353L, 388L, 79L, 3L, 
296L, 268L, 155L, 111L, 63L, 130L, 96L, 281L, 337L, 146L, 304L, 
45L, 72L, 255L, 22L, 212L, 53L, 379L, 326L, 274L, 174L, 164L, 
9L, 361L, 124L, 319L), score = c(0.92, 0.89, 0.87, 0.87, 0.86, 
0.86, 0.85, 0.85, 0.84, 0.84, 0.79, 0.79, 0.78, 0.78, 0.77, 0.76, 
0.75, 0.75, 0.75, 0.75, 0.74, 0.74, 0.73, 0.72, 0.72, 0.71, 0.68, 
0.67, 0.65, 0.63)), .Names = c("rank", "sequence", "start_position", 
"score"), row.names = c(NA, -30L), class = c("tbl_df", "tbl", 
"data.frame"))

我想做的是以下。查看input$rank，我想将input$score 下的分数相加，直到保持input$rank 下的序列。

例如，考虑第一个序列，从 rows 1:36 开始（input$rank 下的第 37 个值是 1 - 注意input$rank 下有重复值），我的总和是 26.76 - 我是通过 sum(input$score[1:36]) 获得的。

我想过在for或while循环中插入break或next，尽管我对循环中的这些参数不太熟悉。

【问题讨论】：

你能添加你的预期输出吗？
可以是向量，也可以是其他列，只要保持顺序即可。然后我将输出绑定到另一个数据框
“直到 input$rank 下的序列被保持”是什么意思？另请注意，在您发布的示例数据中，只有 30 行。
是的，我刚刚粘贴了部分数据框，因为它有超过 4k 行。保持序列我的意思是序列没有中断，即是连续的。在我的示例中，序列从 1:23 开始不间断，然后排名下的下一个数字是 1，这意味着序列被中断。
对，我意识到你不能粘贴整个 4k 行，我想暗示的是你发布的数据不允许我们重现你要求的输出，因为它只有 30 行长（不是 36 行）。难怪费德里科想不通。如果您可以重新处理示例数据以使其反映您在实际数据中的特征，那将会很有帮助。你可以编一个等级，比如1, 1, 2, 3, 4, 1, 2, 3, 3, 5, 5, 1, 2, 4 或其他东西，作为更简单、更简洁的真实mcve 的基础

标签： r for-loop while-loop break next

【解决方案1】：

aggregate(input[c(1,4)],by=list(input$rank),sum)

这是你想要的吗？

【讨论】：

不是真的，很遗憾。从我编辑的答案中，您可以看到，考虑到第一个序列，我将有 26.76 作为输出，如果我使用您的代码，则此值不存在。除此之外，我总共有 65 个序列，使用您的代码的输出给了我 47 个序列。

【解决方案2】：

希望这更接近您正在寻找的内容。我所做的是针对小于零的值测试秩向量的差异版本，并对其进行累积和。然后将结果用作对aggregate() 的调用中的分组向量。

set.seed(1)
rank <- c(1, 2, 3, 5, 5, 1, 2, 2, 3, 1, 2, 4, 4, 5)
score <- round(runif(length(rank)), 2)
input <- data.frame(rank, score)
input <- cbind(group=cumsum(c(-1, diff(input$rank)) < 0), input)
input
#    group rank score
# 1      1    1  0.27
# 2      1    2  0.37
# 3      1    3  0.57
# 4      1    5  0.91
# 5      1    5  0.20
# 6      2    1  0.90
# 7      2    2  0.94
# 8      2    2  0.66
# 9      2    3  0.63
# 10     3    1  0.06
# 11     3    2  0.21
# 12     3    4  0.18
# 13     3    4  0.69
# 14     3    5  0.38

aggregate(score ~ group, data=input, sum)
#   group score
# 1     1  2.32
# 2     2  3.13
# 3     3  1.52

【讨论】：