The Challenges of Surveying Heavy-tail Distributions for Use in Infectious Disease Dynamics
===========================================================================================

* Michael E. DeWitt
* Nicholas Kortessis
* John W. Sanders
* Candice J. McNeil

## Abstract

Sexual networks often have heavy-tails, where a small number of exceptional individuals in a population have many more sexual partners than the average (e.g., more than five standard deviations). Heavy-tails pose challenges when surveying this group, as these exceptional individuals are uncommon in the population (and so hard to detect), but have disproportionate impact on epidemiological questions, such as those related to the spread of sexually transmitted diseases. In essence, omitting these individuals is a severe error. In this modeling study, we use prior estimates of the distribution of sexual partners amongst men who have sex with men to explore the implication of different sample sizes on survey estimates. We find that even large surveys consistently fail to capture the variance of the sexual network. Surveys of heavy-tailed sexual networks should be designed with this high variance in mind so as not to underestimate the disease dynamics. The failure to adequately capture the variance within a heavy-tailed network has strong implications for infectious disease dynamics and modeling as disease dynamics are often driven by the heavy-tail.

## Introduction

In early 2022, mpox (formerly known as monkeypox) emerged out of traditionally endemic areas of eastern and central Africa and was identified across the world [1]. Monkeypox virus, the virus which causes mpox, was traditionally thought to have lower epidemic potential due to low rates of human-to-human transmission [2]. Human-to-human transmission typically occurs by close and prolonged skin-to-skin contact, in addition to fomites and respiratory droplets [3]. Prior epidemiological analysis and modeling studies had estimated basic reproduction numbers less than 1, indicating low potential for epidemic growth [2, 4]. However, many gay, bisexual, and other men who have sex with men were being diagnosed after reported sexual contact [5], suggesting sustained mpox spread. Mathematical modeling studies were conducted to explain mpox spread despite the thought that the disease should stutter to a stop. These studies indicated that the highly diffuse sexual networks among men who have sex with men may allow for the virus to achieve a basic reproduction number greater than one [6]. Furthermore, Endo and colleagues have referred to their findings as evidence of “heavy-tailed sexual contact networks”, a concept related to probability distributions in which probability density for larger values is higher than that of the exponential distribution. This phenomena relates to the superspreader phenomena which has been observed in other infectious disease processes where a few infections generate a large number of infections due to large variations in individual infectiousness[7, 8]. In the context of sexual networks, there is extremely high variability between reported sexual partners with a very small proportion of individuals having accounting for most sexual partners (Fig. 1). These few, exceptional individuals, greatly increase the population-level basic reproductive number. The outbreak of mpox has shown that understanding the role of these sexual networks is vital in modeling the potential for infectious disease transmission.

![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/07/06/2023.07.05.23292248/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/F1)

Figure 1: Distribution of sexual partners drawn from a truncated Weibull distribution.

The role of sexual networks on the transmission dynamics of infections and sexually transmitted transmitted infections more specifically has been well understood [9–12]. As more contacts are made, there are more opportunities for infected hosts to transmit the pathogen to their contacts. These insights have been integrated into many mathematical models estimating how infections move within sexual networks [13, 14]. However, many of the data used to generate the network models are limited in the number of participants due to the private nature of the data and the vast heterogeneity in sexual practices. As shown by Endo et al., the heavy-tail of the number of partners complicates understanding these unseen networks [6]. Furthermore, the role that sample size plays in appropriately capturing the heterogeneity of these heavy-tailed networks remains largely unexplored. Achieving proper estimates for modeling sexual networks is vital in developing strategies for disease control and outbreak prevention. In this paper we explore the implications of different sample sizes given existing literature on sexual networks.

## Material and Methods

### Data sources

