Summarize allele frequencies

Usage

summarize_allele_freqs(
  mcmc_results,
  lower_quantile = 0.025,
  upper_quantile = 0.975,
  merge_chains = TRUE
)

Arguments

mcmc_results: Result of calling run_mcmc()
lower_quantile: The lower quantile of the posterior distribution to return
upper_quantile: The upper quantile of the posterior distribution to return
merge_chains: boolean indicating that all chain results should be merged

Value

Data frame with one row per allele (per chain if merge_chains is FALSE) giving posterior quantiles and mean of its frequency, along with locus and allele identifiers.

Details

Summarize individual allele frequencies from the posterior distribution of sampled allele frequencies

Examples

allele_freq_summary <- summarize_allele_freqs(mcmc_results)
head(allele_freq_summary)
#>      post_allele_freqs_lower post_allele_freqs_med post_allele_freqs_upper
#> L1.1             0.019395342            0.03279833              0.05071323
#> L1.2             0.058118432            0.08353058              0.11538900
#> L1.3             0.092535655            0.13331380              0.16875155
#> L1.4             0.061060573            0.08807241              0.12006423
#> L1.5             0.608060258            0.66022103              0.70981749
#> L2.1             0.005734225            0.01441477              0.02644067
#>      post_allele_freqs_mean locus allele
#> L1.1             0.03324264    L1   L1_1
#> L1.2             0.08465452    L1   L1_2
#> L1.3             0.13340537    L1   L1_3
#> L1.4             0.08897809    L1   L1_4
#> L1.5             0.65971939    L1   L1_5
#> L2.1             0.01469600    L2   L2_1