🚀 go-pugleaf

Well I came up with a novel way to test heterogeneity of gene frequencies
between samples of haploid individuals from several populations (sounds
like Karl Pearson, so what? :-)). The question is how to compare other
methods (like X^2, Fisher's exact, conditional X^2 etc) with it. I
simulated genetic drift by the coalescence process, then sampled from
diverged populations, repeated that process many times to get empirical
powers for the different degrees of divergence.. Therefore, under the
neutral scenario I have an idea how things behave.

Now, how people implement selection? Motoo Kimura's treatment of random
intensity selection when there is no dominance, gives U-shaped allele
frequency distribution. Is it theoretically valid to generate steady-state
distribution of allele frequencies for each population from a Beta(a<1,
b<1) density to approximate Kimura's results?  What shape should allele
frequency profile follow for a non-equilibrium situation? I could not find
if there is any theory done for that case. I'd appreciate any references,
and hints.


Ivan Zimogorov