X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=unweighted.cpp;h=3694d263fb1039ea5208926301b7bd18a38875a3;hb=503b23b959230453d95dc812fbbe088b3be6910c;hp=405dfd2b3feaf4c1a5327526936a7baaaa4aba0a;hpb=8fd4a602269974e6548f4a387dddd6bd80a999ab;p=mothur.git diff --git a/unweighted.cpp b/unweighted.cpp index 405dfd2..3694d26 100644 --- a/unweighted.cpp +++ b/unweighted.cpp @@ -64,11 +64,11 @@ EstOutput Unweighted::getValues(Tree* t) { //if i's children are from the same group then i's pcount size will be 1 //if copyIpcount.size() = 0 they are from a branch that is entirely from a group the user doesn't want if (copyIpcount.size() == 0) { } - else if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += t->tree[i].getBranchLength(); } + else if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += abs(t->tree[i].getBranchLength()); } //add i's BL to total if it is from the groups the user wants if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() != 0)) { - totalBL += t->tree[i].getBranchLength(); + totalBL += abs(t->tree[i].getBranchLength()); } /**********************************************************************/ @@ -78,16 +78,16 @@ EstOutput Unweighted::getValues(Tree* t) { if (t->tree[rc].getRChild() == -1) { //if rc is a valid group and rc has a BL if ((inUsersGroups(t->tree[rc].getGroup(), groups) == true) && (t->tree[rc].getBranchLength() != -1)) { - UniqueBL += t->tree[rc].getBranchLength(); - totalBL += t->tree[rc].getBranchLength(); + UniqueBL += abs(t->tree[rc].getBranchLength()); + totalBL += abs(t->tree[rc].getBranchLength()); } } if (t->tree[lc].getLChild() == -1) { //if lc is a valid group and lc has a BL if ((inUsersGroups(t->tree[lc].getGroup(), groups) == true) && (t->tree[lc].getBranchLength() != -1)) { - UniqueBL += t->tree[lc].getBranchLength(); - totalBL += t->tree[lc].getBranchLength(); + UniqueBL += abs(t->tree[lc].getBranchLength()); + totalBL += abs(t->tree[lc].getBranchLength()); } } @@ -138,11 +138,11 @@ EstOutput Unweighted::getValues(Tree* t) { //if i's children are from the same group then i's pcount size will be 1 //if copyIpcount.size() = 0 they are from a branch that is entirely from a group the user doesn't want if (copyIpcount.size() == 0) { } - else if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += t->tree[i].getBranchLength(); } + else if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += abs(t->tree[i].getBranchLength()); } //add i's BL to total if it is from the groups the user wants if ((t->tree[i].getBranchLength() != -1) && (copyIpcount.size() != 0)) { - totalBL += t->tree[i].getBranchLength(); + totalBL += abs(t->tree[i].getBranchLength()); } /**********************************************************************/ @@ -152,16 +152,16 @@ EstOutput Unweighted::getValues(Tree* t) { if (t->tree[rc].getRChild() == -1) { //if rc is a valid group and rc has a BL if ((inUsersGroups(t->tree[rc].getGroup(), groups) == true) && (t->tree[rc].getBranchLength() != -1)) { - UniqueBL += t->tree[rc].getBranchLength(); - totalBL += t->tree[rc].getBranchLength(); + UniqueBL += abs(t->tree[rc].getBranchLength()); + totalBL += abs(t->tree[rc].getBranchLength()); } } if (t->tree[lc].getLChild() == -1) { //if lc is a valid group and lc has a BL if ((inUsersGroups(t->tree[lc].getGroup(), groups) == true) && (t->tree[lc].getBranchLength() != -1)) { - UniqueBL += t->tree[lc].getBranchLength(); - totalBL += t->tree[lc].getBranchLength(); + UniqueBL += abs(t->tree[lc].getBranchLength()); + totalBL += abs(t->tree[lc].getBranchLength()); } } @@ -252,11 +252,11 @@ EstOutput Unweighted::getValues(Tree* t, string groupA, string groupB) { //if i's children are from the same group then i's pcount size will be 1 //if copyIpcount.size() = 0 they are from a branch that is entirely from a group the user doesn't want if (copyIpcount.size() == 0) { } - else if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += copyTree->tree[i].getBranchLength(); } + else if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += abs(copyTree->tree[i].getBranchLength()); } //add i's BL to total if it is from the groups the user wants if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() != 0)) { - totalBL += copyTree->tree[i].getBranchLength(); + totalBL += abs(copyTree->tree[i].getBranchLength()); } /**********************************************************************/ @@ -266,16 +266,16 @@ EstOutput Unweighted::getValues(Tree* t, string groupA, string groupB) { if (copyTree->tree[rc].getRChild() == -1) { //if rc is a valid group and rc has a BL if ((inUsersGroups(copyTree->tree[rc].getGroup(), groups) == true) && (copyTree->tree[rc].getBranchLength() != -1)) { - UniqueBL += copyTree->tree[rc].getBranchLength(); - totalBL += copyTree->tree[rc].getBranchLength(); + UniqueBL += abs(copyTree->tree[rc].getBranchLength()); + totalBL += abs(copyTree->tree[rc].getBranchLength()); } } if (copyTree->tree[lc].getLChild() == -1) { //if lc is a valid group and lc has a BL if ((inUsersGroups(copyTree->tree[lc].getGroup(), groups) == true) && (copyTree->tree[lc].getBranchLength() != -1)) { - UniqueBL += copyTree->tree[lc].getBranchLength(); - totalBL += copyTree->tree[lc].getBranchLength(); + UniqueBL += abs(copyTree->tree[lc].getBranchLength()); + totalBL += abs(copyTree->tree[lc].getBranchLength()); } } @@ -332,11 +332,11 @@ EstOutput Unweighted::getValues(Tree* t, string groupA, string groupB) { //if i's children are from the same group then i's pcount size will be 1 //if copyIpcount.size() = 0 they are from a branch that is entirely from a group the user doesn't want if (copyIpcount.size() == 0) { } - else if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { UniqueBL += copyTree->tree[i].getBranchLength(); } + else if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() == 1)) { abs(UniqueBL += copyTree->tree[i].getBranchLength()); } //add i's BL to total if it is from the groups the user wants if ((copyTree->tree[i].getBranchLength() != -1) && (copyIpcount.size() != 0)) { - totalBL += copyTree->tree[i].getBranchLength(); + totalBL += abs(copyTree->tree[i].getBranchLength()); } /**********************************************************************/ @@ -346,16 +346,16 @@ EstOutput Unweighted::getValues(Tree* t, string groupA, string groupB) { if (copyTree->tree[rc].getRChild() == -1) { //if rc is a valid group and rc has a BL if ((inUsersGroups(copyTree->tree[rc].getGroup(), groups) == true) && (copyTree->tree[rc].getBranchLength() != -1)) { - UniqueBL += copyTree->tree[rc].getBranchLength(); - totalBL += copyTree->tree[rc].getBranchLength(); + UniqueBL += abs(copyTree->tree[rc].getBranchLength()); + totalBL += abs(copyTree->tree[rc].getBranchLength()); } } if (copyTree->tree[lc].getLChild() == -1) { //if lc is a valid group and lc has a BL if ((inUsersGroups(copyTree->tree[lc].getGroup(), groups) == true) && (copyTree->tree[lc].getBranchLength() != -1)) { - UniqueBL += copyTree->tree[lc].getBranchLength(); - totalBL += copyTree->tree[lc].getBranchLength(); + UniqueBL += abs(copyTree->tree[lc].getBranchLength()); + totalBL += abs(copyTree->tree[lc].getBranchLength()); } }