R函数实现氨基酸三字母代码转单字母的错误排查与解决

你原来的代码出错是因为gsub的第二个参数直接把所有单字母代码拼接成了一个长字符串——不管匹配到哪个三字母代码，都会替换成这个完整长串，这就导致了ARNDCEQGHILKMFPSTWYV3=ARNDCEQGHILKMFPSTWYV20这种错误结果。

最佳解决方案（依赖stringr包） #

stringr的str_replace_all函数原生支持用命名向量作为替换字典，写法简洁且高效：

library(stringr)

substitute_codes <- function(data, col_name) {
  # 氨基酸三字母→单字母映射字典
  code_dict <- c("ALA" = "A", "ARG" = "R", "ASN" = "N", "ASP" = "D",
                 "CYS" = "C", "GLU" = "E", "GLN" = "Q", "GLY" = "G",
                 "HIS" = "H", "ILE" = "I", "LEU" = "L", "LYS" = "K",
                 "MET" = "M", "PHE" = "F", "PRO" = "P", "SER" = "S",
                 "THR" = "T", "TRP" = "W", "TYR" = "Y", "VAL" = "V")
  
  # 批量替换
  data[[col_name]] <- str_replace_all(data[[col_name]], code_dict)
  return(data)
}

测试验证：

# 创建测试数据框
test_df <- data.frame(mutation = c("ALA3=VAL20", "ARG10=GLU5", "PHE15=TYR25"))
# 执行替换
test_df <- substitute_codes(test_df, "mutation")
# 查看结果
print(test_df$mutation)
# 输出：[1] "A3=V20"   "R10=E5"   "F15=Y25"

不依赖第三方包的Base R方案 #

如果不想引入stringr，可以通过循环逐个替换实现：

substitute_codes_base <- function(data, col_name) {
  code_dict <- c("ALA" = "A", "ARG" = "R", "ASN" = "N", "ASP" = "D",
                 "CYS" = "C", "GLU" = "E", "GLN" = "Q", "GLY" = "G",
                 "HIS" = "H", "ILE" = "I", "LEU" = "L", "LYS" = "K",
                 "MET" = "M", "PHE" = "F", "PRO" = "P", "SER" = "S",
                 "THR" = "T", "TRP" = "W", "TYR" = "Y", "VAL" = "V")
  
  # 遍历字典逐个替换
  for (code in names(code_dict)) {
    data[[col_name]] <- gsub(code, code_dict[code], data[[col_name]])
  }
  return(data)
}

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