+ mothurOut("Calculating observed distance... "); cout.flush();
+ for (int i = lines[0]->start; i < lines[0]->end; i++) {
+ vector<float> obsi = decalc->calcObserved(querySeqs[i], bestfit[i], windows[i], windowSizes[i]);
+ obsDistance[i] = obsi;
+ }
+ mothurOut("Done."); mothurOutEndLine();
+
+
+
+ mothurOut("Finding variability... "); cout.flush();
+ for (int i = lines[0]->start; i < lines[0]->end; i++) {
+ vector<float> q = decalc->findQav(windows[i], windowSizes[i], probabilityProfile, h[i]);
+ Qav[i] = q;
+ }
+ mothurOut("Done."); mothurOutEndLine();
+
+
+
+ mothurOut("Calculating alpha... "); cout.flush();
+ for (int i = lines[0]->start; i < lines[0]->end; i++) {
+ float alpha = decalc->getCoef(obsDistance[i], Qav[i]);
+ seqCoef.push_back(alpha);
+ }
+ mothurOut("Done."); mothurOutEndLine();
+
+
+
+ mothurOut("Calculating expected distance... "); cout.flush();
+ for (int i = lines[0]->start; i < lines[0]->end; i++) {
+ vector<float> exp = decalc->calcExpected(Qav[i], seqCoef[i]);
+ expectedDistance[i] = exp;
+ }
+ mothurOut("Done."); mothurOutEndLine();
+
+
+
+ mothurOut("Finding deviation... "); cout.flush();
+ for (int i = lines[0]->start; i < lines[0]->end; i++) {
+ float de = decalc->calcDE(obsDistance[i], expectedDistance[i]);
+ DE[i] = de;
+
+ it = trimmed[i].begin();
+ float dist = decalc->calcDist(querySeqs[i], bestfit[i], it->first, it->second);
+ deviation[i] = dist;
+ }
+ mothurOut("Done."); mothurOutEndLine();