LR_table: Compute pairwise likelihood ratios

Description Usage Arguments Value Author(s) Examples

View source: R/LR_table.R

Description

For each pair of a specified vector of profiles the likelihood ratios are computed. The list can include all populations in the data or only a subset. We may for inferral purposes restrict to ratios including at least one "accepted" population.

Usage

1
2
LR_table(result_df, lr_populations = NULL, only_accepted = TRUE,
  CI = 0.95, digits = NULL, keep_logP = FALSE)

Arguments

result_df

The output from genogeo

lr_populations

A vector of population names (pop in result_df). If NULL all populations are used.

only_accepted

Restrict the ratios to include minimum one accepted population.

CI

The level of confidence interval to be computed

digits

If rounding of the output should be performed.

keep_logP

Logical. Should the logP's be returned in output

Value

A tibble with numerator and denominator populations with their log10 LR and uncertainty.

Author(s)

Torben Tvedebrink tvede@math.aau.dk

Examples

1
2
3
4
5
6
df_ <- simulate_pops(pop_n = 4, 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)
LR_table(result)

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