如何用nchar函数从蛋白质组学数据集中提取长度4-10的肽段

你可以通过以下步骤实现需求：

1. 先定义肽段向量：

x <- c("DFGVEVDNNSYK","EFHIR","VGQEYQPR","AYQVEGDWSGAAFMLVAGAIAGSIEVENLNTR","SHQADR","AILDALK","LR")

2. 利用nchar()获取每个肽段的字符长度，结合逻辑条件筛选出长度在4-10之间的肽段：

# 筛选长度为4到10的肽段（包含4和10）
filtered_peptides <- x[nchar(x) >= 4 & nchar(x) <= 10]

3. 查看筛选结果：

print(filtered_peptides)
# 输出结果：
# [1] "EFHIR"    "VGQEYQPR" "SHQADR"   "AILDALK"

关键说明： #

• nchar(x)会返回向量中每个元素的字符长度，得到一个和原向量长度一致的数值向量；
• 逻辑表达式nchar(x) >=4 & nchar(x) <=10会生成布尔向量，TRUE对应符合长度要求的肽段；
• 用该布尔向量对原向量索引，即可提取出目标肽段。

内容的提问来源于stack exchange，提问作者Bhagyashree Swarge