turpin@cs.utexas.edu (Russell Turpin) (01/25/91)
----- Cross-posted to sci.bio. ----- In article <1991Jan24.220906.16423@usenet.ins.cwru.edu> rpetsche@mrg.CWRU.EDU (rolfe g petschek) writes: > Hmmm. As a statistical physicist I can not help but notice that the > period of time required required for a single lineage to overtake all > other lineages is a random variable with calculable distribution, > dependent on the model. ... Absolutely. You were not the first to notice this. I was not the first to notice this. Ah, well, sometimes we get beat to the punch. > ... Well, I tried lots of models (admittedly crudely) and I did > not get sense from this statement except when (a) teh size of the > breeding population is *very* small (need N^2 generations assuming > neither advantage or growth of size N) or (b) their is significant > advantage passed *only* through the female line ... For biological reasons, forget (b). But for the most obvious model reflecting biological constraints, the number of expected generations is (surprise!) *linear* in N. This discussion came up some months ago in another newsgroup. Not being a statistical physicists, I wrote a model in C rather than attempting an analytic solution. When I posted the results, someone else informed me that the results were well-known, and are known as the Wright-Fisher model. (He did not send me a reference.) Shown below are the number of females in a fixed-size population, the range in the number of generations my simulation required for one female line to push out all others, and the expected number that is analytically calculated. #generations (N) Wright-Fisher F least most average ---- ----- ---- ------- 8 5 37 16 128 172 625 256 1024 2392 4163 2048 (For F=1024, I only ran my model twice, because it took so long on my MacPlus. The other values of F were run several times.) The assumptions in the model are (1) a fixed-size female population, and (2) every woman has a fixed probability of not having daughters. (The number of men and number of sons produced just plain doesn't matter.) In my simulation, I assumed that each "excess" daughter had equal chance of coming from any of the women who did have daughters in the previous generation. A slightly different assumption is made in the Wright-Fisher model. > ... Did you have some model in mind which would answer these > confusions? 10,000 generations is not long (its square root is > 100, and recently, at least, there have been more than 100 > persons in the breeding population.) 2000 generations of early man is only 30,000 years, and that is how long the Wright-Fischer model predicts it took for a single female line to emerge in a population of 1024 women. Playing around with my simulation, I was left with the intuition (which is what a good model provides) that allowing the population to increase and decrease, so that it was only sporadically as low as 1024 women, would actually *increase* the chance of a single female line emerging, especially if one takes into account that small groups of related people, are likely to survive or die together. Russell