HyPhy message board - Implementing a categorical frequency matrix

	Welcome, Guest. Please Login

Home

Help

HyPhy message board › HYPHY Package › Feature Additions › Implementing a categorical frequency matrix

(Moderators: Sergei, Simon)

‹ Previous Topic | Next Topic ›

Pages: 1

Send Topic

Implementing a categorical frequency matrix (Read 4838 times)

LeavittJM YaBB Newbies Offline Feed your monkey! Posts: 3	Implementing a categorical frequency matrix Mar 31^st, 2011 at 11:02am I am implementing a categorical frequency matrix in HYPHY using the ConstructCategoryMatrix command. However I am having some trouble interpreting the resulting matrix. Do the values for a particular codon still need to be normalized? Thanks!
Back to top	IP Logged

Sergei YaBB Administrator Offline Datamonkeys are forever... Posts: 1658 UCSD Gender:	Re: Implementing a categorical frequency matrix Reply #1 - Mar 31^st, 2011 at 11:13am Hi there, ConstructCategoryMatrix takes a lot of options which control what it generates; one of the options will output normalized posterior probabilities. I'll post an example batch file later today. Sergei
Back to top	Associate Professor Division of Infectious Diseases Division of Biomedical Informatics School of Medicine University of California San Diego WWW IP Logged

LeavittJM

YaBB Newbies

Offline

Feed your monkey!

Posts: 3

Re: Implementing a categorical frequency matrix
Reply #2 - Mar 31^st, 2011 at 12:41pm

Thank you for the quick reply.

I will check out your batch file when you put it up. However, I have one quick question to make sure I understand the output from the ConstructCategoryMatrix command.

For a CategoryMatrix of {2,1}, if it gives you {0.1} {0.3}, then there is a distribution of 25% and 75% between the two categories.

Thanks again.

IP Logged

Sergei

YaBB Administrator

Offline

Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male

Re: Implementing a categorical frequency matrix
Reply #3 - Apr 1^st, 2011 at 12:38am

Hi there,

By default, the function returns conditional likelihoods: Pr (site | rate class). In your example, the probability of observing a site given the model in class 1 is 0.1, and given the model in class 2 is 0.3

Most of the time one wants Pr (rate class | site), which can be obtained by a simple application of the Bayes rule:

Code:

Pr (rate class | site) = Pr (site|rate class) * Pr (rate class) / Pr (site)

Pr (site) is simply the sum of all Pr (site|rate class) * Pr (rate class) (the normalization factor), and Pr (rate class) can be obtained by calling

Code:

ConstructCategoryMatrix (priors, lfID, WEIGHTS);

HTH,
Sergei

PS I could probably give you more pointers if I had a better idea of what the final objective was.

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

WWW

IP Logged

LeavittJM

YaBB Newbies

Offline

Feed your monkey!

Posts: 3

Re: Implementing a categorical frequency matrix
Reply #4 - Apr 6^th, 2011 at 2:41pm

Hi Sergei,

The following code functions correctly, though I had to use for loops instead of matrix functions.

Is there a way to do it faster?

Thanks for the help,
John

Code:

ConstructCategoryMatrix(catmat,theLnLik,COMPLETE);
ConstructCategoryMatrix (priors, theLnLik, WEIGHTS);

result = {Rows(catmat),Columns(catmat)};

for( i=0; i<Rows(catmat); i=i+1)
{
  for(j=0; j<Columns(catmat); j=j+1)
  {
    result[i][j] = catmat[i][j]*priors[i][0];
  }
}

for( j=0; j<Columns(catmat); j=j+1)
{
  sum = 0;
  for( i=0; i<Rows(catmat); i=i+1)
  {
    sum = sum+result[i][j];
  }
  for( i=0; i<Rows(catmat); i=i+1)
  {
    result[i][j] = result[i][j]/sum;
  }
}
fprintf(stdout,"\nresult:",result,"\n");

IP Logged

Sergei YaBB Administrator Offline Datamonkeys are forever... Posts: 1658 UCSD Gender:	Re: Implementing a categorical frequency matrix Reply #5 - Apr 13^th, 2011 at 1:49pm Hi John, There really is not a built-in way to do it much faster (in terms of lines of code). HyPhy template batch files contain some HBL functions to do the same, but your code is perfectly adequate. You can define it as a function and include it in other files for later use. Sergei
Back to top	Associate Professor Division of Infectious Diseases Division of Biomedical Informatics School of Medicine University of California San Diego WWW IP Logged

Pages: 1

Send Topic

‹ Previous Topic | Next Topic ›

« Home

‹ Board

Top of this page