ape (version 4.1)

as.alignment: Conversion Among DNA Sequence Internal Formats

Description

These functions transform a set of DNA sequences among various internal formats.

Usage

as.alignment(x)
as.DNAbin(x, ...)

# S3 method for character as.DNAbin(x, ...)

# S3 method for list as.DNAbin(x, ...)

# S3 method for alignment as.DNAbin(x, ...)

# S3 method for DNAString as.DNAbin(x, ...)

# S3 method for DNAStringSet as.DNAbin(x, ...)

# S3 method for PairwiseAlignmentsSingleSubject as.DNAbin(x, ...)

# S3 method for DNAMultipleAlignment as.DNAbin(x, ...)

# S3 method for DNAbin as.character(x, ...)

Arguments

x

a matrix or a list containing the DNA sequences, or an object of class "alignment".

further arguments to be passed to or from other methods.

Value

an object of class "alignment" in the case of "as.alignment"; an object of class "DNAbin" in the case of "as.DNAbin"; a matrix of mode character or a list containing vectors of mode character in the case of "as.character".

Details

For as.alignment, the sequences given as argument should be stored as matrices or lists of single-character strings (the format used in ape before version 1.10). The returned object is in the format used in the package seqinr to store aligned sequences.

as.DNAbin is a generic function with methods so that it works with sequences stored into vectors, matrices, or lists. It can convert some S4 classes from the package Biostrings in BioConductor. For consistency within ape, this uses an S3-style syntax. To convert objects of class "DNAStringSetList", see the examples.

as.character is a generic function: the present method converts objects of class "DNAbin" into the format used before ape 1.10 (matrix of single characters, or list of vectors of single characters). This function must be used first to convert objects of class "DNAbin" into the class "alignment".

See Also

DNAbin, read.dna, read.GenBank, write.dna

Examples

Run this code
# NOT RUN {
data(woodmouse)
x <- as.character(woodmouse)
x[, 1:20]
str(as.alignment(x))
identical(as.DNAbin(x), woodmouse)
### conversion from BioConductor:
# }
# NOT RUN {
if (require(Biostrings)) {
data(phiX174Phage)
X <- as.DNAbin(phiX174Phage)
## base frequencies:
base.freq(X) # from ape
alphabetFrequency(phiX174Phage) # from Biostrings
### for objects of class "DNAStringSetList"
X <- lapply(x, as.DNAbin) # a list of lists
### to put all sequences in a single list:
X <- unlist(X, recursive = FALSE)
class(X) <- "DNAbin"
}
# }

Run the code above in your browser using DataCamp Workspace