if(qvalues[i] < threshold){
m->mothurOut("Feature " + toString((i+1)) + " is significant, q = ");
cout << qvalues[i];
- m->mothurOutJustToLog(toString(pvalues[i])); m->mothurOutEndLine();
+ m->mothurOutJustToLog(toString(qvalues[i])); m->mothurOutEndLine();
}
}
vector<double> MothurMetastats::calc_qvalues(vector<double>& pValues) {
try {
+ /* cout << "x <- c(" << pValues[0];
+ for (int l = 1; l < pValues.size(); l++){
+ cout << ", " << pValues[l];
+ }
+ cout << ")\n";*/
+
int numRows = pValues.size();
vector<double> qvalues(numRows, 0.0);