5 * Created by Thomas Ryabin on 5/18/09.
6 * Copyright 2009 __MyCompanyName__. All rights reserved.
14 /***********************************************************************/
15 EstOutput Shen::getValues(SAbundVector* rank){
20 double n = (double)rank->getNumSeqs();
21 double f1 = (double)rank->get(1);
23 double D_rare = 0; //I didn't know what this was. Simply replace the '0' with the appropriate expression.
24 double C_rare = 1; //I didn't know what this was. Simply replace the '1' with the appropriate expression.
25 double Y_rare = 1; //I didn't know what this was. Simply replace the '1' with the appropriate expression.
27 double f0 = D_rare/C_rare + f1/C_rare * Y_rare*Y_rare - D_rare;
29 data[0] = f0 * (1 - pow(1 - f1/n/f0, m));
34 cout << "Standard Error: " << e.what() << " has occurred in the Coverage class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
38 cout << "An unknown error has occurred in the Coverage class function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
44 /***********************************************************************/