residuecount: Total number of residues in an ACNUC list
In seqinr: Biological Sequences Retrieval and Analysis

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/residuecount.R

Computes the total number of residues (nucleotides or aminoacids) in all sequences of the list of specified rank.

1	residuecount(lrank, socket = autosocket())

`lrank`	the list rank on the ACNUC server
`socket`	an object of class `sockconn` connecting to a remote ACNUC database (default is a socket to the last opened database).

A single numeric value corresponding to the total number of residues or NA in case of problem.

J.R. Lobry

http://doua.prabi.fr/databases/acnuc.html

citation("seqinr")

choosebank, query, glr

## Not run:  
  ### Need internet connection
  choosebank("emblTP")
  mylist <- query("mylist", "t=CDS", virtual = TRUE)
  stopifnot(residuecount(glr("mylist")) == 1611439240)
  stopifnot(is.na(residuecount(glr("unknowlist")))) # A warning is issued
  
## End(Not run)