嗨,我是一名生物学家,也是R编码的新手。我正在尝试为设计师引物编写一个小程序;你为DNA输入了一串字符ATGC,它反转顺序并交换DNA碱基(A变为T,T变为A,G变为C,C变为G)。 我可以让它将相反的DNA返回给字符串中的第一个字符,但它不会遍历整个字符串。任何帮助非常感谢!
primer <- function(input_string){
last_position <- nchar(input_string)
output <- ""
for (position in 1:last_position) {
one_character <- substr(input_string, position, position)
if(one_character=="A"){
new_character <- "T" }
if (one_character=="T"){
new_character <- "A"}
if (one_character=="G"){
new_character <- "C"}
if (one_character=="C"){
new_character <- "G"}
output <- paste(new_character, output, sep="")
return(output)
}
}