1 package jalview.ext.archaeopteryx;
3 import jalview.analysis.TreeBuilder;
4 import jalview.datamodel.SequenceI;
5 import jalview.ext.forester.ForesterConversions;
6 import jalview.ext.forester.ForesterMatrix;
7 import jalview.util.MessageManager;
9 import java.util.HashMap;
12 import org.forester.evoinference.matrix.distance.DistanceMatrix;
13 import org.forester.phylogeny.Phylogeny;
14 import org.forester.phylogeny.PhylogenyNode;
15 import org.forester.phylogeny.data.NodeData;
16 import org.forester.phylogeny.data.Sequence;
19 * Class for converting trees made in Jalview (through TreeBuilder) to trees
20 * compatible with Forester (Phylogeny objects).
22 * Note that this currently demands a 1:1 relationship between tree nodes and
23 * the sequences used for generating them.
25 * @author kjvanderheide
28 public class ArchaeopteryxTreeBuilder
29 implements ExternalTreeBuilderI<Phylogeny, PhylogenyNode>
31 protected final SequenceI[] sequences;
33 protected final DistanceMatrix distances;
35 protected final TreeBuilder jalviewTree;
37 public String treeTitle;
39 private final Phylogeny aptxTree;
41 private PhylogenyNode rootNode;
43 private final Map<SequenceI, PhylogenyNode> alignmentBoundNodes;
45 private final Map<PhylogenyNode, SequenceI> nodesBoundAlignment;
47 public ArchaeopteryxTreeBuilder(final TreeBuilder calculatedTree)
49 jalviewTree = calculatedTree;
50 sequences = jalviewTree.getSequences();
51 distances = ForesterMatrix.convertJalviewToForester(
52 jalviewTree.getDistances(), sequences);
53 aptxTree = new Phylogeny();
54 rootNode = new PhylogenyNode();
55 alignmentBoundNodes = new HashMap<>(distances.getSize());
56 nodesBoundAlignment = new HashMap<>(distances.getSize());
62 public Phylogeny buildTree(final PhylogenyNode treeRoot)
78 public Phylogeny buildTree()
81 for (SequenceI sequence : sequences)
83 Sequence seq = ForesterConversions
84 .createForesterSequence(sequence, true);
85 PhylogenyNode sequenceNode = new PhylogenyNode(sequence.getName());
86 NodeData nodeData = sequenceNode.getNodeData();
87 nodeData.setSequence(seq);
89 nodesBoundAlignment.put(sequenceNode, sequence);
90 alignmentBoundNodes.put(sequence, sequenceNode);
91 rootNode.addAsChild(sequenceNode);
95 aptxTree.setRoot(rootNode);
97 treeTitle = generateTreeName();
98 aptxTree.setName(treeTitle);
105 public Map<SequenceI, PhylogenyNode> getAlignmentBoundNodes()
107 return alignmentBoundNodes;
111 public Map<PhylogenyNode, SequenceI> getNodesBoundAlignment()
113 return nodesBoundAlignment;
116 private Phylogeny clusterNodes()
122 * Formats a localised title for the tree panel, like
124 * Neighbour Joining Using BLOSUM62
126 * For a tree loaded from file, just uses the file name
131 public String generateTreeName() // Move this and add selection region to the
132 // title when applicable
134 if (treeTitle != null) // will currently never happen, loaded tree file will
135 // take a different path
142 * i18n description of Neighbour Joining or Average Distance method
144 String treecalcnm = MessageManager
145 .getString("label.tree_calc_" + jalviewTree.getClass()
146 .getSimpleName().substring(0, 2).toLowerCase());
148 * short score model name (long description can be too long)
150 String smn = jalviewTree.getScoreModel().getName();
153 * put them together as <method> Using <model>
155 final String ttl = MessageManager
156 .formatMessage("label.treecalc_title", treecalcnm, smn);