ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/gclib/cuffcompare/cuffcompare.cpp
(Generate patch)
# Line 1451 | Line 1451
1451      fsn=(100.0*(double)ps->locusATP)/ps->total_rloci; //(ps->locusATP+ps->locusAFN);
1452      fprintf(fout, "       Locus level: \t%5.1f\t%5.1f\t%5.1f\t%5.1f\n",sn, sp, fsn, fsp);
1453      //fprintf(fout, "                   (locus TP=%d, total ref loci=%d)\n",ps->locusTP, ps->total_rloci);
1454 +    fprintf(fout, "\nMatching intron chains: %7d\n",ps->ichainTP);
1455 +    fprintf(fout, "         Matching loci: %7d\n",ps->locusTP);
1456 +    fprintf(fout, "\n");
1457      sn=(100.0*(double)ps->m_exons)/(ps->total_rexons);
1458 <    fprintf(fout, "   Missed exons:\t%d/%d (%5.1f%%)\n",ps->m_exons, ps->total_rexons, sn);
1458 >    fprintf(fout, "          Missed exons: %7d/%d\t(%5.1f%%)\n",ps->m_exons, ps->total_rexons, sn);
1459      sn=(100.0*(double)ps->w_exons)/(ps->total_qexons);
1460 <    fprintf(fout, "    Wrong exons:\t%d/%d (%5.1f%%)\n",ps->w_exons, ps->total_qexons,sn);
1460 >    fprintf(fout, "           Novel exons: %7d/%d\t(%5.1f%%)\n",ps->w_exons, ps->total_qexons,sn);
1461      if (ps->total_rintrons>0) {
1462 <      sn=(100.0*(double)ps->m_introns)/(ps->total_rintrons);
1463 <      fprintf(fout, " Missed introns:\t%d/%d (%5.1f%%)\n",ps->m_introns, ps->total_rintrons, sn);
1464 <      }
1462 >    sn=(100.0*(double)ps->m_introns)/(ps->total_rintrons);
1463 >    fprintf(fout, "        Missed introns: %7d/%d\t(%5.1f%%)\n",ps->m_introns, ps->total_rintrons, sn);
1464 >    }
1465      if (ps->total_qintrons>0) {
1466 <      sn=(100.0*(double)ps->w_introns)/(ps->total_qintrons);
1467 <      fprintf(fout, "  Wrong introns:\t%d/%d (%5.1f%%)\n",ps->w_introns, ps->total_qintrons,sn);
1468 <      }
1466 >    sn=(100.0*(double)ps->w_introns)/(ps->total_qintrons);
1467 >    fprintf(fout, "         Novel introns: %7d/%d\t(%5.1f%%)\n",ps->w_introns, ps->total_qintrons,sn);
1468 >    }
1469      if (ps->total_rloci>0) {
1470 <      sn=(100.0*(double)ps->m_loci)/(ps->total_rloci);
1471 <      fprintf(fout, " Missed loci:\t%d/%d (%5.1f%%)\n",ps->m_loci, ps->total_rloci, sn);
1472 <      }
1470 >    sn=(100.0*(double)ps->m_loci)/(ps->total_rloci);
1471 >    fprintf(fout, "           Missed loci: %7d/%d\t(%5.1f%%)\n",ps->m_loci, ps->total_rloci, sn);
1472 >    }
1473      if (ps->total_qloci>0) {
1474 <      sn=(100.0*(double)ps->w_loci)/(ps->total_qloci);
1475 <      fprintf(fout, "  Wrong loci:\t%d/%d (%5.1f%%)\n",ps->w_loci, ps->total_qloci,sn);
1476 <      }
1474 >    sn=(100.0*(double)ps->w_loci)/(ps->total_qloci);
1475 >    fprintf(fout, "            Novel loci: %7d/%d\t(%5.1f%%)\n",ps->w_loci, ps->total_qloci,sn);
1476 >    }
1477  
1478    }
1479   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines