Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
different rate matrix for every codon (Read 3347 times)
YuanlingRui
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 3
different rate matrix for every codon
Sep 3rd, 2013 at 1:44am
 
Dear Hyphy users,
   I am doing some work about coding sequences evolution that giving different rate matrix (Q matrix) for every codon. I want to know, 1) whether do I have to use the function LikelihoodFunction id = (datafilterid_0,tree_0,...,<datafilterid_n>,<tree_n>,), in which the n is codon number; 2) if I have to use the function, how about the time-consuming.
Thank you very much!

Regards!
Yuanling Xia
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: different rate matrix for every codon
Reply #1 - Sep 5th, 2013 at 1:59pm
 
Hi Yuanling,

In the most general case, you do need to use the form you describe. You probably also want to constrain branch lengths between trees on each codon, otherwise you'll have a separate set of branch lengths (e.g. the t parameter) for every codon. Do you have prototype code? I can look at it and help you ensure that everything works smoothly and correctly.

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
 
YuanlingRui
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 3
Re: different rate matrix for every codon
Reply #2 - Sep 9th, 2013 at 1:39am
 
Hi Sergei,
    Thank you for answer, when I write out the prototype code, I will ask for you help again.

Regards
Yuanling
Back to top
 
 
IP Logged
 
YuanlingRui
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 3
Re: different rate matrix for every codon
Reply #3 - Oct 26th, 2013 at 7:44pm
 
Hi Sergei,

    I have send the prototype code to you. The time-consuming will increase greatly when I give many partitions using function LikelihoodFunction. Thank you very much!

Regards
Yuanling
Back to top
 
 
IP Logged