genogeo: Likelihood ratio tests for AIMs

In genogeographer: Methods for Analysing Forensic Ancestry Informative Markers

Description

Computes the likelihood ratio test statistics for each population in a database of reference populations.

Usage

genogeo(profile, df, CI = 0.95, min_n = 75, grouping = "pop",
  tilt = FALSE, ...)

Arguments

profile	The AIMs profile encoded as returned by the profile_AA_x0 function.
df	The database of reference populations as returned by the pops_to_DB function.
CI	The confidence level used to reject or accept the various hypotheses (between 0 and 1).
min_n	Minimum number of individuals in each database sample
grouping	should "pop" (the default) or "meta" be used for aggregating the results. Can also be "cluster" if this variable is defined in the input database.
tilt	Should exponential titling be used to obtain more accurate $p$-values in the distribution's tail (currently not implemented)
...	Further arguments that are passed to other functions

Value

A tibble containing the $z$-scores, $p$-values etc for each population.

Examples

df_ <- simulate_pops(pop_n = 20, aims_n = 50)
df_db <- pops_to_DB(df_)
profile <- random_AIMs_profile(df_db, keep_pop = TRUE)
profile$pop[1] # The true population
result <- genogeo(profile[,c("locus","x0")], df = df_db)

genogeographer documentation built on Sept. 27, 2019, 5:03 p.m.