SNP datasets generated by DArT have an index, RepAvg, generated by reproducing the data independently for 30 of alleles that give a repeatable result, averaged over both alleles for each locus. In the case of fragment presence/absence data (SilicoDArT), repeatability is the percentage of scores that are repeated in the technical replicate dataset.
gl.report.reproducibility(
x,
plot.display = TRUE,
plot.theme = theme_dartR(),
plot.colors = NULL,
plot.dir = NULL,
plot.file = NULL,
verbose = NULL
)
An unaltered genlight object
Name of the genlight object containing the SNP or presence/absence (SilicoDArT) data [required].
Specify if plot is to be produced [default TRUE].
Theme for the plot. See Details for options [default theme_dartR()].
Vector with two color names for the borders and fill [default c("#2171B5", "#6BAED6")].
Directory to save the plot RDS files [default as specified by the global working directory or tempdir()]
Filename (minus extension) for the RDS plot file [Required for plot save]
Verbosity: 0, silent or fatal errors; 1, begin and end; 2, progress log; 3, progress and results summary; 5, full report [default 2, unless specified using gl.set.verbosity].
Custodian: Arthur Georges -- Post to https://groups.google.com/d/forum/dartr
The function displays a table of minimum, maximum, mean and quantiles for
repeatbility against possible thresholds that might subsequently be
specified in gl.filter.reproducibility
.
If plot.display=TRUE, display also includes a boxplot and a histogram to guide
in the selection of a threshold for filtering on repeatability.
If plot.file is specified, plots are saved to the directory specified by the user, or the global
default working directory set by gl.set.wd() or to the tempdir().
For examples of themes, see:
gl.filter.reproducibility
Other matched report:
gl.filter.excess.het()
,
gl.report.allna()
,
gl.report.callrate()
,
gl.report.hamming()
,
gl.report.locmetric()
,
gl.report.maf()
,
gl.report.overshoot()
,
gl.report.pa()
,
gl.report.rdepth()
,
gl.report.secondaries()
,
gl.report.taglength()
# \donttest{
# SNP data
out <- gl.report.reproducibility(testset.gl)
# }
# Tag P/A data
out <- gl.report.reproducibility(testset.gs)
Run the code above in your browser using DataLab