Bioinformatics Toolbox™

graphpred2path - Convert predecessor indices to paths

Syntax

path = graphpred2path(pred, D)

Arguments

`pred`	Row vector or matrix of predecessor node indices. The value of the root (or source) node in `pred` must be 0.
`D`	Destination node in `pred`.

Description

Tip For introductory information on graph theory functions, see Graph Theory Functions in the Bioinformatics Toolbox User's Guide.

path = graphpred2path(pred, D) traces back a path by following the predecessor list in pred starting at destination node D.

The value of the root (or source) node in pred must be 0. If a NaN is found when following the predecessor nodes, graphpred2path returns an empty path.

If `pred` is a ...	And `D` is a ...	Then `path` is a ...
row vector of predecessor node indices	scalar	row vector listing the nodes from the root (or source) to `D`.
row vector of predecessor node indices	row vector	row cell array with every column containing the path to the destination for every element in `D`.
matrix	scalar	column cell array with every row containing the path for every row in `pred`.
matrix	row vector	matrix cell array with every row containing the paths for the respective row in `pred`, and every column containing the paths to the respective destination in `D`.

Note If D is omitted, the paths to all the destinations are calculated for every predecessor listed in pred.

Examples

Create a phytree object from the phylogenetic tree file for the GLR_HUMAN protein.

tr = phytreeread('pf00002.tree')
    Phylogenetic tree object with 33 leaves (32 branches)

View the phytree object.
```
view(tr)
```
From the phytree object, create a connection matrix to represent the phylogenetic tree.
```
[CM,labels,dist] = getmatrix(tr);
```

Find the nodes from the root to one leaf in the phylogenetic tree created from the phylogenetic tree file for the GLR_HUMAN protein.

root_loc = size(CM,1)

root_loc =

    65

glr_loc = strmatch('GLR',labels)

glr_loc =

    28

[T,PRED]=graphminspantree(CM,root_loc);
PATH = graphpred2path(PRED,glr_loc)

PATH =

    65    64    53    52    46    45    44    43    28

References