Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Perinatol. Author manuscript; available in PMC 2017 September 23.
Published in final edited form as:
PMCID: PMC5446293

Network Analysis: A Novel Method for Mapping Neonatal Acute Transport Patterns in California

Sarah N. Kunz, MD MPH,1,2 John A. F. Zupancic, MD ScD,1,2 Joseph Rigdon, PhD,3 Ciaran S. Phibbs, PhD,4,5 Henry C. Lee, MD MS,4,6 Jeffrey B. Gould, MD MPH,4,6 Jure Leskovec, PhD,7 and Jochen Profit, MD MPH4,6



To use network analysis to describe the pattern of neonatal transfers in California, to compare empirical sub-networks with established referral regions, and to determine factors associated with transport outside the originating sub-network.

Study Design

This cross-sectional database study included 6546 infants <28 days old transported within California in 2012. After generating a graph representing acute transfers between hospitals (n=6696), we used community detection techniques to identify more tightly connected sub-networks. These empirically-derived sub-networks were compared to state-defined regional referral networks. Reasons for transfer between empirical sub-networks were assessed using logistic regression.


Empirical sub-networks showed significant overlap with regulatory regions (p <0.001). Transfer outside the empirical sub-network was associated with major congenital anomalies (p<0.001), need for surgery (p=0.01), and insurance as the reason for transfer (p<0.001).


Network analysis accurately reflected empirical neonatal transfer patterns, potentially facilitating quantitative, rather than qualitative, analysis of regionalized health care delivery systems.


Admission rates to neonatal intensive care units (NICUs) in the United States have been increasing.1 Care and outcomes for these newborns vary widely, which is not fully explained by differences in underlying clinical risk.2 Differences in how neonatal care is organized across regionalized health care delivery networks may account for some of this variation.3-8

In order to optimally match patient need with neonatal care delivery capabilities, a program of regionalized care was implemented in the United States, beginning in the 1970s.9 Regionalization involves a linking of hospitals in a coordinated system of communication, learning and response. Over the past several decades, these networks have been shaped by regulatory policy, financial incentives, and market forces.4,5,8 For example, the Affordable Care Act has led to a wave of consolidation of care that is expected to have a significant effect on the shape of regionalized care delivery for newborns.10

A rich literature has documented the benefits of maternal transport and the adverse effects of acute inter-hospital transport on neonatal morbidity and mortality.4,6,11-16 However, findings of a detrimental effect of neonatal transport have not been uniform.17,18 This suggests that reasons for transport may vary, as do their consequences on outcomes. Previous attempts to assess the effect of regionalization on outcomes have been limited by their reliance on hospital-level characteristics or before-after study designs, rather than an understanding of the flow of patients through a care delivery system.3,6,7,9,12,19,20 Accounting for the characteristics of the patients and their movement through the care delivery network would provide for a more granular understanding of the way care networks function for infants requiring acute transport.

Network analysis, which stems from graph theory in mathematics, analyzes the structure of relationships; it has been utilized in many fields, such as sociology, biology, and economics.21,22 Network analysis has previously been used to examine the structure of adult intensive care unit transports in Texas,23 but has not yet been applied to the field of neonatal regionalization. The goal of this study was to test the construct validity of network analysis as a tool to quantify the linkages between hospitals providing neonatal care in California.



We undertook a cross-sectional database study including all neonatal transfers (up to the age of 28 days) conducted by California Perinatal Quality Care Collaborative (CPQCC) member hospitals that occurred between January 1, 2012 and December 31, 2012.24 Data are collected on infants by CPQCC member hospitals using any of the following inclusion criteria: birth weight 401 to 1500 grams, gestational age 22 0/7 weeks to 29 6/7 weeks, or for infants >1500 grams either death, surgery, intubation or positive pressure support for more than 4 hours, readmitted for total bilirubin ≥25 and/or exchange transfusion, early bacterial sepsis, or acute transfer. Acute transfer was defined as an infant who requires acute resolution of medical problems and who is transferred in order to obtain care that is not provided, or that cannot be effectively provided at the referring institution (e.g., staffing/census issues, insurance).25 Within the CPQCC dataset, the California Perinatal Transport System, a system of over 100 specialized NICUs that participate in the transport of critically ill infants in California, collects data on acute neonatal transports.26 Transports originating from or traveling to hospitals not included in the CPQCC network, from outside the state of California, or from outside the United States, and infants admitted from home, were excluded from the cohort.


