## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Interlocking directorates in Irish companies using a latent space model for bipartite networks

Contributed by Adrian E. Raftery, April 19, 2016 (sent for review April 10, 2015; reviewed by David L. Banks and Tian Zheng)

## Significance

We develop a statistical model for the evolution of the network of leading Irish company directorates over 11 years, before and after the financial crisis of 2008. We focus on company interlocks, whereby a director simultaneously sits on more than one company board. Our analysis indicates that the level of director interlockingness increased from 2003, reaching a peak in 2009, and has generally stabilized since then.

## Abstract

We analyze the temporal bipartite network of the leading Irish companies and their directors from 2003 to 2013, encompassing the end of the Celtic Tiger boom and the ensuing financial crisis in 2008. We focus on the evolution of company interlocks, whereby a company director simultaneously sits on two or more boards. We develop a statistical model for this dataset by embedding the positions of companies and directors in a latent space. The temporal evolution of the network is modeled through three levels of Markovian dependence: one on the model parameters, one on the companies’ latent positions, and one on the edges themselves. The model is estimated using Bayesian inference. Our analysis reveals that the level of interlocking, as measured by a contraction of the latent space, increased before and during the crisis, reaching a peak in 2009, and has generally stabilized since then.

In the past two decades, Ireland has seen big changes to its economic prosperity. From the mid-1990s to 2007, there was a period of rapid economic growth, often referred to as the Celtic Tiger. However, in 2008, there was a major financial and banking crisis that had a devastating effect on the economy and was associated with a series of banking scandals. Since then, there has been considerable interest in ascertaining the roots of the crisis.

Interlocking directorships have been identified as one possible factor in the Irish economic crash (1), possibly contributing to “group think” in the economy and lack of autonomy of individual companies. This possibility is in line with more general literature, which suggests that interlocking directorships may be associated with poorer performance and lower value of companies (2⇓⇓–5), social embeddedness that limits their effectiveness (6), excessive remuneration of directors (7, 8) and conflicts of interest, lack of commitment of directors, and lack of diversity (9, 10).

In this article, we use bipartite networks to investigate interlocking directorships in companies listed on the Irish Stock Exchange (ISE) between 2003 and 2013 inclusive. Bipartite networks have two node types such that two nodes can be linked only if they are of different types. Here, one set of nodes consists of the directors who fill board positions and the other consists of the boards themselves, each of which corresponds to a unique company or organization.

We develop a dynamic statistical model which captures key features of the evolution of interlocking directorships over time. Our data cover the period from 2003 to 2013, allowing an investigation of network behavior in the years before and after the crisis. A key aspect of the model that we develop is that it allows for a statistically principled visual representation of the network’s evolution over time. Clancy et al. (1) examined Irish interlocking directorates, but over a shorter period from 2005 to 2007, which did not cover the financial crash, and they did not carry out any formal statistical modeling.

The modeling task for the type of data we consider has two primary goals. First, we aim to capture patterns of interaction between directors and companies within a given year. Second, we model the temporal evolution of the network and so would like to capture persistence in the interlocking board behavior over time.

Latent space models have been developed for static one-mode networks (11⇓⇓⇓–15). By one-mode, we mean that network contains one node type, with directed or undirected links between different nodes. Statistical models have also been developed for dynamic social networks evolving in time (16⇓⇓–19). Among these models (16) extended the latent space model of ref. 11 to allow for networks that evolve over time by allowing the latent positions of the nodes to change according to a random walk.

Several models have been proposed for static bipartite networks and cooccurrence data (20). Models for dynamic data of this type have also been proposed (10, 21, 22).

The latent space model we develop is for dynamic bipartite networks. The approach most similar to the one we adopt here is that of Sarkar et al. (23), which is a dynamic adaptation of the latent space model for cooccurrence data of Globerson et al. (24). Sarkar et al. (23) examined cooccurrences of words and authors in text corpora. The model used attaches a latent position in a low-dimensional real space to each word and author. It models the probability of cooccurrence of words and authors using the empirical distribution of word and author occurrences and a decreasing function of the distance between their latent positions. An advantage of this model is that it does not make any assumption about the distributions of ties directly, allowing natural modeling of cooccurrence counts through only the empirical distribution of author/word occurrences.

We model evolution over time of the ties between nodes using a Markov model, extending the model of Sarkar and Moore (16), who assumed the nodes to be independent given the latent positions. We found that the conditional independence assumption misses the persistence in links observed in our director board membership network, due to directors retaining a board position for a number of consecutive years. The model we propose does account for this persistence.

## Methods

### Data.

