Unweighted estimate of autosomal theta from subpopulations of individual ids

Estimates unweighted autosomal theta for a number of subpopulations given a list of pids (individual ids).

estimate_autotheta_subpops_unweighted_pids(population, subpops, assume_HWE)

Arguments

population

Population obtain from simulation

subpops

List of individual pids

assume_HWE

if the alleles themselves are used instead of genotypes

Value

Estimate of autosomal theta

Details

Assumes that pedigrees_all_populate_autosomal() was used first to populate autosomal genotypes.

Based on Weir and Goudet, Genetics 2017: http://www.genetics.org/content/early/2017/05/26/genetics.116.198424