Following the work of Endo et al. we examine the truncated discrete Weibull distributions which were fitted to empirical data from sexual partnerships from the National Surveys of Sexual Attitudes and Lifestyles (NATSAL) from three survey waves: Natsal-2 (1999–2000; 12,110 participants) (26), Natsal-3 (2010–2012; 15,162 participants) [6, 15]. The NATSAL data includes information on the reported number of sexual partners over a given year for those aged 18-44. In this analysis we considered only the distribution of partners amongst men who have sex with other men (MSM) which were estimated from 409 participants in the combined data sets. The Weibull distribution, a continuous probability distribution with support for positive real numbers, has two parameters: a scale parameter *α* and shape parameter *β*. The shape of the distribution is primarily governed by the scale parameter with values less than one resulting in more mass towards zero while scale values greater than one have more mass in the tails. However, the mass in the tails is greater than that of the exponential distribution. When the scale parameter is one, the Weibull distribution reduces to the exponential distribution, a distribution which is commonly used to model the number of contacts in social networks. In order to account for the discrete nature of sexual partners, we passed the continuous estimates generated from the Weibull distribution through a Poisson distribution in order to account for the discrete nature of sexual contacts. We utilized the parameter estimates for MSM sexual networks from Endo et al. in (Table 1). These parameters result in a mean number of annual sexual partners of 9.24 amongst MSM, however, nearly 70% of the population has fewer than 5 partners, but nearly 5% have more than 30 (Fig. 1). The estimates for this distribution are broadly consistent with the findings from the 4,904 participant study by Weiss and colleagues who found a mean number of annual active partners of 8.5 [13].

View this table:
[Table 1:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/T1)

Table 1: Weibull distribution parameters estimates from Endo et al.

### Estimation of survey representativeness

To understand the role of sample size on estimates of the number of sexual partners within an MSM network, we generated a synthetic population of 100,000 participants using a truncated Weibull distribution to represent the population of MSM and their number of sexual partners in a given geography. Similar methods have been described previously [6, 16] to approximate the distribution of contacts. We then sampled this synthetic population with increasing sample sizes and calculated summary statistics. This process was repeated for 1,000 different synthetic populations.

For each sample size, we calculated several summary statistics including the statistical coverage of the true mean number of sexual partners, the mean absolute error (MAE), and mean absolute percent error (MAPE). In order to assess the statistical coverage, the sampled values were fit using a Poisson model and the estimated rate parameter was recorded. The use of the Poisson distribution for contacts approximates an Erdos-Renyi model for the distribution of contacts in a random network [17]. Erdos-Renyi models are often used to parameterize network models for random network connections an have been used extensively in modeling sexual partnerships [13, 18]. To examine the results in the context of epidemiological estimates, we then used the properties of the sample to estimate the basic reproduction number, *R*, given a recover rate, *γ*, of 10 days and a contact rate, *β*, of 0.1%. Assuming that sexual partnerships are formed at random and in proportion to their count, we utilized the methodology described in Keeling and Rohani to estimate the contact rate using the mean, **M**, and variance, **V**, of the number of contacts in the network as shown in equation 1 [19].

![Formula][1]</img> 
When generating synthetic populations, the heavy-tailed distribution was observed in the quantity of individuals who have a greater number of partners(Figure 1).

## Results

We found that poor coverage of the true mean number of partners was observed under all sampling strategies (Figure 2).

![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/07/06/2023.07.05.23292248/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/F2)

Figure 2: Statistical coverage of the true mean number of sexual partners in the synthetic population based on survey sample size.

When examining the estimates, we find that surveys with smaller sample sizes consistently underestimated the true number of sexual partners (Figure 5, panel A). Similarly, because of this underestimation of the mean number of partners as associated smaller estimated variance, the true basic reproduction number was generally underestimated (Figure 5, panel B). In essence, it is much more likely to sample individuals with below average number of sexual contacts than above average sexual contacts, which consistently biases estimates of the mean to be lower than the true value.

When investigating the average expected bias of the estimate across all iterations of a given sample size, we found that there was no clear convergence to a single value for bias (Figure 4, panel A). Similarly when examining the average proportion of the true variance captured compared to the true variance, we find that the variance was consistently underestimated (Figure 4, panel B).

![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/07/06/2023.07.05.23292248/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/F3)

Figure 3: The estimated and true values of the mean number of sexual partners, 95% confidence intervals from a Poisson regression shown (A). The calculated and true contact rate given synthetic epidemiological parameters (B).

![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/07/06/2023.07.05.23292248/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/F4)

Figure 4: The mean value of the bias for the average number of sexual partners based on the sample size (A). The mean proportion of the variance in the number of sexual partners recovered as a function of the sample size (B).

However as the sample size increased both the mean absolute error and the mean absolute percent error began to converge after samples of 2,000 individuals were reached (Figure 3, panels A and B).

## Discussion

In our study of different surveying approaches of heavy-tailed networks, we found poor statistical coverage of the mean number of partners. We found that the average error, both on absolute (MAE) and proportional scales (MAPE) decreased rapidly after 2,000 individuals were sampled. Given the poor statistical coverage and consistent underestimation of the true variance of the network, estimates of the basic reproduction number for contact networks expected to have heavy tails are bound to be biased far below their true value. And the heavier the tail, the larger the bias in estimation of the basic reproductive number specifically and transmission rates more generally. In essence, the individuals most influential on transmission are scarce in the population and are unlikely to be sampled.

These findings suggest that when surveying heavy-tailed social networks, many more individuals need to be sampled than is standard in most approaches assuming distributions closer to normal. As such, care must be taken to sample a sufficient number of the population in order to properly characterize these networks. These findings have strong implications for future field surveillance and modeling efforts.

### Challenges of heavy-tailed distributions

A key finding of this analysis is that because of extreme events, here sexual partners, the mean degree of the network will be poorly measured if insufficient samples of the population are taken. As demonstrated in our modeling study, this is further complicated by the fact that there did not appear to be convergence to a consistent value of bias as the sample size increase. Additionally, because of the large samples required to capture comparatively rare, extreme events, the variance estimates in smaller samples will fail to capture the true, larger variance of the population. As the proportion of the population in the heavy-tail of the distribution increases, the total variance will increase. Thus for the same estimator, more samples will be required in heavy-tailed distributions than distributions with thinner tails. Use of the Poisson distribution to model these effects will result in overconfidence in the uncertainty surrounding the estimates of the rate parameter due to the high variance of the true parameter despite growing sample sizes.

### Infectious disease dynamics

Early studies of the dynamics of sexually transmitted infections in social networks resulted in discovery of reservoirs of infectious diseases due to heterogeneous mixing patterns of partners [20, 21]. Dense, heavy-tailed sexual networks allow for higher persistence of infection due to their scale free nature [22–24]. This is in contrast to so called “small-world” networks which are best displayed as ring lattice where local connections tends to die out because of local pockets of immunity [16]. Additionally, sexual networks have been shown to be highly diffuse with partners traveling across domestic and international boundaries [25, 26]. Providing reliable estimates for these networks is vital in forming interventions to reduce the likelihood of outbreaks and controlling infection. Importantly, the role these networks play extends to understanding the dynamics of multi-drug resistance organisms such as *Neisseria gonorrhoeae* [27] and sexually transmitted *Shigella flexneri* [28]. It is because of the long tails of these distributions that there is a higher potential risk for explosive outbreaks and persistence with the human population.

Better characterization of these networks could be performed to supplement the information regarding the number of partners. Information about demographics, travel history, and other potential features of a given population could be used to supplement partner information to better capture the true network dynamics. For example, knowledge of participation in group sexual activities could allow modelers to better capture information about the long tail of sexual activities, estimate this population size and thus inform not only survey approaches but also better understand how disease may move within the network.

### Implications for modelers