A list of all Irish-based companies listed on the ISE from 2003 to 2013 was obtained from the Exchange. The company registration office number and address of each of the remaining companies was obtained from the Irish Companies Registration Office website. Hence, for each company, a list of their directors was found. The data collected were rearranged into a bipartite dynamic graph, whereby for each year edges between directors and companies were created according to boards’ compositions. The network was made up of 1,009 directors, 91 companies, and 4,952 edges. Because the main focus of our analysis is the study of director interlocks, we removed those companies that did not involve any interlocking links over the study period, in other words, those companies that did not share any of their directors with any other board over the 11 years of the study. The resulting dataset (corresponding to 761 directors, 59 companies, and 3,855 edges) was used to carry out our analysis. Only 22 of the companies were quoted on the ISE for the entire time period considered, whereas all other companies either joined at a later stage or left early or both. There was also a significant number of directors who sat only on one board, sometimes for only a year. These directors can be seen as peripheral to the Irish network of boards and directors. Table 1 gives an overview of the number of directors and boards in each year of the study.

### Exploratory Analysis.

An interlocked director in any given year is a director who sits on two or more boards in that same year. It can be seen from Table 1 that the proportion of interlocked directors increased from 2006 to 2010 and declined from 2011 to 2013.

The number of interlocking links for a node is related to the excess degree, defined as follows: for a director with degree *d*, it is equal to

### A Dynamic Latent Space Model for Interlocking Boards.

The network is described by the adjacency cube

for

We use a latent space approach to model the positions of companies and their directors over time. This approach uses the idea that the positions of companies and their directors lie in a *D*-dimensional latent Euclidean space, so that a director which lies close to a company will tend to sit on that company board. We denote latent director and company positions at time *t*, it may persist over time. We similarly include a nonedge persistence parameter, *t* also typically persists over time. Separating the edge and nonedge persistences is an important facet of our model, because the observed temporal bipartite graph is sparse and the effect of persistence of edges and nonedges is markedly different, as can be seen in *Results and Discussion*. Conditional on these latent positions, the log-odds of director *i* sitting on board *j* in year *t* is

In fact, we found for the Irish boards data that modeling changes in director positions over time gave little extra information but was considerably more computationally demanding; hence, these are kept fixed. The persistence feature is included through the combination of the intercepts *i* and *j* at time *t*. If the state of the edge in the previous time frame is not known, the quantity *T*-dimensional vectors

We account for the missing data by considering the contribution of boards only in those years where they are quoted in the ISE. As a matter of fact, from 2003 to 2013, no boards joined or left the ISE more than once. Hence, for each board the period of activity is made up of consecutive years only. A similar reasoning is applied to directors, whereas their contribution to the likelihood is considered only in the years where they sit on one board at least. The rationale leading to this simplification is that the likelihood information used to estimate positions is mainly carried by edges that are present rather than those that are missing. Hence, we speculate that neglecting the contribution of inactive boards and directors does not cause relevant changes in our conclusions.

Therefore, we express the likelihood of the network as

where *t*.

### Bayesian Model.

We formulate our dynamic latent space model as a Bayesian hierarchical model. We denote by *s* and with mean μ and variance ν, and by *s* with shape *u* and rate *v*. We use a random walk model for the boards’ latent positions, specified by

The directors’ positions are assumed to arise independently from a *D*-variate zero mean Gaussian prior with spherical covariance, hence

Note that the parameters corresponding to inactive boards and directors are included in the prior modeling. The intercept parameters follow a random walk prior:

The precision parameters *a* and rate *b*. (The precision is the reciprocal of the variance.) The remaining parameters

Markov Chain Monte Carlo (MCMC) sampling was used to sample from the joint posterior distribution of all of the model parameters. One issue is that the likelihood is invariant to rotations, reflections, and translations of the latent positions, because the likelihood depends on the latent positions through the distances *SI Appendix*.

### Quantifying the Extent of Interlockingness.

Quantifying the extent of interlockingness in a network of companies and directors may be of interest to policy makers and the public, especially in the aftermath of the economic crises in Ireland and beyond. A useful byproduct of the output of our statistical model is that the distribution of latent positions can be used to derive a metric of interlockingness in each year. In fact, we expect that directors located near the center of the space will more likely sit on more than one board, whereas those having a peripheral position will tend to have fewer connections. Equivalently, the centrality of boards affects the interlocking level of the network in the same way. For this reason, we assess the variation of interlocking through the empirical variance of boards, which may be computed for each iteration of the MCMC algorithm at each time

where

### Model Choice.

We also compare our model with two simpler models for the same data. The three models are as follows:

Model 1: the model described in this section.

Model 2: a special case of Model 1, where the latent distances are all set to zero. Eq.

**1**therefore simplifies to

Model 3: a special case of Model 1, where the persistence feature is omitted. Therefore, the log-odds of Eq.

**1**is expressed as

Model choice is carried out using the deviance information criteria (DIC) (25), which is estimated from the output of the MCMC algorithm.