Objective 1: Representation of neonatal transfers using a network graph

We employed techniques of network analysis to characterize the movement of patients across the system of hospitals.21,22 This approach utilizes a network graph to represent mathematically the structure, direction, and intensity of relationships between distinct entities using a set of connected nodes. In this network graph, hospitals were represented by nodes; neonatal transfers between hospitals were represented as connections, or edges, between them. Edges between nodes can be represented either as arrows, representing the direction of the transfer between two hospitals (known as a “directed” network), or as simple lines, in which case the direction of the transfer is not denoted (an “undirected” or “symmetric” network). For clarity of presentation, we collapsed all transfers between each pair of hospitals into a single edge, and the network was symmetrized by removing the directionality of transfers.

There are several methods of representing the volume of transfers between each pair of hospitals, most commonly to proportionally widen the line as number of transfers increases. However, given the large size of this network and its relative density, we represented volume of transfers into each hospital by the size of the node. The American Academy of Pediatrics (AAP) level designation of each hospital, as derived in a 2012 survey of NICU directors,27 was denoted by the shape of the node.

After establishing the network structure, we characterized the most interactive communities of the network – that is, those groups of hospitals between which transfers occur most frequently, or regional referral networks – using network community detection techniques.28 Community detection separates the nodes of a network into groups (that is, communities) where there are many connections, while simultaneously attempting to minimize the number of connections between groups. In this way, community detection techniques break a network down along its most natural divisions. In our analysis, communities were identified using a hierarchical agglomeration algorithm for detecting community structure, a bottom-up approach that starts by considering each node as its own community, and successively merges pairs of communities.29 We did not pre-specify a desired number of communities (henceforth termed “sub-networks”), thereby allowing the algorithm to determine the optimal number of clusters to maximize within-cluster transfers and minimize between-cluster transfers.

Objective 2: Comparison of empirical sub-networks with state-defined regulatory regions

In order to establish construct validity of the network analysis approach, we assessed whether the networks derived from the analysis would be consistent with our knowledge of existing transfer patterns across California. To this end, we compared the network analysis derived sub-networks (henceforth termed “empirical”) to perinatal referral regions defined by the state of California (henceforth termed “regulatory”). The Regional Perinatal Programs of California (RPPC), established by the California Department of Public Health, has divided the state into perinatal referral regions.30 Perinatal centers are tasked with conducting oversight and outreach on quality improvement to other hospitals in the region.31,32 There are 9 RPPC-designated referral regions, as well as 2 regions defined by Kaiser Permanente (splitting the state between north and south, and here included under the designation “regulatory”).30 We generated a heat map, in which the row-by-column transports between empirical and regulatory regions were specified, to visualize the degree of alignment between these two types of networks. We quantified agreement between empirical sub-networks and regulatory regions using a chi-square test.

Objective 3: Identifying factors associated with transfer outside empirical sub-networks

In order to determine factors associated with cross-sub-network transfers, all transfers were coded according to whether they occurred between two hospitals within the same empirical sub-network or crossed sub-network boundaries. Variables included gestational age, birth weight, sex, race/ethnicity, major congenital anomaly, and documented reason for transfer (medical/diagnostic services, surgery, insurance, or bed availability). Congenital anomalies comprised central nervous system defects, congenital heart defects, gastrointestinal defects, genitourinary defects, chromosomal anomalies, pulmonary defects, vascular/lymphatic defects, and other conditions such as skeletal dysplasia, congenital diaphragmatic hernia, and hydrops fetalis.33 We used multivariable logistic regression to assess cross-sub-network transfer, using p<0.05 as the threshold for significance. Network visualization and statistical analyses were undertaken with R using the igraph package (R Development Core Team, Vienna, Austria).34 The study was approved by the institutional review boards of the investigators.


Description of patient population

