X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sharedace.cpp;fp=sharedace.cpp;h=c8ba4fa2347735dc27640321c286384960ddb134;hb=e4c80376cc4533f66c8dfc18f3e1a86a60ac17fe;hp=24fe5ef00e369a7a221a9bcabc4e416492a29e14;hpb=8f7164a88df624fd0a8d1eddbb3d744463cc9ecb;p=mothur.git diff --git a/sharedace.cpp b/sharedace.cpp index 24fe5ef..c8ba4fa 100644 --- a/sharedace.cpp +++ b/sharedace.cpp @@ -18,10 +18,10 @@ EstOutput SharedAce::getValues(vector shared) { string label; label = shared[0]->getLabel(); - int fARare, fBRare, S12Rare, S12Abund, S12, f11, tempA, tempB, t10, t01, t11, t21, t12, t22, C12Numerator; + double fARare, fBRare, S12Rare, S12Abund, S12, f11, tempA, tempB, t10, t01, t11, t21, t12, t22, C12Numerator; fARare = 0; fBRare = 0; S12Rare = 0; S12Abund = 0; S12 = 0; f11 = 0; t10 = 0; t01 = 0; t11= 0; t21= 0; t12= 0; t22= 0; C12Numerator = 0; - float Sharedace, C12, part1, part2, part3, part4, part5, Gamma1, Gamma2, Gamma3; + double Sharedace, C12, part1, part2, part3, part4, part5, Gamma1, Gamma2, Gamma3; /*fARare = number of OTUs with one individual found in A and less than or equal to 10 in B. fBRare = number of OTUs with one individual found in B and less than or equal to 10 in A. @@ -89,7 +89,7 @@ EstOutput SharedAce::getValues(vector shared) { if (isnan(Gamma3) || isinf(Gamma3)) { Gamma3 = 0; } if (isnan(part1) || isinf(part1)) { part1 = 0; } if (isnan(part2) || isinf(part2)) { part2 = 0; } - + part3 = fARare * Gamma1; part4 = fBRare * Gamma2; part5 = f11 * Gamma3;