1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
/*
* whittaker.cpp
* Mothur
*
* Created by Pat Schloss on 4/23/09.
* Copyright 2009 Patrick D. Schloss. All rights reserved.
*
*/
#include "whittaker.h"
/***********************************************************************/
EstOutput Whittaker::getValues(vector<SharedRAbundVector*> shared){
try{
data.resize(1);
int countA = 0;
int countB = 0;
int sTotal = shared[0]->getNumBins();
for(int i=0;i<sTotal;i++){
if(shared[0]->getAbundance(i) != 0){ countA++; }
if(shared[1]->getAbundance(i) != 0){ countB++; }
}
data[0] = 2-2*sTotal/(float)(countA+countB);
return data;
}
catch(exception& e) {
m->errorOut(e, "Whittaker", "getValues");
exit(1);
}
}
/***********************************************************************/
|
