# AnalyzeAlignment: Analyze an alignment In GraphAlignment: GraphAlignment

## Description

Analyze an alignment between two networks.

## Usage

 `1` ```AnalyzeAlignment(A, B, R, P, lookupNode, epsilon, clamp) ```

## Arguments

 `A` adjacency matrix for network A `B` adjacency matrix for network B `R` node similarity matrix `P` permutation vector `lookupNode` node bin lookup vector `epsilon` node similarity threshold `clamp` clamp values to range when performing bin lookups

## Details

This function analyzes an alignment and returns various characteristics.

## Value

The return value is a list containing the results. Defined values are:

na - number of aligned node pairs

nb - number of aligned node pairs where neither partner has appreciable sequence similarity with any node in the other network. Appreciable sequence similarity means that r>epsilon. Formally, this is the number of aligned node pairs (ia, ib), where no jb exists such that R[ia, jb] > epsilon and no ja such that R[ja, ib] > epsilon.

nc - number of aligned node pairs, where the partners have no appreciable sequence similarity, but one or both of them has appreciable sequence similarity with some other node in the other network. Formally, this is the number of aligned node pairs (ia, ib) with R[ia, ib] < epsilon but jb or ja exists, such that R[ia, jb] > epsilon or R[ja, ib] > epsilon.

## Author(s)

Joern P. Meier, Michal Kolar, Ville Mustonen, Michael Laessig, and Johannes Berg

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22``` ``` ex<-GenerateExample(dimA=22, dimB=22, filling=.5, covariance=.6, symmetric=TRUE, numOrths=10, correlated=seq(1,18)) pinitial<-InitialAlignment(psize=34, r=ex\$r, mode="reciprocal") lookupLink<-seq(-2,2,.5) linkParams<-ComputeLinkParameters(ex\$a, ex\$b, pinitial, lookupLink) lookupNode<-c(-.5,.5,1.5) nodeParams<-ComputeNodeParameters(dimA=22, dimB=22, ex\$r, pinitial, lookupNode) al<-AlignNetworks(A=ex\$a, B=ex\$b, R=ex\$r, P=pinitial, linkScore=linkParams\$ls, selfLinkScore=linkParams\$ls, nodeScore1=nodeParams\$s1, nodeScore0=nodeParams\$s0, lookupLink=lookupLink, lookupNode=lookupNode, bStart=.1, bEnd=30, maxNumSteps=50) AnalyzeAlignment(A=ex\$a, B=ex\$b, R=ex\$r, P=al, lookupNode, epsilon=.5) ```

### Example output

```\$na
 18

\$nb
 8

\$nc
 0
```

GraphAlignment documentation built on Nov. 8, 2020, 6:56 p.m.