working on slayer bug

[mothur.git] / slayer.cpp

diff --git a/slayer.cpp b/slayer.cpp
index ba1bf337d4084e58416b06c48cb6a359bb5b05b1..8575faf4312450c78cc8db249259b036edb91e0d 100644 (file)
--- a/slayer.cpp
+++ b/slayer.cpp
@@ -297,8 +297,8 @@ int Slayer::bootstrapSNPS(vector<snps> left, vector<snps> right, float& BSA, flo
                int count_A = 0; // sceneario QLA,QRB supported
                int count_B = 0; // sceneario QLB,QRA supported
        
-               int numLeft = max(1, int(left.size() * (percentSNPSample/(float)100) + 0.5));
-               int numRight = max(1, int(right.size() * (percentSNPSample/(float)100) + 0.5));
+               int numLeft = max(1, int(left.size() * percentSNPSample/(float)100 + 0.5));
+               int numRight = max(1, int(right.size() * percentSNPSample/(float)100 + 0.5));
 
                for (int i = 0; i < iters; i++) {
                        //random sampling with replacement.
@@ -365,8 +365,8 @@ int Slayer::bootstrapSNPS(vector<snps> left, vector<snps> right, float& BSA, flo
 
 
 
-               BSA = ((float) count_A / (float) iters) * 100;
-               BSB = ((float) count_B / (float) iters) * 100;
+               BSA = (float) count_A / (float) iters * 100;
+               BSB = (float) count_B / (float) iters * 100;
 //cout << "bsa = " << BSA << " bsb = " << BSB << endl;
 
                return 0;