The sample included 6546 infants, representing 6696 inter-hospital transports. Characteristics of these patients and transfers are shown in Table 1. Of note, 14% (n=915) of infants were very low birth weight, 24% (n=1572) of the infants had a major congenital anomaly, and 2% (n=150) underwent more than one transfer. Of 296 hospitals, 167 were designated as Level 1, 24 as Level 2, 86 as Level 3, and 19 as Level 4 units. In 91% (n=6122) of transfers, the destination hospital was of a higher level than the originating hospital, while 6% (n=422) transferred to the same level of care (Table 2).

Table 1
Sample characteristics of acutely transported infants in California, 2012
Table 2
Originating hospital AAP1 level versus receiving hospital AAP level

Objective 1: Representation of neonatal transfers using a network graph

The network graph of the neonatal transport network in California, overlaid onto a map of the RPPC regulatory regions, is shown in Figure 1.Node shape denotes AAP level of care, and node size reflects the number of transfers into that hospital.

Figure 1
Network graph of acute neonatal transports in California

The clustering algorithm empirically identified 11 sub-networks, the same number as the regulatory regions. Each sub-network is represented by a different node color. The largest of these empirical sub-networks included 46 hospitals, while the smallest included 7 hospitals. Nine of the 11 empirical sub-networks included at least one Level 4 NICU.

Objective 2: Comparison of empirical sub-networks with state-defined regulatory regions

Table 3 shows the number of hospitals coexisting in each pair of empirical sub-network and RPPC regulatory regions. The degree of overlap between the empirical sub-networks and the regulatory regions is significant (p<0.001). The minimum amount of mismatch (that is, the number of hospitals not falling into the largest overlapping cell per sub-network/regulatory region pair) was 0%, denoting perfect alignment of one RPPC regulatory region with one empirical sub-network. The greatest amount of mismatch for hospitals within the empirical sub-networks was 53% (that is, only 47% of that sub-network's hospitals fell into a single regulatory region), while the maximum mismatch was 43% when the analysis was based on regulatory regions. Of note, Kaiser North was the only regulatory region that did not map onto its own empirical sub-network; the empirical sub-networks with which it overlapped were more highly connected to other regulatory regions. Conversely, there were several empirical sub-networks that overlapped with the Los Angeles regulatory region; only the empirical sub-network with the largest number of overlapping hospitals (cyan) was considered as Los Angeles' best-matching sub-network.

Table 3
Overlap between empirical sub-networks and RPPC1 regulatory regions

Objective 3: Identifying factors associated with transfer outside empirical sub-networks

Overall, 91% (n=6086) of neonatal transports remained within their empirical sub-network, while 83% (n=5534) remained within their regulatory region (Table 4a and and4b).4b). Comparison of characteristics between transports that stayed within an empirical sub-network compared to those that crossed empirical sub-network boundaries is shown in Table 5.Significant predictors of cross-sub-network transport included the presence of a major congenital anomaly (OR 1.79 [95% CI, 1.47-2.19], p<0.001), need for surgery (OR 1.42 [95% CI, 1.08-1.88], p=0.01), and insurance as the documented reason for transfer (OR 5.13 [95% CI, 3.89-6.76], p<0.001) (Table 5). Hispanic infants were significantly less likely to cross sub-network boundaries during transport (OR 0.81 [95% CI, 0.67-0.98], p=0.03).

Table 4a
RPPC region of originating versus receiving hospital for all transports
Table 4b
Empirical sub-network of originating versus receiving hospital for all transports
Table 5
Infant characteristics associated with transfer within versus out of empirical sub-networks


This paper describes the use of network analysis as a novel methodology for accurately visualizing regionalized neonatal health care delivery systems. By conducting the analysis in a highly regulated state in which the existing referral relationships are well understood, we have demonstrated that network analysis can accurately depict known transport patterns. Furthermore, we showed that infant transports out of an empirical sub-network are associated with patient factors, including the presence of a major congenital anomaly, need for surgery, and insurance status; again aligning with our clinical expectations. Our findings are important, because they support the validity of network analysis as a tool to empirically quantify the degree of regionalization of neonatal care networks. Network analysis could thus be a powerful tool to define, analyze, and improve care at the network level in regions that lack a strong regulatory structure.

Taking an empirical network approach to neonatal regionalization addresses several limitations of prior studies of regionalization. First, definitions of NICU level of care vary considerably across the United States, confounding the definition of regionalization and, subsequently, studies based on these levels.18,35 While network analysis can incorporate these levels of care, the visualization of the networks and the calculation of network metrics does not depend on them. Moreover, before-after studies of the impact of regionalization programs on neonatal outcomes fall prey to confounding by secular trends, inadvertently reflecting concurrent changes in the field.9 Network analysis elucidates how the flow of patients changes over time, clarifying the impact of regionalization itself, regardless of other changes. Having a method to quantify neonatal care at a systems level represents a more rigorous and complete approach to assessing regionalization.

The finding of significant overlap between empirical and regulatory care networks suggests that hospitals are generally, but not always, transporting neonates in accordance with state-defined neonatal referral regions. Regulatory regions may mirror current practice because the regions are already optimally drawn, or hospitals may simply be working within their state-defined boundaries. While the two types of networks do overlap, empirical networks may represent a better unit of regional analysis, as they have fewer out-of-network transfers than the regulatory networks.

Although more than 90% of transports stay within their empirical sub-network, the patients who do cross sub-networks may illustrate factors shaping current neonatal care networks and may also identify suboptimal patterns of care delivery. Infants with major congenital anomalies, who comprised nearly a quarter of our patient population, typically require specialized care after delivery and thus frequently require transport to a regional referral center to receive subspecialty or surgical services. With prenatal ultrasound screening, most of these infants could be identified prenatally and their postnatal needs anticipated, allowing for delivery at an appropriate center. Likewise, patients who require surgery are more likely to be transported beyond their empirical sub-network, for related reasons of appropriate care availability. The high rate of transport across sub-network boundaries for these infants may indicate highly successful regionalization, as such subspecialized care cannot realistically be delivered by a larger number of centers. On the other hand, this could signify a breakdown in regionalization, if providers are bypassing closer, equally capable centers. For example, the clustering analysis detected real historical relationships between hospitals, showing that certain hospitals that are geographically disconnected from a regional referral center are a part of that center's empirical sub-network.

It is important to recognize that these findings do not directly address quality of care; however, they hold promise for optimizing care networks for efficiency and quality. Network analysis techniques that allow comparisons between sub-networks have the potential to quantify the association between care delivery network structure and quality of care. For example, network analysis techniques can measure how tightly connected a network is; if such cohesiveness measures are correlated with neonatal outcomes, there may be opportunities to improve the quality of neonatal care by optimizing the networks in which care is delivered. Moreover, network analysis allows for visualization of care patterns that appear suboptimal, such as when infants are transported long distances, bypassing closer centers that could provide appropriate care. While many of these instances can be explained by historical relationships between hospitals or providers, visualizing these cases would allow policy makers to identify areas for improvement in the flow of neonatal patients through care networks.

Our findings thus have important health policy implications. Quantifying the degree of regionalization of neonatal care networks will allow between-network comparisons as well as measurement of changes over time. For example, while there have been descriptions of deregionalization in the United States, longitudinal comparison of network graphs will allow quantification of this deregionalization.5 Furthermore, policy-makers will be able to measure the effects of new legislation such as the Affordable Care Act on care delivery systems. Also, as hospital network mergers and acquisitions occur, such quantitative approaches hold promise not only in evaluating the effects of shifting hospital systems, but also in informing such decisions prior to their execution.

This study must be viewed within the context of its design. Not all neonatal transports in California were included, as we excluded non-acute transports, as well as interstate and international transfers. We are currently working to incorporate non-acute transports for convalescent care and maternal transports into our analyses. Furthermore, because the data are based on self-report, some data may be missing or inaccurately classified. However, as the CPQCC database captures over 90% of neonates cared for in California NICUs, and the results reflect our expectations of the transport networks in California, this sample appears adequate. In addition, while we anticipate that network analysis techniques will be generalizable to less regulated states based on the adaptability of the underlying mathematical principles, these methods need to be empirically tested in other regions, particularly those with weaker regulatory oversight, differing geography, or known reliance on interstate transfer. Finally, these data do not assess the degree to which care is provided to meet specific goals, such as parental preference for a certain institution.

