popkin {popkin} | R Documentation |
Given the biallelic genotypes of n individuals, this function returns the n-by-n kinship matrix such that the kinship estimate between the most distant subpopulations is zero on average (this sets the ancestral population T to the most recent common ancestor population).
popkin(X, subpops = NULL, n = NA, loci_on_cols = FALSE, mem_lim = NA, lociOnCols = FALSE, memLim = NA)
X |
Genotype matrix, BEDMatrix object, or a function X(m) that returns the genotypes of all individuals at m successive locus blocks each time it is called, and NULL when no loci are left. |
subpops |
The length-n vector of subpopulation assignments for each individual. If missing, every individual is effectively treated as a different population. |
n |
Number of individuals (required only when X is a function, ignored otherwise).
If n is missing but |
loci_on_cols |
If |
mem_lim |
Memory limit in GB, used to break up genotype data into chunks for very large datasets. Note memory usage is somewhat underestimated and is not controlled strictly. Default in Linux and Windows is 70 % of the free system memory, otherwise it is 1GB (OSX and other systems). DEPRECATED PARAMETER NAMES. These generate a warning that they are deprecated, and will be removed in a future release. |
lociOnCols |
Same as |
memLim |
Same as |
The subpopulation assignments are only used to estimate the baseline kinship (the zero value).
If the user wants to re-estimate the kinship matrix using different subpopulation labels,
it suffices to rescale it using rescale_popkin
(as opposed to starting from the genotypes again, which gives the same answer less efficiently).
The matrix X must have values only in c(0,1,2,NA)
, encoded to count the number of reference alleles at the locus, or NA
for missing data.
The estimated n-by-n kinship matrix. If X has names for the individuals, they will be copied to the rows and columns of this kinship matrix.
# Construct toy data X <- matrix(c(0,1,2,1,0,1,1,0,2), nrow = 3, byrow = TRUE) # genotype matrix subpops <- c(1,1,2) # subpopulation assignments for individuals # NOTE: for BED-formatted input, use BEDMatrix! # "file" is path to BED file (excluding .bed extension) ## library(BEDMatrix) ## X <- BEDMatrix(file) # load genotype matrix object kinship <- popkin(X, subpops) # calculate kinship from genotypes and subpopulation labels