R/make_otu_names.R

In jfq3/RDPutils: R Utilities for Processing RDPTool Output

make_otu_names <-
function(otu_nums, otu_format="R") {
  new.clstr.name <- matrix("",length(otu_nums),1)
  if(otu_format=="R") {
    n <- nchar(as.character(max(otu_nums)))
    for (i in 1:length(otu_nums)) {
      k <- nchar(as.character(otu_nums[i]))
      k <- n-k
      if (k==0) {pad<-""} else {
        pad <- ""
        for (j in 1:k) {
          pad <- ""
          for (j in 1:k){
            pad <- paste(pad, "0", sep="")
          }
        }
      }
      new.clstr.name[i,1] <- paste("OTU_", pad, as.character(otu_nums[i]), sep="")
    }
  } else if (otu_format == "biom") {
    for (i in 1:length(otu_nums)) {
      new.clstr.name[i,1] <- paste("cluster_", as.character(otu_nums[i]), sep="")
    }
  } else {
    return("otu_format not recognized")
  }

  return(new.clstr.name)
}