Function to compute a theoretical z-score given source split times and admixture proportions.

Source: R/theory.R

This function computes the theoretical z-score of an f4-statistic of the form f4(PO,P2,PX,P1) as a function of the admixture time a, source (P1 & P2) split time s, and admixture proportion alpha. It assumes a single lineage is sampled from each population.

theoretical_zscore(t, a, s, alpha)

Arguments

t

Twigstats cutoff time in units of 2Ne generations

a

Admixture time in units of 2Ne generations

s

Source split time in units of 2Ne generations

alpha

Admixture proportions. Proportion of P2 in PX.

Value

Returns the theoretical z-score value.