Dr. Lulz
YaBB Newbies
Offline
Curious HyPhy user
Posts: 9

If I have a category matrix with 4 classes (0,1,2,3) and I use ConstructCategoryMatrix, it will return the value for each position.
I am interested in cases where I have 2 categories each with 4 classes and the sitespecific class refinement for each category. How do I find this information from the number like "12" when the original categories are 03 classes each?
I assume that it is category1*4 + category2, but how do I distinguish between category1 and category2? is category1 always the one that has the "alpha" variables in the category statement? as the reference line? and not the "beta"?
Ie:
category category1 = (4, EQUAL, MEAN, GammaDist(_x_,alpha,alpha), CGammaDist(_x_,alpha,alpha), 0, 1e25, CGammaDist(_x_, alpha+1, alpha));
category category2 = (4, EQUAL, MEAN, GammaDist(_x_,alpha,beta), CGammaDist(_x_,alpha,beta), 0, 1e25, CGammaDist(_x_, alpha+1, beta)*alpha/beta);
