🚀 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 do 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