GENEALOGY-DNA-L Archives

Archiver > GENEALOGY-DNA > 2007-07 > 1183279645


From: "Bob Bootle" <>
Subject: Re: [DNA] Help please with marker 458
Date: Sun, 1 Jul 2007 09:47:25 +0100
In-Reply-To: <a81622ac0706291814v778c80ccl9882806072865762@mail.gmail.com>


Joe,

I am a newby observer of your discussions, but have become
interested in your model run, and your table of 4984 haplotypes.

I don't understand the relationship between diversity and average
mutation rate.

Would you kindly assist my education by either, listing the mutation
rate used on each of the 37 markers, so that I may make a comparison
with the diversity of the results in your table;
or perhaps give me a sentence that will put me right on this.

Regards
Bob Bootle

....................................................................
Joe wrote :-
OK, I just did a run of a population size of 10,000 people for 10,000
years, with different mutation rates ranging from 0.002-0.008/gen for
37 single-copy markers. During the run 2,117,190 sims were born and
there were 180,345 yDNA mutations. So for each male birth (~1,059,868)
there were 180,345/1,059,868 = 0.170 mutations. For 37 markers, that's
an average rate per marker of about 0.00460. Taking the actual average
of the mutation rates modeled for the 37 markers yields 0.00461, so it
looks like the drand48() function did pretty well.

Here's a tab-separated file of the haplotypes of 4,984 living males at
the end of the run:

http://coolohio.com/dna/hapsim.tab

Hopefully you can use that as input to your model.

BTW, as another check on the data, here's a haplotype distance
histogram of the 4,894 males, similar to the one in your paper:

http://coolohio.com/dna/hapdist.jpg

Joe

-------------------------------
To unsubscribe from the list, please send an email to
with the word 'unsubscribe' without the
quotes in the subject and the body of the message


--
No virus found in this incoming message.
Checked by AVG Free Edition.
Version: 7.5.476 / Virus Database: 269.9.12/878 - Release Date: 28/06/2007
17:57


This thread: