Subsample n loci from a genlight object and return as a genlight object
This is a support script, to subsample a genlight {adegenet} object based on loci. Two methods are used to subsample, random and based on information content (avgPIC).
gl.subsample.loci(x, n, method = "random", mono.rm = FALSE, verbose = NULL)
x |
– name of the genlight object containing the SNP data [required] |
n |
– number of loci to include in the subsample [required] |
method |
– "random", in which case the loci are sampled at random; or avgPIC, in which case the top n loci ranked on information content (AvgPIC) are chosen [default 'random'] |
mono.rm |
– delete monomorphic loci before sampling [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 or as specified using gl.set.verbosity] |
A genlight object with n loci
Arthur Georges (Post to https://groups.google.com/d/forum/dartr)
# SNP data gl2 <- gl.subsample.loci(testset.gl, n=200, method="pic") # Tag P/A data gl2 <- gl.subsample.loci(testset.gl, n=100, method="random")
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.