X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=structchord.cpp;fp=structchord.cpp;h=0000000000000000000000000000000000000000;hb=4a877efa127e56e81a21f53cfdbbfd3bfbe8c4ff;hp=e85846f9d657466a2609837f5db4fe897f9dc472;hpb=a6cf29fa4dac0909c7582cb1094151d34093ee76;p=mothur.git diff --git a/structchord.cpp b/structchord.cpp deleted file mode 100644 index e85846f..0000000 --- a/structchord.cpp +++ /dev/null @@ -1,55 +0,0 @@ -/* - * structchord.cpp - * Mothur - * - * Created by westcott on 12/15/10. - * Copyright 2010 Schloss Lab. All rights reserved. - * - */ - -#include "structchord.h" - -/***********************************************************************/ -EstOutput StructChord::getValues(vector shared) { - try { - data.resize(1,0); - - double sumAj2 = 0.0; - double sumBj2 = 0.0; - - //calc the 2 denominators - for (int i = 0; i < shared[0]->getNumBins(); i++) { - - int Aij = shared[0]->getAbundance(i); - int Bij = shared[1]->getAbundance(i); - - //(Aij) ^ 2 - sumAj2 += (Aij * Aij); - sumBj2 += (Bij * Bij); - } - - sumAj2 = sqrt(sumAj2); - sumBj2 = sqrt(sumBj2); - - //calc sum - double sum = 0.0; - for (int i = 0; i < shared[0]->getNumBins(); i++) { - - int Aij = shared[0]->getAbundance(i); - int Bij = shared[1]->getAbundance(i); - - sum += (((Aij / sumAj2) - (Bij / sumBj2)) * ((Aij / sumAj2) - (Bij / sumBj2))); - } - - data[0] = sqrt(sum); - - if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; } - - return data; - } - catch(exception& e) { - m->errorOut(e, "StructChord", "getValues"); - exit(1); - } -} -/***********************************************************************/