Even though Kaiser provides tightly integrated care, Kaiser North was not identified by our algorithm as a separate sub-network. We believe that this is likely a function of the cases we considered and how Kaiser is organized. Because Kaiser puts great emphasis on having mothers deliver at appropriate settings, many cases that would be neonatal transfers in other systems are antenatal referrals/transports and thus are not included in our data. Further, this exclusion shifts the balance of types of infants transferred by Kaiser to more complex cases that need subspecialty care provided at a small number of hospitals (e.g., cardiac surgery). Since Kaiser contracts with non-Kaiser hospitals for many of these types of cases, these appropriate referrals are, by definition, out of network, while transfers from many of the non-Kaiser hospitals who contract with the same tertiary hospitals will be in-network in our analyses. This is supported by additional analyses (not shown), showing that infants transported outside of Kaiser North disproportionately required subspecialty care. Finally, although we have employed a well-established type of network community detection algorithm29 that appears to accurately map transport sub-networks, it is possible that our network analysis algorithm was not a good tool to compare networks that are distributed by membership rather than location.

In conclusion, applying network analysis techniques to the field of neonatal regionalization offers an unprecedented view of how care delivery networks are shaped. In the future, such quantification of neonatal referral networks could be used to correlate network structure with patient outcomes. Network analysis has the potential to significantly influence regionalized neonatal care delivery systems.


Dr. Kunz's effort was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development [5T32HD075727-02; PI – Finkelstein] and through a Marshall Klaus Perinatal Research Award from the American Academy of Pediatrics. Dr. Profit's effort was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development [R01 HD083368-01, PI - Profit] and the Stanford Child Health Research Institute [1111239-285-JHACT; PI - Profit].

Funding source: Dr. Kunz's effort was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development [5T32HD075727-02, PI – Finkelstein] and through a Marshall Klaus Perinatal Research Award from the American Academy of Pediatrics. Dr. Profit's effort was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development [R01 HD083368-01, PI - Profit] and the Stanford Child Health Research Institute [1111239-285-JHACT; PI - Profit].


Conflict of Interest: The authors have no conflicts of interest relevant to this article to disclose.


