Demerelate (version 0.8-1)

queller: Queller and Goodnight estimator of relatedness

Description

Calculation of symetric allele sharing using the RE-RAT method and allele frequencies for correction based on a reference population.

Usage

queller(pop1, pop2, allele.column, row.position,
             col.position, matrix.share, ref.pop)

Arguments

pop1

Population data containing allele information. All individual names are found in row names of pairwise relatedness matrix.

pop2

Population data containing allele information. All individual names are found in column names of pairwise relatedness matrix.

allele.column

Column number where loci information is taken from i.e. allele.column and allele.column+1.

row.position

Row position of the pairwise relatedness matrix which contains the name of one individual from population data pop1.

col.position

Column position of the pairwise relatedness matrix which contains the name of one individual from population data pop2.

matrix.share

Pairwise relatedness matrix containing individual names of pop1 and pop2.

ref.pop

Reference population to calculate allele frequencies necessary for the calculation of the estimator.

Value

• share
  value of $r_{xy}$. NA is returned if estimator could not be calculated.

Details

Be aware that unlike in some other programs single individuals are not removed from the calculation to correct for any bias. If you want to use custom reference population to calculate overall allele frequencies, which may be corrected for single individuals or parts of the population, these population need to be constructed manually and indicated by ref.pop. Otherwise it may be feasible to use the whole population information to calculate allele frequency distributions.

References

Oliehoek, P. A. et al. (2006) Estimating relatedness between individuals in general populations with a focus on their use in conservation programs. Genetics, 173, 483-496. Queller, D.C. and Goodnight, K.F. (1989) Estimating relatedness using genetic markers. Evolution, 43, 258-275.

See Also

random.pairs

Examples

## internal function not intended for direct usage