Remove specified populations from a genelight {adegenet} object
Individuals are assigned to populations based on the specimen metadata data file (csv) used with gl.read.dart() or gs.read.dart().
gl.drop.pop( x, pop.list, as.pop = NULL, recalc = FALSE, mono.rm = FALSE, verbose = NULL )
x |
– name of the genlight object containing SNP genotypes or Tag P/A data (SilicoDArT) [required] |
pop.list |
– a list of populations to be removed [required] |
as.pop |
– temporarily assign another metric to represent population for the purposes of deletions [default NULL] |
recalc |
– Recalculate the locus metadata statistics [default FALSE] |
mono.rm |
– Remove monomorphic loci [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] |
The script, having deleted populations, optionally identifies resultant monomorphic loci or loci with all values missing and deletes them (using gl.filter.monomorphs.r). The script also optionally recalculates statistics made redundant by the deletion of individuals from the dataset.
The script returns a genlight object with the new population assignments and the recalculated locus metadata.
A genlight object with the reduced data
Arthur Georges (Post to https://groups.google.com/d/forum/dartr)
# SNP data gl2 <- gl.drop.pop(testset.gl, pop.list=c("EmsubRopeMata","EmvicVictJasp")) gl2 <- gl.drop.pop(testset.gl, pop.list=c("EmsubRopeMata","EmvicVictJasp"), mono.rm=TRUE,recalc=TRUE) gl2 <- gl.drop.pop(testset.gl, pop.list=c("Male","Unknown"),as.pop="sex") # Tag P/A data gs2 <- gl.drop.pop(testset.gs, pop.list=c("EmsubRopeMata","EmvicVictJasp"))
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.