+ m->errorOut(e, "Ccode", "getT");
+ exit(1);
+ }
+}
+/**************************************************************************************************/
+float Ccode::getF(int numseq) {
+ try {
+
+ float fvalue = 0;
+
+ /* F-Snedecor critical values for v1=1 and different degrees of freedom v2 and alpha 0.05 */
+ if (numseq > 120) fvalue = 3.84;
+ else if (numseq > 60) fvalue = 3.92;
+ else if (numseq > 40) fvalue = 4.00;
+ else if (numseq > 30) fvalue = 4.08;
+ else if (numseq > 29) fvalue = 4.17;
+ else if (numseq > 28) fvalue = 4.18;
+ else if (numseq > 27) fvalue = 4.20;
+ else if (numseq > 26) fvalue = 4.21;
+ else if (numseq > 25) fvalue = 4.23;
+ else if (numseq > 24) fvalue = 4.24;
+ else if (numseq > 23) fvalue = 4.26;
+ else if (numseq > 22) fvalue = 4.28;
+ else if (numseq > 21) fvalue = 4.30;
+ else if (numseq > 20) fvalue = 4.32;
+ else if (numseq > 19) fvalue = 4.35;
+ else if (numseq > 18) fvalue = 4.38;
+ else if (numseq > 17) fvalue = 4.41;
+ else if (numseq > 16) fvalue = 4.45;
+ else if (numseq > 15) fvalue = 4.49;
+ else if (numseq > 14) fvalue = 4.54;
+ else if (numseq > 13) fvalue = 4.60;
+ else if (numseq > 12) fvalue = 4.67;
+ else if (numseq > 11) fvalue = 4.75;
+ else if (numseq > 10) fvalue = 4.84;
+ else if (numseq > 9) fvalue = 4.96;
+ else if (numseq > 8) fvalue = 5.12;
+ else if (numseq > 7) fvalue = 5.32;
+ else if (numseq > 6) fvalue = 5.59;
+ else if (numseq > 5) fvalue = 5.99;
+ else if (numseq > 4) fvalue = 6.61;
+ else if (numseq > 3) fvalue = 7.71;
+ else if (numseq > 2) fvalue = 10.1;
+ else if (numseq > 1) fvalue = 18.5;
+ else if (numseq > 0) fvalue = 161;
+ else if (numseq <= 0) {
+ m->mothurOut("Two or more reference sequences are required, your data will be flawed.\n"); m->mothurOutEndLine();
+ }
+
+ return fvalue;
+ }
+ catch(exception& e) {
+ m->errorOut(e, "Ccode", "getF");