Delete Alignment Gaps in DNA Sequences
These functions remove gaps ("-"
) in a sample of DNA sequences.
del.gaps(x) del.colgapsonly(x, threshold = 1, freq.only = FALSE) del.rowgapsonly(x, threshold = 1, freq.only = FALSE)
x |
a matrix, a list, or a vector containing the DNA
sequences; only matrices for |
threshold |
the largest gap proportion to delete the column or row. |
freq.only |
if |
del.gaps
remove all gaps, so the returned sequences may not
have all the same lengths and are therefore returned in a list.
del.colgapsonly
removes the columns with a proportion at least
threshold
of gaps. Thus by default, only the columns with gaps
only are removed (useful when a small matrix is extracted from a large
alignment). del.rowgapsonly
does the same for the rows.
The sequences can be either in "DNAbin"
or in another format,
but the returned object is always of class "DNAbin"
.
del.gaps
returns a vector (if there is only one input sequence)
or a list of class "DNAbin"
; del.colgapsonly
and
del.rowgapsonly
return a matrix of class "DNAbin"
or a
numeric vector (with names for the second function) if freq.only
= TRUE
.
Emmanuel Paradis
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.