1. Harrison W, Goodman D. Epidemiologic Trends in Neonatal Intensive Care, 2007-2012. JAMA Pediatr. 2015;169:855–862. [PubMed]
2. Rogowski JA, Staiger DO, Horbar JD. Variations in the quality of care for very-low-birthweight infants: implications for policy. Health Aff (Millwood) 2004;23:88–97. [PubMed]
3. Paneth N, et al. Newborn intensive care and neonatal mortality in low-birth-weight infants: a population study. N Engl J Med. 1982;307:149–155. [PubMed]
4. Phibbs CS, et al. Level and volume of neonatal intensive care and mortality in very-low-birth-weight infants. N Engl J Med. 2007;356:2165–2175. [PubMed]
5. Holmstrom ST, Phibbs CS. Regionalization and mortality in neonatal intensive care. Pediatr Clin North Am. 2009;56:617–630. [PubMed]
6. Lasswell SM, Barfield WD, Rochat RW, Blackmon L. Perinatal regionalization for very low-birth-weight and very preterm infants: a meta-analysis. JAMA. 2010;304:992–1000. [PubMed]
7. Lorch SA, Baiocchi M, Ahlberg CE, Small DS. The differential impact of delivery hospital on the outcomes of premature infants. Pediatrics. 2012;130:270–278. [PMC free article] [PubMed]
8. Kastenberg ZJ, Lee HC, Profit J, Gould JB, Sylvester KG. Effect of deregionalized care on mortality in very low-birth-weight infants with necrotizing enterocolitis. JAMA Pediatr. 2015;169:26–32. [PubMed]
9. Rashidian A, et al. The effectiveness of regionalization of perinatal care services--a systematic review. Public Health. 2014;128:872–885. [PubMed]
10. Profit J, Wise PH, Lee HC. Consequences of the Affordable Care Act for sick newborns. Pediatrics. 2014;134:e1284–1286. [PubMed]
11. Phibbs CS, Bronstein JM, Buxton E, Phibbs RH. The effects of patient volume and level of care at the hospital of birth on neonatal mortality. JAMA. 1996;276:1054–1059. [PubMed]
12. Gortmaker S, Sobol A, Clark C, Walker DK, Geronimus A. The survival of very low-birth weight infants by level of hospital of birth: a population study of perinatal systems in four states. Am J Obstet Gynecol. 1985;152:517–524. [PubMed]
13. Bowman E, Doyle LW, Murton LJ, Roy RN, Kitchen WH. Increased mortality of preterm infants transferred between tertiary perinatal centres. BMJ. 1988;297:1098–1100. [PMC free article] [PubMed]
14. Palmer KG, et al. Effect of inborn versus outborn delivery on clinical outcomes in ventilated preterm neonates: secondary results from the NEOPAIN trial. J Perinatol. 2005;25:270–275. [PubMed]
15. Mohamed MA, Aly H. Transport of premature infants is associated with increased risk for intraventricular haemorrhage. Arch Dis Child Fetal Neonatal Ed. 2010;95:F403–407. [PubMed]
16. Arora P, et al. Impact of interhospital transport on the physiologic status of very low-birth-weight infants. Am J Perinatol. 2014;31:237–244. [PubMed]
17. Cifuentes J, et al. Mortality in low birth weight infants according to level of neonatal care at hospital of birth. Pediatrics. 2002;109:745–751. [PubMed]
18. Profit J, et al. The Association of Level of Care With NICU Quality. Pediatrics. 2016;137:1–9. [PMC free article] [PubMed]
19. McCormick MC, Shapiro S, Starfield BH. The regionalization of perinatal services. Summary of the evaluation of a national demonstration program. JAMA. 1985;253:799–804. [PubMed]
20. Neto MT. Perinatal care in Portugal: effects of 15 years of a regionalized system. Acta Paediatr. 2006;95:1349–1352. [PubMed]
21. Jackson MO. Social and economic networks. Vol. 3. Princeton university press; 2008.
22. Easley D, Kleinberg J. Networks, Crowds, and Markets: Reasoning About a Highly Connected World. Cambridge University Press; 2010.
23. Iwashyna TJ, Christie JD, Moody J, Kahn JM, Asch DA. The structure of critical care transfer networks. Med Care. 2009;47:787–793. [PMC free article] [PubMed]
24. Gould JB. The role of regional collaboratives: the California Perinatal Quality Care Collaborative model. Clin Perinatol. 2010;37:71–86. [PubMed]
25. CPQCC Network Database. Manual of Definitions: For Infants Born in 2012. 2011;Version 12.1:28–33.
26. Gould JB, Danielsen BH, Bollman L, Hackel A, Murphy B. Estimating the quality of neonatal transport in California. J Perinatol. 2013;33:964–970. [PubMed]
27. American Academy of Pediatrics Committee on, F & Newborn. Levels of neonatal care. Pediatrics. 2012;130:587–597. [PubMed]
28. Fortunato S. Community detection in graphs. Physics Reports. 2010;486:75–174.
29. Clauset A, Newman ME, Moore C. Finding community structure in very large networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004;70:066111. [PubMed]
30. Health CDoP. Regional Perinatal Programs of California (RPPC) 2015
31. Kattwinkel J, Nowacek GA, Cook LJ, Hurt H, Short JG. Perinatal outreach education. A continuation strategy for a basic program. Am J Perinatol. 1984;1:335–340. [PubMed]
32. Nowacek GA, Cook LJ, Kattwinkel J. Assessment of transportability of a perinatal education program. South Med J. 1983;76:1490–1492. [PubMed]
33. CPQCC Network Database. Manual of Definitions: For Infants Born in 2016. 2015;Version 15.1
34. Kolaczyk ED, Cs´rdi GB. Statistical analysis of network data with R. Springer; 2014.
35. Blackmon LR, Barfield WD, Stark AR. Hospital neonatal services in the United States: variation in definitions, criteria, and regulatory status, 2008. J Perinatol. 2009;29:788–794. [PubMed]