These findings suggest that higher rates of survey participants may be needed in order to more fully capture the dynamics of these heavy-tailed network. Because the infectious disease dynamics of these networks is largely driven by the heavy-tail of the distribution is it vital that care is taken to adequately capture the bias and variance. As we have shown in our simulation, depending on the sample size, the estimated basic reproduction number from a sample for a pathogen could be 5 times lower than the true value (Figure 5, panel A). This has large implications for the anticipated number of cases, potential number of hospitals, and depending on the virulence, the number of deaths. Furthermore, different control strategies might be adopted with information regarding the basic reproduction number [29].

![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/07/06/2023.07.05.23292248/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2023/07/06/2023.07.05.23292248/F5)

Figure 5: The estimated and true values of the mean number of sexual partners, 95% confidence intervals from a Poisson regression shown (A). The calculated and true contact rate given synthetic epidemiological parameters (B).

Furthermore, modeling studies have shown that when risk-based models are used, the binning strategy for risk can have an impact on the inferences [30]. Furthermore, there are likely additional spatial, behavioral, and psycho-social phenomena that may need to be captured in order to capture how these social networks are connected. We did not assess the role of assortative mixing or homophily where those who are more like one another associate more strongly. Studies have shown that outbreak dynamics and effectiveness of different intervention strategies may differ depending on the overall contact structure and the degree of assortative mixing [31, 32]. As more data become available, these demographic, social, and spatial features could be user to better capture the tail behavior in models considering dense sexual networks. In the absence of data, modeler should be sensitive to the fact that existing data may categorically underestimate the true variance in networks and conduct sensitivity analyse to capture the fuller range of possibilities.

## Financial support

MD and CM receive support administered through the US Centers for Disease Control and Prevention’s Epidemiology and Laboratory Capacity for the Prevention and Control of Infectious Diseases Cooperative Agreement (CK19-1904) for Strengthening the US Response to Resistant Gonorrhea.

## Conflicts of interest

All authors have no relevant conflicts of interest to disclose.

## Data availability

