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
36
37
38
39
40
41
42
|
/*
* invsimpson.cpp
* Mothur
*
* Created by Pat Schloss on 8/20/10.
* Copyright 2010 Schloss Lab. All rights reserved.
*
*/
#include "invsimpson.h"
#include "simpson.h"
/***********************************************************************/
EstOutput InvSimpson::getValues(SAbundVector* rank){
try {
//vector<double> simpsonData(3,0);
data.resize(3,0);
vector<double> simpData(3,0);
Simpson* simp = new Simpson();
simpData = simp->getValues(rank);
if(!util.isEqual(simpData[0], 0)){
data[0] = 1/simpData[0];
data[1] = 1/simpData[2];
data[2] = 1/simpData[1];
}
else{
data.assign(3,1);
}
delete simp;
return data;
}
catch(exception& e) {
m->errorOut(e, "InvSimpson", "getValues");
exit(1);
}
}
/***********************************************************************/
|
