https://github.com/cran/ape
Raw File
Tip revision: c44657c7122de6f1ff126170a82fdac9211b53d5 authored by Emmanuel Paradis on 01 August 2012, 00:00:00 UTC
version 3.0-5
Tip revision: c44657c
seg.sites.Rd
\name{seg.sites}
\alias{seg.sites}
\title{
  Find Segregating Sites in DNA Sequences
}
\usage{
seg.sites(x)
}
\arguments{
  \item{x}{a matrix or a list which contains the DNA sequences.}
}
\description{
  This function gives the indices of segregating (polymorphic) sites in
  a sample of DNA sequences.
}
\details{
  If the sequences are in a list, all the sequences must be of the same
  length. Ambiguous nucleotides are ignored.
}
\value{
  A numeric (integer) vector giving the indices of the segregating
  sites.
}
\author{Emmanuel Paradis}
\note{
  The present version looks for the sites which are ``variable'' in the
  data in terms of different \emph{letters}. This may give unexpected
  results if there are ambiguous bases in the data.
}
\seealso{
  \code{\link{base.freq}}, \code{\link[pegas]{theta.s}},
  \code{\link[pegas]{nuc.div}}
}
\examples{
data(woodmouse)
y <- seg.sites(woodmouse)
y
length(y)
}
\keyword{univar}
back to top