string noOutliers, outliers;
if ((!filter) && (seqMask == "")) {
- noOutliers = templateFileName + "pintail.quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.quan";
}else if ((!filter) && (seqMask != "")) {
- noOutliers =templateFileName + "pintail.masked.quan";
+ noOutliers =m->getRootName(m->getSimpleName(templateFileName)) + "pintail.masked.quan";
}else if ((filter) && (seqMask != "")) {
- noOutliers = templateFileName + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "masked.quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "masked.quan";
}else if ((filter) && (seqMask == "")) {
- noOutliers = templateFileName + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "quan";
+ noOutliers = m->getRootName(m->getSimpleName(templateFileName)) + "pintail.filtered." + m->getSimpleName(m->getRootName(fastafile)) + "quan";
}
decalc->removeObviousOutliers(quantilesMembers, templateSeqs.size());
}
}
//***************************************************************************************************************
-int Pintail::print(ostream& out, ostream& outAcc) {
+Sequence Pintail::print(ostream& out, ostream& outAcc) {
try {
+
int index = ceil(deviation);
//is your DE value higher than the 95%
for (int m = 0; m < expectedDistance.size(); m++) { out << expectedDistance[m] << '\t'; }
out << endl;
- return 0;
+ return *querySeq;
}
catch(exception& e) {
}
#ifdef USE_MPI
//***************************************************************************************************************
-int Pintail::print(MPI_File& out, MPI_File& outAcc) {
+Sequence Pintail::print(MPI_File& out, MPI_File& outAcc) {
try {
- bool results = false;
+
string outputString = "";
int index = ceil(deviation);
MPI_File_write_shared(outAcc, buf, length, MPI_CHAR, &statusAcc);
delete buf;
- results = true;
+ return *querySeq;
}
outputString += "Observed\t";
MPI_File_write_shared(out, buf2, length, MPI_CHAR, &status);
delete buf2;
- return results;
+ return *querySeq;
}
catch(exception& e) {
m->errorOut(e, "Pintail", "print");
}
in.close();
- remove(s.c_str());
+ m->mothurRemove(s);
}
#else