X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;ds=sidebyside;f=sharedanderbergs.cpp;fp=sharedanderbergs.cpp;h=0000000000000000000000000000000000000000;hb=4a877efa127e56e81a21f53cfdbbfd3bfbe8c4ff;hp=cbb9d3003ed84631e5251f9baf6738a165c83c7e;hpb=a6cf29fa4dac0909c7582cb1094151d34093ee76;p=mothur.git diff --git a/sharedanderbergs.cpp b/sharedanderbergs.cpp deleted file mode 100644 index cbb9d30..0000000 --- a/sharedanderbergs.cpp +++ /dev/null @@ -1,48 +0,0 @@ -/* - * sharedanderberg.cpp - * Mothur - * - * Created by Sarah Westcott on 3/23/09. - * Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved. - * - */ - -#include "sharedanderbergs.h" - -/***********************************************************************/ - -EstOutput Anderberg::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] = 1.0 - S12 / ((float)((2 * S1) + (2 * S2) - (3 * S12))); - - if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; } - - return data; - } - catch(exception& e) { - m->errorOut(e, "Anderberg", "getValues"); - exit(1); - } -} - -/***********************************************************************/