Serveur d'exploration Covid (26 mars)

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models

Identifieur interne : 000984 ( Pmc/Corpus ); précédent : 000983; suivant : 000985

A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models

Auteurs : David A. Rolls ; Peng Wang ; Emma Mcbryde ; Philippa Pattison ; Garry Robins

Source :

RBID : PMC:4640514

Abstract

We compare two broad types of empirically grounded random network models in terms of their abilities to capture both network features and simulated Susceptible-Infected-Recovered (SIR) epidemic dynamics. The types of network models are exponential random graph models (ERGMs) and extensions of the configuration model. We use three kinds of empirical contact networks, chosen to provide both variety and realistic patterns of human contact: a highly clustered network, a bipartite network and a snowball sampled network of a “hidden population”. In the case of the snowball sampled network we present a novel method for fitting an edge-triangle model. In our results, ERGMs consistently capture clustering as well or better than configuration-type models, but the latter models better capture the node degree distribution. Despite the additional computational requirements to fit ERGMs to empirical networks, the use of ERGMs provides only a slight improvement in the ability of the models to recreate epidemic features of the empirical network in simulated SIR epidemics. Generally, SIR epidemic results from using configuration-type models fall between those from a random network model (i.e., an Erdős-Rényi model) and an ERGM. The addition of subgraphs of size four to edge-triangle type models does improve agreement with the empirical network for smaller densities in clustered networks. Additional subgraphs do not make a noticeable difference in our example, although we would expect the ability to model cliques to be helpful for contact networks exhibiting household structure.


Url:
DOI: 10.1371/journal.pone.0142181
PubMed: 26555701
PubMed Central: 4640514

Links to Exploration step

