R/103-getDNAGenbank.R

In BioMedR: Generating Various Molecular Representations for Chemicals, Proteins, DNAs/RNAs and Their Interactions

#' Get DNA/RNA Sequences from Genbank by GI ID
#'
#' Get DNA/RNA Sequences from Genbank by GI ID
#'
#' This function get DNA/RNA sequences from Genbank by GI ID(s).
#'
#' @param id A character vector, as the GI ID(s).
#'
#' @return A list, each component contains one of the DNA/RNA sequences.
#'
#' @keywords Genbank
#'
#' @aliases BMgetDNAGenBank
#'
#' @author Min-feng Zhu <\email{wind2zhu@@163.com}>
#'
#' @seealso See \code{\link{readFASTA}} for reading FASTA format files.
#'
#' @export BMgetDNAGenBank
#' 
#' @importFrom RCurl getURL
#'
#' @examples
#' \donttest{
#' # Network latency may slow down this example
#' # Only test this when your connection is fast enough
#' 
#' ids = c(2, 11)
#' BMgetDNAGenBank(ids)}


BMgetDNAGenBank = function (id) {
  id = as.character(id)
  
  n = length(id)
  
  dna = vector('list', n)
  
  for (i in 1:n) {
    url = paste('https://www.ncbi.nlm.nih.gov/sviewer/?db=nuccore&val=', id[i], 
              '&page_size=5&fmt_mask=0&report=fasta&retmode=text&page=1&page_size=5', sep = '')
    genb = RCurl::getURL(url) 
    sequences = strsplit(genb[[1]], split = "\n")[[1]]
    start = 2
    end = length(sequences)
    dna[[i]]=  paste(sequences[start:end], collapse = "")
  }
  
  gi_name = lapply(1:n, function(i) paste('gi', id[i], sep = "_"))
 
  names(dna) = gi_name
  return(dna)
}