genogeo {genogeographer} | R Documentation |
Likelihood ratio tests for AIMs
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 |
df |
The database of reference populations as returned by the |
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 |
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)