Skip to Main Content Skip to Search
Product Documentation

seqlinkage - Construct phylogenetic tree from pairwise distances

Syntax

PhyloTree = seqlinkage(Distances)
PhyloTree = seqlinkage(Distances, Method)
PhyloTree = seqlinkage(Distances, Method, Names)

Arguments

Distances

Matrix or vector of pairwise distances, such as returned by the seqpdist function.

Method

String that specifies a distance method. Choices are:

  • 'single'

  • 'complete'

  • 'average' (default)

  • 'weighted'

  • 'centroid'

  • 'median'

Names

Specifies alternative labels for leaf nodes. Choices are:

  • Vector of structures, each with a Header or Name field

  • Cell array of strings

The elements must be unique. The number of elements must comply with the number of samples used to generate the pairwise distances in Dist.

Description

PhyloTree = seqlinkage(Distances) returns a phylogenetic tree object from the pairwise distances, Distances, between the species or products. Distances is a matrix or vector of pairwise distances, such as returned by the seqpdist function.

PhyloTree = seqlinkage(Distances, Method) creates a phylogenetic tree object using a specified patristic distance method. The available methods are:

'single'

Nearest distance (single linkage method)

'complete'

Furthest distance (complete linkage method)

'average' (default)

Unweighted Pair Group Method Average (UPGMA, group average).

'weighted'

Weighted Pair Group Method Average (WPGMA)

'centroid'

Unweighted Pair Group Method Centroid (UPGMC)

'median'

Weighted Pair Group Method Centroid (WPGMC)

PhyloTree = seqlinkage(Distances, Method, Names) passes a list of unique names to label the leaf nodes (for example, species or products) in a phylogenetic tree object.

Examples

Build Phylogenetic Tree from Pairwise Distances

Build a phylogenetic tree from pairwise distances, specifying both a distance-computing method and leaf names.

Create an array of structures representing a multiple alignment of amino acids:

seqs = fastaread('pf00002.fa');

Measure the Jukes-Cantor pairwise distances between sequences:

distances = seqpdist(seqs,'method','jukes-cantor','indels','pair');
You will use the output argument distances,
a vector containing biological distances between each pair of sequences,
as an input argument to seqlinkage.
Build the phylogenetic tree for the multiple sequence
alignment from pairwise distances. Specify the method to compute the
distances of the new nodes to all other nodes. Provide leaf names:
phylotree = seqlinkage(distances,'single',seqs)
Phylogenetic tree object with 32 leaves (31 branches)
View the phylogenetic tree:
view(phylotree)


See Also
cluster | phytree | phytreewrite | plot | seqneighjoin | seqpdist | view




  

 



  

  









    
    
  
  
  
  
  
    
        
    


        
        




  

  

    
 Free Computational Biology Interactive Kit

    
See how to analyze, visualize, and model biological data and systems using MathWorks products.

    Get free kit 

  

    
Trials Available

    
Try the latest computational biology products.

    Get trial software 



     
     
     
    
    
   

      

         © 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS