#include "parsimony.h"
/**************************************************************************************************/
+
EstOutput Parsimony::getValues(Tree* t) {
try {
+ globaldata = GlobalData::getInstance();
+
data.resize(1,0);
int score = 0;
-
+
+ //create pgroups that reflect the groups the user want to use
+ for(int i=t->getNumLeaves();i<t->getNumNodes();i++){
+ t->tree[i].pGroups = (t->mergeUserGroups(i));
+ }
+
for(int i=t->getNumLeaves();i<t->getNumNodes();i++){
int lc = t->tree[i].getLChild();
int rc = t->tree[i].getRChild();
- if(t->tree[i].pGroups.size() > t->tree[rc].pGroups.size() || t->tree[i].pGroups.size() > t->tree[lc].pGroups.size()){
+ int iSize = t->tree[i].pGroups.size();
+ int rcSize = t->tree[rc].pGroups.size();
+ int lcSize = t->tree[lc].pGroups.size();
+
+ //if isize are 0 then that branch is to be ignored
+ if (iSize == 0) { }
+ else if ((rcSize == 0) || (lcSize == 0)) { }
+ //if you have more groups than either of your kids then theres been a change.
+ else if(iSize > rcSize || iSize > lcSize){
score++;
}
}
+
data[0] = score;
return data;
}
}
+
/**************************************************************************************************/
+