Estimate autosomal theta from genotypes

Source: R/RcppExports.R

Estimate autosomal theta for one subpopulation given a sample of genotypes.

estimate_autotheta_1subpop_genotypes(genotypes, return_estimation_info = FALSE)

Arguments

genotypes

Matrix of genotypes: two columns (allele1 and allele2) and a row per individual

return_estimation_info

Whether to return the quantities used to estimate theta

Value

List:

  • theta

    • estimate: Vector of length 1 containing estimate of theta or NA if it could not be estimated

    • error: true if an error happened, false otherwise

    • details: contains description if an error happened

    • estimation_info: If return_estimation_info = true: a list with information used to estimate theta. Else NULL.

Details

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