2
set.seed (12345)

data <- paste(sample(c("A","C","G","T"),100000,replace=TRUE,prob=rep(0.25,4)))
data <- ifelse(data=="A",1,0)

假设我将数据转换为 1(所需字符)和 0(其他)。然后在每个位置求和。如果一个位置的总和与下一个位置的总和匹配,我们停止;否则我们进行总和并将总和存储在位置上。然后对应位置的最大和将给出序列的最大长度。

我得到了算法,但无法编码。请帮忙。

4

1 回答 1

5

rle功能是您想要的:

set.seed(12345)
data = sample(c('A', 'C', 'G', 'T'), 100000, replace = TRUE, prob = rep(0.25, 4))

run_lengths = rle(data == 'A')
(result = max(run_lengths$lengths[run_lengths$values]))
# [1] 10

获得最长跑的位置有点困难。您可以使用which.max它,但我们之前已经过滤掉了所有非A结果。相反,我们也可以将所有非A运行设置为0 ——这样,它们仍然存在,但不会是最大值:

only_a = ifelse(run_lengths$values, run_lengths$lengths, 0)
longest_run_index = which.max(only_a)

现在我们需要计算从longest_run_index到 里面的索引A。我们通过将这个索引之前的所有运行的长度相加来做到这一点:

index = sum(run_lengths$lengths[seq_len(longest_run_index - 1)]) + 1
data[index : (index + result - 1)]
# [1] "A" "A" "A" "A" "A" "A" "A" "A" "A" "A"
于 2015-09-24T19:14:37.403 回答