New constructors for turning NewickFile objects into renderable trees
authorjprocter <Jim Procter>
Fri, 18 Mar 2005 16:53:46 +0000 (16:53 +0000)
committerjprocter <Jim Procter>
Fri, 18 Mar 2005 16:53:46 +0000 (16:53 +0000)
mapped onto the current set of sequences by matching leafnames to
sequence names (where one string must contain the other, rather than
direct equivalence).

src/jalview/analysis/NJTree.java

index 7324b6c..fb0ec33 100755 (executable)
@@ -5,6 +5,8 @@ import jalview.util.*;
 import jalview.schemes.ResidueProperties;\r
 import java.util.*;\r
 \r
+import jalview.io.NewickFile;\r
+\r
 public class NJTree {\r
 \r
   Vector cluster;\r
@@ -44,6 +46,73 @@ public class NJTree {
   public NJTree(SequenceNode node) {\r
     top = node;\r
     maxheight = findHeight(top);\r
+\r
+  }\r
+  // Private SequenceID class to do fuzzy .equals() method for Hashtable.\r
+\r
+  private class SeqIdname {\r
+    String id;\r
+\r
+    SeqIdname(String s) {\r
+      id = new String(s);\r
+    }\r
+    public int hashCode() {\r
+      return (id.substring(0,4).hashCode());\r
+    }\r
+    public boolean equals(Object s) {\r
+      if (s instanceof SeqIdname) {\r
+        return this.equals((SeqIdname) s);\r
+      } else {\r
+        if (s instanceof String) {\r
+          return this.equals((String) s);\r
+        }\r
+      }\r
+      return false;\r
+    }\r
+\r
+\r
+    public boolean equals(SeqIdname s) {\r
+      if (id.startsWith(s.id) || s.id.startsWith(id)) {\r
+        return true;\r
+      }\r
+      return false;\r
+    }\r
+\r
+    public boolean equals(String s) {\r
+      if (id.startsWith(s) || s.startsWith(id)) {\r
+        return true;\r
+      }\r
+      return false;\r
+    }\r
+  }\r
+\r
+  public NJTree(SequenceI[] seqs, NewickFile treefile) {\r
+    top = treefile.getTree();\r
+    maxheight = findHeight(top);\r
+    Hashtable names = new Hashtable();\r
+    for (int i = 0; i < seqs.length; i++)\r
+    {\r
+      names.put(new SeqIdname(seqs[i].getDisplayId()), seqs[i]);\r
+    }\r
+    Vector leaves = new Vector();\r
+    findLeaves(top, leaves);\r
+    int i = 0;\r
+    int namesleft = seqs.length;\r
+    SequenceNode j;\r
+    SeqIdname nam;\r
+    while (i < leaves.size())\r
+    {\r
+      j = (SequenceNode) leaves.elementAt(i++);\r
+      nam = new SeqIdname(j.getName());\r
+      if ((namesleft>-1)\r
+          && names.containsKey(nam))\r
+      {\r
+        j.setElement(names.get(nam));\r
+        namesleft--;\r
+      } else {\r
+        j.setElement(new Sequence(nam.id, "THISISAPLACEHLDER"));\r
+      }\r
+    }\r
   }\r
 \r
   public NJTree(SequenceI[] sequence,int start, int end) {\r