Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Mean substitution rates (Read 2192 times)
sspielman
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 2
Mean substitution rates
Apr 13th, 2010 at 3:15pm
 
Hello,
I am trying to find the mean substitution rate for particular regions in a gene using HKY85 model, but thus far have only been able to find per-site values for substitution rates using siterates.bf, and inputing my data, defining a partition, providing HyPhy with a tree, etc. returns only a tv/ts value.
Any suggestions would be much appreciated!
-Stephanie
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Mean substitution rates
Reply #1 - Apr 13th, 2010 at 5:10pm
 
Hi Stephanie,

What do mean by an average rate? HyPhy will estimate branch lengths for you (measured as expected subs/unit time), but and if you want the actual rate, you need to estimate the evolutionary time by other means (fossil record, molecular clock etc). Could you elaborate on the problem you are trying to solve?

Sergei

PS You may find section 2.2. of Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login useful as background reading.
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged
 
sspielman
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 2
Re: Mean substitution rates
Reply #2 - Apr 14th, 2010 at 12:04pm
 
Essentially, I have a gene with several regions which are under different functional constraints, so sites within a region should have similar substitution rates. For instance, if sites 1-50 are in region 1 and 51-75 are in region 2, I am looking to find the mean of the per-site substitution rates within region 1 and again for region 2.
Is there a way to find the confidence interval or standard deviation for the individual site rates returned by SiteRates.bf that would allow me to calculate the mean with error values?
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Mean substitution rates
Reply #3 - Apr 15th, 2010 at 9:26pm
 
Hi Stephanie,

Do you have time information, i.e. are you talking about actual rates (e.g. subs/site/MY), or simply genetic distances (subs/site/unit time)?

In case it is the latter, then all you need is to perform a relative-ratio type analysis for your partitions. Averaging site-by-site estimates is not a good idea (mostly because they will have vastly differing variances), I would recommend simply obtaining a mean (per partition) estimate of substitution rates. Take a look at Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login

This will allow you to estimate the ratio between substitution rates in your multiple regions.

To convert this to actual rates, you need more information (e.g. divergence time estimates)

HTH,
Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged