public String getOutputRecord(int tagIndex, double[][] theResults) {
   int GNum = 0;
   for (int i = 0; i < theResults.length; i++) {
     if (theResults[i] != null) GNum++;
   }
   StringBuilder sb = new StringBuilder();
   sb.append(BaseEncoder.getSequenceFromLong(tbt.getTag(tagIndex))).append("\t");
   sb.append(String.valueOf(tbt.getNumberOfTaxaWithTag(tagIndex))).append("\t");
   sb.append(String.valueOf(GNum)).append("\t");
   // chr, bestSite, refAlignment.getPositionInLocus(bestSite), bestP, countSig, (double)-1,
   // bestNewTestFamilies.length, (double)this.getFamilyCode(bestNewTestFamilies)}
   for (int i = 0; i < theResults.length; i++) {
     if (theResults[i] == null) continue;
     sb.append(String.valueOf((int) theResults[i][0])).append("\t");
     sb.append(String.valueOf((int) theResults[i][1])).append("\t");
     sb.append(String.valueOf((int) theResults[i][2])).append("\t");
     sb.append(String.valueOf(theResults[i][3])).append("\t");
     sb.append(String.valueOf((int) theResults[i][4])).append("\t");
     sb.append(String.valueOf(theResults[i][5])).append("\t");
     sb.append(String.valueOf((int) theResults[i][6])).append("\t");
     sb.append(String.valueOf((long) theResults[i][7])).append("\t");
   }
   sb.append("\n");
   return sb.toString();
 }
コード例 #2
0
 public void run() {
   long starttime = System.currentTimeMillis();
   ArrayList<String> resultList = new ArrayList();
   for (int i = 0; i < tbt.getTagCount(); i++) {
     if (tbt.getNumberOfTaxaWithTag(i) < 30) continue;
     double[] bestR = {-1, -1, -1, 1, -1};
     int blastChr = -1, blastPos = -1, bestAlignment = -1, refDiv = -1;
     long[] testTag = tbt.getTag(i);
     long[] testTagDist;
     double[][] theResults = new double[theAnchorChrInBits.length][];
     ScanChromosomeMTR[] scanOnChr = new ScanChromosomeMTR[theAnchorChrInBits.length];
     for (int cn = 0; cn < theAnchorChrInBits.length; cn++) {
       testTagDist =
           getTagsInBits(
               tbt, i, tbt2anchRedirect[cn], theAnchorChrInBits[cn].getSequenceCount());
       scanOnChr[cn] =
           new ScanChromosomeMTR(
               theAnchorChrInBits[cn], testTagDist, cn, theResults, 0.000001, 1, blastPos);
       scanOnChr[cn].scan();
     }
     int countRealSig = 0;
     double[] pRank = new double[theAnchorChrInBits.length];
     for (int cn = 0; cn < theAnchorChrInBits.length; cn++) {
       double[] r = theResults[cn];
       pRank[cn] = r[3];
       if (r[3] < bestR[3]) {
         bestR = r.clone();
         bestAlignment = cn;
       }
       if (r[3] < 0.000001) countRealSig++;
     }
     Arrays.sort(pRank);
     double[][] bestResWithNewThreshold = new double[1][];
     testTagDist =
         getTagsInBits(
             tbt,
             i,
             tbt2anchRedirect[bestAlignment],
             theAnchorChrInBits[bestAlignment].getSequenceCount());
     ScanChromosomeMTR bestChrNewThres =
         new ScanChromosomeMTR(
             theAnchorChrInBits[bestAlignment],
             testTagDist,
             0,
             bestResWithNewThreshold,
             pRank[1],
             1,
             blastPos);
     bestChrNewThres.scan();
     int countOfSitesBetterThanNextBestChr = (int) bestResWithNewThreshold[0][4];
     String s =
         String.format(
             "%s %d %d %d %d %d %d %d %g %d %d %d %g %g %d %d %d %n",
             BaseEncoder.getSequenceFromLong(testTag),
             i,
             blastChr,
             blastPos,
             refDiv,
             (int) bestR[0],
             (int) bestR[1],
             (int) bestR[2],
             bestR[3],
             (int) bestR[4],
             tbt.getNumberOfTaxaWithTag(i),
             countRealSig,
             Math.log10(pRank[1] / pRank[0]),
             Math.log10(pRank[theAnchorChrInBits.length / 2] / pRank[0]),
             countOfSitesBetterThanNextBestChr,
             bestChrNewThres.minSigPos,
             bestChrNewThres.maxSigPos);
     resultList.add(s);
   }
   result = resultList.toArray(new String[resultList.size()]);
   double secs = (double) (System.currentTimeMillis() - starttime) * 0.001;
   System.out.println(
       "Task "
           + taskID
           + " with "
           + tbt.getTagCount()
           + " tags is finished in "
           + secs
           + " seconds");
 }