如何在R中替换向量字符并实现DNA序列反向互补转换功能？

完整实现代码 #

首先补全你已有的两个辅助函数（你之前的代码缺少闭合的大括号）：

string.to.vec <- function(s) {
  strsplit(s, "")[[1]]
}

vec.to.string <- function(v) {
  paste(v, collapse = "")
}

核心的reverse.dna函数用命名向量做碱基映射即可实现替换，逻辑非常简单：

reverse.dna <- function(s) {
  # 定义ATCG的互补碱基映射表
  comp_base <- c("A" = "T", "T" = "A", "G" = "C", "C" = "G")
  # 拆分字符串为单个碱基的向量
  base_vec <- string.to.vec(s)
  # 批量替换为互补碱基
  comp_vec <- comp_base[base_vec]
  # 反转互补序列后拼接为字符串返回
  return(vec.to.string(rev(comp_vec)))
}

运行测试 #

> s <- "CCCTTAG"
> reverse.dna(s)
[1] "CTAAGGG"

极简版本 #

如果不需要单独保留两个辅助函数，可以把逻辑合并为更短的版本：

reverse.dna <- function(s) {
  comp_base <- c("A"="T", "T"="A", "G"="C", "C"="G")
  paste(rev(comp_base[strsplit(toupper(s), "")[[1]]]), collapse = "")
}

注：函数内增加了toupper()处理，兼容输入小写碱基的场景，避免映射出错。

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