Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
What is the amino acid states order in a Hyphy matrix 20x20? (Read 1299 times)
Miguel
Junior Member
**
Offline


Hi Hyphy!

Posts: 53
CBMSO, CSIC (Spain)
Gender: male
What is the amino acid states order in a Hyphy matrix 20x20?
Apr 27th, 2014 at 7:49am
 
Hi There,

This is a very easy question. I am wondering what is the amino acid order in the 20x20 matrices of Hyphy (i.e., WAG.mdl), I mean to which amino acid state belongs each number?. For example, ModelMatrixName[0][0] "might" correspond to A -> A, but what is the list of the corresponding states for all numbers 0 - 19?

Many thanks!  Smiley

Best,
Miguel
Back to top
 
WWW WWW  
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: What is the amino acid states order in a Hyphy matrix 20x20?
Reply #1 - Apr 27th, 2014 at 9:00am
 
Hi Miguel,

The model itself does not assign any characters to the indices; the mapping happens at the level of a data filter -- the same 20-state model can produce arbitrary sets of 20 characters. The default amino-acid mapping is alphabetical, i.e. ACDEFGHIKLMNPQRSTVWY; this is what amino-acid filters read from alignment files adhere to.

You can use, for example,

Code:
	GetDataInfo (chars, filter, "CHARACTERS");
 



to store the 0:N-1 character mapping of the N-state "filter" in the "chars" variable.

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
 
Miguel
Junior Member
**
Offline


Hi Hyphy!

Posts: 53
CBMSO, CSIC (Spain)
Gender: male
Re: What is the amino acid states order in a Hyphy matrix 20x20?
Reply #2 - Apr 28th, 2014 at 5:44pm
 
Hi Sergei,

Thanks!
So, do you mean something like?
Code:
chars = "ACDEFGHIKLMNPQRSTVWY";
GetDataInfo (chars, filteredData, "CHARACTERS");
 



Anyway, given by default it is the alphabetic order, if I just give the values in the matrix according to such an order I guess I do not need to use "GetDataInfo", correct?

Thanks again!  Smiley
B,
Miguel
Back to top
 
WWW WWW  
IP Logged