Question: How To Calculate Genome Clock Rate From Snp'S Clock Rate
gravatar for polangxin
7.9 years ago by
United States
polangxin80 wrote:

a lot of sample mapping to 1 genome(as reference), call snps, then join them as a Multiple Sequenc Aligenment file, use BEAST calculate the clock.rate is 9e-4.

the ref is about 40000bp, and SNPs is 6000.

i think 9e-4 is the snps' clock rate, am i right? then how i calcute the specie's clock rate?


ADD COMMENTlink modified 7.9 years ago by Lee Katz3.1k • written 7.9 years ago by polangxin80
gravatar for Lee Katz
7.9 years ago by
Lee Katz3.1k
Atlanta, GA
Lee Katz3.1k wrote:

You are close. I am assuming that you got the clock rate from Tracer or something comparable. Next you'll have to consider the search space. If you start with 6kb out of 40kb, then I believe that you will need to multiply your clock rate by 6/40. 9E-4*6/40=0.000135 snps/genome/unit time.

Please anyone, correct me if I am wrong. I haven't had to do this yet but it makes sense to me, conceptually.

ADD COMMENTlink written 7.9 years ago by Lee Katz3.1k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1485 users visited in the last hour