X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Futil%2FComparison.java;h=e48c0cac61ade1262c60ca13e062d3615977d9be;hb=9807ea8cc4e867378c52d51b832168d4bc611d29;hp=8041d8074a4364df0ec23245f167e6bf3366924d;hpb=588042b69abf8e60bcc950b24c283933c7dd422f;p=jalview.git

diff --git a/src/jalview/util/Comparison.java b/src/jalview/util/Comparison.java
index 8041d80..e48c0ca 100755
--- a/src/jalview/util/Comparison.java
+++ b/src/jalview/util/Comparison.java
@@ -21,160 +21,254 @@ package jalview.util;
 import jalview.datamodel.*;
 
 
-public class Comparison {
-    public static String GapChars = " .-";
-
-    public static float compare(SequenceI ii, SequenceI jj) {
-        return Comparison.compare(ii, jj, 0, ii.getLength() - 1);
+/**
+ * DOCUMENT ME!
+ *
+ * @author $author$
+ * @version $Revision$
+ */
+public class Comparison
+{
+  /** DOCUMENT ME!! */
+  public static String GapChars = " .-";
+
+  /**
+   * DOCUMENT ME!
+   *
+   * @param ii DOCUMENT ME!
+   * @param jj DOCUMENT ME!
+   *
+   * @return DOCUMENT ME!
+   */
+  public static float compare(SequenceI ii, SequenceI jj)
+  {
+    return Comparison.compare(ii, jj, 0, ii.getLength() - 1);
+  }
+
+  /**
+   * 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.getSequence();
+    String sj = jj.getSequence();
+
+    int ilen = si.length() - 1;
+    int jlen = sj.length() - 1;
+
+    while (jalview.util.Comparison.isGap(si.charAt(start + ilen)))
+    {
+      ilen--;
     }
 
-    /**
-     * 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.getSequence();
-        String sj = jj.getSequence();
-
-        int ilen = si.length() - 1;
-        int jlen = sj.length() - 1;
-
-        while (jalview.util.Comparison.isGap(si.charAt(start + ilen))) {
-            ilen--;
-        }
+    while (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++;
         }
 
-        int count = 0;
-        int match = 0;
-        float pid = -1;
+        count++;
+      }
 
-        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++;
-                }
+      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++;
-            }
+        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++;
-                }
+      pid = (float) match / (float) jlen * 100;
+    }
 
-                count++;
-            }
+    return pid;
+  }
+
+  /**
+   * this is a gapped PID calculation
+   *
+   * @param s1 SequenceI
+   * @param s2 SequenceI
+   * @return float
+   */
+  public static float PID(SequenceI s1, SequenceI s2)
+  {
+    int len;
+
+    if (s1.getSequence().length() > s2.getSequence().length())
+    {
+      len = s1.getSequence().length();
+    }
+    else
+    {
+      len = s2.getSequence().length();
+    }
 
-            pid = (float) match / (float) jlen * 100;
+    int bad = 0;
+
+    for (int i = 0; i < len; i++)
+    {
+      char chr1;
+      char chr2;
+
+      if (i < s1.getSequence().length())
+      {
+        chr1 = Character.toUpperCase(s1.getSequence().charAt(i));
+      }
+      else
+      {
+        chr1 = '.';
+      }
+
+      if (i < s2.getSequence().length())
+      {
+        chr2 = Character.toUpperCase(s2.getSequence().charAt(i));
+      }
+      else
+      {
+        chr2 = '.';
+      }
+
+      if (! (jalview.util.Comparison.isGap(chr1)) &&
+          ! (jalview.util.Comparison.isGap(chr2)))
+      {
+        if (chr1 != chr2)
+        {
+          bad++;
         }
-
-        return pid;
+      }
     }
 
-    /**
-     * this is a gapped PID calculation
-     *
-     * @param s1 SequenceI
-     * @param s2 SequenceI
-     * @return float
-     */
-    public static float PID(SequenceI s1, SequenceI s2) {
-        int len;
-
-        if (s1.getSequence().length() > s2.getSequence().length()) {
-            len = s1.getSequence().length();
-        } else {
-            len = s2.getSequence().length();
-        }
+    return ( (float) 100 * (len - bad)) / len;
+  }
 
-        int bad = 0;
-
-        for (int i = 0; i < len; i++) {
-            char chr1;
-            char chr2;
-
-            if (i < s1.getSequence().length()) {
-                chr1 = s1.getSequence().charAt(i);
-            } else {
-                chr1 = '.';
-            }
-
-            if (i < s2.getSequence().length()) {
-                chr2 = s2.getSequence().charAt(i);
-            } else {
-                chr2 = '.';
-            }
-
-            if (!(jalview.util.Comparison.isGap(chr1)) &&
-                    !(jalview.util.Comparison.isGap(chr2))) {
-                if (chr1 != chr2) {
-                    bad++;
-                }
-            }
-        }
+  // Another pid with region specification
+  public static float PID(SequenceI s1, SequenceI s2, int start, int end)
+  {
+    int len;
 
-        return ((float) 100 * (len - bad)) / len;
+    if (s1.getSequence().length() > s2.getSequence().length())
+    {
+      len = s1.getSequence().length();
+    }
+    else
+    {
+      len = s2.getSequence().length();
     }
 
-    // Another pid with region specification
-    public static float PID(SequenceI s1, SequenceI s2, int start, int end) {
-        int len;
-
-        if (s1.getSequence().length() > s2.getSequence().length()) {
-            len = s1.getSequence().length();
-        } else {
-            len = s2.getSequence().length();
-        }
+    if (end < len)
+    {
+      len = end;
+    }
 
-        if (end < len) {
-            len = end;
-        }
+    if (len < start)
+    {
+      start = len - 1; // we just use a single residue for the difference
+    }
 
-        if (len < start) {
-            start = len - 1; // we just use a single residue for the difference
+    int bad = 0;
+
+    for (int i = start; i < len; i++)
+    {
+      char chr1;
+      char chr2;
+
+      if (i < s1.getSequence().length())
+      {
+        chr1 = Character.toUpperCase(s1.getSequence().charAt(i));
+      }
+      else
+      {
+        chr1 = '.';
+      }
+
+      if (i < s2.getSequence().length())
+      {
+        chr2 = Character.toUpperCase(s2.getSequence().charAt(i));
+      }
+      else
+      {
+        chr2 = '.';
+      }
+
+      if (! (jalview.util.Comparison.isGap(chr1)) &&
+          ! (jalview.util.Comparison.isGap(chr2)))
+      {
+        if (chr1 != chr2)
+        {
+          bad++;
         }
+      }
+    }
 
-        int bad = 0;
-
-        for (int i = start; i < len; i++) {
-            char chr1;
-            char chr2;
-
-            if (i < s1.getSequence().length()) {
-                chr1 = s1.getSequence().charAt(i);
-            } else {
-                chr1 = '.';
-            }
-
-            if (i < s2.getSequence().length()) {
-                chr2 = s2.getSequence().charAt(i);
-            } else {
-                chr2 = '.';
-            }
-
-            if (!(jalview.util.Comparison.isGap(chr1)) &&
-                    !(jalview.util.Comparison.isGap(chr2))) {
-                if (chr1 != chr2) {
-                    bad++;
-                }
-            }
+    return ( (float) 100 * (len - bad)) / len;
+  }
+
+  /**
+   * DOCUMENT ME!
+   *
+   * @param c DOCUMENT ME!
+   *
+   * @return DOCUMENT ME!
+   */
+  public static boolean isGap(char c)
+  {
+    return (c == '-' || c == '.' || c == ' ') ? true : false;
+  }
+
+  public static 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++;
         }
-
-        return ((float) 100 * (len - bad)) / len;
+      }
+      i++;
     }
 
-    public static boolean isGap(char c) {
-        return ((c != '.') && (c != '-') && (c != ' ')) ? false : true;
-    }
+    if ( (nt / (nt + aa)) > 0.85f)
+      return true;
+    else
+      return false;
+
+  }
 }