seqinr (version 3.3-3)

residuecount: Total number of residues in an ACNUC list

Description

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

Usage

residuecount(lrank, socket = autosocket())

Arguments

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).

Value

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

References

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

citation("seqinr")

See Also

choosebank, query, glr

Examples

## 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)