November 12, 2015
By

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

I spent a decent chunk of my morning trying to figure out how to construct a sparse adjacency matrix for use with `graph.adjacency()`. I’d have thought that this would be rather straight forward, but I tripped over a few subtle issues with the Matrix package. My biggest problem (which in retrospect seems rather trivial) was that elements in my adjacency matrix were occupied by the pipe symbol.

```> adjacency[1:10,1:10]
10 x 10 sparse Matrix of class 'ngCMatrix'

[1,] . . . . . | . . . .
[2,] . . . . . . . | . .
[3,] . . . . . . . . . .
[4,] . . . . . . . . . .
[5,] . . . . | . . . . .
[6,] . . . . . . . . . .
[7,] . . . . . . . . . .
[8,] . . . . . . . . . .
[9,] . . . . . . . . . .
[10,] . | . . . . . . . .
```

Of course, the error message I was encountering didn’t point me to this fact. No, that would have been far too simple! The solution is highlighted in the sample code below: you need to specify the symbol used for the occupied sites in the sparse matrix.

```> library(Matrix)
>
> set.seed(1)
>
> edges = data.frame(i = 1:20, j = sample(1:20, 20, replace = TRUE))
>
> adjacency = sparseMatrix(i = as.integer(edges\$i),
+                          j = as.integer(edges\$j),
+                          x = 1,
+                          dims = rep(20, 2),
+                          use.last.ij = TRUE
+ )
```

The resulting adjacency matrix then looks like this:

```> adjacency[1:10,1:10]
10 x 10 sparse Matrix of class 'dgCMatrix'

[1,] . . . . . 1 . . . .
[2,] . . . . . . . 1 . .
[3,] . . . . . . . . . .
[4,] . . . . . . . . . .
[5,] . . . . 1 . . . . .
[6,] . . . . . . . . . .
[7,] . . . . . . . . . .
[8,] . . . . . . . . . .
[9,] . . . . . . . . . .
[10,] . 1 . . . . . . . .
```

And can be passed into `graph.adjacency()` without any further issues.

```> library(igraph)
```

The post Graph from Sparse Adjacency Matrix appeared first on Exegetic Analytics.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.