This contains Pat's bug fixes and updates. It represents mothur v.1.3.0

[mothur.git] / shen.cpp

diff --git a/shen.cpp b/shen.cpp
index 61291a62554a3ccfb0a1b35bd761fb5bfcf1e881..a2f61889c7a88e1b43c3646470dd96224a298519 100644 (file)
--- a/shen.cpp
+++ b/shen.cpp
@@ -28,7 +28,7 @@ EstOutput Shen::getValues(SAbundVector* rank){
                Ace* calc = new Ace(abund);
                EstOutput ace = calc->getValues(rank);
                
-               double f0 = n - ace[0];
+               double f0 = ace[0]-rank->getNumBins();
                
                data[0] = f0 * (1 - pow(1 - f1/n/f0, m));