## Results and Discussion

The DIC values for the three competing models are provided in Table 2, along with the number of effective parameters and the actual number of parameters. Our model achieves the lowest DIC and, therefore, we focus on it for the remainder of the analysis.

Estimates of the intercept and precision parameters are shown in Tables 3–5. The precision parameters

Fig. 3 illustrates the fact that the persistence feature introduced in our model plays an important role because the log-odds for the probability of an edge can be markedly different when the same edge was present in the previous year. The difference in the posterior SDs between the two types of the intercept parameters is a consequence of the network being sparse: more data are available for the estimation of the

The posterior means of the latent positions of both directors and companies are shown in Fig. 4 by a time-collapsed snapshot. For the boards, a clustering effect can be noted due to the little movement over time, whereas for directors, clusters of points arise due to boards’ compositions remaining relatively unchanged over time. The cluster of companies and directors just below the center of the space corresponds to a higher density region of the graph that includes most of the leading companies of the Irish economy, whereas isolated boards typically are companies with a low market capitalization or a short longevity in the ISE.

The empirical variance, outlined in Eq. **6**, and shown in Fig. 5, *Left* clearly exhibits a downward trend from 2003 to 2009, signaling a contraction of the latent space. This result is in agreement with the increase of interlocking over this period, as presented in Fig. 1. From 2009 to 2013, the empirical variance index appears to stabilize, coinciding with an upswing in the economy. The motivation behind the contraction of the latent space can be traced back to two main causes: the movement of boards toward the center of the space and the variations induced by boards joining the ISE after 2003 or leaving before 2013 (recall that the empirical variance is evaluated using only active boards). The plot in Fig. 5, *Right* shows the average distance between the center and board *j* at the times of its first and last appearance on the ISE, for every

The combination of the evolution of boards’ positions and the changes in the composition of the ISE basket give a reasonable justification for the decrease of the empirical variance and hence for the contraction of the latent space. This finding yields a grounded theoretical support for the increase of interlockingness argued in the paper.

## Conclusions

We have proposed a latent space model for dynamic bipartite networks, motivated by Irish corporate board interlocks before and after the 2008 financial crash. The model can capture heterogeneity in the network by means of the latent positions as well as the persistence of links through a modeling structure for intercept parameters. In our application, we have mainly focused on the relation between the variation of interlocking level in the Irish directorates network and the compression/expansion of the latent space. Our model has captured the dynamic evolution of this interlockingness, exhibiting a contraction of the latent space up until 2009 and a stabilization thereafter. This result may provide some support for the view that there are disadvantages to interlocking, and that increasing interlocking may be a warning sign of, and possible contributing factor to future financial instability.

## Acknowledgments

We are grateful to the reviewers and Katherine Stovel for helpful comments. This research was supported by the Insight Centre for Data Analytics through Science Foundation Ireland Grant SFI/12/RC/2289 (to N.F.); by Science Foundation Ireland Grants 12/IP/1424 (to N.F. and R.R.), 08/IN.1/I1879 (to J.W.), and 11/W.1/I2079 (to A.E.R.); and by National Institutes of Health Grants R01 HD054511 and R01 HD070936 (to A.E.R.).

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: raftery{at}uw.edu.

Author contributions: N.F. and A.E.R. designed research; N.F., R.R., J.W., and A.E.R. performed research; N.F., R.R., J.W., and A.E.R. analyzed data; and N.F., J.W., and A.E.R. wrote the paper.

Reviewers: D.L.B., Duke University; and T.Z., Columbia University.

The authors declare no conflict of interest.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1606295113/-/DCSupplemental.

Freely available online through the PNAS open access option.

## References

- ↵.
- Clancy P,
- O’Connor N,
- Dillon K

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Raftery AE,
- Niu X,
- Hoff PD,
- Yeung KY

- ↵
- ↵
- ↵.
- Melia M,
- Shen X

- Snijders TAB

- ↵.
- Foulds J,
- DuBois C,
- Asuncion AU,
- Butts CT,
- Smyth P

*JMLR Workshop Conf Proc*15:287–295. - ↵
- ↵
- ↵.
- Tang L,
- Liu H,
- Zhang J,
- Nazeri Z

- ↵.
- Sun Y,
- Tang J,
- Han J,
- Chen C,
- Gupta M

- ↵.
- Sarkar P,
- Siddiqi SM,
- Gordon GJ

*JMLR Workshop Conf Proc*2:420–427. - ↵.
- Globerson A,
- Chechik G,
- Pereira F,
- Tishby N

- ↵

## Citation Manager Formats

### More Articles of This Classification

### Physical Sciences

### Statistics

### Related Content

- No related articles found.

### Cited by...

- No citing articles found.