r - 使用 gsub 的互补序列

Question

我正在尝试将 DNA 链的互补序列存储在一个向量中。它应该改变“T”的“A”和“G”的“C”，反之亦然，问题是，我需要这发生在第一个向量上并正确打印互补序列。这是我尝试但卡住的方法:

pilot_sequence <- c("C","G","A","T","C","C","T","A","T")

complement_sequence_display <- function(pilot_sequence){
  complement_chain_Incom <- gsub("A", "T", pilot_sequence)
  complement_chain <- paste(complement_chain_Incom, collapse = "")
  cat("Complement sequence: ", complement_chain, "\n")
} 
complement_chain_Incom <- gsub("A","T", pilot_sequence)
complement_chain <- paste(complement_chain_Incom, collapse= "")
complement_sequence_display(pilot_sequence)

我得到的答案是:CGTTCCTTT，只有第二个和倒数第二个 T 是正确的，我如何解决其余字母？

pilot_sequence向量为字符类型，函数没有执行错误。

Answer 1

这是 chartr 函数的理想用例:

chartr("ATGC","TACG",pilot_sequence)

输出:

[1] "G" "C" "T" "A" "G" "G" "A" "T" "A"