miga-base 1.2.17.0 → 1.2.17.2

data/utils/enveomics/enveomics.R/man/enve.recplot2.__peakHist.Rd

@@ -0,0 +1,24 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.__peakHist}
+\alias{enve.recplot2.__peakHist}
+\title{Enveomics: Recruitment Plot (2) Peak S4 Class - Internal Ancillary Function}
+\usage{
+enve.recplot2.__peakHist(x, mids, counts = TRUE)
+}
+\arguments{
+\item{x}{\code{\link{enve.RecPlot2.Peak}} object}
+
+\item{mids}{Midpoints}
+
+\item{counts}{Counts}
+}
+\value{
+A numeric vector of counts (histogram)
+}
+\description{
+Internal ancillary function (see \code{\link{enve.RecPlot2.Peak}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.__whichClosestPeak.Rd

@@ -0,0 +1,22 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.__whichClosestPeak}
+\alias{enve.recplot2.__whichClosestPeak}
+\title{Enveomics: Recruitment Plot (2) Peak Finder - Internal Ancillary Function}
+\usage{
+enve.recplot2.__whichClosestPeak(peak, peaks)
+}
+\arguments{
+\item{peak}{Query \code{\link{enve.RecPlot2.Peak}} object}
+
+\item{peaks}{list of \code{\link{enve.RecPlot2.Peak}} objects}
+}
+\value{
+A numeric index out of \code{peaks}.
+}
+\description{
+Internal ancillary function (see \code{\link{enve.recplot2.findPeaks}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.changeCutoff.Rd

@@ -0,0 +1,22 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.changeCutoff}
+\alias{enve.recplot2.changeCutoff}
+\title{Enveomics: Recruitment Plot (2) Change Cutoff}
+\usage{
+enve.recplot2.changeCutoff(rp, new.cutoff = 98)
+}
+\arguments{
+\item{rp}{\code{\link{enve.RecPlot2}} object.}
+
+\item{new.cutoff}{New cutoff to use.}
+}
+\value{
+The modified \code{\link{enve.RecPlot2}} object.
+}
+\description{
+Change the intra-species cutoff of an existing recruitment plot.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.compareIdentities.Rd

@@ -0,0 +1,52 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.compareIdentities}
+\alias{enve.recplot2.compareIdentities}
+\title{Enveomics: Recruitment Plot (2) Compare Identities}
+\usage{
+enve.recplot2.compareIdentities(
+  x,
+  y,
+  method = "hellinger",
+  smooth.par = NULL,
+  pseudocounts = 0,
+  max.deviation = 0.75
+)
+}
+\arguments{
+\item{x}{First \code{\link{enve.RecPlot2}} object.}
+
+\item{y}{Second \code{\link{enve.RecPlot2}} object.}
+
+\item{method}{Distance method to use. This should be (an unambiguous abbreviation of)
+one of:
+\itemize{
+   \item{"hellinger"
+         (\emph{Hellinger, 1090, doi:10.1515/crll.1909.136.210}),}
+   \item{"bhattacharyya"
+         (\emph{Bhattacharyya, 1943, Bull. Calcutta Math. Soc. 35}),}
+   \item{"kl" or "kullback-leibler"
+         (\emph{Kullback & Leibler, 1951, doi:10.1214/aoms/1177729694}), or}
+   \item{"euclidean"}
+}}
+
+\item{smooth.par}{Smoothing parameter for cubic spline smoothing. Use 0 for no smoothing.
+Use \code{NULL} to automatically determine this value using leave-one-out
+cross-validation (see \code{smooth.spline} parameter \code{spar}).}
+
+\item{pseudocounts}{Smoothing parameter for Laplace smoothing. Use 0 for no smoothing, or
+1 for add-one smoothing.}
+
+\item{max.deviation}{Maximum mean deviation between identity breaks tolerated (as percent
+identity). Difference in number of \code{id.breaks} is never tolerated.}
+}
+\value{
+A \strong{numeric} indicating the distance between the objects.
+}
+\description{
+Compare the distribution of identities between two
+\code{\link{enve.RecPlot2}} objects.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.coordinates.Rd

@@ -0,0 +1,29 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.coordinates}
+\alias{enve.recplot2.coordinates}
+\title{Enveomics: Recruitment Plot (2) Coordinates}
+\usage{
+enve.recplot2.coordinates(x, bins)
+}
+\arguments{
+\item{x}{\code{\link{enve.RecPlot2}} object.}
+
+\item{bins}{Vector of selected bins to return. It can be a vector of logical values
+with the same length as \code{x$pos.breaks-1} or a vector of integers. If
+missing, returns the coordinates of all windows.}
+}
+\value{
+Returns a data.frame with five columns: \code{name.from} (character),
+\code{pos.from} (numeric), \code{name.to} (character), \code{pos.to}
+(numeric), and \code{seq.name} (character).
+The first two correspond to sequence and position of the start point of the
+bin. The next two correspond to the sequence and position of the end point of
+the bin. The last one indicates the name of the sequence (if defined).
+}
+\description{
+Returns the sequence name and coordinates of the requested position bins.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.corePeak.Rd

@@ -0,0 +1,21 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.corePeak}
+\alias{enve.recplot2.corePeak}
+\title{Enveomics: Recruitment Plot (2) Core Peak Finder}
+\usage{
+enve.recplot2.corePeak(x)
+}
+\arguments{
+\item{x}{\code{list} of \code{\link{enve.RecPlot2.Peak}} objects.}
+}
+\value{
+A \code{\link{enve.RecPlot2.Peak}} object.
+}
+\description{
+Finds the peak in a list of peaks that is most likely to represent the
+"core genome" of a population.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.extractWindows.Rd

@@ -0,0 +1,45 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.extractWindows}
+\alias{enve.recplot2.extractWindows}
+\title{Enveomics: Recruitment Plot (2) Extract Windows}
+\usage{
+enve.recplot2.extractWindows(
+  rp,
+  peak,
+  lower.tail = TRUE,
+  significance = 0.05,
+  seq.names = FALSE
+)
+}
+\arguments{
+\item{rp}{Recruitment plot, a \code{\link{enve.RecPlot2}} object.}
+
+\item{peak}{Peak, an \code{\link{enve.RecPlot2.Peak}} object. If list, it is assumed to
+be a list of \code{\link{enve.RecPlot2.Peak}} objects, in which case the core
+peak is used (see \code{\link{enve.recplot2.corePeak}}).}
+
+\item{lower.tail}{If \code{FALSE}, it returns windows significantly above the peak in
+sequencing depth.}
+
+\item{significance}{Significance threshold (alpha) to select windows.}
+
+\item{seq.names}{Returns subject sequence names instead of a vector of Booleans. If
+the recruitment plot was generated with named position bins (e.g, using
+\code{pos.breaks=0} or a two-column \code{pos.breaks.tsv}), it returns a
+vector of characters (the sequence identifiers), otherwise it returns a
+data.frame with a name column and two columns of coordinates.}
+}
+\value{
+Returns a vector of logicals if \code{seq.names = FALSE}.
+If \code{seq.names = TRUE}, it returns a data.frame with five columns:
+\code{name.from}, \code{name.to}, \code{pos.from}, \code{pos.to}, and
+\code{seq.name} (see \code{\link{enve.recplot2.coordinates}}).
+}
+\description{
+Extract windows significantly below (or above) the peak in sequencing
+depth.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.Rd

@@ -0,0 +1,34 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks}
+\alias{enve.recplot2.findPeaks}
+\title{Enveomics: Recruitment Plot (2) Peak Finder}
+\usage{
+enve.recplot2.findPeaks(x, method = "emauto", ...)
+}
+\arguments{
+\item{x}{An \code{\link{enve.RecPlot2}} object.}
+
+\item{method}{Peak-finder method. This should be one of:
+\itemize{
+   \item \strong{emauto}
+   (Expectation-Maximization with auto-selection of components)
+   \item \strong{em}
+   (Expectation-Maximization)
+   \item \strong{mower}
+   (Custom distribution-mowing method)
+}}
+
+\item{...}{Any additional parameters supported by
+\code{\link{enve.recplot2.findPeaks}}.}
+}
+\value{
+Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
+}
+\description{
+Identifies peaks in the population histogram potentially indicating
+sub-population mixtures.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_e.Rd

@@ -0,0 +1,23 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.__em_e}
+\alias{enve.recplot2.findPeaks.__em_e}
+\title{Enveomics: Recruitment Plot (2) EM Peak Finder - Internal Ancillary Function Expectation}
+\usage{
+enve.recplot2.findPeaks.__em_e(x, theta)
+}
+\arguments{
+\item{x}{Vector of log-transformed sequencing depths}
+
+\item{theta}{Parameters list}
+}
+\value{
+A list with components \code{ll} (numeric) the log-likelihood, and
+\code{posterior} (numeric) the posterior probability.
+}
+\description{
+Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.em}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_m.Rd

@@ -0,0 +1,24 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.__em_m}
+\alias{enve.recplot2.findPeaks.__em_m}
+\title{Enveomics: Recruitment Plot (2) Em Peak Finder - Internal Ancillary Function Maximization}
+\usage{
+enve.recplot2.findPeaks.__em_m(x, posterior)
+}
+\arguments{
+\item{x}{Vector of log-transformed sequencing depths}
+
+\item{posterior}{Posterior probability}
+}
+\value{
+A list with components \code{mu} (numeric) the estimated mean,
+\code{sd} (numeric) the estimated standard deviation, and \code{alpha}
+(numeric) the estimated alpha parameter.
+}
+\description{
+Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.em}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__emauto_one.Rd

@@ -0,0 +1,31 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.__emauto_one}
+\alias{enve.recplot2.findPeaks.__emauto_one}
+\title{Enveomics: Recruitment Plot (2) EMauto Peak Finder - Internal Ancillary Function}
+\usage{
+enve.recplot2.findPeaks.__emauto_one(x, comp, do_crit, best, verbose, ...)
+}
+\arguments{
+\item{x}{\code{\link{enve.RecPlot2}} object.}
+
+\item{comp}{Components.}
+
+\item{do_crit}{Function estimating the criterion.}
+
+\item{best}{Best solution thus far.}
+
+\item{verbose}{If verbose.}
+
+\item{...}{Additional parameters for \code{\link{enve.recplot2.findPeaks.em}}.}
+}
+\value{
+Updated solution with the same structure as \code{best}.
+}
+\description{
+Internal ancillary function (see
+\code{\link{enve.recplot2.findPeaks.emauto}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mow_one.Rd

@@ -0,0 +1,56 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.__mow_one}
+\alias{enve.recplot2.findPeaks.__mow_one}
+\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder - Internal Ancillary Function 1}
+\usage{
+enve.recplot2.findPeaks.__mow_one(
+  lsd1,
+  min.points,
+  quant.est,
+  mlv.opts,
+  fitdist.opts,
+  with.skewness,
+  optim.rounds,
+  optim.epsilon,
+  n.total,
+  merge.logdist,
+  verbose,
+  log
+)
+}
+\arguments{
+\item{lsd1}{Vector of log-transformed sequencing depths}
+
+\item{min.points}{Minimum number of points}
+
+\item{quant.est}{Quantile estimate}
+
+\item{mlv.opts}{List of options for \code{mlv}}
+
+\item{fitdist.opts}{List of options for \code{fitdist}}
+
+\item{with.skewness}{If skewed-normal should be used}
+
+\item{optim.rounds}{Maximum number of optimization rounds}
+
+\item{optim.epsilon}{Minimum difference considered negligible}
+
+\item{n.total}{Global number of windows}
+
+\item{merge.logdist}{Attempted \code{merge.logdist} parameter}
+
+\item{verbose}{If verbose}
+
+\item{log}{If log-transformed depths}
+}
+\value{
+Return an \code{enve.RecPlot2.Peak} object.
+}
+\description{
+Internal ancillary function (see
+\code{\link{enve.recplot2.findPeaks.mower}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mower.Rd

@@ -0,0 +1,20 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.__mower}
+\alias{enve.recplot2.findPeaks.__mower}
+\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder - Internal Ancillary Function 2}
+\usage{
+enve.recplot2.findPeaks.__mower(peaks.opts)
+}
+\arguments{
+\item{peaks.opts}{List of options for \code{\link{enve.recplot2.findPeaks.__mow_one}}}
+}
+\value{
+A list of \code{enve.RecPlot2.Peak} objects.
+}
+\description{
+Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.mower}}).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.em.Rd

@@ -0,0 +1,51 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.em}
+\alias{enve.recplot2.findPeaks.em}
+\title{Enveomics: Recruitment Plot (2) Em Peak Finder}
+\usage{
+enve.recplot2.findPeaks.em(
+  x,
+  max.iter = 1000,
+  ll.diff.res = 1e-08,
+  components = 2,
+  rm.top = 0.05,
+  verbose = FALSE,
+  init,
+  log = TRUE
+)
+}
+\arguments{
+\item{x}{An \code{\link{enve.RecPlot2}} object.}
+
+\item{max.iter}{Maximum number of EM iterations.}
+
+\item{ll.diff.res}{Maximum Log-Likelihood difference to be considered as convergent.}
+
+\item{components}{Number of distributions assumed in the mixture.}
+
+\item{rm.top}{Top-values to remove before finding peaks, as a quantile probability.
+This step is useful to remove highly conserved regions, but can be
+turned off by setting \code{rm.top=0}. The quantile is determined
+\strong{after} removing zero-coverage windows.}
+
+\item{verbose}{Display (mostly debugging) information.}
+
+\item{init}{Initialization parameters. By default, these are derived from k-means
+clustering. A named list with vectors for \code{mu}, \code{sd}, and
+\code{alpha}, each of length \code{components}.}
+
+\item{log}{Logical value indicating if the estimations should be performed in
+natural logarithm units. Do not change unless you know what you're
+doing.}
+}
+\value{
+Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
+}
+\description{
+Identifies peaks in the population histogram using a Gaussian Mixture
+Model Expectation Maximization (GMM-EM) method.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.emauto.Rd

@@ -0,0 +1,43 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.emauto}
+\alias{enve.recplot2.findPeaks.emauto}
+\title{Enveomics: Recruitment Plot (2) Emauto Peak Finder}
+\usage{
+enve.recplot2.findPeaks.emauto(
+  x,
+  components = seq(1, 5),
+  criterion = "aic",
+  merge.tol = 2L,
+  verbose = FALSE,
+  ...
+)
+}
+\arguments{
+\item{x}{An \code{\link{enve.RecPlot2}} object.}
+
+\item{components}{A vector of number of components to evaluate.}
+
+\item{criterion}{Criterion to use for components selection. Must be one of:
+\code{aic} (Akaike Information Criterion), \code{bic} or \code{sbc}
+(Bayesian Information Criterion or Schwarz Criterion).}
+
+\item{merge.tol}{When attempting to merge peaks with very similar sequencing depth, use
+this number of significant digits (in log-scale).}
+
+\item{verbose}{Display (mostly debugging) information.}
+
+\item{...}{Any additional parameters supported by
+\code{\link{enve.recplot2.findPeaks.em}}.}
+}
+\value{
+Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
+}
+\description{
+Identifies peaks in the population histogram using a Gaussian Mixture
+Model Expectation Maximization (GMM-EM) method with number of components
+automatically detected.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.mower.Rd

@@ -0,0 +1,82 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.findPeaks.mower}
+\alias{enve.recplot2.findPeaks.mower}
+\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder}
+\usage{
+enve.recplot2.findPeaks.mower(
+  x,
+  min.points = 10,
+  quant.est = c(0.002, 0.998),
+  mlv.opts = list(method = "parzen"),
+  fitdist.opts.sn = list(distr = "sn", method = "qme", probs = c(0.1, 0.5, 0.8), start
+    = list(omega = 1, alpha = -1), lower = c(0, -Inf, -Inf)),
+  fitdist.opts.norm = list(distr = "norm", method = "qme", probs = c(0.4, 0.6), start =
+    list(sd = 1), lower = c(0, -Inf)),
+  rm.top = 0.05,
+  with.skewness = TRUE,
+  optim.rounds = 200,
+  optim.epsilon = 1e-04,
+  merge.logdist = log(1.75),
+  verbose = FALSE,
+  log = TRUE
+)
+}
+\arguments{
+\item{x}{An \code{\link{enve.RecPlot2}} object.}
+
+\item{min.points}{Minimum number of points in the quantile-estimation-range
+\code{(quant.est)} to estimate a peak.}
+
+\item{quant.est}{Range of quantiles to be used in the estimation of a peak's
+parameters.}
+
+\item{mlv.opts}{Ignored. For backwards compatibility.}
+
+\item{fitdist.opts.sn}{Options passed to \code{fitdist} to estimate the standard deviation if
+\code{with.skewness=TRUE}. Note that the \code{start} parameter will be
+ammended with \code{xi=estimated} mode for each peak.}
+
+\item{fitdist.opts.norm}{Options passed to \code{fitdist} to estimate the standard deviation if
+\code{with.skewness=FALSE}. Note that the \code{start} parameter will be
+ammended with \code{mean=estimated} mode for each peak.}
+
+\item{rm.top}{Top-values to remove before finding peaks, as a quantile probability.
+This step is useful to remove highly conserved regions, but can be
+turned off by setting \code{rm.top=0}. The quantile is determined
+\strong{after} removing zero-coverage windows.}
+
+\item{with.skewness}{Allow skewness correction of the peaks. Typically, the
+sequencing-depth distribution for a single peak is left-skewed, due
+partly (but not exclusively) to fragmentation and mapping sensitivity.
+See \emph{Lindner et al 2013, Bioinformatics 29(10):1260-7} for an
+alternative solution for the first problem (fragmentation) called
+"tail distribution".}
+
+\item{optim.rounds}{Maximum rounds of peak optimization.}
+
+\item{optim.epsilon}{Trace change at which optimization stops (unless \code{optim.rounds} is
+reached first). The trace change is estimated as the sum of square
+differences between parameters in one round and those from two rounds
+earlier (to avoid infinite loops from approximation).}
+
+\item{merge.logdist}{Maximum value of \code{|log-ratio|} between centrality parameters in peaks
+to attempt merging. The default of ~0.22 corresponds to a maximum
+difference of 25\%.}
+
+\item{verbose}{Display (mostly debugging) information.}
+
+\item{log}{Logical value indicating if the estimations should be performed in
+natural logarithm units. Do not change unless you know what you're
+doing.}
+}
+\value{
+Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
+}
+\description{
+Identifies peaks in the population histogram potentially indicating
+sub-population mixtures, using a custom distribution-mowing method.
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.peak-class.Rd

@@ -0,0 +1,59 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\docType{class}
+\name{enve.RecPlot2.Peak-class}
+\alias{enve.RecPlot2.Peak-class}
+\alias{enve.RecPlot2.Peak}
+\title{Enveomics: Recruitment Plot (2) Peak - S4 Class}
+\description{
+Enve-omics representation of a peak in the sequencing depth histogram
+of a Recruitment plot (see \code{\link{enve.recplot2.findPeaks}}).
+}
+\section{Slots}{
+
+\describe{
+\item{\code{dist}}{\code{(character)}
+Distribution of the peak. Currently supported: \code{norm} (normal) and \code{sn}
+(skew-normal).}
+
+\item{\code{values}}{\code{(numeric)}
+Sequencing depth values predicted to conform the peak.}
+
+\item{\code{values.res}}{\code{(numeric)}
+Sequencing depth values not explained by this or previously identified
+peaks.}
+
+\item{\code{mode}}{\code{(numeric)}
+Seed-value of mode anchoring the peak.}
+
+\item{\code{param.hat}}{\code{(list)}
+Parameters of the distribution. A list of two values if dist=\code{norm} (sd
+and mean), or three values if dist=\code{sn}(omega=scale, alpha=shape, and
+xi=location). Note that the "dispersion" parameter is always first and
+the "location" parameter is always last.}
+
+\item{\code{n.hat}}{\code{(numeric)}
+Number of bins estimated to be explained by this peak. This should
+ideally be equal to the length of  \code{values}, but it's not an integer.}
+
+\item{\code{n.total}}{\code{(numeric)}
+Total number of bins from which the peak was extracted. I.e., total
+number of position bins with non-zero sequencing depth in the recruitment
+plot (regardless of peak count).}
+
+\item{\code{err.res}}{\code{(numeric)}
+Error left after adding the peak (mower) or log-likelihood (em or emauto).}
+
+\item{\code{merge.logdist}}{\code{(numeric)}
+Attempted \code{merge.logdist} parameter.}
+
+\item{\code{seq.depth}}{\code{(numeric)}
+Best estimate available for the sequencing depth of the peak (centrality).}
+
+\item{\code{log}}{\code{(logical)}
+Indicates if the estimation was performed in natural logarithm space.}
+}}
+
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}

data/utils/enveomics/enveomics.R/man/enve.recplot2.seqdepth.Rd

@@ -0,0 +1,27 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/recplot2.R
+\name{enve.recplot2.seqdepth}
+\alias{enve.recplot2.seqdepth}
+\title{Enveomics: Recruitment Plot (2) Sequencing Depth}
+\usage{
+enve.recplot2.seqdepth(x, sel, low.identity = FALSE)
+}
+\arguments{
+\item{x}{\code{\link{enve.RecPlot2}} object.}
+
+\item{sel}{Window(s) for which the sequencing depth is to be calculated. If not
+passed, it returns the sequencing depth of all windows.}
+
+\item{low.identity}{A logical indicating if the sequencing depth is to be estimated only
+with low-identity matches. By default, only high-identity matches are
+used.}
+}
+\value{
+Returns a numeric vector of sequencing depths (in bp/bp).
+}
+\description{
+Calculate the sequencing depth of the given window(s).
+}
+\author{
+Luis M. Rodriguez-R [aut, cre]
+}