modified some calculators and other stuff - pds

[mothur.git] / bergerparker.cpp

diff --git a/bergerparker.cpp b/bergerparker.cpp
index 893f2a95dad68e6f9ba1638eb29be43071a711aa..c9f3565f985c52f1e744f4e320207243a6bba11d 100644 (file)
--- a/bergerparker.cpp
+++ b/bergerparker.cpp
@@ -15,8 +15,7 @@ EstOutput BergerParker::getValues(SAbundVector* rank){
        try {
                data.resize(1,0);
                //Berger-Parker index
-               double BP = (double)rank->getNumSeqs()/(double)rank->getMaxRank();
-               //cout << "BP index = " << 1/BP << "\n\n";
+               double BP = (double)rank->getMaxRank() / (double)rank->getNumSeqs();
                
                data[0] = BP;
                if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; }