X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=shen.cpp;fp=shen.cpp;h=a2f61889c7a88e1b43c3646470dd96224a298519;hb=cdcf99b7760701e6869ca8e4e6e91c8e8c4ae186;hp=61291a62554a3ccfb0a1b35bd761fb5bfcf1e881;hpb=edec660ee2be477ac1539109d949a60cb7e16243;p=mothur.git

diff --git a/shen.cpp b/shen.cpp
index 61291a6..a2f6188 100644
--- 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));