genogeographer (version 0.1.19)

genogeo: 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 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

# NOT RUN {
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)
# }