/*
- * Jalview - A Sequence Alignment Editor and Viewer (Version 2.4)
- * Copyright (C) 2008 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle
+ * Jalview - A Sequence Alignment Editor and Viewer (Version 2.7)
+ * Copyright (C) 2011 J Procter, AM Waterhouse, J Engelhardt, LM Lui, 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 file is part of Jalview.
*
- * 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.
+ * Jalview 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 3 of the License, or (at your option) any later version.
*
- * 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
+ * Jalview 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 Jalview. If not, see <http://www.gnu.org/licenses/>.
*/
package jalview.analysis;
* and original alignment data represented by Cigar strings.
*
* @param seqs
- * SequenceI[]
+ * SequenceI[]
* @param odata
- * Cigar[]
+ * Cigar[]
* @param treefile
- * NewickFile
+ * NewickFile
*/
public NJTree(SequenceI[] seqs, AlignmentView odata, NewickFile treefile)
{
* Creates a new NJTree object from a tree from an external source
*
* @param seqs
- * SequenceI which should be associated with leafs of treefile
+ * SequenceI which should be associated with leafs of treefile
* @param treefile
- * A parsed tree
+ * A parsed tree
*/
public NJTree(SequenceI[] seqs, NewickFile treefile)
{
* Creates a new NJTree object.
*
* @param sequence
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param type
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param pwtype
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param start
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param end
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public NJTree(SequenceI[] sequence, AlignmentView seqData, String type,
String pwtype, int start, int end)
sdata.addOperation(CigarArray.M, end - start + 1);
this.seqData = new AlignmentView(sdata, start);
}
-
+ // System.err.println("Made seqData");// dbg
if (!(type.equals("NJ")))
{
type = "AV";
{
if (ResidueProperties.getScoreMatrix(pwtype) == null)
{
- type = "BLOSUM62";
+ pwtype = "BLOSUM62";
}
}
distance = findDistances(this.seqData
.getSequenceStrings(Comparison.GapChars.charAt(0)));
-
+ // System.err.println("Made distances");// dbg
makeLeaves();
+ // System.err.println("Made leaves");// dbg
noClus = cluster.size();
cluster();
+ // System.err.println("Made clusters");// dbg
+
}
/**
- * DOCUMENT ME!
+ * Generate a string representation of the Tree
*
- * @return DOCUMENT ME!
+ * @return Newick File with all tree data available
*/
public String toString()
{
jalview.io.NewickFile fout = new jalview.io.NewickFile(getTopNode());
- return fout.print(false, true); // distances only
+ return fout.print(isHasBootstrap(),
+ isHasDistances(), isHasRootDistance()); // output all data available for tree
}
/**
* used when the alignment associated to a tree has changed.
*
* @param alignment
- * Vector
+ * Vector
*/
public void UpdatePlaceHolders(Vector alignment)
{
}
}
}
+ /**
+ * rename any nodes according to their associated sequence.
+ * This will modify the tree's metadata! (ie the original NewickFile or newly generated BinaryTree's label data)
+ */
+ public void renameAssociatedNodes() {
+ applyToNodes(new NodeTransformI() {
+ @Override
+ public void transform(BinaryNode node)
+ {
+ Object el = node.element();
+ if (el!=null && el instanceof SequenceI)
+ {
+ node.setName(((SequenceI)el).getName());
+ }
+ }
+ });
+ }
/**
* DOCUMENT ME!
*/
* DOCUMENT ME!
*
* @param i
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param j
- * DOCUMENT ME!
+ * DOCUMENT ME!
*
* @return DOCUMENT ME!
*/
* DOCUMENT ME!
*
* @param tmpi
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param tmpj
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param dist
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void findNewNJDistances(SequenceNode tmpi, SequenceNode tmpj,
float dist)
* DOCUMENT ME!
*
* @param tmpi
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param tmpj
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param dist
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void findNewDistances(SequenceNode tmpi, SequenceNode tmpj,
float dist)
* DOCUMENT ME!
*
* @param i
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param j
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void findClusterDistance(int i, int j)
{
* DOCUMENT ME!
*
* @param i
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param j
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void findClusterNJDistance(int i, int j)
{
* DOCUMENT ME!
*
* @param i
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param j
- * DOCUMENT ME!
+ * DOCUMENT ME!
*
* @return DOCUMENT ME!
*/
{
try
{
- score += pwmatrix.getPairwiseScore(sequenceString[i]
- .charAt(k), sequenceString[j].charAt(k));
+ score += pwmatrix.getPairwiseScore(
+ sequenceString[i].charAt(k),
+ sequenceString[j].charAt(k));
} catch (Exception ex)
{
System.err.println("err creating BLOSUM62 tree");
/*
* else if (pwtype.equals("SW")) { float max = -1;
*
- * for (int i = 0; i < (noseqs - 1); i++) { for (int j = i; j < noseqs; j++) {
- * AlignSeq as = new AlignSeq(sequence[i], sequence[j], "pep");
+ * for (int i = 0; i < (noseqs - 1); i++) { for (int j = i; j < noseqs; j++)
+ * { AlignSeq as = new AlignSeq(sequence[i], sequence[j], "pep");
* as.calcScoreMatrix(); as.traceAlignment(); as.printAlignment(System.out);
* distance[i][j] = (float) as.maxscore;
*
* if (max < distance[i][j]) { max = distance[i][j]; } } }
*
- * for (int i = 0; i < (noseqs - 1); i++) { for (int j = i; j < noseqs; j++) {
- * distance[i][j] = max - distance[i][j]; distance[j][i] = distance[i][j]; } } }/
+ * for (int i = 0; i < (noseqs - 1); i++) { for (int j = i; j < noseqs; j++)
+ * { distance[i][j] = max - distance[i][j]; distance[j][i] = distance[i][j];
+ * } } }/
*/
}
* Search for leaf nodes.
*
* @param node
- * root node to search from
+ * root node to search from
* @param leaves
- * Vector of leaves to add leaf node objects too.
+ * Vector of leaves to add leaf node objects too.
*
* @return Vector of leaf nodes on binary tree
*/
}
if ((node.left() == null) && (node.right() == null)) // Interior node
- // detection
+ // detection
{
leaves.addElement(node);
* Find the leaf node with a particular ycount
*
* @param node
- * initial point on tree to search from
+ * initial point on tree to search from
* @param count
- * value to search for
+ * value to search for
*
* @return null or the node with ycound=count
*/
/*
* #see findLeaf(SequenceNode node, count)
- *
*/
public Object _findLeaf(SequenceNode node, int count)
{
* printNode is mainly for debugging purposes.
*
* @param node
- * SequenceNode
+ * SequenceNode
*/
public void printNode(SequenceNode node)
{
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void findMaxDist(SequenceNode node)
{
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param threshold
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void groupNodes(SequenceNode node, float threshold)
{
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*
* @return DOCUMENT ME!
*/
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void printN(SequenceNode node)
{
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void reCount(SequenceNode node)
{
ycount = 0;
+ _lycount = 0;
+ // _lylimit = this.node.size();
_reCount(node);
}
+ private long _lycount = 0, _lylimit = 0;
+
/**
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void _reCount(SequenceNode node)
{
+ // if (_lycount<_lylimit)
+ // {
+ // System.err.println("Warning: depth of _recount greater than number of nodes.");
+ // }
if (node == null)
{
return;
}
+ _lycount++;
if ((node.left() != null) && (node.right() != null))
{
+
_reCount((SequenceNode) node.left());
_reCount((SequenceNode) node.right());
((SequenceNode) node).count = 1;
((SequenceNode) node).ycount = ycount++;
}
+ _lycount--;
}
/**
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void swapNodes(SequenceNode node)
{
* DOCUMENT ME!
*
* @param node
- * DOCUMENT ME!
+ * DOCUMENT ME!
* @param dir
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public void changeDirection(SequenceNode node, SequenceNode dir)
{
{
return hasRootDistance;
}
+
/**
* apply the given transform to all the nodes in the tree.
+ *
* @param nodeTransformI
*/
public void applyToNodes(NodeTransformI nodeTransformI)
{
- for (Enumeration nodes = node.elements(); nodes.hasMoreElements();
- nodeTransformI.transform((BinaryNode)nodes.nextElement()))
+ for (Enumeration nodes = node.elements(); nodes.hasMoreElements(); nodeTransformI
+ .transform((BinaryNode) nodes.nextElement()))
;
}
}
* Creates a new Cluster object.
*
* @param value
- * DOCUMENT ME!
+ * DOCUMENT ME!
*/
public Cluster(int[] value)
{
git://source.jalview.org / jalview.git / blobdiff