X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=sharedochiai.cpp;fp=sharedochiai.cpp;h=0000000000000000000000000000000000000000;hp=b49fa4a3cc5730b0a97b0eefe30f952a7671f932;hb=4a877efa127e56e81a21f53cfdbbfd3bfbe8c4ff;hpb=a6cf29fa4dac0909c7582cb1094151d34093ee76 diff --git a/sharedochiai.cpp b/sharedochiai.cpp deleted file mode 100644 index b49fa4a..0000000 --- a/sharedochiai.cpp +++ /dev/null @@ -1,48 +0,0 @@ -/* - * sharedochiai.cpp - * Mothur - * - * Created by Sarah Westcott on 3/23/09. - * Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved. - * - */ - -#include "sharedochiai.h" - -/***********************************************************************/ - -EstOutput Ochiai::getValues(vector shared) { - try { - double S1, S2, S12, tempA, tempB; - S1 = 0; S2 = 0; S12 = 0; tempA = 0; tempB = 0; - - /*S1, S2 = number of OTUs observed or estimated in A and B - S12=number of OTUs shared between A and B */ - - data.resize(1,0); - - for (int i = 0; i < shared[0]->getNumBins(); i++) { - //store in temps to avoid multiple repetitive function calls - tempA = shared[0]->getAbundance(i); - tempB = shared[1]->getAbundance(i); - - if (tempA != 0) { S1++; } - if (tempB != 0) { S2++; } - - //they are shared - if ((tempA != 0) && (tempB != 0)) { S12++; } - } - - data[0] = S12 / ((float)pow((S1 * S2), 0.5)); - - if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; } - - return data; - } - catch(exception& e) { - m->errorOut(e, "Ochiai", "getValues"); - exit(1); - } -} - -/***********************************************************************/