\
Section author: Catherine Lozupone
An example of how to calculate the pairwise distances for a set of sequences.
NOTE: UPGMA should not be used for phylogenetic reconstruction.
>>> from cogent import LoadSeqs
>>> from cogent.phylo import distance
>>> from cogent.cluster.UPGMA import upgma
Import a substitution model (or create your own)
>>> from cogent.evolve.models import HKY85
Load the alignment.
>>> al = LoadSeqs("data/test.paml")
Create a pairwise distances object calculator for the alignment, providing a substitution model instance.
>>> d = distance.EstimateDistances(al, submodel= HKY85())
>>> d.run()
Now use this matrix to build a UPGMA cluster.
>>> mycluster = upgma(d.getPairwiseDistances())
>>> print mycluster.asciiArt()
/-NineBande
/edge.1--|
| | /-HowlerMon
/edge.0--| \edge.2--|
| | \-Human
-root----| |
| \-DogFaced
|
\-Mouse
We demonstrate saving this UPGMA cluster to a file.
>>> mycluster.writeToFile('test_upgma.tree')