5 * Created by Thomas Ryabin on 5/18/09.
6 * Copyright 2009Schloss Lab UMASS Amherst. All rights reserved.
14 /***********************************************************************/
15 EstOutput Shen::getValues(SAbundVector* rank){
21 double n = (double)rank->getNumSeqs();
22 double f1 = (double)rank->get(1);
24 Ace* calc = new Ace(abund);
25 EstOutput ace = calc->getValues(rank);
27 double f0 = ace[0]-rank->getNumBins();
29 data[0] = f0 * (1 - pow(1 - f1/n/f0, m));
36 errorOut(e, "Shen", "getValues");
42 /***********************************************************************/