Filter loci in a genlight {adegenet} object based on average repeatability of alleles at a locus
SNP datasets generated by DArT have an index, RepAvg, generated by reproducing the data independently for 30 RepAvg is the proportion of alleles that give a repeatable result, averaged over both alleles for each locus.
gl.filter.RepAvg(x, threshold = 0.99, verbose = NULL)
x |
– name of the genlight object containing the SNP data [required] |
threshold |
– threshold value below which loci will be removed [default 0.99] |
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] |
SilicoDArT datasets generated by DArT have a smilar index, Reproducibility. For these fragment presence/absence data, repeatability is the percentage of scores that are repeated in the technical replicate dataset.
Returns a genlight object retaining loci with repeatability (Repavg or Reproducibility) greater than the specified threshold.
Arthur Georges (Post to https://groups.google.com/d/forum/dartr)
# SNP data gl.report.reproducibility(testset.gl) result <- gl.filter.reproducibility(testset.gl, threshold=0.99, verbose=3) # Tag P/A data gl.report.reproducibility(testset.gs) result <- gl.filter.reproducibility(testset.gs, threshold=0.99)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.