ENc.values {vhcub}

Effective Number of Codons (ENc).

Description

Measure the Effective Number of Codons (ENc) of DNA sequence. Using its modified version (Novembre, 2002).

Usage

ENc.values(df.fasta,genetic.code = "1",threshold=0)

Arguments

df.fasta

a data frame with seq_name and its DNA sequence.

genetic.code

a single string that uniquely identifies a genetic code to use.

threshold

optional numeric, specifying sequence length, in codons, used for filtering.

Details

For more information about ENc Novembre, 2002.

Value

A data.frame containing the computed ENc values for each DNA sequences within df.fasta.

Author(s)

Ali Mostafa Anwar [e-mail] and Mohmed Soudy [e-mail]

Examples

# read DNA from fasta file
fasta <- fasta.read("virus.fasta", "host.fasta")
fasta.v <- fasta[[1]]
fasta.h <- fasta[[2]]


# Calculate ENc
enc.df.v <- ENc.values(fasta.v)

enc.df.h <- ENc.values(fasta.h)