working on pam

[mothur.git] / homovacommand.cpp

diff --git a/homovacommand.cpp b/homovacommand.cpp
index 49b8b64a0ebea8d5c6dae20118a19e0b814806b1..488b68e0b6b69a8eb8e351d39be9ff09193fb36f 100644 (file)
--- a/homovacommand.cpp
+++ b/homovacommand.cpp
@@ -308,7 +308,7 @@ double HomovaCommand::runHOMOVA(ofstream& HOMOVAFile, map<string, vector<int> >
                        vector<double> ssWithinRandVector;
                        map<string, vector<int> > randomizedGroup = getRandomizedGroups(groupSampleMap);
                        double bValueRand = calcBValue(randomizedGroup, ssWithinRandVector);
-                       if(bValueRand > bValueOrig){    counter++;      }
+                       if(bValueRand >= bValueOrig){   counter++;      }
                }
                
                double pValue = (double) counter / (double) iters;