Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Ancestral sequence reconstruction (Read 7342 times)
sergios-orestis
YaBB Newbies
*
Offline


osx/centos

Posts: 10
New York City
Gender: male
Ancestral sequence reconstruction
Mar 11th, 2010 at 12:15pm
 
Hi Sergei,

I'd like to estimate ancestral sequences given a topology, an alignment (AA in my case) and a model (WAG, JTT, REV, maybe LG). Any chance this routine can be isolated from SLAC etc? I thought a bf was there, but I don't see it. The point is to compare with the PAML (aaml) results.

Many thanks,
sergios
Back to top
 
WWW WWW  
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Ancestral sequence reconstruction
Reply #1 - Mar 11th, 2010 at 10:13pm
 
Hi Sergios,

There are three different ways to reconstruct ancestors in HyPhy: joint likelihood (default), marginal likelihood and sampled likelihood. If you are happy with the joint likelihood approach, then the BF you are looking for is a post-processor 'Ancestral Sequences': it will become available after you run AnalyzeNucProtData.bf, or any other BF that builds a likelihood function. In the GUI version you can access it via Analysis->Results->Ancestral Sequences; in the command line version use the -p command line argument. The BF itself is in TemplateBatchFiles/post_ancestors.bf.

Cheers,
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
 
sergios-orestis
YaBB Newbies
*
Offline


osx/centos

Posts: 10
New York City
Gender: male
Re: Ancestral sequence reconstruction
Reply #2 - Mar 12th, 2010 at 1:13pm
 
Thanks, Sergei. It worked fine. How could I try the other methods, too, and what is the joint lik reference, please?

Best,
sergios
Back to top
 
WWW WWW  
IP Logged