Reports summary of RepAvg (repeatability averaged over both alleles for each locus) or reproducibility (repeatability of the scores for fragment presence/absence)
Source:R/gl.report.reproducibility.r
gl.report.reproducibility.RdSNP 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.
Usage
gl.report.reproducibility(
x,
plot.out = TRUE,
plot_theme = theme_dartR(),
plot_colors = two_colors,
save2tmp = FALSE,
verbose = NULL
)Arguments
- x
Name of the genlight object containing the SNP or presence/absence (SilicoDArT) data [required].
- plot.out
If TRUE, displays a plot to guide the decision on a filter threshold [default TRUE].
- plot_theme
Theme for the plot. See Details for options [default theme_dartR()].
- plot_colors
List of two color names for the borders and fill of the plots [default two_colors].
- save2tmp
If TRUE, saves any ggplots and listings to the session temporary directory (tempdir) [default FALSE].
- verbose
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].
Details
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.out=TRUE, display also includes a boxplot and a histogram to guide in the selection of a threshold for filtering on repeatability.
If save2tmp=TRUE, ggplots and relevant tabulations are saved to the session's temp directory (tempdir)
For examples of themes, see:
See also
Other report functions:
gl.report.bases(),
gl.report.callrate(),
gl.report.diversity(),
gl.report.hamming(),
gl.report.heterozygosity(),
gl.report.hwe(),
gl.report.ld.map(),
gl.report.locmetric(),
gl.report.maf(),
gl.report.monomorphs(),
gl.report.overshoot(),
gl.report.pa(),
gl.report.parent.offspring(),
gl.report.rdepth(),
gl.report.secondaries(),
gl.report.sexlinked(),
gl.report.taglength()
Author
Custodian: Arthur Georges – Post to https://groups.google.com/d/forum/dartr
Examples
# \donttest{
# SNP data
out <- gl.report.reproducibility(testset.gl)
#> Starting gl.report.reproducibility
#> Processing genlight object with SNP data
#> Reporting Repeatability by Locus
#> No. of loci = 255
#> No. of individuals = 250
#> Minimum : 0.959459
#> 1st quartile : 1
#> Median : 1
#> Mean : 0.9981525
#> 3r quartile : 1
#> Maximum : 1
#> Missing Rate Overall: 0.12
#>
#> Quantile Threshold Retained Percent Filtered Percent
#> 1 100% 1.000000 214 83.9 41 16.1
#> 2 95% 1.000000 214 83.9 41 16.1
#> 3 90% 1.000000 214 83.9 41 16.1
#> 4 85% 1.000000 214 83.9 41 16.1
#> 5 80% 1.000000 214 83.9 41 16.1
#> 6 75% 1.000000 214 83.9 41 16.1
#> 7 70% 1.000000 214 83.9 41 16.1
#> 8 65% 1.000000 214 83.9 41 16.1
#> 9 60% 1.000000 214 83.9 41 16.1
#> 10 55% 1.000000 214 83.9 41 16.1
#> 11 50% 1.000000 214 83.9 41 16.1
#> 12 45% 1.000000 214 83.9 41 16.1
#> 13 40% 1.000000 214 83.9 41 16.1
#> 14 35% 1.000000 214 83.9 41 16.1
#> 15 30% 1.000000 214 83.9 41 16.1
#> 16 25% 1.000000 214 83.9 41 16.1
#> 17 20% 1.000000 214 83.9 41 16.1
#> 18 15% 0.997674 217 85.1 38 14.9
#> 19 10% 0.994536 230 90.2 25 9.8
#> 20 5% 0.984694 243 95.3 12 4.7
#> 21 0% 0.959459 255 100.0 0 0.0
#> Completed: gl.report.reproducibility
#>
# }
# Tag P/A data
out <- gl.report.reproducibility(testset.gs)
#> Starting gl.report.reproducibility
#> Processing genlight object with Presence/Absence (SilicoDArT) data
#> Reporting Repeatability by Locus
#> No. of loci = 255
#> No. of individuals = 218
#> Minimum : 0.952941
#> 1st quartile : 1
#> Median : 1
#> Mean : 0.9983001
#> 3r quartile : 1
#> Maximum : 1
#> Missing Rate Overall: 0.04
#>
#> Quantile Threshold Retained Percent Filtered Percent
#> 1 100% 1.000000 228 89.4 27 10.6
#> 2 95% 1.000000 228 89.4 27 10.6
#> 3 90% 1.000000 228 89.4 27 10.6
#> 4 85% 1.000000 228 89.4 27 10.6
#> 5 80% 1.000000 228 89.4 27 10.6
#> 6 75% 1.000000 228 89.4 27 10.6
#> 7 70% 1.000000 228 89.4 27 10.6
#> 8 65% 1.000000 228 89.4 27 10.6
#> 9 60% 1.000000 228 89.4 27 10.6
#> 10 55% 1.000000 228 89.4 27 10.6
#> 11 50% 1.000000 228 89.4 27 10.6
#> 12 45% 1.000000 228 89.4 27 10.6
#> 13 40% 1.000000 228 89.4 27 10.6
#> 14 35% 1.000000 228 89.4 27 10.6
#> 15 30% 1.000000 228 89.4 27 10.6
#> 16 25% 1.000000 228 89.4 27 10.6
#> 17 20% 1.000000 228 89.4 27 10.6
#> 18 15% 1.000000 228 89.4 27 10.6
#> 19 10% 0.995745 230 90.2 25 9.8
#> 20 5% 0.986607 243 95.3 12 4.7
#> 21 0% 0.952941 255 100.0 0 0.0
#> Completed: gl.report.reproducibility
#>