multidna2alignment: Convert from multidna into alignment (seqinr)

In apex: Phylogenetic Methods for Multiple Gene Data

Description

The functions multidna2alignment and multiphyDat2alignment concatenates separate sequences and return an alignment object of the seqinr package.

Usage

multidna2alignment(x, genes = TRUE)

multiphyDat2alignment(x, genes = TRUE)

Arguments

x	a multidna or multiphyDat object.
genes	an optional vector indicating the genes to retain for the concatenation; any way to subset the list in x@dna or x@seq is acceptable; by default, all genes are used.

Value

a alignment object

Author(s)

Thibaut Jombart [e-mail], Zhian N. Kamvar, Klaus Schliep

See Also

• concatenate

• as.alignment to convert single DNAbin objects.

Examples

## Not run: 
## simple conversion with nicely ordered output
data(woodmouse)
genes <- list(gene1=woodmouse[,1:500], gene2=woodmouse[,501:965])
x <- new("multidna", genes)
x
y <- multidna2alignment(x)
y
x2 <- multidna2multiphyDat(x)
z <- multiphyDat2alignment(x2)

## End(Not run)

apex documentation built on April 14, 2020, 5:32 p.m.