PMC:4640514

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models</title>
<author>
<name sortKey="Rolls, David A" sort="Rolls, David A" uniqKey="Rolls D" first="David A." last="Rolls">David A. Rolls</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Wang, Peng" sort="Wang, Peng" uniqKey="Wang P" first="Peng" last="Wang">Peng Wang</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Mcbryde, Emma" sort="Mcbryde, Emma" uniqKey="Mcbryde E" first="Emma" last="Mcbryde">Emma Mcbryde</name>
<affiliation>
<nlm:aff id="aff002">
<addr-line>Department of Medicine-RMH, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Pattison, Philippa" sort="Pattison, Philippa" uniqKey="Pattison P" first="Philippa" last="Pattison">Philippa Pattison</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Robins, Garry" sort="Robins, Garry" uniqKey="Robins G" first="Garry" last="Robins">Garry Robins</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">26555701</idno>
<idno type="pmc">4640514</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4640514</idno>
<idno type="RBID">PMC:4640514</idno>
<idno type="doi">10.1371/journal.pone.0142181</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000984</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000984</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models</title>
<author>
<name sortKey="Rolls, David A" sort="Rolls, David A" uniqKey="Rolls D" first="David A." last="Rolls">David A. Rolls</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Wang, Peng" sort="Wang, Peng" uniqKey="Wang P" first="Peng" last="Wang">Peng Wang</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Mcbryde, Emma" sort="Mcbryde, Emma" uniqKey="Mcbryde E" first="Emma" last="Mcbryde">Emma Mcbryde</name>
<affiliation>
<nlm:aff id="aff002">
<addr-line>Department of Medicine-RMH, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Pattison, Philippa" sort="Pattison, Philippa" uniqKey="Pattison P" first="Philippa" last="Pattison">Philippa Pattison</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Robins, Garry" sort="Robins, Garry" uniqKey="Robins G" first="Garry" last="Robins">Garry Robins</name>
<affiliation>
<nlm:aff id="aff001">
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint>
<date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>We compare two broad types of empirically grounded random network models in terms of their abilities to capture both network features and simulated Susceptible-Infected-Recovered (SIR) epidemic dynamics. The types of network models are exponential random graph models (ERGMs) and extensions of the configuration model. We use three kinds of empirical contact networks, chosen to provide both variety and realistic patterns of human contact: a highly clustered network, a bipartite network and a snowball sampled network of a “hidden population”. In the case of the snowball sampled network we present a novel method for fitting an edge-triangle model. In our results, ERGMs consistently capture clustering as well or better than configuration-type models, but the latter models better capture the node degree distribution. Despite the additional computational requirements to fit ERGMs to empirical networks, the use of ERGMs provides only a slight improvement in the ability of the models to recreate epidemic features of the empirical network in simulated SIR epidemics. Generally, SIR epidemic results from using configuration-type models fall between those from a random network model (i.e., an Erdős-Rényi model) and an ERGM. The addition of subgraphs of size four to edge-triangle type models does improve agreement with the empirical network for smaller densities in clustered networks. Additional subgraphs do not make a noticeable difference in our example, although we would expect the ability to model cliques to be helpful for contact networks exhibiting household structure.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct>
<analytic>
<author>
<name sortKey="Welch, D" uniqKey="Welch D">D Welch</name>
</author>
<author>
<name sortKey="Bansal, S" uniqKey="Bansal S">S Bansal</name>
</author>
<author>
<name sortKey="Hunter, Dr" uniqKey="Hunter D">DR Hunter</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Newman, Mej" uniqKey="Newman M">MEJ Newman</name>
</author>
<author>
<name sortKey="Strogatz, Sh" uniqKey="Strogatz S">SH Strogatz</name>
</author>
<author>
<name sortKey="Watts, Dj" uniqKey="Watts D">DJ Watts</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Molloy, M" uniqKey="Molloy M">M Molloy</name>
</author>
<author>
<name sortKey="Reed, B" uniqKey="Reed B">B Reed</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Frank, O" uniqKey="Frank O">O Frank</name>
</author>
<author>
<name sortKey="Strauss, D" uniqKey="Strauss D">D Strauss</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Newman, Mej" uniqKey="Newman M">MEJ Newman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kenah, E" uniqKey="Kenah E">E Kenah</name>
</author>
<author>
<name sortKey="Robins, Jm" uniqKey="Robins J">JM Robins</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Meyers, La" uniqKey="Meyers L">LA Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Newman, Mej" uniqKey="Newman M">MEJ Newman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Newman, Mej" uniqKey="Newman M">MEJ Newman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Volz, Em" uniqKey="Volz E">EM Volz</name>
</author>
<author>
<name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
<author>
<name sortKey="Galvani, A" uniqKey="Galvani A">A Galvani</name>
</author>
<author>
<name sortKey="Ancel Meyers, L" uniqKey="Ancel Meyers L">L Ancel Meyers</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="House, T" uniqKey="House T">T House</name>
</author>
<author>
<name sortKey="Keeling, Mj" uniqKey="Keeling M">MJ Keeling</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Smieszek, T" uniqKey="Smieszek T">T Smieszek</name>
</author>
<author>
<name sortKey="Fiebig, L" uniqKey="Fiebig L">L Fiebig</name>
</author>
<author>
<name sortKey="Scholz, R" uniqKey="Scholz R">R Scholz</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Karrer, B" uniqKey="Karrer B">B Karrer</name>
</author>
<author>
<name sortKey="Newman, Mej" uniqKey="Newman M">MEJ Newman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
<author>
<name sortKey="Kalish, Y" uniqKey="Kalish Y">Y Kalish</name>
</author>
<author>
<name sortKey="Lusher, D" uniqKey="Lusher D">D Lusher</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
<author>
<name sortKey="Snijders, T" uniqKey="Snijders T">T Snijders</name>
</author>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
<author>
<name sortKey="Handcock, M" uniqKey="Handcock M">M Handcock</name>
</author>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Snijders, Tab" uniqKey="Snijders T">TAB Snijders</name>
</author>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
<author>
<name sortKey="Handcock, Ms" uniqKey="Handcock M">MS Handcock</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pattison, Pe" uniqKey="Pattison P">PE Pattison</name>
</author>
<author>
<name sortKey="Robins, Gl" uniqKey="Robins G">GL Robins</name>
</author>
<author>
<name sortKey="Snijders, Tab" uniqKey="Snijders T">TAB Snijders</name>
</author>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Goodman, La" uniqKey="Goodman L">LA Goodman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Handcock, Ms" uniqKey="Handcock M">MS Handcock</name>
</author>
<author>
<name sortKey="Gile, Kj" uniqKey="Gile K">KJ Gile</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Rolls, Da" uniqKey="Rolls D">DA Rolls</name>
</author>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
<author>
<name sortKey="Jenkinson, R" uniqKey="Jenkinson R">R Jenkinson</name>
</author>
<author>
<name sortKey="Pattison, Pe" uniqKey="Pattison P">PE Pattison</name>
</author>
<author>
<name sortKey="Robins, Gl" uniqKey="Robins G">GL Robins</name>
</author>
<author>
<name sortKey="Sacks Davis, R" uniqKey="Sacks Davis R">R Sacks-Davis</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Salganik, Mj" uniqKey="Salganik M">MJ Salganik</name>
</author>
<author>
<name sortKey="Heckathorn, Dd" uniqKey="Heckathorn D">DD Heckathorn</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Snijders, Tab" uniqKey="Snijders T">TAB Snijders</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gile, Kj" uniqKey="Gile K">KJ Gile</name>
</author>
<author>
<name sortKey="Handcock, Ms" uniqKey="Handcock M">MS Handcock</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Illenberger, J" uniqKey="Illenberger J">J Illenberger</name>
</author>
<author>
<name sortKey="Flotterod, G" uniqKey="Flotterod G">G Flötteröd</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Salathe, M" uniqKey="Salathe M">M Salathé</name>
</author>
<author>
<name sortKey="Kazandjieva, M" uniqKey="Kazandjieva M">M Kazandjieva</name>
</author>
<author>
<name sortKey="Lee, Jw" uniqKey="Lee J">JW Lee</name>
</author>
<author>
<name sortKey="Levis, P" uniqKey="Levis P">P Levis</name>
</author>
<author>
<name sortKey="Feldman, Mw" uniqKey="Feldman M">MW Feldman</name>
</author>
<author>
<name sortKey="Jones, Jh" uniqKey="Jones J">JH Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bearman, Ps" uniqKey="Bearman P">PS Bearman</name>
</author>
<author>
<name sortKey="Moody, J" uniqKey="Moody J">J Moody</name>
</author>
<author>
<name sortKey="Stovel, K" uniqKey="Stovel K">K Stovel</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gilbert, En" uniqKey="Gilbert E">EN Gilbert</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hunter, Dr" uniqKey="Hunter D">DR Hunter</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
<author>
<name sortKey="Volz, Em" uniqKey="Volz E">EM Volz</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Skvoretz, J" uniqKey="Skvoretz J">J Skvoretz</name>
</author>
<author>
<name sortKey="Faust, K" uniqKey="Faust K">K Faust</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
<author>
<name sortKey="Sharpe, K" uniqKey="Sharpe K">K Sharpe</name>
</author>
<author>
<name sortKey="Robins, Gl" uniqKey="Robins G">GL Robins</name>
</author>
<author>
<name sortKey="Pattison, Pe" uniqKey="Pattison P">PE Pattison</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Wang, P" uniqKey="Wang P">P Wang</name>
</author>
<author>
<name sortKey="Robins, G" uniqKey="Robins G">G Robins</name>
</author>
<author>
<name sortKey="Pattison, P" uniqKey="Pattison P">P Pattison</name>
</author>
<author>
<name sortKey="Lazega, E" uniqKey="Lazega E">E Lazega</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Erd S, P" uniqKey="Erd S P">P Erdős</name>
</author>
<author>
<name sortKey="Renyi, A" uniqKey="Renyi A">A Rényi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
<author>
<name sortKey="Slim, Ac" uniqKey="Slim A">AC Slim</name>
</author>
<author>
<name sortKey="Volz, Em" uniqKey="Volz E">EM Volz</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Holme, P" uniqKey="Holme P">P Holme</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Machens, A" uniqKey="Machens A">A Machens</name>
</author>
<author>
<name sortKey="Gesualdo, F" uniqKey="Gesualdo F">F Gesualdo</name>
</author>
<author>
<name sortKey="Rizzo, C" uniqKey="Rizzo C">C Rizzo</name>
</author>
<author>
<name sortKey="Tozzi, A" uniqKey="Tozzi A">A Tozzi</name>
</author>
<author>
<name sortKey="Barrat, A" uniqKey="Barrat A">A Barrat</name>
</author>
<author>
<name sortKey="Cattuto, C" uniqKey="Cattuto C">C Cattuto</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Isella, L" uniqKey="Isella L">L Isella</name>
</author>
<author>
<name sortKey="Stehle, J" uniqKey="Stehle J">J Stehlé</name>
</author>
<author>
<name sortKey="Barrat, A" uniqKey="Barrat A">A Barrat</name>
</author>
<author>
<name sortKey="Cattuto, C" uniqKey="Cattuto C">C Cattuto</name>
</author>
<author>
<name sortKey="Pinton, J" uniqKey="Pinton J">J Pinton</name>
</author>
<author>
<name sortKey="Van Den Broeck, W" uniqKey="Van Den Broeck W">W Van den Broeck</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Stehle, J" uniqKey="Stehle J">J Stehlé</name>
</author>
<author>
<name sortKey="Voirin, N" uniqKey="Voirin N">N Voirin</name>
</author>
<author>
<name sortKey="Barrat, A" uniqKey="Barrat A">A Barrat</name>
</author>
<author>
<name sortKey="Cattuto, C" uniqKey="Cattuto C">C Cattuto</name>
</author>
<author>
<name sortKey="Isella, L" uniqKey="Isella L">L Isella</name>
</author>
<author>
<name sortKey="Pinton, J" uniqKey="Pinton J">J Pinton</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Stehle, J" uniqKey="Stehle J">J Stehlé</name>
</author>
<author>
<name sortKey="Voirin, N" uniqKey="Voirin N">N Voirin</name>
</author>
<author>
<name sortKey="Barrat, A" uniqKey="Barrat A">A Barrat</name>
</author>
<author>
<name sortKey="Cattuto, C" uniqKey="Cattuto C">C Cattuto</name>
</author>
<author>
<name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author>
<name sortKey="Isella, L" uniqKey="Isella L">L Isella</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Salathe, M" uniqKey="Salathe M">M Salathé</name>
</author>
<author>
<name sortKey="Bengtsson, L" uniqKey="Bengtsson L">L Bengtsson</name>
</author>
<author>
<name sortKey="Bodnar, Tj" uniqKey="Bodnar T">TJ Bodnar</name>
</author>
<author>
<name sortKey="Brewer, Dd" uniqKey="Brewer D">DD Brewer</name>
</author>
<author>
<name sortKey="Brownstein, Js" uniqKey="Brownstein J">JS Brownstein</name>
</author>
<author>
<name sortKey="Buckee, C" uniqKey="Buckee C">C Buckee</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Barabasi, Al" uniqKey="Barabasi A">AL Barabási</name>
</author>
<author>
<name sortKey="Albert, R" uniqKey="Albert R">R Albert</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Watts, Dj" uniqKey="Watts D">DJ Watts</name>
</author>
<author>
<name sortKey="Strogatz, Sh" uniqKey="Strogatz S">SH Strogatz</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pastor Satorras, R" uniqKey="Pastor Satorras R">R Pastor-Satorras</name>
</author>
<author>
<name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, CA USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">26555701</article-id>
<article-id pub-id-type="pmc">4640514</article-id>
<article-id pub-id-type="publisher-id">PONE-D-14-27406</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0142181</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models</article-title>
<alt-title alt-title-type="running-head">Comparing Epidemics on ERGMs and Configuration Type Network Models</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Rolls</surname>
<given-names>David A.</given-names>
</name>
<xref ref-type="aff" rid="aff001">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="currentaff001">
<sup>¤</sup>
</xref>
<xref ref-type="corresp" rid="cor001">*</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Peng</given-names>
</name>
<xref ref-type="aff" rid="aff001">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>McBryde</surname>
<given-names>Emma</given-names>
</name>
<xref ref-type="aff" rid="aff002">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pattison</surname>
<given-names>Philippa</given-names>
</name>
<xref ref-type="aff" rid="aff001">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Robins</surname>
<given-names>Garry</given-names>
</name>
<xref ref-type="aff" rid="aff001">
<sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff001">
<label>1</label>
<addr-line>Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</aff>
<aff id="aff002">
<label>2</label>
<addr-line>Department of Medicine-RMH, University of Melbourne, Melbourne, VIC 3010, Australia</addr-line>
</aff>
<contrib-group>
<contrib contrib-type="editor">
<name>
<surname>Moreno</surname>
<given-names>Yamir</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1">
<addr-line>University of Zaragoza, SPAIN</addr-line>
</aff>
<author-notes>
<fn fn-type="COI-statement" id="coi001">
<p>
<bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con" id="contrib001">
<p>Conceived and designed the experiments: DAR PW EM PP GR. Performed the experiments: DAR. Analyzed the data: DAR. Contributed reagents/materials/analysis tools: DAR PW. Wrote the paper: DAR EM GR.</p>
</fn>
<fn fn-type="current-aff" id="currentaff001">
<label>¤</label>
<p>Current address: Melbourne School of Psychological Sciences, University of Melbourne, Melbourne, VIC 3010, Australia</p>
</fn>
<corresp id="cor001">* E-mail:
<email>drolls@unimelb.edu.au</email>
</corresp>
</author-notes>
<pub-date pub-type="collection">
<year>2015</year>
</pub-date>
<pub-date pub-type="epub">
<day>10</day>
<month>11</month>
<year>2015</year>
</pub-date>
<volume>10</volume>
<issue>11</issue>
<elocation-id>e0142181</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>6</month>
<year>2014</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>10</month>
<year>2015</year>
</date>
</history>
<permissions>
<copyright-statement>© 2015 Rolls et al</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Rolls et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:type="simple" xlink:href="pone.0142181.pdf"></self-uri>
<abstract>
<p>We compare two broad types of empirically grounded random network models in terms of their abilities to capture both network features and simulated Susceptible-Infected-Recovered (SIR) epidemic dynamics. The types of network models are exponential random graph models (ERGMs) and extensions of the configuration model. We use three kinds of empirical contact networks, chosen to provide both variety and realistic patterns of human contact: a highly clustered network, a bipartite network and a snowball sampled network of a “hidden population”. In the case of the snowball sampled network we present a novel method for fitting an edge-triangle model. In our results, ERGMs consistently capture clustering as well or better than configuration-type models, but the latter models better capture the node degree distribution. Despite the additional computational requirements to fit ERGMs to empirical networks, the use of ERGMs provides only a slight improvement in the ability of the models to recreate epidemic features of the empirical network in simulated SIR epidemics. Generally, SIR epidemic results from using configuration-type models fall between those from a random network model (i.e., an Erdős-Rényi model) and an ERGM. The addition of subgraphs of size four to edge-triangle type models does improve agreement with the empirical network for smaller densities in clustered networks. Additional subgraphs do not make a noticeable difference in our example, although we would expect the ability to model cliques to be helpful for contact networks exhibiting household structure.</p>
</abstract>
<funding-group>
<funding-statement>DAR was partly funded by the Australian Research Council (Discovery project grant 0987730,
<ext-link ext-link-type="uri" xlink:href="http://www.arc.gov.au/ncgp/dp/dp_default.htm">http://www.arc.gov.au/ncgp/dp/dp_default.htm</ext-link>
) and partly by the University of Melbourne. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts>
<fig-count count="14"></fig-count>
<table-count count="7"></table-count>
<page-count count="36"></page-count>
</counts>
<custom-meta-group>
<custom-meta id="data-availability">
<meta-name>Data Availability</meta-name>
<meta-value>The Supporting Information files contain MATLAB datasets with the node degrees, graph statistics, and SIR simulation results from this study along with the node degrees of snowball sampled empirical and simulated networks using estimators one to six. All empirical networks are held by third parties. Data for the high school contact network is publicly available with that paper. Our depiction of the relationships network is new, based on manually recoding the network depiction in Bearman et al. For the original network data, interested parties should contact the authors of that paper or the AddHealth study team. Data for the PWID network is held by the Burnet Institute. Researchers with serious enquiries should contact Dr. Margaret Hellard (
<email>hellard@burnet.edu.au</email>
).</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
<notes>
<title>Data Availability</title>
<p>The Supporting Information files contain MATLAB datasets with the node degrees, graph statistics, and SIR simulation results from this study along with the node degrees of snowball sampled empirical and simulated networks using estimators one to six. All empirical networks are held by third parties. Data for the high school contact network is publicly available with that paper. Our depiction of the relationships network is new, based on manually recoding the network depiction in Bearman et al. For the original network data, interested parties should contact the authors of that paper or the AddHealth study team. Data for the PWID network is held by the Burnet Institute. Researchers with serious enquiries should contact Dr. Margaret Hellard (
<email>hellard@burnet.edu.au</email>
).</p>
</notes>
</front>
<body>
<sec sec-type="intro" id="sec001">
<title>Introduction</title>
<p>It is recognised that contact networks can play an important role in studying and understanding epidemic dynamics [
<xref rid="pone.0142181.ref001" ref-type="bibr">1</xref>
]. For these purposes, network models are important for at least two reasons. First, researchers may have an empirical network, but want to simulate additional, similar networks to avoid “over-fitting” to the particular observed network. Second, researchers may have no data on empirical contact patterns, so use a model, possibly informed by an empirical or textbook distribution, for example. One family of network models arises from a collection of network research expanding on the standard “configuration model” [
<xref rid="pone.0142181.ref002" ref-type="bibr">2</xref>
,
<xref rid="pone.0142181.ref003" ref-type="bibr">3</xref>
]. An alternative approach uses the family of Exponential Random Graph models [
<xref rid="pone.0142181.ref004" ref-type="bibr">4</xref>
], combined with simulation of disease transmission. In this paper we try to make some progress in bridging these two network modelling domains in the context of network epidemiology. We do not try to investigate the effects of node degree variation or clustering on epidemic dynamics, per se. (Welch et al. [
<xref rid="pone.0142181.ref001" ref-type="bibr">1</xref>
] has a nice summary of current results.) Nor do we study the suitability of theoretical solutions, that are asymptotically correct, to finite empirical network sizes. (We recognise some researchers use variations of the configuration model, but then simulate a disease transmission process.) Rather, we are trying to make some progress in understanding which subgraphs are relevant to the study of epidemic dynamics and how the strengths and shortcomings of each modelling domain affect conclusions on epidemic dynamics. A key feature of our work is a focus on empirical contact networks, which therefore have realistic patterns of human contact.</p>
<p>We can list several desirable properties of network models for the study of epidemic dynamics. 1) Key network features (e.g., giant component size, location of phase transitions) can be found without Monte Carlo simulation of an ensemble of networks. 2) They capture “relevant” features of networks (e.g., clustering, homophily), although a sufficient list of such features to study epidemic dynamics is unknown, and which network features are relevant probably varies depending on disease features such as infectivity. 3) Epidemic dynamics are solvable. Key quantities of interest like final size, epidemic duration, and the infectivity threshold for epidemics, are solvable without explicit stochastic simulation of the epidemic (e.g., using a system of ODEs or bond percolation results). 4) Variability in quantities of interest are also solvable. While deterministic models of the infection process provide average (in some sense) values, probabilistic models can provide the variability, and maybe even a distribution.</p>
<sec id="sec002">
<title>Configuration-Type Models</title>
<p>Research based on the configuration model [
<xref rid="pone.0142181.ref002" ref-type="bibr">2</xref>
,
<xref rid="pone.0142181.ref003" ref-type="bibr">3</xref>
] attempts to create a model with these desirable properties. The “standard configuration model” creates networks with a specified degree sequence. In short, a node
<italic>i</italic>
with degree
<italic>d</italic>
<sub>
<italic>i</italic>
</sub>
is allocated
<italic>d</italic>
<sub>
<italic>i</italic>
</sub>
“stubs”. Network ties are formed by connecting pairs of stubs at random, thus achieving a network with desired degree sequence. The standard configuration model has proven to be amenable to the study of Susceptible-Infected-Recovered (SIR) epidemics (e.g, [
<xref rid="pone.0142181.ref005" ref-type="bibr">5</xref>
<xref rid="pone.0142181.ref007" ref-type="bibr">7</xref>
]). For example, a reproduction ratio,
<italic>R</italic>
<sub>0</sub>
, and epidemic threshold analogous to those of compartmental models can be expressed. Define transmissibility [
<xref rid="pone.0142181.ref005" ref-type="bibr">5</xref>
],
<italic>T</italic>
, in a homogeneous population with a constant infectious period to be the probability an infected node infects a susceptible neighbour along a single edge. In the early stages of an outbreak on a configuration model network, the expected number of infections caused by a newly infected single node is
<disp-formula id="pone.0142181.e001">
<alternatives>
<graphic xlink:href="pone.0142181.e001.jpg" id="pone.0142181.e001g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M1">
<mml:mtable displaystyle="true">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>T</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</alternatives>
<label>(1)</label>
</disp-formula>
where
<italic>E</italic>
(
<italic>d</italic>
) is the mean node degree and
<italic>E</italic>
(
<italic>d</italic>
<sup>2</sup>
) is the mean of the squared node degree (closely related to the variance). For heterogeneous populations and non-constant infectious periods the expression can be generalized to
<disp-formula id="pone.0142181.e002">
<alternatives>
<graphic xlink:href="pone.0142181.e002.jpg" id="pone.0142181.e002g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M2">
<mml:mtable displaystyle="true">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</alternatives>
</disp-formula>
where
<italic>E</italic>
(
<italic>T</italic>
) is the average transmissibility over all infectious and susceptible nodes [
<xref rid="pone.0142181.ref006" ref-type="bibr">6</xref>
,
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
]. Using bond percolation arguments, an epidemic is possible only when
<italic>R</italic>
<sub>0</sub>
> 1 (as for compartmental models), or equivalently the transmissibility exceeds a critical transmissibility [
<xref rid="pone.0142181.ref009" ref-type="bibr">9</xref>
],
<italic>T</italic>
<sub>
<italic>c</italic>
</sub>
, given by
<disp-formula id="pone.0142181.e003">
<alternatives>
<graphic xlink:href="pone.0142181.e003.jpg" id="pone.0142181.e003g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M3">
<mml:mtable displaystyle="true">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</alternatives>
</disp-formula>
</p>
<p>In more general configuration-type networks with triangles and other small-cycle subgraphs,
<xref ref-type="disp-formula" rid="pone.0142181.e001">Eq (1)</xref>
no longer applies. Although the degree distribution will be the same, fewer susceptible nodes will be available for an infection in the early stages of an outbreak, and we would expect the critical transmissibility to be higher, although the extent to which
<xref ref-type="disp-formula" rid="pone.0142181.e001">Eq (1)</xref>
must be modified is the subject of research. Miller [
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
] considers approximations of
<italic>R</italic>
<sub>0</sub>
using additional subgraphs, and suggests clustering affects the final epidemic size only if the typical node degrees are small or clustering is very high. The expression for unclustered networks in
<xref ref-type="disp-formula" rid="pone.0142181.e001">Eq (1)</xref>
also serves as a useful guideline and bound.</p>
<p>The “edge-triangle” model is a generalisation of the configuration model that includes clustering [
<xref rid="pone.0142181.ref010" ref-type="bibr">10</xref>
,
<xref rid="pone.0142181.ref011" ref-type="bibr">11</xref>
]. In short, single edges are counted separately from triangles. Each node
<italic>i</italic>
has
<italic>s</italic>
<sub>
<italic>i</italic>
</sub>
single edges (“stubs”) and
<italic>t</italic>
<sub>
<italic>i</italic>
</sub>
triangle vertices (“corners”). Networks are formed by connecting pairs of stubs at random and triples of corners at random. The joint distribution of stubs and corners is the analogue to the degree distribution in this model. Volz et al. [
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
] derive a deterministic solution for an SIR epidemic process on edge-triangle configuration networks. They also derive a solution for final epidemic size of clustered networks that is asymptotically exact for large network sizes using bond percolation approximation methods. Importantly they also extend their results on epidemic dynamics beyond edge-triangle configuration networks to a further generalisation that includes cliques (i.e., fully connected subgraphs) of size 4 and above to model household structure. Cliques of size
<italic>n</italic>
are included in the network model by analogy with triangles- size
<italic>n</italic>
clique corners are counted in addition to stubs and triangle corners. We note two of their key results. 1) The global clustering coefficient is not sufficient to determine the full epidemiological impact of clustering. Comparing with a different model [
<xref rid="pone.0142181.ref013" ref-type="bibr">13</xref>
] such that the two models are calibrated to have the same degree distribution, the same global clustering coefficient and the same disease parameters, epidemic dynamics can be quite different, especially when clustering is extensive. Thus, the
<italic>specific nature of the clustering</italic>
appears relevant for epidemic dynamics. 2) Clustering slows epidemics and reduces final size, at least for a particular subclass of edge-triangle networks considered (negative binomial degree distribution, tunable fraction of edges in triangles, mean and variance of the degree distribution held constant).</p>
<p>Separately, a simulation study [
<xref rid="pone.0142181.ref014" ref-type="bibr">14</xref>
] has considered the situations when contact network models, not random mixing, are warranted. When the number of contacts per day is large, or the transmissibility of the infection is high (e.g. mumps, measles, pertussis, chickenpox), random mixing models produce adequate results. On the other hand, when numbers of contacts and transmissibility are low (e.g. MRSA, possibly Ebola), network contact structure is important.</p>
<p>A key weakness of the edge-triangle networks [
<xref rid="pone.0142181.ref010" ref-type="bibr">10</xref>
,
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
] is that triangles (and larger cliques) share a node but not an edge (i.e., no overlapping “motifs”, referred to here as subgraphs). A generalisation that allows arbitrary subgraphs beyond cliques has been proposed [
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
] for which network features (and importantly, percolation results for epidemic dynamics) remain solvable. For example, the necessary generating functions are provided for edges, triangles, squares, diamonds (a square with diagonal, also known as a 2-triangle), and 4-cliques. The joint distribution for node “roles” is the key distribution of interest. These five subgraphs require six node roles (i.e., two for the diamond). Analogous to the degree sequence is the “role sequence” which is a vector for each node, counting the number of each role in which that node participates. More general subgraphs could be treated similarly, although the necessary generating function may be complicated or difficult to express. A key question is
<italic>which subgraphs should be included in the network model for the study of epidemics?</italic>
</p>
</sec>
<sec id="sec003">
<title>Exponential Random Graph Models</title>
<p>Exponential Random Graph models (ERGMs) [
<xref rid="pone.0142181.ref004" ref-type="bibr">4</xref>
,
<xref rid="pone.0142181.ref016" ref-type="bibr">16</xref>
] are a particular class of network models that have proven useful in modelling social networks. As a class of network models, ERGMs offer several advantages. 1) They allow general subgraphs, called “configurations”, without additional machinery. 2) Configurations can share edges as well as nodes. 3) Importantly, they are grounded in hypotheses about social processes underlying network formation. 4) They are parameterised to capture network features relevant to human interaction (e.g., clustering, homophily, social circuit dependence [
<xref rid="pone.0142181.ref017" ref-type="bibr">17</xref>
,
<xref rid="pone.0142181.ref018" ref-type="bibr">18</xref>
]) in a manner that aids interpretation. 5) Very parsimonious models can be shown to capture a large number of network features well. Under a homogeneity assumption whereby all structurally identical subgraphs are equally probable, these models have the form
<disp-formula id="pone.0142181.e004">
<alternatives>
<graphic xlink:href="pone.0142181.e004.jpg" id="pone.0142181.e004g" mimetype="image" position="anchor" orientation="portrait"></graphic>
<mml:math id="M4">
<mml:mtable displaystyle="true">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:mrow>
<mml:mtext>Pr</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold">Y</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi mathvariant="bold">y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mo form="prefix">exp</mml:mo>
<mml:mo>(</mml:mo>
<mml:munder>
<mml:mo></mml:mo>
<mml:mi>C</mml:mi>
</mml:munder>
<mml:msub>
<mml:mi>θ</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="bold">y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
<mml:mo>/</mml:mo>
<mml:mi>κ</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</alternatives>
<label>(2)</label>
</disp-formula>
where
<bold>Y</bold>
= [
<italic>Y</italic>
<sub>
<italic>ij</italic>
</sub>
] is an
<italic>n</italic>
×
<italic>n</italic>
binary matrix of variables denoting whether a tie is present (1) or absent (0) and
<bold>y</bold>
= [
<italic>y</italic>
<sub>ij</sub>
] denotes a realisation of
<bold>Y</bold>
. Summation is over all possible subgraph configurations
<italic>C</italic>
,
<italic>θ</italic>
<sub>
<italic>C</italic>
</sub>
is a parameter corresponding to the configuration
<italic>C</italic>
and is non-zero only if all pairs of networks variables in
<italic>C</italic>
are conditionally dependent,
<italic>z</italic>
<sub>
<italic>C</italic>
</sub>
(
<bold>y</bold>
) = ∏
<sub>(
<italic>i</italic>
,
<italic>j</italic>
)∈
<italic>C</italic>
</sub>
<italic>y</italic>
<sub>
<italic>ij</italic>
</sub>
is the network statistic corresponding to the configuration
<italic>C</italic>
(indicating whether all ties in
<italic>C</italic>
are observed in
<bold>y</bold>
), and
<italic>κ</italic>
is a normalizing constant which ensures
<xref ref-type="disp-formula" rid="pone.0142181.e004">Eq (2)</xref>
describes a proper distribution. The homogeneity assumption means there can be one parameter for each configuration, thus reducing the number of sufficient statistics to specify the model.</p>
<p>As has been pointed out (e.g., [
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
]), in practice ERGMs are difficult to analyse and can display pathological behaviours that limit their usefulness. The former concern is certainly true in the general case, although there has been some progress [
<xref rid="pone.0142181.ref019" ref-type="bibr">19</xref>
]. Realistically, the study of both ERGM networks and epidemic dynamics on ERGM networks must rely on stochastic simulation of both the networks and epidemics on the networks. Pathological behaviours, which often make model fitting difficult for the earlier Markov models, are less serious with new social circuit specifications for ERGMs [
<xref rid="pone.0142181.ref016" ref-type="bibr">16</xref>
,
<xref rid="pone.0142181.ref020" ref-type="bibr">20</xref>
]. In addition, several other concerns might be raised. 1) Model fitting currently requires Markov Chain Monte Carlo maximum likelihood estimation (MCMCMLE) techniques that become computationally intensive as the number of nodes and edges grows. This imposes a practical constraint on the size of networks that can be modelled with ERGMs (although the new snowball methods—see below—open possibilities for parameter estimation of large networks). 2) Model fitting has generally assumed the entire network has been observed, which is unrealistic for large networks. However, recent results [
<xref rid="pone.0142181.ref021" ref-type="bibr">21</xref>
] have relaxed this “whole network” requirement such that a “snowball sample” [
<xref rid="pone.0142181.ref022" ref-type="bibr">22</xref>
] can be used. An alternative approach for less than whole-network data has also been proposed [
<xref rid="pone.0142181.ref023" ref-type="bibr">23</xref>
]. 3) Model fitting can be tricky, requiring the “right” collection of configurations to be included to adequately capture network features.</p>
<p>One of the empirical networks used here is actually from a chain-referral sampling technique called “snowball sampling” [
<xref rid="pone.0142181.ref022" ref-type="bibr">22</xref>
]. It is a network of people who co-inject in Melbourne, Australia [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
], an example of a “hidden population” for which snowball sampling (and other chain-referral techniques like “respondent-driven sampling” (RDS) [
<xref rid="pone.0142181.ref025" ref-type="bibr">25</xref>
]) are often used. There are additional difficulties when fitting a network model to a network sample rather than a network census [
<xref rid="pone.0142181.ref026" ref-type="bibr">26</xref>
<xref rid="pone.0142181.ref028" ref-type="bibr">28</xref>
]. 1) High degree nodes are more likely to appear in the sample, so their probability will be over-estimated. 2) The samples are sensitive to the initial seed nodes chosen, which in practice may be a convenience sample rather than a random sample. 3) Homophily within the population can lead to some sub-groups being over- or under-represented in the sample, particularly if the set of seed nodes poorly represents the population. While there are techniques to fit ERGMs to snowball sampled data [
<xref rid="pone.0142181.ref021" ref-type="bibr">21</xref>
], we are not aware of a technique to estimate a network’s node role distribution from a snowball sample. For the purposes of this paper we propose a method to estimate the joint distribution of stubs and corners, thus fitting an edge-triangle model.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="sec004">
<title>Methods</title>
<p>Our basic approach is to start with an empirical contact network and fit edge-triangle type and ERGM type network models from which we simulate a number of networks. We then simulate SIR epidemics on both the empirical network and simulated networks. We compare simulated networks with the empirical network, and SIR results from using the simulated networks and the empirical network. A schematic representation of the study method is shown in
<xref ref-type="fig" rid="pone.0142181.g001">Fig 1</xref>
.</p>
<fig id="pone.0142181.g001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g001</object-id>
<label>Fig 1</label>
<caption>
<title>Schematic representation of the study method.</title>
<p>For each empirical network we fit a number of network models. For each network model we simulate 100 networks and then simulate 1000 Susceptible-Infected-Recovered (SIR) epidemics per simulated network. Simulated networks are compared with the empirical network in terms of network features and simulated epidemic features.</p>
</caption>
<graphic xlink:href="pone.0142181.g001"></graphic>
</fig>
<sec id="sec005">
<title>Empirical Networks</title>
<p>We use three kinds of empirical contact networks which we further describe below. All networks are assumed symmetric.</p>
<sec id="sec006">
<title>High School Contact Networks</title>
<p>The high school data [
<xref rid="pone.0142181.ref029" ref-type="bibr">29</xref>
] describes the cumulative time in close proximity between pairs of people over one day in a U.S. high school. Networks derived from this data provide examples of contact networks relevant to infections where close proximity, but not necessarily physical contact, is important. Such SIR-type infections include corona virus, influenza, norovirus, rhinovirus, varicella and measles. Pertussis (i.e., whooping cough) is an example of a bacterial SIR-type infection (at least in time scales up to years) for which close proximity is relevant.</p>
<p>In total, in this data there are 789 people, comprised of 656 students, 73 teachers, 55 staff, and 5 other individuals. One person has no contacts. By specifying a lower bound
<italic>c</italic>
for contact duration we can define a network in which a tie represents total contact of at least
<italic>c</italic>
time units. For this study we consider networks given by 75, 60 and 6 minute durations (referred to as HS75, HS60 and HS6 respectively), which give a range of densities. Clearly the network ties from larger
<italic>c</italic>
are a subset of the ties from smaller
<italic>c</italic>
. We use these three networks to illustrate the effects of increasing network density and connectivity. Details of the networks are provided in
<xref rid="pone.0142181.t001" ref-type="table">Table 1</xref>
. (To confine comparisons of simulated epidemics to similarly sized networks, isolates have been removed from HS75 and HS60 networks.) Mean node degree increases from 2.5 to 57.5 as shorter duration contacts are included. We also show maximal clique sizes. A maximal clique is a clique that cannot be made larger by including an adjacent node (i.e., not a subset of a larger clique). Maximal clique sizes range from 1 to 4 nodes for HS75, to between 1 and 25 nodes for HS6. Loosely, HS75 is a sparse network with small cliques while HS6 is dense with high mean numbers of contacts and large cliques.</p>
<table-wrap id="pone.0142181.t001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t001</object-id>
<label>Table 1</label>
<caption>
<title>Summary of High School Contact Networks.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t001g" xlink:href="pone.0142181.t001"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Network</th>
<th align="center" rowspan="1" colspan="1">Num. Nodes</th>
<th align="center" rowspan="1" colspan="1">Num. Isolates</th>
<th align="center" rowspan="1" colspan="1">Tie Duration (min.)</th>
<th align="center" rowspan="1" colspan="1">Density</th>
<th align="center" rowspan="1" colspan="1">Mean Node Degree</th>
<th align="center" rowspan="1" colspan="1">Max. Clique Sizes</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">HS75</td>
<td align="center" rowspan="1" colspan="1">538</td>
<td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">75</td>
<td align="char" char="." rowspan="1" colspan="1">0.005</td>
<td align="char" char="." rowspan="1" colspan="1">2.48</td>
<td align="center" rowspan="1" colspan="1">1–4</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">HS60</td>
<td align="center" rowspan="1" colspan="1">688</td>
<td align="center" rowspan="1" colspan="1">0</td>
<td align="center" rowspan="1" colspan="1">60</td>
<td align="char" char="." rowspan="1" colspan="1">0.008</td>
<td align="char" char="." rowspan="1" colspan="1">5.42</td>
<td align="center" rowspan="1" colspan="1">1–6</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">HS6</td>
<td align="center" rowspan="1" colspan="1">789</td>
<td align="center" rowspan="1" colspan="1">5</td>
<td align="center" rowspan="1" colspan="1">6</td>
<td align="char" char="." rowspan="1" colspan="1">0.073</td>
<td align="char" char="." rowspan="1" colspan="1">57.48</td>
<td align="center" rowspan="1" colspan="1">1–25</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t001fn001">
<p>Property summary of empirical networks from high school data with contact durations at least 75 minutes (HS75), 60 minutes (HS60) and 6 minutes (HS6).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec007">
<title>Relationships Network</title>
<p>The relationships network of Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
] is based on data collected from a U.S. high school in 1994, as part of wave 1 of the Add Health study. If a student reported having a special romantic relationship in the last 18 months, they were asked to describe their three most recent relationships (including any current ones) and up to three individuals with whom they had a non-romantic sexual relationship in the previous 18 months. A network tie represents a romantic, or a non-romantic sexual, relationship between students of the high school or between a student of the high school and another of the feeder middle school. Our version of the network, shown in
<xref ref-type="fig" rid="pone.0142181.g002">Fig 2</xref>
, was manually recoded starting from Fig 2 of Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
] so we have an attribute for gender, but not the additional attributes collected in the study. This contact network is relevant for infections where transmission involves intimate or sexual contact. Such SIR-type infections include cytomegalovirus (CMV) and Epstein-Barr virus (EBV).</p>
<fig id="pone.0142181.g002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g002</object-id>
<label>Fig 2</label>
<caption>
<title>High school relationships network based on the relationships network of Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
].</title>
<p>A network tie indicates a romantic or non-romantic sexual relationship was reported by one of the incident nodes. Gender is denoted by node colour (blue-male, pink-female). This version of the network was manually re-coded starting from Fig 2 of Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
].</p>
</caption>
<graphic xlink:href="pone.0142181.g002"></graphic>
</fig>
<p>In the relationships network there are 573 nodes and 477 edges, including one male-male edge and one female-female edge. If one ignored those two same-sex edges the network would be bipartite, although we illustrate more general models here. Notably, the network has only one triangle, one 9-star and three overlapping 4-cycles. In particular, clustering is not a significant feature of the network.</p>
</sec>
<sec id="sec008">
<title>PWID Contact Network</title>
<p>For this study we use the same empirical PWID network described previously [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] and shown in
<xref ref-type="fig" rid="pone.0142181.g003">Fig 3</xref>
. A tie represents two people have engaged in injecting behaviour at the same place and time within three months prior to interview. The population of interest are in three urban neighbourhoods in Melbourne, Australia. A number of personal details are known for each respondent. Location, age (<25 years or >25 years), gender and injecting frequency (less than daily, at least daily) have previously been included in a contact network model [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] for this data. This contact network is relevant for infections involving blood-to-blood contact. Hepatitis B is an example of an infection that can be approximated as SIR-type.</p>
<fig id="pone.0142181.g003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g003</object-id>
<label>Fig 3</label>
<caption>
<title>Snowball sample of PWID contact network used for model estimation.</title>
<p>Waves are indicated by shade, from wave 0 (black) to wave 2 (light gray).</p>
</caption>
<graphic xlink:href="pone.0142181.g003"></graphic>
</fig>
<p>Importantly, although the data collection used network-based methods, a complete network census was not performed. Rather, we can think of the empirical network as a snowball sample, consisting of seed nodes (“wave 0”, “zone 0”) and nodes in waves
<italic>i</italic>
,
<italic>i</italic>
≥ 1. Assuming every node “remembers” to nominate all its network neighbours, nodes in wave
<italic>i</italic>
,
<italic>i</italic>
≥ 1, are all the nodes nominated by nodes in wave
<italic>i</italic>
− 1 that are not in wave
<italic>i</italic>
− 1 themselves. (See Rolls et al. [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] for more details on addressing the effect of missing nominations on the zone structure.)</p>
</sec>
</sec>
<sec id="sec009">
<title>Network Models for the High School Networks</title>
<p>For the high school networks we consider several configuration-type models and an exponential random graph model (ERGM). For comparison purposes we also consider an Erdős-Rényi model [
<xref rid="pone.0142181.ref031" ref-type="bibr">31</xref>
] (sometimes called the
<italic>G</italic>
(
<italic>n</italic>
,
<italic>p</italic>
) model for
<italic>n</italic>
nodes and edge probability
<italic>p</italic>
). For each network model we generate a random sample of 100 networks.</p>
<sec id="sec010">
<title>Configuration-Type Models</title>
<p>We consider six configuration-type models: the standard configuration model (CM), the edge-triangle model (ET), a variation of the edge-triangle model that adds subgraphs of size four (+sub4), a variation that then adds 3-triangles and 4-triangles to the first variation (+34tri), a variation that then adds a subgraph with three connected triangles to the second variation (+truss), and a fourth variation that then adds maximal cliques of size five and above to the third variation (+clqs5+). Some of the subgraphs have been proposed elsewhere [
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
,
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
], others help to mitigate the restriction that modelled triangles cannot share edges, and others are motivated by hypotheses on human tie formation.
<xref ref-type="fig" rid="pone.0142181.g004">Fig 4</xref>
shows the subgraphs, with node colour used to denote the various node roles. The number of node roles is shown in column three. Columns 4–9 show which subgraphs are included in each model variation by the presence of an “X”.</p>
<fig id="pone.0142181.g004" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g004</object-id>
<label>Fig 4</label>
<caption>
<title>Subgraphs for edge-triangle models.</title>
<p>For each subgraph, node roles are distinguished by node colour. The number of node roles for each subgraph is shown in column 3. Inclusion of a subgraph is shown (X) for each of the model variations: standard configuration model (CM), edge-triangle model (ET), variation one with subgraphs of four nodes (+sub4), variation two with 3-triangles and 4-triangles (+34tri), variation three with a “truss”, and variation four with maximal cliques of five or more nodes (+clqs5+).</p>
</caption>
<graphic xlink:href="pone.0142181.g004"></graphic>
</fig>
<p>The standard configuration model [
<xref rid="pone.0142181.ref002" ref-type="bibr">2</xref>
,
<xref rid="pone.0142181.ref003" ref-type="bibr">3</xref>
] and the edge-triangle model [
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
,
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
] have been proposed previously. Our first variation of the edge-triangle model adds the subgraphs of four nodes (square, diamond, 4-clique) are described in [
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
] which also provides the relevant generating functions. The node role sequence is a 6-tuple for each node.</p>
<p>The second variation adds to our first variation the 3-triangles and 4-triangles, which are generalisations of the diamond. They are specific cases of the more general
<italic>k</italic>
-triangles that have proven useful for ERGM modelling [
<xref rid="pone.0142181.ref020" ref-type="bibr">20</xref>
]. Theoretically they can arise in connection with
<italic>social circuit dependence</italic>
[
<xref rid="pone.0142181.ref017" ref-type="bibr">17</xref>
]. Loosely, social circuit dependence captures the idea that people whose contacts are connected are themselves more likely to be connected.</p>
<p>Our third variation adds to our second variation a subgraph we call a “truss”. This is a five node subgraph with three triangles, a maximum of two sharing any given edge. This experimental subgraph is intended to capture additional clustering where several triangles share edges in a manner not captured by the other subgraphs composed of multiple triangles (i.e., the 3-triangle, the 4-triangle, and the 4-clique). Note that since the diamond is a subgraph of the truss, the node roles for truss membership are assigned with priority over the node roles for diamond membership.</p>
<p>Our fourth variation adds to our third variation maximal cliques of size five and above. Cliques are useful for modelling household structure (e.g., [
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
]) where close contact between all individuals is assumed. The general assumption of homogeneous mixing between
<italic>n</italic>
people can be regarded as contact happening in a size
<italic>n</italic>
clique. Of course, as network density increases, large maximal cliques arise as a consequence of the network having many edges. Since the edge-triangle type configuration models do not allow subgraphs to share edges, high density graphs without explicitly modelled cliques would count many edges as stubs, and a typical network from the model would lack much of this structure.</p>
<p>To model cliques, the maximal cliques of each size
<italic>n</italic>
are counted in the empirical network. Associated with each size clique is a single node role. That is, in a maximal clique of size
<italic>n</italic>
, each and every node is connected to
<italic>n</italic>
− 1 other nodes to form the clique. So if the network has only one occurrence of a size
<italic>n</italic>
clique, using the node role sequence the same nodes will appear in a size
<italic>n</italic>
clique in the simulated network. But if there is more than one occurrence of a size
<italic>n</italic>
clique, size
<italic>n</italic>
random samples are taken from the collection of all nodes in size
<italic>n</italic>
cliques.</p>
<p>Earlier work on the edge-triangle [
<xref rid="pone.0142181.ref010" ref-type="bibr">10</xref>
] and subgraph [
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
] models describe how to generate networks starting from the node role sequence (or its distribution). They do not say how to count the subgraphs to create the role sequence or estimate its distribution. To model an empirical network this is most critical. Here we describe our method to form the node role sequence by decomposing the empirical network. For example, to fit an edge-triangle model we first count triangle corners at each node, subject to the condition that two triangles cannot share an edge. Thus, once an edge is identified for a triangle, it is removed from use in future triangles. In this sense we “decompose” the empirical network by removing edges as node roles are assigned. Edges not counted as part of a triangle are counted as stubs at their incidence nodes and account for all remaining edges. This approach can be generalised when additional subgraphs are included in the model, as described in
<xref ref-type="supplementary-material" rid="pone.0142181.s001">S1 Appendix</xref>
for the decomposition for the model with the largest number of subgraphs (i.e., our fourth variation). Finally we note that how one decomposes the empirical network into non-overlapping subgraphs is not unique. In this study the decomposition generally occurs in decreasing size/complexity, starting with the cliques of at least five nodes (most nodes to fewest nodes) and ending with stubs, as shown in
<xref ref-type="supplementary-material" rid="pone.0142181.s001">S1 Appendix</xref>
. Note that node roles for triangles are assigned before node roles for 4-cycles to give the most favourable opportunity to capture clustering. Node roles for the truss are assigned after 3-triangles and 4-cliques, but before diamonds, which are a subgraph of the truss.</p>
<p>To simulate an edge-triangle network, for example, we take the node role sequence and randomly join pairs of stubs and triples of triangle corners as described elsewhere [
<xref rid="pone.0142181.ref010" ref-type="bibr">10</xref>
,
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
]. In particular, we use the observed role sequence instead of fitting a bivariate probability distribution to the number of stubs and corners. Since the quality of that fit would affect our results, using the observed node role sequence provides more controlled comparisons. Joining random stubs and corners can create self-loops and multiple edges. It has been suggested to maintain these ties [
<xref rid="pone.0142181.ref015" ref-type="bibr">15</xref>
] since excluding them makes theoretical calculations more difficult. Asymptotically, as the number of nodes grows the density of self-loops and multiple edges goes to zero, so they are not such a concern for larger networks. As in [
<xref rid="pone.0142181.ref012" ref-type="bibr">12</xref>
], we choose to delete them. Self-loops do not add to epidemic dynamics. Multiple edges do not fit within the framework of a binary incidence matrix, and are inconsistent with decomposing the empirical network into subgraphs that do not share edges.</p>
</sec>
<sec id="sec011">
<title>Exponential Random Graph Models</title>
<p>To specify an ERGM one provides a parameter for each configuration included in the model. When combined in a model, ERGM parameters represent a balance of competing effects whereby positive (negative) values indicate a tendency to increase (decrease) occurrences of that configuration in the model. The general approach to fitting an ERGM is to select a collection of configurations and use MCMCMLE to estimate parameters while MCMC simulation is used to evaluate model fit and simulate from models. All ERGMs for the high school networks in this study contain five configurations for structure and one configuration for homophily based on person type. The
<italic>edge</italic>
configuration controls network density. An “isolates” term controls the number of isolates produced by the model. Where the number of isolates in the empirical network is zero, this parameter helps keep the number of isolates small. Note that in the special case where there are no isolates in an empirical network, more negative values of the isolates parameter would work equally well.</p>
<p>Three additional parameters are based on a parsimonious specification of Snijders et al. [
<xref rid="pone.0142181.ref020" ref-type="bibr">20</xref>
]. These are the
<italic>alternating k-star</italic>
,
<italic>alternating two-path</italic>
and
<italic>alternating k-triangle</italic>
statistics, respectively (see
<xref ref-type="supplementary-material" rid="pone.0142181.s002">S2 Appendix</xref>
for definitions). In an equivalent formulation [
<xref rid="pone.0142181.ref032" ref-type="bibr">32</xref>
], analogous statistics are called the “geometrically weighted degree”, the “geometrically weighted dyadic shared partner” and the “geometrically weighted edgewise shared partner” statistics, respectively. Alternating
<italic>k</italic>
-stars provide a way to model the node degree distribution. When occurring together in a model, alternating
<italic>k</italic>
-triangles model triadic closure and alternating two-paths model the prerequisite for triadic closure. They are also useful for modelling social circuit dependence [
<xref rid="pone.0142181.ref017" ref-type="bibr">17</xref>
,
<xref rid="pone.0142181.ref018" ref-type="bibr">18</xref>
]. While it is not clear
<italic>a priori</italic>
that the simplifying assumptions underlying these statistics should hold for any network, these statistics have proven extremely effective in providing good model fits for many empirical networks.</p>
</sec>
</sec>
<sec id="sec012">
<title>Network Models for the Relationships Network</title>
<p>Because the relationships network has 475 male-female edges and just one male-male and one female-female edge, the male-female bipartite network is an important part of the network structure. We consider a bipartite version of the configuration model that maintains the node degree distribution and gender structure of ties (i.e., male-female, male-male, female-female). We also consider an ERGM and a multilevel ERGM which both explicitly capture the gender structure of ties. For comparison purposes only, we also consider two models that ignore the gender structure of ties: the edge-triangle model and a random network model that randomly rearranges the network ties. For each network model we generate a random sample of 100 networks.</p>
<sec id="sec013">
<title>Configuration-Type Models</title>
<p>For a configuration-type model that considers the gender structure of ties we form the node role sequence for the male-female network and form new networks by taking pairs of stubs with the restriction that each edge is formed by one female and one male stub. This fixes the degree sequence for both the male and female groups, and is analogous to the approach of extending the configuration model to include assortativity by type [
<xref rid="pone.0142181.ref033" ref-type="bibr">33</xref>
]. Note there are no triangles in the bipartite network so triangle corners are not needed. For the purposes of comparing with other models we explicitly put the two same-sex edges into the simulated networks, at the same incident nodes as in the observed network, but acknowledge these two edges will have a negligible impact on the SIR disease simulations.</p>
</sec>
<sec id="sec014">
<title>Exponential Random Graph Models</title>
<p>For the ERGM, instead of using an homophily parameter to force most edges into the bipartite structure, we explicitly fix the presence and absence of same-sex ties so that the bipartite network is the network that is modelled. The latter approach is more computationally efficient and recognises the lack of same-sex ties makes model fitting more difficult. To obtain a model with satisfactory goodness-of-fit, the edges of the 9-star were explicitly fixed also. For this network, the 9-star is not typical of the network and violates the homogeneity assumption that is key to ERGM modelling.</p>
<p>An alternative approach to modelling this network would start with a bipartite ERGM [
<xref rid="pone.0142181.ref034" ref-type="bibr">34</xref>
<xref rid="pone.0142181.ref036" ref-type="bibr">36</xref>
] and explicitly include (or ignore) the two same-sex edges. Instead, we treat the network as an example of a
<italic>multilevel</italic>
network and use a multilevel ERGM [
<xref rid="pone.0142181.ref037" ref-type="bibr">37</xref>
] to model the male-female, female-female and male-male networks together with interactions between the three networks. (Note that the same-sex networks are characterised not just by the presence of one edge each, but also by the
<italic>absence</italic>
of all other possible same-sex edges.) To simplify the modelling we fix the two same-sex edges and model just the male-female bipartite network and any cross-level interactions between the networks. To obtain a model with satisfactory goodness-of-fit, edges of the 9-star, three overlapping 4-cycles and the one triangle were fixed. (Thus, 17 of 475 edges were fixed.) Finally, we flatten the network to create a unipartite network with nodes that have a binary gender attribute.</p>
</sec>
</sec>
<sec id="sec015">
<title>Network Models for the PWID Network</title>
<p>For the PWID network we consider an edge-triangle model, an ERGM, and an Erdős-Rényi model (strictly for comparison). For each network model we generate a random sample of 100 networks. Because the empirical PWID network is really a snowball sample, fitting a network requires some care. Techniques that assume the entire network was observed will be biased. Since a whole network census may be infeasible, fitting to snowball sample data is not an uncommon problem.</p>
<p>Note that when starting from a snowball sample, the size of the population is unknown so this must be determined before networks can be simulated. Previously, an ERGM-based, model-dependent estimate of the size was reported as 524 people [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] and shown to be plausible using independent sources of data. To facilitate comparisons we simulate networks from all three models using 524 nodes.</p>
<sec id="sec016">
<title>Configuration-Type Models</title>
<p>The edge-triangle model is fit to the data using the technique described in
<xref ref-type="supplementary-material" rid="pone.0142181.s003">S3 Appendix</xref>
. In short, we consider six estimators for the node degree distribution, form an estimate for a node’s number of triangle corners conditional on it’s degree, simulate networks and form snowball samples to identify which estimated node degree distribution produces snowball samples most like the observed sample. The choice of an estimator for the node degree distribution is not straightforward and balances effects such as bias from the higher likelihood of including nodes of high degree in a snowball sample, reduced data from excluding particular waves of the snowball, and deviations from simplifying assumptions and approximations.
<xref ref-type="supplementary-material" rid="pone.0142181.s003">S3 Appendix</xref>
includes results of a simulation sub-study to assess the node degree estimators.</p>
</sec>
<sec id="sec017">
<title>Exponential Random Graph Models</title>
<p>Details on fitting an ERGM to this snowball sample network have appeared elsewhere [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] and use a technique called “conditional estimation” [
<xref rid="pone.0142181.ref021" ref-type="bibr">21</xref>
]. Note that for the PWID network, only waves 0—2 were used in model fitting. In addition, for technical reasons related to hypotheses about how people limit their risk of infection, only components of size three and above in the empirical network were used in Rolls et al. [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] but we ignore that restriction in simulated networks here to simplify the comparisons. Fig 6 of Rolls et al. [
<xref rid="pone.0142181.ref024" ref-type="bibr">24</xref>
] shows the largest component from one simulated network using the ERGM.</p>
</sec>
</sec>
<sec id="sec018">
<title>SIR Simulations</title>
<p>All outbreak simulations proceed in the same way and are based on nodes being in one of three states: susceptible, infected, or recovered. Parameters for the simulations are shown in
<xref rid="pone.0142181.t002" ref-type="table">Table 2</xref>
. At time 0 some initial number of seed nodes are explicitly infected. At each time epoch, transmission only occurs along network edges, independently for each infected-susceptible pair, according to a Bernoulli distribution with probability of infection that is fixed for each simulation. Infections are one time unit in duration. Simulations continue until there are no infectious nodes remaining. For each network, 1000 outbreak simulations are performed. As is usual, we choose a minimum outbreak size to define an “epidemic”. Across the three high school networks the minimum outbreak size is about 2.2%, up to rounding. We report the final size, epidemic durations and proportion of simulations yielding epidemics.</p>
<table-wrap id="pone.0142181.t002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t002</object-id>
<label>Table 2</label>
<caption>
<title>Parameters for SIR Simulations.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t002g" xlink:href="pone.0142181.t002"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1"></th>
<th align="center" rowspan="1" colspan="1">HS75</th>
<th align="center" rowspan="1" colspan="1">HS60</th>
<th align="center" rowspan="1" colspan="1">HS6</th>
<th align="center" rowspan="1" colspan="1">Relationships</th>
<th align="center" rowspan="1" colspan="1">PWID</th>
</tr>
<tr>
<th align="center" rowspan="1" colspan="1">Probability of Transmission</th>
<th align="center" colspan="5" rowspan="1">experimental variable</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">Number of seeds</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">10</td>
<td align="center" rowspan="1" colspan="1">10</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Duration of infection</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
<td align="center" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">“Epidemic” min. size (
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
)</td>
<td align="center" rowspan="1" colspan="1">12</td>
<td align="center" rowspan="1" colspan="1">15</td>
<td align="center" rowspan="1" colspan="1">18</td>
<td align="center" rowspan="1" colspan="1">40</td>
<td align="center" rowspan="1" colspan="1">20</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t002fn001">
<p>Parameters for SIR simulations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec019">
<title>Analysis</title>
<p>ERGM parameter estimates, goodness-of-fit and network simulation were performed using PNet [
<xref rid="pone.0142181.ref038" ref-type="bibr">38</xref>
] (high school networks) or MPNet [
<xref rid="pone.0142181.ref035" ref-type="bibr">35</xref>
<xref rid="pone.0142181.ref037" ref-type="bibr">37</xref>
] (relationships network). Estimation and simulation for configuration-model type networks was performed using Matlab [
<xref rid="pone.0142181.ref039" ref-type="bibr">39</xref>
]. Maximal cliques were found using the
<monospace>clique.census</monospace>
command of the SNA package of statnet [
<xref rid="pone.0142181.ref040" ref-type="bibr">40</xref>
]. All SIR simulations and boxplots were completed using Matlab [
<xref rid="pone.0142181.ref039" ref-type="bibr">39</xref>
]. Confidence intervals use a Gaussian approximation.</p>
</sec>
</sec>
<sec sec-type="results" id="sec020">
<title>Results</title>
<sec id="sec021">
<title>High School Contact Networks</title>
<p>Configuration-type models and ERGMs were successfully fit to the HS75 network. Edge-triangle type models here use the node role sequence, which is not parsimonious, so we do not list the counts for the various node roles of each node here. On the other hand, ERGMs often have only a few parameters so the model is easily expressed.
<xref rid="pone.0142181.t003" ref-type="table">Table 3</xref>
shows the ERGM specification for a model of the HS75 network.
<xref rid="pone.0142181.t004" ref-type="table">Table 4</xref>
shows the goodness-of-fit results (burn-in 10
<sup>8</sup>
iterations, 20 000 iterations per network, 4000 networks). Column 1 lists the various statistics that are observed. Column 2 shows the value of the statistic for the empirical network (“Sample”). Columns 3 and 4 show the mean (“Mean”) and sample standard deviation (“Std. Dev.”), respectively, of the statistic over simulated networks. Column 5 shows the “
<italic>t</italic>
-ratio”, defined as (
<italic>X</italic>
<italic>m</italic>
)/
<italic>s</italic>
where
<italic>X</italic>
is the observed value, and
<italic>m</italic>
and
<italic>s</italic>
are the sample mean and sample standard deviation across the simulated networks, respectively. For each parameter, using standard methodology, the model is considered to be a good fit for the empirical network if the
<italic>t</italic>
-ratio has absolute value less than 0.1 for parameters explicitly in the model (shown in boldface) and less than 2 otherwise. Under a Gaussian distribution approximation these roughly correspond to the middle 8% and 95% of a parameter’s distribution, respectively. Column six shows this requirement expressed as an interval (“Goodness-of-Fit Interval”) given by (
<italic>m</italic>
<italic>cs</italic>
,
<italic>m</italic>
+
<italic>cs</italic>
) where
<italic>c</italic>
= 0.1 for statistics explicitly fit in the model and
<italic>c</italic>
= 2 otherwise. For each parameter, the model is a good fit if the observed statistic lies within it’s goodness-of-fit interval. The model itself is a good fit if all the parameters have good fit. The ability of this model specification to closely capture all these features of the empirical network, including ones not explicitly modelled, is strong evidence that a useful model has been specified.</p>
<table-wrap id="pone.0142181.t003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t003</object-id>
<label>Table 3</label>
<caption>
<title>HS75 Network ERGM Specification.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t003g" xlink:href="pone.0142181.t003"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Statistic</th>
<th align="center" rowspan="1" colspan="1">Estimates</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">Edge</td>
<td align="char" char="." rowspan="1" colspan="1">-7.954149</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Isolates</td>
<td align="char" char="." rowspan="1" colspan="1">-10.994985</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-star (
<italic>λ</italic>
= 1.50)</td>
<td align="char" char="." rowspan="1" colspan="1">0.828329</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-triangle (
<italic>λ</italic>
= 3.00)</td>
<td align="char" char="." rowspan="1" colspan="1">1.392038</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-2-path (
<italic>λ</italic>
= 3.00)</td>
<td align="char" char="." rowspan="1" colspan="1">-0.099633</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Homophily: person type</td>
<td align="char" char="." rowspan="1" colspan="1">0.259939</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t003fn001">
<p>ERGM specification for a model of the HS75 network.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="pone.0142181.t004" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t004</object-id>
<label>Table 4</label>
<caption>
<title>HS75 Network ERGM Goodness-of-fit.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t004g" xlink:href="pone.0142181.t004"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Observed Statistic</th>
<th align="center" rowspan="1" colspan="1">Sample</th>
<th align="center" rowspan="1" colspan="1">Mean</th>
<th align="center" rowspan="1" colspan="1">Std. Dev.</th>
<th align="center" rowspan="1" colspan="1">
<italic>t</italic>
-ratio</th>
<th align="center" rowspan="1" colspan="1">Goodness-of-Fit Interval</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Edge</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>666</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>665.497</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>25.841</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.019</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(662.91, 668.08)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">2-star</td>
<td align="center" rowspan="1" colspan="1">1615</td>
<td align="char" char="." rowspan="1" colspan="1">1598.32</td>
<td align="char" char="." rowspan="1" colspan="1">161.321</td>
<td align="char" char="." rowspan="1" colspan="1">0.103</td>
<td align="center" rowspan="1" colspan="1">(1275.68, 1920.96)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">3-star</td>
<td align="center" rowspan="1" colspan="1">1483</td>
<td align="char" char="." rowspan="1" colspan="1">1420.114</td>
<td align="char" char="." rowspan="1" colspan="1">262.421</td>
<td align="char" char="." rowspan="1" colspan="1">0.24</td>
<td align="center" rowspan="1" colspan="1">(895.27, 1944.96)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">4-star</td>
<td align="center" rowspan="1" colspan="1">1093</td>
<td align="char" char="." rowspan="1" colspan="1">990.192</td>
<td align="char" char="." rowspan="1" colspan="1">313.791</td>
<td align="char" char="." rowspan="1" colspan="1">0.328</td>
<td align="center" rowspan="1" colspan="1">(362.61, 1617.77)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">5-star</td>
<td align="center" rowspan="1" colspan="1">649</td>
<td align="char" char="." rowspan="1" colspan="1">566.198</td>
<td align="char" char="." rowspan="1" colspan="1">311.516</td>
<td align="char" char="." rowspan="1" colspan="1">0.266</td>
<td align="center" rowspan="1" colspan="1">(-56.83, 1189.23)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Triangles</td>
<td align="center" rowspan="1" colspan="1">80</td>
<td align="char" char="." rowspan="1" colspan="1">78.913</td>
<td align="char" char="." rowspan="1" colspan="1">14.606</td>
<td align="char" char="." rowspan="1" colspan="1">0.074</td>
<td align="center" rowspan="1" colspan="1">(49.70, 108.13)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Isolates</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>0</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.005</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.069</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>-0.069</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(0.00, 0.01)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">2-Triangle</td>
<td align="center" rowspan="1" colspan="1">65</td>
<td align="char" char="." rowspan="1" colspan="1">55.432</td>
<td align="char" char="." rowspan="1" colspan="1">29.701</td>
<td align="char" char="." rowspan="1" colspan="1">0.322</td>
<td align="center" rowspan="1" colspan="1">(-3.97, 114.83)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Bow tie</td>
<td align="center" rowspan="1" colspan="1">57</td>
<td align="char" char="." rowspan="1" colspan="1">62.922</td>
<td align="char" char="." rowspan="1" colspan="1">42.339</td>
<td align="char" char="." rowspan="1" colspan="1">-0.14</td>
<td align="center" rowspan="1" colspan="1">(-21.76, 147.60)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">4-Cycle</td>
<td align="center" rowspan="1" colspan="1">64</td>
<td align="char" char="." rowspan="1" colspan="1">43.681</td>
<td align="char" char="." rowspan="1" colspan="1">20.232</td>
<td align="char" char="." rowspan="1" colspan="1">1.004</td>
<td align="center" rowspan="1" colspan="1">(3.22, 84.15)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Alt.
<italic>k</italic>
-star (
<italic>λ</italic>
= 1.50)</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>967.347</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>965.926</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>69.811</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.02</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(958.94, 972.91)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Alt.
<italic>k</italic>
-triangle (
<italic>λ</italic>
= 3.00)</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>219.741</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>219.328</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>36.8</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.011</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(215.65, 223.01)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Alt.
<italic>k</italic>
-2-path (
<italic>λ</italic>
= 3.00)</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>1574.074</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>1570.435</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>153.481</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.024</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(1555.09, 1585.78)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">
<bold>Homophily: person type</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>507</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>506.628</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>23.016</bold>
</td>
<td align="char" char="." rowspan="1" colspan="1">
<bold>0.016</bold>
</td>
<td align="center" rowspan="1" colspan="1">
<bold>(504.33, 508.93)</bold>
</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Std. Dev. degree dist.</td>
<td align="center" rowspan="1" colspan="1">1.534</td>
<td align="char" char="." rowspan="1" colspan="1">1.511</td>
<td align="char" char="." rowspan="1" colspan="1">0.082</td>
<td align="char" char="." rowspan="1" colspan="1">0.282</td>
<td align="center" rowspan="1" colspan="1">(1.35, 1.68)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Skewness degree dist.</td>
<td align="center" rowspan="1" colspan="1">1.217</td>
<td align="char" char="." rowspan="1" colspan="1">1.114</td>
<td align="char" char="." rowspan="1" colspan="1">0.141</td>
<td align="char" char="." rowspan="1" colspan="1">0.731</td>
<td align="center" rowspan="1" colspan="1">(0.83, 1.40)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Global Clustering Coefficient</td>
<td align="center" rowspan="1" colspan="1">0.149</td>
<td align="char" char="." rowspan="1" colspan="1">0.147</td>
<td align="char" char="." rowspan="1" colspan="1">0.017</td>
<td align="char" char="." rowspan="1" colspan="1">0.074</td>
<td align="center" rowspan="1" colspan="1">(0.11, 0.18)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Mean Local Clustering</td>
<td align="center" rowspan="1" colspan="1">0.11</td>
<td align="char" char="." rowspan="1" colspan="1">0.1</td>
<td align="char" char="." rowspan="1" colspan="1">0.014</td>
<td align="char" char="." rowspan="1" colspan="1">0.723</td>
<td align="center" rowspan="1" colspan="1">(0.07, 0.13)</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Variance Local Clustering</td>
<td align="center" rowspan="1" colspan="1">0.06</td>
<td align="char" char="." rowspan="1" colspan="1">0.051</td>
<td align="char" char="." rowspan="1" colspan="1">0.008</td>
<td align="char" char="." rowspan="1" colspan="1">1.116</td>
<td align="center" rowspan="1" colspan="1">(0.04, 0.07)</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t004fn001">
<p>Goodness-of-fit result for the HS75 network ERGM. Each row shows a different network feature considered. Statistics explicitly modelled are shown in boldface. Empirical values are shown along with the sample mean (“Mean”) and sample standard deviation (“Std. Dev.”) across 4000 simulations. The sample mean and standard deviation are expressed as a
<italic>t</italic>
-ratio in column five and a goodness-of-fit interval in column six. For all parameters the absolute
<italic>t</italic>
-ratio is always smaller than 0.1 for parameters explicitly modelled, and smaller than 2 otherwise. Equivalently, the observed statistic always lies within it’s goodness-of-fit interval. Based on standard methodology, the network model re-creates all these features of the empirical network to an acceptable standard.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>
<xref ref-type="fig" rid="pone.0142181.g005">Fig 5</xref>
illustrates how well the various models capture various features of the empirical network as boxplots across the simulated networks for each model. As is common in social network analysis, we allow subgraphs to overlap when counting, recognising the idea that a person can be a member of several different groups at once. Boxes show the interquartile range. The central line denotes the median, the whiskers show the range of data not considered outliers, and outliers are shown individually. Where there is no variability the box collapses down to the median. The corresponding value for the empirical network is shown by the horizontal line. Results are reported for an Erdős-Rényi model (“ER”) with the same mean node degree as the empirical network, the standard configuration model (“CM”), an edge-triangle model (“ET”), and the four variations which include size four subgraphs (“+sub4”), 3- and 4-triangles (“+34 tri”), trusses (“+truss”), and cliques of size 5 and above (“+clqs5+”). Each variation listed includes the subgraphs of the models listed before it. Since the largest maximal clique in the HS75 network has only four nodes, the final variation should not show any differences. Finally, we also show results for the ERGM (“ERGM”). (Results for a a random network model that simply randomly rearranges edges [
<xref rid="pone.0142181.ref041" ref-type="bibr">41</xref>
], sometimes called the
<italic>G</italic>
(
<italic>n</italic>
,
<italic>m</italic>
) model for
<italic>n</italic>
nodes and
<italic>m</italic>
edges and referred to hereafter as “ER fix” are indistinguishable from the Erdős-Rényi model and not shown for brevity.)</p>
<fig id="pone.0142181.g005" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g005</object-id>
<label>Fig 5</label>
<caption>
<title>Graph statistics for HS75 network models.</title>
<p>Various network statistics shown as boxplots from 100 simulated networks. AKS is alternating
<italic>k</italic>
-star, GCC is global clustering coefficient, AKT is alternating
<italic>k</italic>
-triangle, and “node max. geo. dist.” is the mean over all nodes in the largest component of the nodewise maximum geodesic distance. Results are reported for an Erdős-Rényi model (“ER”), the configuration model (“CM”), the edge-triangle model (“ET”) and the four variations which include size four subgraphs (“+sub4”), 3- and 4-triangles (“+34 tri”), trusses (“+truss”), and cliques of size 5 and above (“+clqs5+”), respectively, and an ERGM (“ERGM”). Values from the observed network shown by horizontal dotted lines.</p>
</caption>
<graphic xlink:href="pone.0142181.g005"></graphic>
</fig>
<p>The three statistics of the first column (2-star, 3-star, alternating
<italic>k</italic>
-star (AKS)) show how the models are able to capture aspects of the node degree distribution. A different way to assess the agreement of node degree distributions with the empirical network is to use a
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit test with 5% significance. With that test, the degree distributions of 38 simulated ERGM networks are rejected. Similarly, 66 networks from the “ER” model and 78 networks from the “ER fix” model are rejected. On the other hand, the configuration-type models maintain the degree sequence. As expected, none of the simulated networks from any of these models are rejected.</p>
<p>Clustering is an important aspect of social networks. The second and third columns of
<xref ref-type="fig" rid="pone.0142181.g005">Fig 5</xref>
show different aspects of clustering. The second column shows the number of triangles (the key feature of clustering), along with the global clustering coefficient and the alternating
<italic>k</italic>
-triangle statistic. The third column shows graph statistics for 2-triangles (i.e., diamonds), 3-triangles, and 4-triangles.</p>
<p>The three statistics of the fourth column (number of isolates, size of the largest component, diameter of the largest component (i.e., the maximum geodesic length between any two nodes in the largest component) show aspects of the connectivity of the networks. The
<italic>diameter</italic>
of the large component is relevant because when the probability of disease transmission is near one, the disease will propagate one hop per time unit. Thus, the large component diameter is an upper bound on the epidemic duration when the probability of disease transmission is near one and nodes are infected and infectious for one time unit. (If there were more than one seed node per component the epidemic duration would typically be less, but the large component diameter would still be indicative.) There is a large difference between the diameter of the large component in the observed network (31) and all the network models. This diameter arises from a connected core and two chains branching out of the core with 8 and 14 nodes each. None of the network models include two of these long chains, so their diameters are smaller.</p>
<p>We can understand the difference in epidemic duration more precisely by noting that every node has a maximum geodesic distance to some other node(s). For probabilities of transmission near one, the distribution of those maxima will govern the distribution of the epidemic duration. Indeed, if nodes are infected and infectious for one time unit the epidemic duration from a randomly chosen seed node will be its maximum geodesic distance, plus one for the seed node infection. For the HS75 network, the mean of the nodewise maximum geodesic distances in the large component is 22.8. For the 468 nodes in the large component, 70 have their maximum to the end of the 8 node chain, and 385 have their maximum to the end of the 14 node chain, 13 are equidistant to both end nodes and three are additional nodes. Thus, these two long chains, especially the 14 node chain, are responsible for the large nodewise maximum geometric distances. Without these large chains, the networks generated by the network models have smaller nodewise maximum geometric distances in general.</p>
<p>The
<italic>size</italic>
of a network’s large component is relevant because theoretical bond percolation results for epidemic final size arise from a related large component in which the probability of edge formation includes disease transmissibility. A tendency for a model to produce a bigger large component suggests bond percolation results will also predict larger epidemic final size.</p>
<p>
<xref ref-type="fig" rid="pone.0142181.g006">Fig 6</xref>
shows results for the estimated final size in the simulated SIR epidemics. For each probability of transmission and network model, results are reported as the mean and 95% confidence interval (across 100 networks) of the mean final size (across 1000 SIR simulations). The similarity in final size across models when the probability of transmission is near one reflects the similarity in large component size.</p>
<fig id="pone.0142181.g006" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g006</object-id>
<label>Fig 6</label>
<caption>
<title>Final size of SIR epidemics using HS75 network models.</title>
<p>For each network model, shown are the mean (with 95% confidence intervals) over 100 simulated networks of the mean over all outbreaks from 1000 SIR simulations. For the observed network, shown are the mean (with 95% confidence interval) over all outbreaks from 1000 SIR simulations.</p>
</caption>
<graphic xlink:href="pone.0142181.g006"></graphic>
</fig>
<p>
<xref ref-type="fig" rid="pone.0142181.g007">Fig 7</xref>
shows similar results for epidemic duration. Note that large (or missing) confidence intervals for small probabilities of transmission arise from small numbers of observed epidemics. (Recall that an epidemic requires a minimum number of infected nodes (
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
) as shown in
<xref rid="pone.0142181.t002" ref-type="table">Table 2</xref>
.) The results for probability of transmission larger than 0.7 (empirical network notably larger than all network models) reflect similar results in large component diameter. More precisely, for the observed network the epidemic duration when the probability of transmission is one is 23.7 (95% CI: 23.5–23.9), which can be understood as 1 time unit for the seed node and then one for each new infection. (Recall the mean of the nodewise maximum geodesic distances in the large component is 22.8.)</p>
<fig id="pone.0142181.g007" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g007</object-id>
<label>Fig 7</label>
<caption>
<title>Epidemic duration of SIR epidemics using HS75 network models.</title>
<p>For each network model, shown are the mean (with 95% confidence intervals) over 100 simulated networks of the mean over all outbreaks from 1000 SIR simulations. For the observed network, shown are the mean (with 95% confidence interval) over all outbreaks from 1000 SIR simulations. Large (or missing) confidence intervals for small probabilities of transmission arise from a small number of observed outbreaks (typically less than five). The large difference between observed network and the network models for probabilities of transmission larger than 0.7 reflects differences in the large component diameters.</p>
</caption>
<graphic xlink:href="pone.0142181.g007"></graphic>
</fig>
<p>Epidemics can vary in how long they take to get started. Miller [
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
] suggests time shifts to account for these differences before temporal averaging. To further investigate differences across network models, we aligned the start of epidemics by accounting for the variable times to achieve
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
. For each simulation we record the time until just before
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
is achieved, and the durations from that time until the peak incidence of the epidemic and until the epidemic is over. We find that across the network models, the mean times until
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
is achieved are indistinguishable from each other and from results for the observed network (not shown for brevity). Further, graphs of the time from
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
to both the epidemic peak and the end of the epidemic resemble
<xref ref-type="fig" rid="pone.0142181.g007">Fig 7</xref>
but rescaled in the vertical direction (not shown for brevity). Importantly, even after accounting for the variability in the time for an epidemic to get started, epidemics still last longer (on average) on the observed network for probabilities of transmission larger than 0.7. Another quantity of interest, the proportion of simulations with epidemic outbreaks is similar to final size across the network models (not shown for brevity). In particular, it illustrates that for configuration-type networks, outbreaks are less likely to occur when additional clustering is present.</p>
<p>A similar study was performed for the HS60 network which has more clustering.
<xref rid="pone.0142181.t005" ref-type="table">Table 5</xref>
shows the model specification for the ERGM. Goodness-of-fit results show excellent agreement with the empirical network, but are not shown for brevity. In
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests with 5% significance, the degree distribution of 52 simulated ERGM networks are rejected. All 100 networks from the “ER” model are rejected. As expected, no degree distributions from configuration-type networks are rejected.</p>
<table-wrap id="pone.0142181.t005" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t005</object-id>
<label>Table 5</label>
<caption>
<title>HS60 Network ERGM Specification.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t005g" xlink:href="pone.0142181.t005"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Statistic</th>
<th align="center" rowspan="1" colspan="1">Estimates</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">Edge</td>
<td align="char" char="." rowspan="1" colspan="1">-8.62622</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Isolates</td>
<td align="char" char="." rowspan="1" colspan="1">-11.891</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-star (
<italic>λ</italic>
= 1.50)</td>
<td align="char" char="." rowspan="1" colspan="1">1.137897</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-triangle (
<italic>λ</italic>
= 3.00)</td>
<td align="char" char="." rowspan="1" colspan="1">1.250777</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-2-path (
<italic>λ</italic>
= 3.00)</td>
<td align="char" char="." rowspan="1" colspan="1">-0.047</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Homophily: person type</td>
<td align="char" char="." rowspan="1" colspan="1">0.092378</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t005fn001">
<p>ERGM specification for a model of the HS60 network.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>
<xref ref-type="fig" rid="pone.0142181.g008">Fig 8</xref>
shows results for various graph statistics, and are qualitatively similar to those for the HS75 network.
<xref ref-type="fig" rid="pone.0142181.g009">Fig 9</xref>
shows results for the final size from the SIR simulations. In contrast to the HS75 results, there is little difference between the configuration-type models, the ERGM and the empirical network. As with the HS75 network, the random networks (“ER”) overestimate the final size for larger probabilities of transmission. The maximum difference in the mean final size is about 8%.</p>
<fig id="pone.0142181.g008" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g008</object-id>
<label>Fig 8</label>
<caption>
<title>Graph statistics for HS60 network models.</title>
<p>Various network statistics shown as boxplots from 100 simulated networks. AKS is alternating
<italic>k</italic>
-star, GCC is global clustering coefficient, AKT is alternating
<italic>k</italic>
-triangle, and “node max. geo. dist.” is the mean over all nodes in the largest component of the nodewise maximum geodesic distance. Results are reported for an Erdős-Rényi model (“ER”), the configuration model (“CM”), the edge-triangle model (“ET”) and the four variations which include size four subgraphs (“+sub4”), 3- and 4-triangles (“+34 tri”), trusses (“+truss”), and cliques of size 5 and above (“+clqs5+”), respectively, and an ERGM (“ERGM”). Values from the observed network shown by horizontal dotted lines.</p>
</caption>
<graphic xlink:href="pone.0142181.g008"></graphic>
</fig>
<fig id="pone.0142181.g009" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g009</object-id>
<label>Fig 9</label>
<caption>
<title>Final size of SIR epidemics using HS60 network models.</title>
<p>For each network model, shown are the mean (with 95% confidence intervals) over 100 simulated networks of the mean over all outbreaks from 1000 SIR simulations. For the observed network, shown are the mean (with 95% confidence interval) over all outbreaks from 1000 SIR simulations.</p>
</caption>
<graphic xlink:href="pone.0142181.g009"></graphic>
</fig>
<p>In results for the epidemic duration (not shown for brevity) the configuration-type models are indistinguishable, but do have shorter epidemics than the observed network for probabilities of transmission above 0.2 (e.g., 20.1 vs. 23.3 when the probability of transmission is 0.25). Results for the ERGM are similar, but the differences with the observed network are smaller (e.g., 20.7 vs. 23.3 when the probability of transmission is 0.25). Results for the “ER” model are quite different qualitatively, and show the largest differences with the observed network over a range of probabilities of transmission. Conclusions on the time to
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
and the times from
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
to epidemic peak and epidemic end are qualitatively similar to those for the HS75 network, with the exception of the “ER” model which shows notable differences with the other models for probabilities of transmission below 0.5.</p>
<p>The HS6 network is the densest network modelled here, and has mean node degree of 57. Because edge-triangle type models cannot have overlapping subgraphs, we would expect difficulty for these models to capture network features. In
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests with 5% significance, the degree distribution of all simulated edge-triangle type networks is rejected. A key reason is surely the number of multiple edges deleted. While the network has 22675 edges, approximately 1100 to 1200 edges are removed because they are multiple edges. A number of graph statistics show little variation across networks. Because of the large density, the edge-triangle type networks are all composed of five isolates and one large component with all the remaining nodes and there is little variation in this large component.</p>
<p>For the ERGM, which is computationally intensive, the large numbers of nodes and edges also make model fitting difficult. In fact, attempts to fit ERGMs were not successful. While fitting a model with a few configurations (e.g., edge + alternating star (
<italic>λ</italic>
= 1.2) + homophily) was successful in capturing some network features, such models did not sufficiently capture a number of network features not explicitly modelled to a standard common in social network analysis. Fitting models with additional configurations, to improve model fit, appears computationally infeasible. Thus, we cannot report results for an ERGM.
<xref ref-type="fig" rid="pone.0142181.g010">Fig 10</xref>
shows results for various graph statistics for the edge-triangle type models. We do not report alternating statistics; in the absence of ERGMs with acceptable goodness-of-fit there are no meaningful values of
<italic>λ</italic>
<sub>
<italic>i</italic>
</sub>
,
<italic>i</italic>
, = 1, 2, 3.</p>
<fig id="pone.0142181.g010" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g010</object-id>
<label>Fig 10</label>
<caption>
<title>Graph statistics for HS6 network models.</title>
<p>Various network statistics shown as boxplots from 100 simulated networks. GCC is global clustering coefficient, and “node max. geo. dist.” is the mean over all nodes in the largest component of the nodewise maximum geodesic distance. Results are reported for an Erdős-Rényi model (“ER”), the configuration model (“CM”), the edge-triangle model (“ET”) and the four variations which include size four subgraphs (“+sub4”), 3- and 4-triangles (“+34 tri”), trusses (“+truss”), and cliques of size 5 and above (“+clqs5+”), respectively. Values from the observed network shown by horizontal dotted lines.</p>
</caption>
<graphic xlink:href="pone.0142181.g010"></graphic>
</fig>
<p>SIR epidemic results for the HS6 models (not shown for brevity) are qualitatively similar to those for HS60. For example, the final size results are indistinguishable between the configuration-type models, and closely follow those of the empirical network. Results from the random network models (“ER” and “ER fix”) are similar to the empirical network for final sizes between 0% and 65%. Beyond that, these two models overestimate the final size by 5–7%. The epidemic duration results are similar. The maximum difference for the mean epidemic duration between the empirical network and the various edge-triangle models is small (about 2 weeks) when the empirical network predicts 16 weeks duration. On the other hand, between the empirical network and the random networks the difference is about 5 weeks when the empirical network predicts 16 weeks. As with the HS75 and HS60 networks, with the exception of the “ER” model, accounting for the differences in the time to
<italic>R</italic>
<sub>
<italic>min</italic>
</sub>
does not reveal any meaningful differences between the models.</p>
</sec>
<sec id="sec022">
<title>Relationships Network</title>
<p>
<xref rid="pone.0142181.t006" ref-type="table">Table 6</xref>
shows the ERGM specification for a model of the relationships network where the 9-star ties and same-sex network structure are fixed. Goodness-of-fit results show excellent agreement with the empirical network, but are not shown for brevity. The positive alternating
<italic>k</italic>
-star parameter shows a tendency to create stars in the network, and helps capture the degree distribution. The positive isolated edges parameter “IsolateEdges” creates the components of size two. The (relatively) large triangle parameter creates triangles around one of the two same-sex edges. As with the HS75 network, there are no isolates in the empirical network so the isolates parameter serves to keep the number of isolates produced by the model near zero.</p>
<table-wrap id="pone.0142181.t006" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t006</object-id>
<label>Table 6</label>
<caption>
<title>Relationships Network ERGM Specification.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t006g" xlink:href="pone.0142181.t006"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Statistic</th>
<th align="center" rowspan="1" colspan="1">Estimates</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">Edge</td>
<td align="char" char="." rowspan="1" colspan="1">-7.4839</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Triangle</td>
<td align="char" char="." rowspan="1" colspan="1">4.5806</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Isolates</td>
<td align="char" char="." rowspan="1" colspan="1">-11.5927</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">IsolateEdges</td>
<td align="char" char="." rowspan="1" colspan="1">0.175</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Alt.
<italic>k</italic>
-star (
<italic>λ</italic>
= 2)</td>
<td align="char" char="." rowspan="1" colspan="1">0.5882</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t006fn001">
<p>ERGM specification for a model of the relationships network.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>
<xref rid="pone.0142181.t007" ref-type="table">Table 7</xref>
shows the multilevel ERGM specification for a model of the relationships network where the 9-star edges, three overlapping 4-cycle edges, same-sex networks and one triangle are fixed. Of the seven parameters, five are used for modelling the bipartite network. The two alternating
<italic>k</italic>
-star parameters assist in modelling the degree distributions of male and female nodes separately, while the isolates parameters serve to keep the numbers of male and female isolates near zero. The final two parameters control the interactions between the bipartite network and either the female-female network (2-star FFM) or the male-male network (2-star MMF). These positive 2-star parameters show a tendency for a node in a same-sex relationship to also appear in male-female relationships. Goodness-of-fit results are presented in
<xref ref-type="supplementary-material" rid="pone.0142181.s004">S4 Appendix</xref>
and illustrate both the additional degree of modelling detail that is possible, and also the amount of similarity that can be demonstrated, between the empirical and simulated networks when taking a multilevel view.</p>
<table-wrap id="pone.0142181.t007" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.t007</object-id>
<label>Table 7</label>
<caption>
<title>Relationships Network Multilevel ERGM Specification.</title>
</caption>
<alternatives>
<graphic id="pone.0142181.t007g" xlink:href="pone.0142181.t007"></graphic>
<table frame="box" rules="all" border="0">
<colgroup span="1">
<col align="left" valign="middle" span="1"></col>
<col align="left" valign="middle" span="1"></col>
</colgroup>
<thead>
<tr>
<th align="center" rowspan="1" colspan="1">Statistic</th>
<th align="center" rowspan="1" colspan="1">Estimates</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="1" colspan="1">Bipartite Edge</td>
<td align="char" char="." rowspan="1" colspan="1">-7.4549</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Bipartite Isolates (female)</td>
<td align="char" char="." rowspan="1" colspan="1">-11.3222</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Bipartite Isolates (male)</td>
<td align="char" char="." rowspan="1" colspan="1">-10.2588</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Bipartite Alt.
<italic>k</italic>
-star (female;
<italic>λ</italic>
= 2)</td>
<td align="char" char="." rowspan="1" colspan="1">0.9782</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">Bipartite Alt.
<italic>k</italic>
-star (male;
<italic>λ</italic>
= 2)</td>
<td align="char" char="." rowspan="1" colspan="1">0.0703</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">2-Star(FFM)</td>
<td align="char" char="." rowspan="1" colspan="1">0.3454</td>
</tr>
<tr>
<td align="center" rowspan="1" colspan="1">2-Star(MMF)</td>
<td align="char" char="." rowspan="1" colspan="1">1.0056</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t007fn001">
<p>Multilevel ERGM specification for a model of the relationships network. The 2-star effects including same-sex edges show a tendency for a node in a same-sex relationship to also appear in male-female relationships.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>
<xref ref-type="fig" rid="pone.0142181.g011">Fig 11</xref>
shows results for various graph statistics obtained from simulations from the five models: random network with fixed number of edges (“ERfix”), edge-triangle (“ET”), bipartite configuration (“bip CM”), ERGM, and multilevel ERGM (“mERGM”). Because the large number of components is a clear feature of the empirical network, and clustering is not, we focus on quantities related to the degree distribution and connectivity. The numbers of 4-cycles and triangles are also shown for the special role they play in the empirical network to illustrate the ability of the network models to capture these features.</p>
<fig id="pone.0142181.g011" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g011</object-id>
<label>Fig 11</label>
<caption>
<title>Graph statistics for the high school relationships network models.</title>
<p>Various network statistics. AKS is alternating
<italic>k</italic>
-star and “node max. geo. dist.” is the mean over all nodes in the largest component of the nodewise maximum geodesic distance. Results are reported for a random network model with fixed number of edges (“ER fix”), the edge-triangle model (“ET”) without respecting gender, a bipartite version of the configuration model (“bip CM”), an ERGM (“ERGM”) and a multilevel ERGM (“mERGM”). The last three models include the gender of both nodes of each edge, so they have only two same-sex edges. For the ERGM, the 9-star and the same-sex network tie structure are fixed, exogenous effects. For the multilevel ERGM, the 9-star, the sole triangle, the three overlapping 4-cycles and the same-sex network structure are all fixed. All models use 573 nodes.</p>
</caption>
<graphic xlink:href="pone.0142181.g011"></graphic>
</fig>
<p>In general, heterogeneity in the genders of relationship pairs (MM,FF,FM) could play a role in the transmission of some diseases, so capturing the numbers of same-sex edges might be important. For the two models that do not respect the gender structure of edges the number of same-sex edges is large (“ERfix”: 239, 95% CI: 237–241; “bip CM”: 235, 95% CI: 232–237). The other three models all have two same-sex edges. In addition, we again used
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests to check the agreement between the degree distributions of the simulated networks and the empirical network. None of the 100 degree distributions from the two stub models (“ET” and “bip CM”) were rejected at 5% significance while some of the ERGM (24) and multilevel ERGM (37) distributions were rejected. Again, this is despite their seemingly close agreement with the alternating
<italic>k</italic>
-star statistic. On the other hand, while the 2-stars play the largest role in that statistic, there are clearly important differences in the 3-star statistics (and presumably higher order stars too). While the differences may be acceptable in terms of the
<italic>t</italic>
-ratio, those differences can still be sufficient for the
<italic>χ</italic>
<sup>2</sup>
test to reject. Unsurprisingly, all 100 degree distributions from the “ERfix” model were rejected.</p>
<p>
<xref ref-type="fig" rid="pone.0142181.g012">Fig 12</xref>
(left panel) shows results for the final size from simulated SIR epidemics across the five models. To maintain general applicability, the SIR epidemics simulated here do not distinguish between the genders of the two incident nodes, or the direction of transmission. Thus, our SIR results understate the role that nodal attribute heterogeneity might play. So, while results for the edge-triangle and bipartite edge-triangle models are indistinguishable, it is not hard to imagine a more complex disease process for which these two models would give different results. The dotted horizontal line shows the fraction of nodes in the large component of the empirical network. As a result of the large numbers of components, the size of the large component effectively limits the size of an epidemic starting from a limited number of seed nodes. An ordering of the empirical network and network models by decreasing large component size (e.g., “ERfix”, the empirical network, the other four) looks essentially the same as an ordering of the final size when the probability of transmission is one (e.g., “ERfix”, the empirical network, the other four).</p>
<fig id="pone.0142181.g012" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g012</object-id>
<label>Fig 12</label>
<caption>
<title>Final size and duration of SIR epidemics using the high school relationships network models.</title>
<p>For each network model, shown are the mean (with 95% confidence intervals) over 100 simulated networks of the mean over all outbreaks from 1000 SIR simulations. For the observed network, shown are the mean (with 95% confidence interval) over all outbreaks from 1000 SIR simulations. The dotted horizontal line in the left panel shows the fraction of nodes in the large component of the empirical network. The large number of smaller components effectively limits the final size of an SIR epidemic starting from a limited number of seed nodes.</p>
</caption>
<graphic xlink:href="pone.0142181.g012"></graphic>
</fig>
<p>
<xref ref-type="fig" rid="pone.0142181.g012">Fig 12</xref>
(right panel) shows results for the epidemic duration from the simulated SIR epidemics across the five models. An ordering of the empirical network and models by decreasing large component size (e.g., empirical network, “ERGM” and “mERGM”, then down to “ERfix”) looks essentially the same as an ordering of the epidemic duration when the probability of transmission is one (e.g., empirical graph, “mERGM”, “ERGM” then down to “ERfix”).</p>
</sec>
<sec id="sec023">
<title>PWID Contact Network</title>
<p>Three networks models were fit to the snowball sample PWID network: an edge-triangle model (“ET”), an ERGM (“ERGM”) and an Erdős-Rényi model (“ER”). All networks are simulated with 524 nodes. The Erdős-Rényi model uses the edge probability (
<italic>p</italic>
= 0.0042) that gives the same mean node degree as the edge-triangle model and empirical network.</p>
<p>
<xref ref-type="fig" rid="pone.0142181.g013">Fig 13</xref>
shows boxplots for various graph statistics from 100 network simulations from each of the models. Unsurprisingly, the edge-triangle networks show levels of clustering between the Erdős-Rényi networks (fewest triangles, least clustering) and the ERGM networks (most triangles, most clustering).
<xref ref-type="fig" rid="pone.0142181.g014">Fig 14</xref>
(left panel) shows results on the final size of SIR epidemics using the PWID network models while
<xref ref-type="fig" rid="pone.0142181.g014">Fig 14</xref>
(right panel) shows similar results for epidemic duration.</p>
<fig id="pone.0142181.g013" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g013</object-id>
<label>Fig 13</label>
<caption>
<title>Graph statistics for the PWID contact network.</title>
<p>Various network statistics. AKS is alternating
<italic>k</italic>
-star, AKT is alternating
<italic>k</italic>
-triangle, GCC is global clustering coefficient, and “node max. geo. dist.” is the mean over all nodes in the largest component of the nodewise maximum geodesic distance. Results are reported for an Erdős-Rényi model (“ER”), the edge-triangle model (“ET”) and an ERGM (“ERGM”). All models use 524 nodes. For snowball sampled network data there is no complete network with which to compare.</p>
</caption>
<graphic xlink:href="pone.0142181.g013"></graphic>
</fig>
<fig id="pone.0142181.g014" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0142181.g014</object-id>
<label>Fig 14</label>
<caption>
<title>Final size and epidemic duration of SIR epidemics using PWID network models.</title>
<p>For each network model, shown are the mean (with 95% confidence intervals) over 100 simulated networks of the mean over all outbreaks from 1000 SIR simulations. For the observed network, shown are the mean (with 95% confidence interval) over all outbreaks from 1000 SIR simulations. For snowball sampled network data there is no complete network with which to compare.</p>
</caption>
<graphic xlink:href="pone.0142181.g014"></graphic>
</fig>
</sec>
</sec>
<sec sec-type="conclusions" id="sec024">
<title>Discussion</title>
<sec id="sec025">
<title>Degree Distribution</title>
<p>All configuration-type models are able to capture the degree distribution well. Since they are based on the empirical node role sequence, they replicate the degree sequence, except for multiple edges and loops. Looking at the graph statistics for 2-stars, 3-stars and the alternating
<italic>k</italic>
-star statistic, the ERGM seems to provide a good fit to features of the degree distribution. Indeed, a small
<italic>t</italic>
-ratio for the alternating
<italic>k</italic>
-star statistic is often used within social network analysis to argue the network model fits the empirical node degree distribution. Results from the
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests for both HS75 and HS60 networks illustrate how a traditional statistical test can lead to different results. Indeed, for HS60 there is little difference between the confidence intervals for the alternating
<italic>k</italic>
-star statistic. Yet, the
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests rejects none of the networks from the configuration-type models, about half from the ERGM, and all the random networks (“ER”).</p>
<p>The results are somewhat different for the highest density network, for which all the simulated edge-triangle type networks are rejected by
<italic>χ</italic>
<sup>2</sup>
goodness-of-fit tests. A key reason is surely the number of multiple edges deleted. While the empirical network has 22675 edges, approximately 1100 to 1200 edges are removed in simulated edge-triangle networks because they are multiple edges. On the other hand, there is very little difference between the models for “2-stars” and “3-stars”. More will be said about these similarities in connection with both clustering and SIR epidemic results.</p>
</sec>
<sec id="sec026">
<title>Clustering</title>
<p>The number of triangles (the key feature of clustering) illustrates important differences between the models. Unsurprisingly, the random network model (“ER”) produces few triangles. More interesting is that the edge-triangle model also does a poor job of capturing all the triangles. Because of the restriction that triangles in the model cannot share an edge, the edge-triangle model captures only a subset of triangles. Thus, triangles are consistently under-represented in the generated networks in comparison to the empirical network. Further, to capture the node degree sequence, stubs must be used, which makes properties of these networks more closely resemble those of Erdős-Rényi networks.</p>
<p>The addition of size four subgraphs makes a noticeable difference in the HS75 and HS60 networks, but not the HS6 network. This suggests the addition of size four subgraphs are an important addition for capturing network features with edge-triangle type models when both the density is not really large and clustering is not negligible. Additional subgraphs make little difference in the examples considered here. However, household structure is not really a part of the high school network, where cliques might be expected to play an important role. More generally, for the high density HS6 network, as with the star configurations, there are few differences in triangles and clustering between the models. This demonstrates that high density is the main explanation for features in that network.</p>
<p>In our examples, ERGMs do a much better job capturing the number of triangles. The median number of triangles is close to the observed empirical values for both HS75 and HS60 networks, and the interquartile ranges cover the empirical values. (In addition, 95% confidence intervals (not shown for brevity) also cover the empirical values.) Results for the global clustering coefficient are qualitatively similar to those for triangles which illustrates the close connection between those two quantities.</p>
<p>Results for 2-triangles (i.e., diamonds), 3-triangles, 4-triangles and alternating
<italic>k</italic>
-triangles reflect different aspects of clustering. Since the ERGMs here are parameterised to capture these subgraphs it is not surprising they do a better job modelling these features. Since the edge-triangle model has a restriction preventing triangles from sharing edges, it is not surprising how few of these are captured by the model. It is also interesting to note that the addition of 3- and 4-triangles, and even the truss, to the model with size 4 subgraphs makes little change to the network fit.</p>
</sec>
<sec id="sec027">
<title>SIR Epidemic Results</title>
<p>Given these differences between configuration-type models and ERGMs (the former models better capture the degree distribution, the latter models better capture clustering) it is unclear which model should best capture epidemic dynamics. For the HS75, HS60 and relationships networks, final size results for ERGMs show closer agreement than the configuration-type models to those of the empirical network. On the other hand, the difference from using a configuration-type model may be small. Results for the HS75 network illustrate some advantage to including subgraphs of size 4, but no advantage for larger subgraphs. After including subgraphs of size 4, differences with using an ERGM may be considered small enough that the configuration-type models are sufficient. As the density increases, the differences in final size from using those subgraphs becomes negligible. We must emphasize again the special role that cliques can play. None of our networks exhibit household structure, whereby a separate network process creates cliques within households which are then loosely connected to each other. We would expect to observe a larger impact from including cliques in those networks.</p>
<p>By comparison, the random network models (“ER”, “ER fix”) show the worst agreement with results from the empirical network. They are indistinguishable from each other in the examples shown here. When the final size is very high (above 60%), final size is consistently overestimated by these models. Below 60% they may either over- or under-estimate, although they under-estimate when the density is not too high. Because these models are not effective in capturing either the degree distribution or clustering, it remains unclear how to attribute differences in epidemic dynamics to these two features.</p>
<p>In the case of high density, highly clustered networks like HS6, there is little variation between the models. Differences from using the empirical network are generally small. Cliques are undoubtedly present, but because the edge density is so high, cliques arise naturally from the large number of edges. This is a different mechanism for clique formation than the one underlying household structure.</p>
<p>Our results are largely consistent with results of Miller [
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
] which argued in the context of configuration models, that 1) for reasonable networks clustering is only important for final size if node degrees are typically small, and 2) as a consequence, for moderately large transmissibility final size will be unaltered by clustering. Our HS75 network has the smallest node degrees and illustrates some effect of clustering on SIR final size, while HS60 and HS6 have larger node degrees and more clustering but reveal no significant differences in SIR final size across network models with varying amounts of clustering.</p>
<p>Miller [
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
] argues that although “clustering” is usually taken to mean triangles, it could also refer to other short cycles such as squares and 2-triangles that can also impact the spread of infectious diseases. Some of the network models considered here (especially the ERGM) are parameterised to capture such network features. However, even in a network with mean node degree of 2.5 we found only small differences arising from the ability of a network model to capture such features. A much larger difference is apparent from what might be considered the
<italic>absence of clustering</italic>
in our HS75 network due to long chains off a highly-connected core. The results on epidemic duration show similarity between edge-triangle models and the ERGM, and all show mean epidemic durations that are too short compared to the observed network. Long chains can slow an epidemic by increasing the number of transmissions required to reach nodes and might also increase the number of cutpoints in the network. Each cutpoint has the potential to interrupt the epidemic by transmission failure. Overall, this suggests there may be advantage in exploring subgraphs that promote longer chains in addition to ones that promote increased clustering.</p>
<p>A common theme in our results is that results for edge-triangle models are between ERGMs and Erdős-Rényi models. This really is not surprising. Edge-triangle models capture some, but not all, of the number of triangles captured by ERGMs, due to overlapping triangles that must become stubs. In addition, edge-triangle models randomly rewire stubs in a manner that removes structure, not unlike the random edges of Erdős-Rényi networks.</p>
<p>Another observation we can make arises from the fact the edge-triangle models have far fewer 3- and 4-triangles than the ERGM, yet show comparable final size results. Our results suggest that while they play an important role in fitting ERGMs, these subgraphs are not playing an important role in determining epidemic dynamics.</p>
</sec>
<sec id="sec028">
<title>Speed vs. Accuracy</title>
<p>In the networks considered here, SIR epidemic results for ERGMs are generally closer than edge-triangle type models to those for an observed network, although the difference becomes insignificant as density increases. (Of course, all networks become clique-like with high enough density.) Clearly this difference must be balanced by the time spent to develop a model, elapsed time to obtain the model and computational effort to work with a model. In the case of the relationships network for example, in developing the ERGM, 31 models were considered (10 with fixed 9-star and same-sex edges) over the course of a week. Some modeller time was saved by performing estimations overnight. Goodness-of-fit results reported here require 2.5 hours to produce and are an essential step to ensure a useful network model has been found. Twelve earlier goodness-of-fit results (2-3 hours each) ruled out earlier models. In a best case scenario, one could anticipate the need to fix features of the network, and get lucky with the choice of configurations in the model. It is still unlikely to develop the model in under two days. On the other hand, fitting an edge-triangle model can be performed with less specialised expertise and in a negligible amount of time. For the purposes of network epidemiology (not network-based inference which is a separate, excellent domain for ERGMs) these issues must be carefully weighed.</p>
<p>For fitting configuration-type models, the most computationally demanding step is decomposing the empirical network into non-overlapping subgraphs to create the node role sequence. Yet, this is not nearly as demanding as fitting an ERGM. As an example, creating the largest node role sequence for the HS6 network (the largest, densest network used here, with many cliques) required less than eight minutes on a PC. After some effort optimizing the counting algorithms, it likely would perform faster. Because an edge is assigned to only one subgraph, a feature of the decomposition is that once an edge is assigned it can be removed from further consideration. Thus, decomposition of the network becomes less demanding as more edges are assigned to subgraphs. All of this suggests the number of nodes and edges for which the the configuration-type models are useful is at least several orders of magnitude larger than for ERGMs, with the usual provisos on memory and storage limitations for big networks, which are independent of the model.</p>
<p>The decomposition into subgraphs can easily accommodate additional subgraphs. An order for the decomposition is clear when subgraphs are nested- do the larger one first. For subgraphs that have some overlap (e.g., a square and triangle with common edge compared with two squares with common edge) the order for the decomposition is not clear and some kind of choice must be made. It is clear one captures triangle clustering better than the other but that may not be the most important consideration.</p>
<p>There may be an upper-bound on the size of subgraphs that are useful for configuration-type models. As subgraphs increase in size by agglomeration of smaller subgraphs sharing common edges, the number of perimeter edges increases. For example, a 7-node hexagon structure formed from six triangles has six perimeter edges. Inclusion of the hexagon means none of those perimeter edges can be a part of another hexagon, k-triangle or other subgraph. The “cost” for including a hexagon might be to reduce the number of triangles in which some nodes participate.</p>
</sec>
<sec id="sec029">
<title>Snowball Sampled Data</title>
<p>The use of snowball sampled data instead of a whole network poses special challenges. A key concern is a method to fit edge-triangle type models to such samples. We considered six estimators for the degree distribution that might form the basis of technique to fit edge-triangle type models to a snowball sample. In our example the sample degree distribution estimator seems to produce an estimate for which simulated networks more closely resemble the network sample. Clearly this is not enough to form a universal rule, but does suggest the naive estimator should be considered in similar cases.</p>
<p>There are interesting similarities between the PWID and relationships networks. For example, both have low mean node degree. ERGMs for both produce networks with large cycles and smaller branching chains. There is a similarity to the SIR results too. That is, for fixed transmission probabilities we can form the ranking ERGM < edge-triangle < random network, for both final size and epidemic duration. While this proves nothing, it may be a more common feature of similar networks (low mean node degree, large cycle, branching chains).</p>
<p>Results from fitting the PWID contact network illustrate a general problem that with sampled networks there is no complete empirical network to serve as “truth”. This makes understanding the biases of the various modelling approaches even more important.</p>
</sec>
<sec id="sec030">
<title>Nodal Attributes</title>
<p>Exponential random graph modelling provides a rich structure to include the role of nodal attributes in network structure and even model both bipartite and multilevel networks. The role of nodal attributes on network structure (i.e., assortative mixing) and the disease transmission process has been developed for the configuration model [
<xref rid="pone.0142181.ref033" ref-type="bibr">33</xref>
,
<xref rid="pone.0142181.ref042" ref-type="bibr">42</xref>
] but we are not aware of similar results for extensions with additional subgraphs. Such results would be a useful development.</p>
</sec>
<sec id="sec031">
<title>An Interesting Observation on the Relationships Network</title>
<p>We note that Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
] suggests a normative rule is in play to explain the small number of 4-cycles in that network. Further they say that when a 4-cycle term was included in an ERGM, after explicitly accounting for the three 4-cycles occurring through the highly-unusual 9-star, the corresponding 4-cycle parameter was negative. In contrast, we find after explicitly accounting for the 4-cycles and 9-star that small numbers of 4-cycles can be generated in a model (without a need for additional nodal attributes) with acceptable goodness-of-fit that does not explicitly prohibit connected structures larger than triangles. Key differences between the approaches are the use of approximate maximum likelihood estimation and the newer alternating
<italic>k</italic>
-star statistic here. Bearman et al. [
<xref rid="pone.0142181.ref030" ref-type="bibr">30</xref>
] used a pseudolikelihood approach and explicitly included 2-stars and 3-stars, which were common approaches at the time.</p>
</sec>
<sec id="sec032">
<title>Limitations and Future Work</title>
<p>The networks considered here are empirical static networks chosen as illustrative of human contact networks. Our goal here is not to dwell on the particular choices made to create those empirical networks from observational data. We note that Holme [
<xref rid="pone.0142181.ref043" ref-type="bibr">43</xref>
] studies optimal methods for creating static networks for epidemiological studies from empirical temporal contact data. We also acknowledge there are modes of contact (e.g., physical non-intimate contact, contact within a household structure) not represented by our examples. Clearly the study described here can be extended to such empirical networks in the future.</p>
<p>Advances in the use of technology (e.g., Wifi, Bluetooth, RFID tags) to collect high frequency contact network data (e.g., [
<xref rid="pone.0142181.ref029" ref-type="bibr">29</xref>
,
<xref rid="pone.0142181.ref044" ref-type="bibr">44</xref>
<xref rid="pone.0142181.ref047" ref-type="bibr">47</xref>
]) are creating new opportunities for network epidemiology (i.e., “digital epidemiology” [
<xref rid="pone.0142181.ref048" ref-type="bibr">48</xref>
]). Part of the opportunity this provides is the exciting ability to study temporal network effects on disease transmission. On the other hand, it is not clear that temporal, rather than static, models will always be necessary. The comparisons by Machens et al. [
<xref rid="pone.0142181.ref044" ref-type="bibr">44</xref>
] illustrate that a static network model that captures contact durations but not temporal ordering can give SEIR results similar to using full empirical data with a time-varying, individual-based contact network with contact start and end events at 20 s resolution. We adopt the view that the timescales for both epidemic dynamics and contact network dynamics determine the extent to which a static network approximation is useful.</p>
<p>Because our focus is on fitting empirical networks we have not included mathematical network models such as scale-free networks [
<xref rid="pone.0142181.ref049" ref-type="bibr">49</xref>
] and small-world networks [
<xref rid="pone.0142181.ref050" ref-type="bibr">50</xref>
]. These models are not parameterised to simultaneously fit the degree distribution and clustering. In the case of scale-free networks, the degree distribution is key to understanding epidemic dynamics [
<xref rid="pone.0142181.ref051" ref-type="bibr">51</xref>
]. Since the edge-triangle type models can use an empirical node degree sequence, or a node degree distribution, we would expect epidemic dynamics of SIR-type infections would be similar. In the case of small-world networks, the key feature is “shortcuts” that make the typical path lengths shorter than an Erdős-Rényi network. Because the edge-triangle type models randomly re-attach stubs, we would not expect networks from these models to have these shortcuts, and so epidemic dynamics should be more similar to an Erdős-Rényi network.</p>
<p>As a matter of model selection, we would expect better model fit when additional subgraphs are used. We have not attempted to formalize a comparison of these models. We use the empirical node role sequences. If multivariate node role distributions were instead fit to the node role sequences, additional modelling parameters and an additional amount of modelling inaccuracy would be introduced. Some trade-off between model complexity and accuracy would be necessary, but this is beyond the scope of this paper. In many cases researchers will have access to the empirical node roles.</p>
<p>We have focussed on SIR-type infections, but recognise that other types of infections (e.g., SI, SIS) are also important. For the relationships network, most sexually-transmitted diseases would not be modelled as SIR-type. For the PWID network, HIV and hepatitis C are very relevant but would not be modelled as SIR-type. It would be interesting to extend our results to more general infection parameters and types of infections. Edge-based compartmental modelling is a new technique that extends epidemiological results for configuration models to dynamic network models and models where the population has additional heterogeneity (demographic features, multiple types of risk behaviour) and the disease has more complicated history [
<xref rid="pone.0142181.ref033" ref-type="bibr">33</xref>
,
<xref rid="pone.0142181.ref042" ref-type="bibr">42</xref>
]. Unfortunately these methods are not appropriate for SIS-type infections, where the risk of reinfection violates the assumptions underlying simplifications in the calculations. The need for direct simulation, possibly with ERGMs, has not been eliminated for SIS-type infections. Our results could also be extended by looking for differences between the models when a vaccination or treatment is administered. These extensions are left for future work.</p>
<p>We have assumed constant infection duration and both homogeneous infectiousness and susceptibility. Further investigation, in particular to the role of triangles and size four subgraphs in configuration-type models for clustered heterogeneous networks, is left for future work. We note that using theoretical arguments (and demonstrated with a single EpiSimS network) Miller [
<xref rid="pone.0142181.ref008" ref-type="bibr">8</xref>
] states that heterogeneous infectivity and susceptibility allows short-cycle clustering to play a more significant role in outbreak probability and final size for infections close to the epidemic threshold, but less important for networks with sufficiently large average degree. Certainly, as average degree increases we would expect epidemic results to more closely resemble non-network mixing results.</p>
<p>Our approach to fitting an edge-triangle type model to a snowball sample is novel but ad hoc. We leave as future work to develop this approach to modelling within a statistical framework. In addition, it would be interesting to extend this idea to subgraphs of four nodes, which may require sampling additional snowball waves.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="sec033">
<title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0142181.s001">
<label>S1 Appendix</label>
<caption>
<title>Configuration-type Model Node Decomposition.</title>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0142181.s001.pdf">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s002">
<label>S2 Appendix</label>
<caption>
<title>ERGM Alternating Graph Statistics.</title>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0142181.s002.pdf">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s003">
<label>S3 Appendix</label>
<caption>
<title>Fitting an edge-triangle model to snowball sampled data.</title>
<p>Describes a novel method to fit an edge-triangle model to snowball sampled data. Six estimators of the node degree distribution are considered. Results from a simulation sub-study are included.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0142181.s003.pdf">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s004">
<label>S4 Appendix</label>
<caption>
<title>Multilevel ERGM Goodness of Fit.</title>
<p>Tables showing the multilevel ERGM bipartite and cross-level goodness-of-fit results.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0142181.s004.pdf">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s005">
<label>S1 Dataset</label>
<caption>
<title>HS75 network files.</title>
<p>Files related to the HS75 network. Contains the nodes degrees, graph statistics, and SIR simulation results, for both the empirical HS75 network and simulated networks from various models.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s005.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s006">
<label>S2 Dataset</label>
<caption>
<title>HS60 network files.</title>
<p>Files related to the HS60 network. Contains the nodes degrees, graph statistics, and SIR simulation results, for both the empirical HS60 network and simulated networks from various models.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s006.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s007">
<label>S3 Dataset</label>
<caption>
<title>HS6 network files.</title>
<p>Files related to the HS6 network. Contains the nodes degrees, graph statistics, and SIR simulation results, for both the empirical HS6 network and simulated networks from various models.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s007.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s008">
<label>S4 Dataset</label>
<caption>
<title>Relationships network files.</title>
<p>Files related to the Relationships network. Contains the nodes degrees, graph statistics, and SIR simulation results, for both the empirical Relationships network and simulated networks from various models.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s008.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s009">
<label>S5 Dataset</label>
<caption>
<title>PWID network files.</title>
<p>Files related to the PWID network. Contains the nodes degrees, graph statistics, and SIR simulation results, for simulated networks from various models fit to the PWID snowball sample.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s009.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0142181.s010">
<label>S6 Dataset</label>
<caption>
<title>Snowball Sampled Node Degrees for Estimators 1 to 6.</title>
<p>Files for evaluating estimators 1 to 6. Node degrees of the snowball-sampled empirical network and snowball-sampled simulated edge-triangle networks fit using estimators one to six in the paper.</p>
<p>(ZIP)</p>
</caption>
<media xlink:href="pone.0142181.s010.zip">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back>
<ack>
<p>The authors would like to thank Jennifer Badham and Joel Miller for helpful comments on an earlier draft of the manuscript.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0142181.ref001">
<label>1</label>
<mixed-citation publication-type="journal">
<name>
<surname>Welch</surname>
<given-names>D</given-names>
</name>
,
<name>
<surname>Bansal</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Hunter</surname>
<given-names>DR</given-names>
</name>
.
<article-title>Statistical inference to advance network models in epidemiology</article-title>
.
<source>Epidemics</source>
.
<year>2011</year>
;
<volume>3</volume>
(
<issue>1</issue>
):
<fpage>38</fpage>
<lpage>45</lpage>
.
<pub-id pub-id-type="doi">10.1016/j.epidem.2011.01.002</pub-id>
<pub-id pub-id-type="pmid">21420658</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref002">
<label>2</label>
<mixed-citation publication-type="journal">
<name>
<surname>Newman</surname>
<given-names>MEJ</given-names>
</name>
,
<name>
<surname>Strogatz</surname>
<given-names>SH</given-names>
</name>
,
<name>
<surname>Watts</surname>
<given-names>DJ</given-names>
</name>
.
<article-title>Random graphs with arbitrary degree distributions and their applications</article-title>
.
<source>Phys Rev E</source>
.
<year>2001</year>
<month>7</month>
;
<volume>64</volume>
:
<fpage>026118</fpage>
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://link.aps.org/doi/10.1103/PhysRevE.64.026118">http://link.aps.org/doi/10.1103/PhysRevE.64.026118</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1103/PhysRevE.64.026118</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref003">
<label>3</label>
<mixed-citation publication-type="journal">
<name>
<surname>Molloy</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Reed</surname>
<given-names>B</given-names>
</name>
.
<article-title>A critical point for random graphs with a given degree sequence</article-title>
.
<source>Random Structures & Algorithms</source>
.
<year>1995</year>
;
<volume>6</volume>
:
<fpage>161</fpage>
<lpage>179</lpage>
.
<pub-id pub-id-type="doi">10.1002/rsa.3240060204</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref004">
<label>4</label>
<mixed-citation publication-type="journal">
<name>
<surname>Frank</surname>
<given-names>O</given-names>
</name>
,
<name>
<surname>Strauss</surname>
<given-names>D</given-names>
</name>
.
<article-title>Markov graphs</article-title>
.
<source>J Am Stat Assoc</source>
.
<year>1986</year>
;
<volume>81</volume>
(
<issue>395</issue>
):
<fpage>832</fpage>
<lpage>842</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.jstor.org/stable/2289017">http://www.jstor.org/stable/2289017</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1080/01621459.1986.10478342</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref005">
<label>5</label>
<mixed-citation publication-type="journal">
<name>
<surname>Newman</surname>
<given-names>MEJ</given-names>
</name>
.
<article-title>Spread of epidemic disease on networks</article-title>
.
<source>Phys Rev E</source>
.
<year>2002</year>
<month>7</month>
;
<volume>66</volume>
(
<issue>1</issue>
):
<fpage>016128</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevE.66.016128</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref006">
<label>6</label>
<mixed-citation publication-type="journal">
<name>
<surname>Kenah</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>JM</given-names>
</name>
.
<article-title>Second look at the spread of epidemics on networks</article-title>
.
<source>Phys Rev E</source>
.
<year>2007</year>
<month>9</month>
;
<volume>76</volume>
:
<fpage>036113</fpage>
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://link.aps.org/doi/10.1103/PhysRevE.76.036113">http://link.aps.org/doi/10.1103/PhysRevE.76.036113</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1103/PhysRevE.76.036113</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref007">
<label>7</label>
<mixed-citation publication-type="journal">
<name>
<surname>Meyers</surname>
<given-names>LA</given-names>
</name>
.
<article-title>Contact network epidemiology: Bond percolation applied to infectious disease prediction and control</article-title>
.
<source>Bull Am Math Soc</source>
.
<year>2006</year>
;
<volume>44</volume>
:
<fpage>63</fpage>
<lpage>87</lpage>
.
<pub-id pub-id-type="doi">10.1090/S0273-0979-06-01148-7</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref008">
<label>8</label>
<mixed-citation publication-type="journal">
<name>
<surname>Miller</surname>
<given-names>JC</given-names>
</name>
.
<article-title>Spread of infectious disease through clustered populations</article-title>
.
<source>J R Soc Interface</source>
.
<year>2009</year>
;
<volume>6</volume>
:
<fpage>1121</fpage>
<lpage>1134</lpage>
.
<pub-id pub-id-type="doi">10.1098/rsif.2008.0524</pub-id>
<pub-id pub-id-type="pmid">19324673</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref009">
<label>9</label>
<mixed-citation publication-type="journal">
<name>
<surname>Newman</surname>
<given-names>MEJ</given-names>
</name>
.
<article-title>Properties of highly clustered networks</article-title>
.
<source>Phys Rev E</source>
.
<year>2003</year>
;
<volume>68</volume>
(
<issue>2</issue>
):
<fpage>026121</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevE.68.026121</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref010">
<label>10</label>
<mixed-citation publication-type="journal">
<name>
<surname>Newman</surname>
<given-names>MEJ</given-names>
</name>
.
<article-title>Random graphs with clustering</article-title>
.
<source>Phys Rev Lett</source>
.
<year>2009</year>
<month>7</month>
;
<volume>103</volume>
:
<fpage>058701</fpage>
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://link.aps.org/doi/10.1103/PhysRevLett.103.058701">http://link.aps.org/doi/10.1103/PhysRevLett.103.058701</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1103/PhysRevLett.103.058701</pub-id>
<pub-id pub-id-type="pmid">19792540</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref011">
<label>11</label>
<mixed-citation publication-type="journal">
<name>
<surname>Miller</surname>
<given-names>JC</given-names>
</name>
.
<article-title>Percolation and epidemics in random clustered networks</article-title>
.
<source>Phys Rev E</source>
.
<year>2009</year>
<month>8</month>
;
<volume>80</volume>
:
<fpage>020901</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevE.80.020901</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref012">
<label>12</label>
<mixed-citation publication-type="journal">
<name>
<surname>Volz</surname>
<given-names>EM</given-names>
</name>
,
<name>
<surname>Miller</surname>
<given-names>JC</given-names>
</name>
,
<name>
<surname>Galvani</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Ancel Meyers</surname>
<given-names>L</given-names>
</name>
.
<article-title>Effects of heterogeneous and clustered contact patterns on infectious disease dynamics</article-title>
.
<source>PLoS Computational Biology</source>
.
<year>2011</year>
<month>6</month>
;
<volume>7</volume>
(
<issue>6</issue>
).
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002042</pub-id>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002042</pub-id>
<pub-id pub-id-type="pmid">21673864</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref013">
<label>13</label>
<mixed-citation publication-type="journal">
<name>
<surname>House</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Keeling</surname>
<given-names>MJ</given-names>
</name>
.
<article-title>Insights from unifying modern approximations to infections on networks</article-title>
.
<source>J R Soc Interface</source>
.
<year>2011</year>
;
<volume>8</volume>
(
<issue>54</issue>
):
<fpage>67</fpage>
<lpage>73</lpage>
.
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1098/rsif.2010.0179</pub-id>
<pub-id pub-id-type="doi">10.1098/rsif.2010.0179</pub-id>
<pub-id pub-id-type="pmid">20538755</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref014">
<label>14</label>
<mixed-citation publication-type="journal">
<name>
<surname>Smieszek</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Fiebig</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Scholz</surname>
<given-names>R</given-names>
</name>
.
<article-title>Models of epidemics: when contact repetition and clustering should be included</article-title>
.
<source>Theoretical Biology and Medical Modelling</source>
.
<year>2009</year>
<month>6</month>
;
<volume>6</volume>
(
<issue>11</issue>
).</mixed-citation>
</ref>
<ref id="pone.0142181.ref015">
<label>15</label>
<mixed-citation publication-type="journal">
<name>
<surname>Karrer</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Newman</surname>
<given-names>MEJ</given-names>
</name>
.
<article-title>Random graphs containing arbitrary distributions of subgraphs</article-title>
.
<source>Phys Rev E</source>
.
<year>2010</year>
<month>12</month>
;
<volume>82</volume>
:
<fpage>066118</fpage>
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://link.aps.org/doi/10.1103/PhysRevE.82.066118">http://link.aps.org/doi/10.1103/PhysRevE.82.066118</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1103/PhysRevE.82.066118</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref016">
<label>16</label>
<mixed-citation publication-type="journal">
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Kalish</surname>
<given-names>Y</given-names>
</name>
,
<name>
<surname>Lusher</surname>
<given-names>D</given-names>
</name>
.
<article-title>An introduction to exponential random graph (p*) models for social networks</article-title>
.
<source>Social Networks</source>
.
<year>2007</year>
;
<volume>29</volume>
:
<fpage>173</fpage>
<lpage>191</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0142181.ref017">
<label>17</label>
<mixed-citation publication-type="journal">
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Snijders</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Handcock</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
.
<article-title>Recent developments in exponential random graph (p*) models for social networks</article-title>
.
<source>Social Networks</source>
.
<year>2007</year>
;
<volume>29</volume>
(
<issue>2</issue>
):
<fpage>192</fpage>
<lpage>215</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://linkinghub.elsevier.com/retrieve/pii/S0378873306000384">http://linkinghub.elsevier.com/retrieve/pii/S0378873306000384</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2006.08.003</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref018">
<label>18</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
.
<article-title>Neighborhood-based models for social networks</article-title>
.
<source>Sociological Methodology</source>
.
<year>2002</year>
;
<volume>32</volume>
(
<issue>1</issue>
):
<fpage>301</fpage>
<lpage>337</lpage>
.
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1111/1467-9531.00119</pub-id>
<pub-id pub-id-type="doi">10.1111/1467-9531.00119</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref019">
<label>19</label>
<mixed-citation publication-type="other">Chatterjee S, Diaconis P. Estimating and understanding exponential random graph models; 2011. ArXiv:1102.2650v3.</mixed-citation>
</ref>
<ref id="pone.0142181.ref020">
<label>20</label>
<mixed-citation publication-type="journal">
<name>
<surname>Snijders</surname>
<given-names>TAB</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Handcock</surname>
<given-names>MS</given-names>
</name>
.
<article-title>New specifications for exponential random graph models</article-title>
.
<source>Sociol Methodol</source>
.
<year>2006</year>
;
<volume>36</volume>
:
<fpage>99</fpage>
<lpage>153</lpage>
.
<pub-id pub-id-type="doi">10.1111/j.1467-9531.2006.00176.x</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref021">
<label>21</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pattison</surname>
<given-names>PE</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>GL</given-names>
</name>
,
<name>
<surname>Snijders</surname>
<given-names>TAB</given-names>
</name>
,
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
.
<article-title>Conditional estimation of exponential random graph models from snowball sampling designs</article-title>
.
<source>Journal of Mathematical Psychology</source>
.
<year>2013</year>
;
<volume>57</volume>
(
<issue>6</issue>
):
<fpage>284</fpage>
<lpage>296</lpage>
. Social Networks.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0022249613000497">http://www.sciencedirect.com/science/article/pii/S0022249613000497</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.jmp.2013.05.004</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref022">
<label>22</label>
<mixed-citation publication-type="journal">
<name>
<surname>Goodman</surname>
<given-names>LA</given-names>
</name>
.
<article-title>Snowball sampling</article-title>
.
<source>Ann Math Statist</source>
.
<year>1961</year>
;
<volume>32</volume>
:
<fpage>148</fpage>
<lpage>170</lpage>
.
<pub-id pub-id-type="doi">10.1214/aoms/1177705148</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref023">
<label>23</label>
<mixed-citation publication-type="journal">
<name>
<surname>Handcock</surname>
<given-names>MS</given-names>
</name>
,
<name>
<surname>Gile</surname>
<given-names>KJ</given-names>
</name>
.
<article-title>Modeling social networks from sampled data</article-title>
.
<source>Ann Appl Stat</source>
.
<year>2010</year>
;
<volume>4</volume>
(
<issue>1</issue>
):
<fpage>5</fpage>
<lpage>25</lpage>
.
<pub-id pub-id-type="doi">10.1214/08-AOAS221</pub-id>
<pub-id pub-id-type="pmid">26561513</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref024">
<label>24</label>
<mixed-citation publication-type="journal">
<name>
<surname>Rolls</surname>
<given-names>DA</given-names>
</name>
,
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Jenkinson</surname>
<given-names>R</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>PE</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>GL</given-names>
</name>
,
<name>
<surname>Sacks-Davis</surname>
<given-names>R</given-names>
</name>
,
<etal>et al</etal>
<article-title>Modelling a disease-relevant contact network of people who inject drugs</article-title>
.
<source>Social Networks</source>
.
<year>2013</year>
;
<volume>35</volume>
(
<issue>4</issue>
):
<fpage>699</fpage>
<lpage>710</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0378873313000543">http://www.sciencedirect.com/science/article/pii/S0378873313000543</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2013.06.003</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref025">
<label>25</label>
<mixed-citation publication-type="journal">
<name>
<surname>Salganik</surname>
<given-names>MJ</given-names>
</name>
,
<name>
<surname>Heckathorn</surname>
<given-names>DD</given-names>
</name>
.
<article-title>Sampling and estimation in hidden populations using respondent-driven sampling</article-title>
.
<source>Sociological Methodology</source>
.
<year>2004</year>
<month>12</month>
;
<volume>34</volume>
(
<issue>1</issue>
):
<fpage>193</fpage>
<lpage>240</lpage>
.
<pub-id pub-id-type="doi">10.1111/j.0081-1750.2004.00152.x</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref026">
<label>26</label>
<mixed-citation publication-type="journal">
<name>
<surname>Snijders</surname>
<given-names>TAB</given-names>
</name>
.
<article-title>Estimation on the basis of snowball samples: how to weight?</article-title>
<source>Bulletin de Méthodologie Sociologique</source>
.
<year>1992</year>
;
<volume>36</volume>
:
<fpage>59</fpage>
<lpage>70</lpage>
.
<pub-id pub-id-type="doi">10.1177/075910639203600104</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref027">
<label>27</label>
<mixed-citation publication-type="journal">
<name>
<surname>Gile</surname>
<given-names>KJ</given-names>
</name>
,
<name>
<surname>Handcock</surname>
<given-names>MS</given-names>
</name>
.
<article-title>Respondent-driven sampling: an assessment of current methodology</article-title>
.
<source>Sociological Methodology</source>
.
<year>2010</year>
;
<volume>40</volume>
(
<issue>1</issue>
):
<fpage>285</fpage>
<lpage>327</lpage>
.
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1111/j.1467-9531.2010.01223.x</pub-id>
<pub-id pub-id-type="doi">10.1111/j.1467-9531.2010.01223.x</pub-id>
<pub-id pub-id-type="pmid">22969167</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref028">
<label>28</label>
<mixed-citation publication-type="journal">
<name>
<surname>Illenberger</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Flötteröd</surname>
<given-names>G</given-names>
</name>
.
<article-title>Estimating network properties from snowball sampled data</article-title>
.
<source>Social Networks</source>
.
<year>2012</year>
;
<volume>34</volume>
(
<issue>4</issue>
):
<fpage>701</fpage>
<lpage>711</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0378873312000548">http://www.sciencedirect.com/science/article/pii/S0378873312000548</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2012.09.001</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref029">
<label>29</label>
<mixed-citation publication-type="journal">
<name>
<surname>Salathé</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Kazandjieva</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Lee</surname>
<given-names>JW</given-names>
</name>
,
<name>
<surname>Levis</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Feldman</surname>
<given-names>MW</given-names>
</name>
,
<name>
<surname>Jones</surname>
<given-names>JH</given-names>
</name>
.
<article-title>A high-resolution human contact network for infectious disease transmission</article-title>
.
<source>P Natl Acad Sci</source>
.
<year>2010</year>
;
<volume>107</volume>
:
<fpage>22020</fpage>
<lpage>22025</lpage>
.
<pub-id pub-id-type="doi">10.1073/pnas.1009094108</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref030">
<label>30</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bearman</surname>
<given-names>PS</given-names>
</name>
,
<name>
<surname>Moody</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Stovel</surname>
<given-names>K</given-names>
</name>
.
<article-title>Chains of affection: The structure of adolescent romantic and sexual networks</article-title>
.
<source>American Journal of Sociology</source>
.
<year>2004</year>
;
<volume>110</volume>
:
<fpage>44</fpage>
<lpage>91</lpage>
.
<pub-id pub-id-type="doi">10.1086/386272</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref031">
<label>31</label>
<mixed-citation publication-type="journal">
<name>
<surname>Gilbert</surname>
<given-names>EN</given-names>
</name>
.
<article-title>Random graphs</article-title>
.
<source>Ann Math Statist</source>
.
<year>1959</year>
;
<volume>30</volume>
:
<fpage>1141</fpage>
<lpage>1144</lpage>
.
<pub-id pub-id-type="doi">10.1214/aoms/1177706098</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref032">
<label>32</label>
<mixed-citation publication-type="journal">
<name>
<surname>Hunter</surname>
<given-names>DR</given-names>
</name>
.
<article-title>Curved exponential family models for social networks</article-title>
.
<source>Social Networks</source>
.
<year>2007</year>
;
<volume>29</volume>
(
<issue>2</issue>
):
<fpage>216</fpage>
<lpage>230</lpage>
.
<pub-id pub-id-type="doi">10.1016/j.socnet.2006.08.005</pub-id>
<pub-id pub-id-type="pmid">18311321</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref033">
<label>33</label>
<mixed-citation publication-type="journal">
<name>
<surname>Miller</surname>
<given-names>JC</given-names>
</name>
,
<name>
<surname>Volz</surname>
<given-names>EM</given-names>
</name>
.
<article-title>Incorporating disease and population structure into models of SIR disease in contact networks</article-title>
.
<source>PLOS ONE</source>
.
<year>2013</year>
<month>8</month>
;
<volume>8</volume>
(
<issue>8</issue>
):
<fpage>e69162</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0069162</pub-id>
<pub-id pub-id-type="pmid">23990880</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref034">
<label>34</label>
<mixed-citation publication-type="journal">
<name>
<surname>Skvoretz</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Faust</surname>
<given-names>K</given-names>
</name>
.
<article-title>Logit models for affiliation networks</article-title>
.
<source>Sociological Methodology</source>
.
<year>1999</year>
;
<volume>29</volume>
(
<issue>1</issue>
):
<fpage>253</fpage>
<lpage>280</lpage>
.
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1111/0081-1750.00066</pub-id>
<pub-id pub-id-type="doi">10.1111/0081-1750.00066</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref035">
<label>35</label>
<mixed-citation publication-type="journal">
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Sharpe</surname>
<given-names>K</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>GL</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>PE</given-names>
</name>
.
<article-title>Exponential random graph (p*) models for affiliation networks</article-title>
.
<source>Social Networks</source>
.
<year>2009</year>
;
<volume>31</volume>
(
<issue>1</issue>
):
<fpage>12</fpage>
<lpage>25</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0378873308000403">http://www.sciencedirect.com/science/article/pii/S0378873308000403</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2008.08.002</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref036">
<label>36</label>
<mixed-citation publication-type="journal">
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
.
<article-title>Exponential random graph model specifications for bipartite networks- A dependence hierarchy</article-title>
.
<source>Social Networks</source>
.
<year>2013</year>
;
<volume>35</volume>
(
<issue>2</issue>
):
<fpage>211</fpage>
<lpage>222</lpage>
. Special Issue on Advances in Two-mode Social Networks.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0378873311000864">http://www.sciencedirect.com/science/article/pii/S0378873311000864</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2011.12.004</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref037">
<label>37</label>
<mixed-citation publication-type="journal">
<name>
<surname>Wang</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Robins</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Pattison</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Lazega</surname>
<given-names>E</given-names>
</name>
.
<article-title>Exponential random graph models for multilevel networks</article-title>
.
<source>Social Networks</source>
.
<year>2013</year>
;
<volume>35</volume>
(
<issue>1</issue>
):
<fpage>96</fpage>
<lpage>115</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/pii/S0378873313000051">http://www.sciencedirect.com/science/article/pii/S0378873313000051</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.socnet.2013.01.004</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref038">
<label>38</label>
<mixed-citation publication-type="other">Wang P, Robins G, Pattison P. PNet: A Program for the Simulation and Estimation of Exponential Random Graph Models; 2006.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://sna.unimelb.edu.au/PNet">http://sna.unimelb.edu.au/PNet</ext-link>
</comment>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref039">
<label>39</label>
<mixed-citation publication-type="other">Mathworks Inc. MATLAB version 7.14.0 (R2012a). Natick, Massachusetts; 2012.</mixed-citation>
</ref>
<ref id="pone.0142181.ref040">
<label>40</label>
<mixed-citation publication-type="other">Handcock MS, Hunter DR, Butts CT, Goodreau SM, Morris M. statnet: Software tools for the Statistical Modeling of Network Data. Seattle, WA; 2003. Version 2.0.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://statnetproject.org">http://statnetproject.org</ext-link>
</comment>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref041">
<label>41</label>
<mixed-citation publication-type="journal">
<name>
<surname>Erdős</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Rényi</surname>
<given-names>A</given-names>
</name>
.
<article-title>On random graphs, I</article-title>
.
<source>Publicationes Mathematicae (Debrecen)</source>
.
<year>1959</year>
;
<volume>6</volume>
:
<fpage>290</fpage>
<lpage>297</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0142181.ref042">
<label>42</label>
<mixed-citation publication-type="journal">
<name>
<surname>Miller</surname>
<given-names>JC</given-names>
</name>
,
<name>
<surname>Slim</surname>
<given-names>AC</given-names>
</name>
,
<name>
<surname>Volz</surname>
<given-names>EM</given-names>
</name>
.
<article-title>Edge-based compartmental modelling for infectious disease spread</article-title>
.
<source>J R Soc Interface</source>
.
<year>2012</year>
;
<volume>9</volume>
:
<fpage>890</fpage>
<lpage>906</lpage>
.
<pub-id pub-id-type="doi">10.1098/rsif.2011.0403</pub-id>
<pub-id pub-id-type="pmid">21976638</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref043">
<label>43</label>
<mixed-citation publication-type="journal">
<name>
<surname>Holme</surname>
<given-names>P</given-names>
</name>
.
<article-title>Epidemiologically optimal static networks from temporal network Data</article-title>
.
<source>PLOS Comput Biol</source>
.
<year>2013</year>
<month>7</month>
;
<volume>9</volume>
(
<issue>7</issue>
):
<fpage>e1003142</fpage>
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371%2Fjournal.pcbi.1003142">http://dx.doi.org/10.1371%2Fjournal.pcbi.1003142</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1003142</pub-id>
<pub-id pub-id-type="pmid">23874184</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref044">
<label>44</label>
<mixed-citation publication-type="journal">
<name>
<surname>Machens</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Gesualdo</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Rizzo</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Tozzi</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Barrat</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Cattuto</surname>
<given-names>C</given-names>
</name>
.
<article-title>An infectious disease model on empirical networks of human contact: bridging the gap between dynamic network data and contact matrices</article-title>
.
<source>BMC Infectious Diseases</source>
.
<year>2013</year>
;
<volume>13</volume>
(
<issue>1</issue>
):
<fpage>185</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2334-13-185</pub-id>
<pub-id pub-id-type="pmid">23618005</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref045">
<label>45</label>
<mixed-citation publication-type="journal">
<name>
<surname>Isella</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Stehlé</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Barrat</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Cattuto</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Pinton</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Van den Broeck</surname>
<given-names>W</given-names>
</name>
.
<article-title>What’s in a crowd? Analysis of face-to-face behavioral networks</article-title>
.
<source>Journal of Theoretical Biology</source>
.
<year>2011</year>
;
<volume>271</volume>
(
<issue>1</issue>
):
<fpage>166</fpage>
<lpage>180</lpage>
.
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.sciencedirect.com/science/article/B6WMD-51M60KS-2/2/cb31bee32b340b3044c724b88779a60e">http://www.sciencedirect.com/science/article/B6WMD-51M60KS-2/2/cb31bee32b340b3044c724b88779a60e</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1016/j.jtbi.2010.11.033</pub-id>
<pub-id pub-id-type="pmid">21130777</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref046">
<label>46</label>
<mixed-citation publication-type="journal">
<name>
<surname>Stehlé</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Voirin</surname>
<given-names>N</given-names>
</name>
,
<name>
<surname>Barrat</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Cattuto</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Isella</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Pinton</surname>
<given-names>J</given-names>
</name>
,
<etal>et al</etal>
<article-title>High-resolution measurements of face-to-face contact patterns in a primary school</article-title>
.
<source>PLOS ONE</source>
.
<year>2011</year>
<month>8</month>
;
<volume>6</volume>
(
<issue>8</issue>
):
<fpage>e23176</fpage>
<comment>Available from: </comment>
<pub-id pub-id-type="doi">10.1371/journal.pone.0023176</pub-id>
<pub-id pub-id-type="doi">10.1371/journal.pone.0023176</pub-id>
<pub-id pub-id-type="pmid">21858018</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref047">
<label>47</label>
<mixed-citation publication-type="journal">
<name>
<surname>Stehlé</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Voirin</surname>
<given-names>N</given-names>
</name>
,
<name>
<surname>Barrat</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Cattuto</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Colizza</surname>
<given-names>V</given-names>
</name>
,
<name>
<surname>Isella</surname>
<given-names>L</given-names>
</name>
,
<etal>et al</etal>
<article-title>Simulation of an SEIR infectious disease model on the dynamic contact network of conference attendees</article-title>
.
<source>BMC Medicine</source>
.
<year>2011</year>
<month>7</month>
;
<volume>9</volume>
(
<issue>87</issue>
).
<comment>Available from:
<ext-link ext-link-type="uri" xlink:href="http://www.biomedcentral.com/1741-7015/9/87">http://www.biomedcentral.com/1741-7015/9/87</ext-link>
</comment>
<pub-id pub-id-type="doi">10.1186/1741-7015-9-87</pub-id>
<pub-id pub-id-type="pmid">21771290</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref048">
<label>48</label>
<mixed-citation publication-type="journal">
<name>
<surname>Salathé</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Bengtsson</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Bodnar</surname>
<given-names>TJ</given-names>
</name>
,
<name>
<surname>Brewer</surname>
<given-names>DD</given-names>
</name>
,
<name>
<surname>Brownstein</surname>
<given-names>JS</given-names>
</name>
,
<name>
<surname>Buckee</surname>
<given-names>C</given-names>
</name>
,
<etal>et al</etal>
<article-title>Digital epidemiology</article-title>
.
<source>PLOS Comput Biol</source>
.
<year>2012</year>
<month>7</month>
;
<volume>8</volume>
(
<issue>7</issue>
):
<fpage>e1002616</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002616</pub-id>
<pub-id pub-id-type="pmid">22844241</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref049">
<label>49</label>
<mixed-citation publication-type="journal">
<name>
<surname>Barabási</surname>
<given-names>AL</given-names>
</name>
,
<name>
<surname>Albert</surname>
<given-names>R</given-names>
</name>
.
<article-title>Emergence of scaling in random networks</article-title>
.
<source>Science</source>
.
<year>1999</year>
;
<volume>286</volume>
(
<issue>5439</issue>
):
<fpage>509</fpage>
<lpage>512</lpage>
.
<pub-id pub-id-type="doi">10.1126/science.286.5439.509</pub-id>
<pub-id pub-id-type="pmid">10521342</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref050">
<label>50</label>
<mixed-citation publication-type="journal">
<name>
<surname>Watts</surname>
<given-names>DJ</given-names>
</name>
,
<name>
<surname>Strogatz</surname>
<given-names>SH</given-names>
</name>
.
<article-title>Collective dynamics of “small-world” networks</article-title>
.
<source>Nature</source>
.
<year>1998</year>
;
<volume>393</volume>
<pub-id pub-id-type="doi">10.1038/30918</pub-id>
<pub-id pub-id-type="pmid">9623998</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0142181.ref051">
<label>51</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pastor-Satorras</surname>
<given-names>R</given-names>
</name>
,
<name>
<surname>Vespignani</surname>
<given-names>A</given-names>
</name>
.
<article-title>Epidemic spreading in scale-free networks</article-title>
.
<source>Phys Rev Lett</source>
.
<year>2001</year>
<month>4</month>
;
<volume>86</volume>
(
<issue>14</issue>
):
<fpage>3200</fpage>
<lpage>3203</lpage>
.
<pub-id pub-id-type="doi">10.1103/PhysRevLett.86.3200</pub-id>
<pub-id pub-id-type="pmid">11290142</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/CovidV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000984 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000984 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Sante
   |area=    CovidV2
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:4640514
   |texte=   A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:26555701" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a CovidV2 

Wicri

This area was generated with Dilib version V0.6.33.
Data generation: Sat Mar 28 17:51:24 2020. Site generation: Sun Jan 31 15:35:48 2021