All data and code are stored in a publicly available GitHub repository ([https://github](https://github). com/wf-id/heavy-tail-sampling).

## Author contribution

All authors contributed equally to the research design, data collection, formal analysis, interpretation of the results, and writing of the paper.

## Data Availability

All data and code are stored in a publicly available GitHub repository [https://github.com/wf-id/heavy-tail-sampling](https://github.com/wf-id/heavy-tail-sampling).

*   Received July 5, 2023.
*   Revision received July 5, 2023.
*   Accepted July 6, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Gessain, A., Nakoune, E. & Yazdanpanah, Y. Monkeypox. New England Journal of Medicine 387 (ed Hardin, C. C.) 1783–1793. issn: 0028-4793, 1533-4406. [http://www.nejm.org/doi/10.1056/NEJMra2208860](http://www.nejm.org/doi/10.1056/NEJMra2208860) (2023) (Nov. 10, 2022).
    
    

2.  2.Jezek, Z. et al. Four Generations of Probable Person-to-Person Transmission of Human Monkeypox. American Journal of Epidemiology 123, 1004–1012. issn: 0002-9262. pmid: 3010703 (June 1986).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxfordjournals.aje.a114328&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3010703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986C417900006&link_type=ISI) 

3.  3.Organization, W. H. Monkeypox [https://www.who.int/news-room/fact-sheets/detail/monkeypox](https://www.who.int/news-room/fact-sheets/detail/monkeypox) (2022).
    
    

4.  4.Fine, P. E. M., JEZEK, Z., Grab, B. & Dixon, H. The Transmission Potential of Monkeypox Virus in Human Populations. International Journal of Epidemiology 17, 643– 650. issn: 0300-5771. https://doi.org/10.1093/ije/17.3.643<x> (2022) (Sept. 1, 1988).
    
    

5.  5.DeWitt, M. E. et al. Global Monkeypox Case Hospitalisation Rates: A Rapid Systematic Review and Meta-Analysis. eClinicalMedicine 54, 101710. issn: 2589-5370. [https://www.sciencedirect.com/science/article/pii/S2589537022004400](https://www.sciencedirect.com/science/article/pii/S2589537022004400) (2022) (Dec. 1, 2022).
    
    

6.  6.Endo, A. et al. Heavy-Tailed Sexual Contact Networks and Monkeypox Epidemiology in the Global Outbreak, 2022. Science 0, eadd4507. [https://www.science.org/doi/full/10.1126/science.add4507](https://www.science.org/doi/full/10.1126/science.add4507) (2022) (Sept. 22, 2022).
    
    

7.  7.Lloyd-Smith, J. O., Schreiber, S. J., Kopp, P. E. & Getz, W. M. Superspreading and the effect of individual variation on disease emergence. en. Nature 438. Citation Key: lloyd-smithSuperspreadingEffectIndividual2005, 355–359. issn: 1476-4687. [https://www.nature.com/articles/nature04153](https://www.nature.com/articles/nature04153) (Nov. 2005).
    
    

8.  8.Endo, A., for the Mathematical Modelling of Infectious Diseases COVID-19 Working Group, C., Abbott, S., Kucharski, A. J. & Funk, S. Estimating the overdispersion in COVID-19 transmission using outbreak sizes outside China. en. Wellcome Open Research 5. Citation Key: endoEstimatingOverdispersionCOVID192020, 67. issn: 2398-502X. [https://wellcomeopenresearch.org/articles/5-67/v1](https://wellcomeopenresearch.org/articles/5-67/v1) (Apr. 2020).
    
    

9.  9.Keeling, M. J. The Effects of Local Spatial Structure on Epidemiological Invasions. Proceedings of the Royal Society B: Biological Sciences 266, 859–867. issn: 0962-8452. pmid: 10343409. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1689913/<x> (2023) (Apr. 22, 1999).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10343409&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

10. 10.Eames, K. T. D. & Keeling, M. J. Modeling Dynamic and Network Heterogeneities in the Spread of Sexually Transmitted Diseases. Proceedings of the National Academy of Sciences 99, 13330–13335. [https://www.pnas.org/doi/abs/10.1073/pnas](https://www.pnas.org/doi/abs/10.1073/pnas). 202244299 (2022) (Oct. 2002).
    
    

11. 11.Keeling, M. J. & Eames, K. T. D. Networks and Epidemic Models. Journal of the Royal Society, Interface 2, 295–307. issn: 1742-5689. pmid: 16849187 (Sept. 22, 2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2005.0051&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16849187&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234342000003&link_type=ISI) 

12. 12.Read, J. M. & Keeling, M. J. Disease Evolution on Networks: The Role of Contact Structure. Proceedings. Biological Sciences 270, 699–708. issn: 0962-8452. pmid: 12713743 (Apr. 7, 2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2002.2305&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12713743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182249200006&link_type=ISI) 

13. 13.Weiss, K. M. et al. Egocentric Sexual Networks of Men Who Have Sex with Men in the United States: Results from the ARTnet Study. Epidemics 30, 100386. issn: 1755-4365. [https://www.sciencedirect.com/science/article/pii/S1755436519301409](https://www.sciencedirect.com/science/article/pii/S1755436519301409) (2023) (Mar. 1, 2020).
    
    

14. 14.Spicknall, I. H. et al. Modeling the Impact of Sexual Networks in the Transmission of *Monkeypox Virus* Among Gay, Bisexual, and Other Men Who Have Sex with Men — United States, 2022. MMWR. Morbidity and Mortality Weekly Report 71, 1131–1135. issn: 0149-2195, 1545-861X. [http://www.cdc.gov/mmwr/volumes/71/wr/mm7135e2.htm?s\_cid=mm7135e2\_w](http://www.cdc.gov/mmwr/volumes/71/wr/mm7135e2.htm?s_cid=mm7135e2_w) (2023) (Sept. 2, 2022).
    
    

15. 15.Johnson, A., of Hygiene and Tropical Medicine. Centre for Sexual and Reproductive Health Research, L. S., Research, N. S. & Mercer, C. National Survey of Sexual Attitudes and Lifestyles, 2010-2012: Secure Access (May 25, 2017). [https://beta.ukdataservice.ac.uk/datacatalogue/doi/?id=8178](https://beta.ukdataservice.ac.uk/datacatalogue/doi/?id=8178) (2023).
    
    

16. 16.Ferrari, M. J., Bansal, S., Meyers, L. A. & Bjørnstad, O. N. Network Frailty and the Geometry of Herd Immunity. Proceedings of the Royal Society B: Biological Sciences 273, 2743–2748. issn: 0962-8452. pmid: 17015324. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1635496/<x> (2023) (Nov. 7, 2006).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17015324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

17. 17.Bollobás, B. Random Graphs 2nd ed. Cambridge Studies in Advanced Mathematics 73. 498 pp. isbn: 978-0-521-80920-7 978-0-521-79722-1 (Cambridge University Press, Cambridge; New York, 2001).
    
    

18. 18.Datta, S., Mercer, C. H. & Keeling, M. J. Capturing sexual contact patterns in modelling the spread of sexually transmitted infections: Evidence using Natsal-3. en. PLoS ONE 13. Citation Key: dattaCapturingSexualContact2018. [https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6211691/](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6211691/) (2018).
    
    

19. 19.Keeling, M. J. & Rohani, P. Modeling Infectious Diseases in Humans and Animals isbn: 978-0-691-11617-4. [https://press.princeton.edu/books/hardcover/9780691116174/modeling-infectious-diseases-in-humans-and-animals](https://press.princeton.edu/books/hardcover/9780691116174/modeling-infectious-diseases-in-humans-and-animals) (2022) (Princeton University Press, Oct. 28, 2007).
    
    

20. 20.Klovdahl, A. S. et al. Social Networks and Infectious Disease: The Colorado Springs Study. Social Science & Medicine (1982) 38, 79–88. issn: 0277-9536. pmid: 8146718 (Jan. 1994).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0277-9536(94)90302-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8146718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994MK27000010&link_type=ISI) 

21. 21.Stigum, H., Falck, W. & Magnus, P. The Core Group Revisited: The Effect of Partner Mixing and Migration on the Spread of Gonorrhea, Chlamydia, and HIV. Mathematical Biosciences 120, 1–23. issn: 0025-5564. https://www.sciencedirect.com/science/article/pii/0025556494900361<x> (2023) (Mar. 1, 1994).
    
    

22. 22.Pastor-Satorras, R. & Vespignani, A. Epidemic Spreading in Scale-Free Networks. Physical Review Letters 86, 3200–3203. [https://link.aps.org/doi/10.1103/PhysRevLett.86.3200](https://link.aps.org/doi/10.1103/PhysRevLett.86.3200) (2023) (Apr. 2, 2001).
    
    

23. 23.Schneeberger, A. et al. Scale-Free Networks and Sexually Transmitted Diseases: A Description of Observed Patterns of Sexual Contacts in Britain and Zimbabwe. Sexually Transmitted Diseases 31, 380. issn: 0148-5717. https://journals.lww.com/stdjournal/Fulltext/2004/06000/Scale\_Free\_Networks\_and\_Sexually_Transmitted.12.aspx<x> (2023) (June 2004).
    
    

24. 24.Ferguson, N. M. & Garnett, G. P. More Realistic Models of Sexually Transmitted Disease Transmission Dynamics: Sexual Partnership Networks, Pair Models, and Moment Closure. en-US. Sexually Transmitted Diseases 27. Citation Key: fergusonMoreRealisticModels2000, 600. issn: 0148-5717. [https://journals.lww.com/stdjournal/Fulltext/2000/11000/More\_Realistic\_Models\_of\_Sexually\_Transmitted.8.aspx](https://journals.lww.com/stdjournal/Fulltext/2000/11000/More\_Realistic\_Models_of_Sexually_Transmitted.8.aspx) (Nov. 2000).
    
    

25. 25.Ibanhes, F. H. et al. Sexuality Transmitted Infections Risk and Interstate Sexual Networks among Men Who Have Sex with Men in New England. Sexually Transmitted Diseases Publish Ahead of Print. issn: 1537-4521, 0148-5717. [https://journals.lww.com/10.1097/OLQ.0000000000001823](https://journals.lww.com/10.1097/OLQ.0000000000001823) (2023) (Apr. 20, 2023).
    
    

26. 26.Endo, A., Jung, S.-m. & Miura, F. Mpox Emergence in Japan: Ongoing Risk of Establishment in Asia. The Lancet 401, 1923–1924. issn: 0140-6736, 1474-547X. pmid: 37269865. [https://www.thelancet.com/journals/lancet/article/PIIS0140-6736(23)00766-3/fulltext?dgcid=raven\_jbs\_etoc\_email](https://www.thelancet.com/journals/lancet/article/PIIS0140-6736(23)00766-3/fulltext?dgcid=raven_jbs_etoc_email) (2023) (June 10, 2023).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37269865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

27. 27.Gieseker, K. E. et al. Demographic and Epidemiological Characteristics Associated With Reduced Antimicrobial Susceptibility to Neisseria Gonorrhoeae in the United States, Strengthening the US Response to Resistant Gonorrhea, 2018 to 2019. Sexually Transmitted Diseases 48, S118–S123. issn: 1537-4521. pmid: 34433798 (12S Suppl 2 Dec. 1, 2021).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34433798&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

28. 28.Thorley, K. et al. Emergence of extensively drug-resistant and multidrug-resistant Shigella flexneri serotype 2a associated with sexual transmission among gay, bisexual, and other men who have sex with men, in England: a descriptive epidemiological study. en. The Lancet Infectious Diseases 23. Citation Key: thorleyEmergenceExtensivelyDrugresistant2023a, 732–739. issn: 1473-3099. [https://www.sciencedirect.com/science/article/pii/S1473309922008076](https://www.sciencedirect.com/science/article/pii/S1473309922008076) (June 2023).
    
    

29. 29.Fraser, C., Riley, S., Anderson, R. M. & Ferguson, N. M. Factors That Make an Infectious Disease Outbreak Controllable. Proceedings of the National Academy of Sciences 101, 6146–6151. issn: 0027-8424, 1091-6490. pmid: 15071187. https://www.pnas.org/content/101/16/6146<x> (2020) (Apr. 20, 2004).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15071187&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

30. 30.Suen, S.-c., Goldhaber-Fiebert, J. D. & Brandeau, M. L. Risk Stratification in Compartmental Epidemic Models: Where to Draw the Line? Journal of theoretical biology 428, 1–17. issn: 0022-5193. pmid: 28606751. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5555315/<x> (2023) (Sept. 7, 2017).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28606751&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F06%2F2023.07.05.23292248.atom) 

31. 31.Gupta, S., Anderson, R. M. & May, R. M. Networks of sexual contacts: implications for the pattern of spread of HIV. eng. *AIDS (London*, England) 3. Citation Key: guptaNetworksSexualContacts1989, 807–817. issn: 0269-9370 (Dec. 1989).
    
    

32. 32.Kiss, I. Z., Green, D. M. & Kao, R. R. The effect of network mixing patterns on epidemic dynamics and the efficacy of disease contact tracing. Journal of the Royal Society Interface 5. Citation Key: kissEffectNetworkMixing2008a, 791–799. issn: 1742-5689. [https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2386895/](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2386895/) (July 2008).

 [1]: /embed/graphic-3.gif