我试图把一条dna链的互补序列存储在一个载体中。它应该把“A”换成“T”,把“C”换成“G”,反之亦然。问题是,我需要在第一个载体上发生这种情况,并正确地打印互补序列。这是我尝试过的,但卡住了:
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向量是字符类型,并且函数显示没有执行错误。
3条答案
按热度按时间db2dz4w81#
这是
chartr
函数的理想用例:输出:
oxalkeyp2#
您可以使用
purrr::map
执行此操作:wgeznvg73#
您可以使用dplyr中的
recode
或者在基R中,创建一个命名向量并使用
match
匹配命名向量中的值位置,然后调用name
获取名称x一个一个一个一个x一个一个二个x