Periodic reordering

Grindrod, Peter and Higham, Desmond J. and Kalna, Gabriela (2010) Periodic reordering. IMA Journal of Numerical Analysis, 30 (1). pp. 195-207. ISSN 0272-4979 (https://doi.org/10.1093/imanum/drp047)

[thumbnail of periodic.pdf]
Preview
 PDF. Filename: periodic.pdf
Accepted Author Manuscript
Download (4MB)| Preview

Abstract

For many networks in nature, science and technology, it is possible to order the nodes so that most links are short-range, connecting near neighbours, and relatively few long-range links, or shortcuts, are present. Given a network as a set of observed links (interactions), the task of finding an ordering of the nodes that reveals such a range dependent structure is closely related to some sparse matrix reordering problems arising in scientific computation. The spectral, or Fiedler vector, approach for sparse matrix reordering has successfully been applied to biological data sets, revealing useful structures and subpatterns. In this work we argue that a periodic analogue of the standard reordering task is also highly relevant. Here, rather than encouraging nonzeros only to lie close to the diagonal of a suitably ordered adjacency matrix, we also allow them to inhabit the off-diagonal corners. Indeed, for the classic small-world model of Watts and Strogatz (Nature, 1998) this type of periodic structure is inherent. We therefore devise and test a new spectral algorithm for periodic reordering. By generalizing the range-dependent random graph class of Grindrod (Phys. Rev. E, 2002) to the periodic case, we can also construct a computable likelihood ratio that suggests whether a given network is inherently linear or periodic. Tests on synthetic data show that the new algorithm can detect periodic structure, even in the presence of noise. Further experiments on real biological data sets then show that some networks are better regarded as periodic than linear. Hence, we find both qualitative (reordered networks plots) and quantitative (likelihood ratios) evidence ofperiodicity in biological networks.

CORE (COnnecting REpositories)