Calculate the geometric median of the posterior distribution of allele frequencies

Calculate the geometric median of the posterior distribution of allele frequencies

Usage

calculate_med_allele_freqs(mcmc_results, merge_chains = TRUE)

Arguments

mcmc_results

Result of calling run_mcmc()

merge_chains

boolean indicating that all chain results should be merged

Value

If merge_chains is TRUE, a named list with one numeric vector of allele frequencies per locus. Otherwise a list with one such list per chain.

Details

Returns the geometric median of the posterior distribution, defined as the point minimizing the L2 distance from each sampled point.

Examples

med_freqs <- calculate_med_allele_freqs(mcmc_results)
med_freqs[[1]]
#> [1] 0.05156748 0.08734610 0.11455374 0.08627017 0.66026251