R/countsubseqs.R

In seqinr: Biological Sequences Retrieval and Analysis

# ==>   countsubseqs&lrank=xx
# <==  code=xx&count=xx
# Returns the number of subsequences in list of rank lrank.
# Code != 0 indicates error.

countsubseqs <- function(lrank, socket = autosocket()){
  #
  # Check argument:
  #
  if(!is.finite(lrank)) stop("wrong lrank argument")
  #
  # Build request:
  #
  request <- paste("countsubseqs&lrank=", lrank, sep = "")
  writeLines(request, socket, sep = "\n")
  answerFromServer <- readLines(socket, n = 1)
  #
  # Check that there is an answer from server:
  #
  if(length(answerFromServer) == 0){
    warning("Empty answer from server")
    return(NA)
  }
  #
  # Check that no error is returned:
  #
  resitem <- parser.socket(answerFromServer)
  if(resitem[1] != "0"){
    warning(paste("error code from server:", answerFromServer))
    return(NA)
  }
  #
  return(as.numeric(resitem[2]))
}

css <- countsubseqs