X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Futil%2FComparison.java;h=67a668a9c66b29d3e4ca6d36a530b8d4900fd901;hb=cbd8cafe5f1d302cf490d34a96f8231659afd402;hp=e6aae969b63c92b25741cc61b79d40459f956f13;hpb=74aa9aa23dce90e649abce71577e56be2c9ac00e;p=jalview.git

diff --git a/src/jalview/util/Comparison.java b/src/jalview/util/Comparison.java
index e6aae96..67a668a 100755
--- a/src/jalview/util/Comparison.java
+++ b/src/jalview/util/Comparison.java
@@ -1,109 +1,226 @@
+/*
+ * Jalview - A Sequence Alignment Editor and Viewer
+ * Copyright (C) 2007 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle
+ *
+ * This program is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU General Public License
+ * as published by the Free Software Foundation; either version 2
+ * of the License, or (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA
+ */
 package jalview.util;
 
 import jalview.datamodel.*;
 
-public class Comparison {
-
-  public static float compare(SequenceI ii, SequenceI jj)
+/**
+ * DOCUMENT ME!
+ *
+ * @author $author$
+ * @version $Revision$
+ */
+public class Comparison
+{
+  /** DOCUMENT ME!! */
+  public static final String GapChars = " .-";
+
+  /**
+   * DOCUMENT ME!
+   *
+   * @param ii DOCUMENT ME!
+   * @param jj DOCUMENT ME!
+   *
+   * @return DOCUMENT ME!
+   */
+  public static final float compare(SequenceI ii, SequenceI jj)
   {
-    return Comparison.compare(ii,jj,0,ii.getLength()-1);
+    return Comparison.compare(ii, jj, 0, ii.getLength() - 1);
   }
-  public static float compare(SequenceI ii, SequenceI jj, int start, int end) {
 
-     String si   = ii.getSequence();
-     String sj   = jj.getSequence();
+  /**
+   * this was supposed to be an ungapped pid calculation
+   * @param ii SequenceI
+   * @param jj SequenceI
+   * @param start int
+   * @param end int
+   * @return float
+   */
+  public static float compare(SequenceI ii, SequenceI jj, int start, int end)
+  {
+    String si = ii.getSequenceAsString();
+    String sj = jj.getSequenceAsString();
 
-     int ilen = end-start;
-     int jlen = end-start;
+    int ilen = si.length() - 1;
+    int jlen = sj.length() - 1;
 
-     if ( jalview.util.Comparison.isGap( si.charAt(start+ilen)) )
+    while (jalview.util.Comparison.isGap(si.charAt(start + ilen)))
+    {
+      ilen--;
+    }
+
+    while (jalview.util.Comparison.isGap(sj.charAt(start + jlen)))
+    {
+      jlen--;
+    }
+
+    int count = 0;
+    int match = 0;
+    float pid = -1;
+
+    if (ilen > jlen)
+    {
+      for (int j = 0; j < jlen; j++)
       {
-       ilen--;
-
-       while (jalview.util.Comparison.isGap( si.charAt(start+ilen)))
-       {
-         ilen--;
-       }
-     }
-
-     if ( jalview.util.Comparison.isGap( sj.charAt(start+jlen)) )
-     {
-       jlen--;
-
-       while (jalview.util.Comparison.isGap( sj.charAt(start+jlen)))
-       {
-         jlen--;
-       }
-     }
-
-     int   count = 0;
-     int   match = 0;
-     float pid   = -1;
-
-     if (ilen > jlen) {
-
-       for (int j = 0; j < jlen; j++) {
-         if (si.substring(start + j,start + j+1).equals(sj.substring(start + j,start + j+1))) {
-           match++;
-         }
-         count++;
-       }
-       pid = (float)match/(float)ilen * 100;
-     } else {
-       for (int j = 0; j < jlen; j++) {
-         if (si.substring(start + j,start + j+1).equals(sj.substring(start + j,start + j+1))) {
-           match++;
-         }
-         count++;
-       }
-       pid = (float)match/(float)jlen * 100;
-     }
+        if (si.substring(start + j, start + j + 1).equals(sj.substring(start +
+            j, start + j + 1)))
+        {
+          match++;
+        }
+
+        count++;
+      }
+
+      pid = (float) match / (float) ilen * 100;
+    }
+    else
+    {
+      for (int j = 0; j < jlen; j++)
+      {
+        if (si.substring(start + j, start + j + 1).equals(sj.substring(start +
+            j, start + j + 1)))
+        {
+          match++;
+        }
+
+        count++;
+      }
+
+      pid = (float) match / (float) jlen * 100;
+    }
 
     return pid;
   }
 
-  /**    */
-  public static float PID(Sequence s1 , Sequence s2)
+  /**
+   * this is a gapped PID calculation
+   *
+   * @param s1 SequenceI
+   * @param s2 SequenceI
+   * @return float
+   */
+  public final static float PID(String seq1, String seq2)
   {
-    int len;
+    return PID(seq1, seq2, 0, seq1.length());
+  }
 
-    if (s1.getSequence().length() > s2.getSequence().length())
-      len = s1.getSequence().length();
-    else
-      len = s2.getSequence().length();
+  static final int caseShift = 'a' - 'A';
 
+  // Another pid with region specification
+  public final static float PID(String seq1, String seq2, int start, int end)
+  {
 
-    int bad = 0;
+    int s1len = seq1.length();
+    int s2len = seq2.length();
 
-    for (int i = 0; i < len; i++)
+    int len = Math.min(s1len, s2len);
+
+    if (end < len)
     {
-      char chr1;
-      char chr2;
+      len = end;
+    }
 
-      if (i < s1.getSequence().length())
-        chr1 = s1.getSequence().charAt(i);
-      else
-        chr1 = '.';
+    if (len < start)
+    {
+      start = len - 1; // we just use a single residue for the difference
+    }
 
+    int bad = 0;
+    char chr1;
+    char chr2;
 
-      if (i < s2.getSequence().length())
-        chr2 = s2.getSequence().charAt(i);
-      else
-        chr2 = '.';
+    for (int i = start; i < len; i++)
+    {
+      chr1 = seq1.charAt(i);
 
+      chr2 = seq2.charAt(i);
 
-      if (!(jalview.util.Comparison.isGap( chr1 ))  &&  !(jalview.util.Comparison.isGap( chr2 )))
+      if ('a' <= chr1 && chr1 <= 'z')
       {
-        if (chr1!=chr2)
-          bad++;
+        // TO UPPERCASE !!!
+        //Faster than toUpperCase
+        chr1 -= caseShift;
+      }
+      if ('a' <= chr2 && chr2 <= 'z')
+      {
+        // TO UPPERCASE !!!
+        //Faster than toUpperCase
+        chr2 -= caseShift;
+      }
+
+      if (chr1 != chr2 && !isGap(chr1) && !isGap(chr2))
+      {
+        bad++;
       }
     }
 
-    return (float)100*(len-bad)/len;
+    return ( (float) 100 * (len - bad)) / len;
   }
 
-  public static boolean isGap(char c)
+  /**
+   * DOCUMENT ME!
+   *
+   * @param c DOCUMENT ME!
+   *
+   * @return DOCUMENT ME!
+   */
+  public static final boolean isGap(char c)
   {
-    return  (c != '.' && c != '-' && c != ' ') ? false : true;
+    return (c == '-' || c == '.' || c == ' ') ? true : false;
+  }
+
+  public static final boolean isNucleotide(SequenceI[] seqs)
+  {
+    int i = 0, iSize = seqs.length, j, jSize;
+    float nt = 0, aa = 0;
+    char c;
+    while (i < iSize)
+    {
+      jSize = seqs[i].getLength();
+      for (j = 0; j < jSize; j++)
+      {
+        c = seqs[i].getCharAt(j);
+        if ('a' <= c && c <= 'z')
+        {
+          c -= ('a' - 'A');
+        }
+
+        if (c == 'A' || c == 'G' || c == 'C' || c == 'T' || c == 'U')
+        {
+          nt++;
+        }
+        else if (!jalview.util.Comparison.isGap(seqs[i].getCharAt(j)))
+        {
+          aa++;
+        }
+      }
+      i++;
+    }
+
+    if ( (nt / (nt + aa)) > 0.85f)
+    {
+      return true;
+    }
+    else
+    {
+      return false;
+    }
+
   }
 }