 Research article
 Open Access
 Published:
Spatial transmission network construction of influenzalike illness using dynamic Bayesian network and vectorautoregressive moving average model
BMC Infectious Diseases volume 21, Article number: 164 (2021)
Abstract
Background
Although vaccination is one of the main countermeasures against influenza epidemic, it is highly essential to make informed prevention decisions to guarantee that limited vaccination resources are allocated to the places where they are most needed. Hence, one of the fundamental steps for decision making in influenza prevention is to characterize its spatiotemporal trend, especially on the key problem about how influenza transmits among adjacent places and how much impact the influenza of one place could have on its neighbors. To solve this problem while avoiding too much additional timeconsuming work on data collection, this study proposed a new concept of spatiotemporal route as well as its estimation methods to construct the influenza transmission network.
Methods
The influenzalike illness (ILI) data of Sichuan province in 21 cities was collected from 2010 to 2016. A joint pattern based on the dynamic Bayesian network (DBN) model and the vector autoregressive moving average (VARMA) model was utilized to estimate the spatiotemporal routes, which were applied to the two stages of learning process respectively, namely structure learning and parameter learning. In structure learning, the firstorder conditional dependencies approximation algorithm was used to generate the DBN, which could visualize the spatiotemporal routes of influenza among adjacent cities and infer which cities have impacts on others in influenza transmission. In parameter learning, the VARMA model was adopted to estimate the strength of these impacts. Finally, all the estimated spatiotemporal routes were put together to form the final influenza transmission network.
Results
The results showed that the period of influenza transmission cycle was longer in Western Sichuan and Chengdu Plain than that in Northeastern Sichuan, and there would be potential spatiotemporal routes of influenza from bordering provinces or municipalities into Sichuan province. Furthermore, this study also pointed out several estimated spatiotemporal routes with relatively high strength of associations, which could serve as clues of hot spot areas detection for influenza surveillance.
Conclusions
This study proposed a new framework for exploring the potentially stable spatiotemporal routes between different places and measuring specific the sizes of transmission effects. It could help making timely and reliable prediction of the spatiotemporal trend of infectious diseases, and further determining the possible key areas of the next epidemic by considering their neighbors’ incidence and the transmission relationships.
Background
Influenza is an acute respiratory infection caused by influenza viruses [1]. According to the globally estimated disease burden attributed to the influenza, it is estimated there are about 1 billion seasonal influenza cases per year on average, among which three million cases are serious, resulting in 250,000 to 500,000 deaths [2].
Although vaccination is one of the main countermeasures against influenza epidemic [3, 4], some drawbacks still exist. For example, various influenza virus strains are prone to undergo antigen drift and antigen conversion, which often makes vaccineinduced immunity wane over the course of a season. Besides, there exists a natural process of decreasing trend of antibody titers after vaccination [5], suggesting people need to be vaccinated periodically, leading to a huge burden to lowandmiddle income countries. Therefore, to compensate for the huge cost by prevention strategies, it is quite necessary to make informed prevention decisions to guarantee that limited resources are allocated to the places where they are most needed [6]. One of the fundamental steps for decision making in influenza prevention is to characterize its spatiotemporal trend, especially on the key problem about how influenza transmits among adjacent places and how much impact the influenza of one place could have on its neighbors. To solve such problems, Fu et al. [7] developed complex networks with population contact data to predict the epidemic trend in a mathematical way. Recently, Pei et al. [8] utilized accessible human mobility data and a metapopulation model for predicting the spatial transmission of influenza in the United States. In addition, another study has utilized birds migration network to predict the trajectory of avian influenza [9]. Those methods require personal contact data including person mobility data, traffic data, avian mobility data and so on. However, there will be some obstacles on the availability of those data, and collecting those data will be too timeconsuming to make rapid strategies of influenza prevention.
To overcome these obstacles of previous researches, this study proposed a new concept called spatiotemporal route to display the potential transmission directions of influenza and to measure the sizes of those transmission effects. This concept was defined as the timelagged association among influenza surveillance data across different places. In addition, since the estimated spatiotemporal route only depends on surveillance data and does not necessarily need personal contact data, it will be convenient to be used for exploring the transmission network with realworld surveillance data within statistical framework. For illustration purpose, this study selected Sichuan province in China as an example, but the concept of spatiotemporal route as well as its estimation methods mentioned below could also be applied to other places.
Materials and methods
Data preparation
The data of this study came from all the sentinels of influenzalike illness (ILI) of Sichuan province from 2010 to 2016. Sentinel surveillance is one of the important measures for infectious diseases surveillance. According to the unified deployment of the National Center for Disease Control and Prevention (CDC) and the realworld situation of Sichuan province, the sentinel surveillance of influenzalike illness is simultaneously conducted in each of the 21 cities in Sichuan province by hospitals, CDCs, and primary health service institutions. According to the report requirement, the medical staffs of monitoring clinics in sentinels recorded the number of ILI and the total number of outpatients in each age group in each department every day, and uploaded the data to the China Influenza Surveillance Information System before midnight every Monday. For this study, the definition of influenzalike illness (ILI) referring to WHO [10] was as follows: a case measured fever of ≥38 C° and cough; with onset within the last 10 days. Besides, to estimate the absolute ILI case number in city j, we collected the data of yearly number of medical outpatients in each city of Sichuan from Sichuan Health Statistics Yearbook 2010 ~ 2012 and Sichuan Health and Family Planning Statistical Yearbook 2013 ~ 2016. We defined the ILI case number in city i and week t as ILI(i, t) so that for any city i (i.e., when i is fixed), ILI(i, t) (t = 1, 2, 3, …) could be regarded as time series, and the task of this study was to model the timelagged correlation between any of the two time series ILI(i, t) and ILI(j, t) (i ≠ j, t = 1, 2, 3,…).
The estimation of spatiotemporal routes
The definition of spatiotemporal routes
Some previous researches engaged in constructing the influenza transmission network by temporal and spatial statistics [11,12,13,14,15,16]. For example, Alonso WJ [16] used Fourier decomposition to find a seasonal southward traveling wave of influenza across Brazil originating from equatorial and low population regions in March–April and moving towards temperate and highly popular regions over a 3month period. In addition, Paul and Held [17] proposed a random effect model (the epidemicendemic model) to consider the transmission effects of neighboring places. From the statistical point of view, the phenomenon that influenza transmitting from place A to place B could be reflected by the timelagged association. Equivalently, the timelagged association could also be visualized by A → B, where the directed arc indicated that node A (i.e., influenza in place A) had a timelagged effect on node B (i.e., influenza in place B), and we defined such directed arc as spatiotemporal route. Furthermore, if the timelagged associations existed in more than two places in the overall study area, then all the arcs would interweave into a network. Such a network could show how historical influenza in one place would influence its neighbors in the near future so as to make some possible inferences on the temporal and spatial transmission features of influenza. To this end, we defined this network as spatial transmission network because it was essentially a set of spatiotemporal routes.
More precisely, the spatiotemporal routes could also be defined in a mathematical way. Let X = {ILI (i,r)}be the set of all the influenza data from different places. Define A the set of arcs between any two places in set X and then the spatiotemporal routes could be defined as network G = (X, A). Specifically, the network G contains two types of information. The first type was the structure information, which was related to arc existence as well as its direction for any pair of two nodes in G. The second type was the parameter information, which measured the strengths of associations among different nodes. Correspondingly, the estimation of spatiotemporal routes consisted of structure learning and parameter learning, which was dedicated to drawing the structure and parameter information respectively from the original data. Specifically, this study used the dynamic Bayesian network (DBN) model for structure learning and the vector autoregressive moving average (VARMA) model for parameter learning. More details were given as below.
The structure learning of spatiotemporal routes by the DBN
The DBN is a dynamic directed acyclic graph (DAG) using nodes and arcs to express the conditional probabilistic dependencies between a set of time series [18]. In the DBN model, an arc is drawn between two variables at successive time points. For example, from ILI(i, t1) to ILI(j, t), which means the ILI cases of city j at time t (e.g., the current week) are conditionally dependent on the ILI cases of city i at time (t1) (e.g., 1 week ago) given the remaining variables at the past time points. Due to its good theoretical properties, the DBN model was used to characterize the gene regulatory network by characterizing the timelag associations of multiple gene expression data [19, 20]. Recently, our previous work had also used simulation studies to prove that the DBN could be very well applied to the infectious disease surveillance data even when confronted with some rigorous challenges such as high noise, nonlinear correlation, small sample and latent variables [21]. Therefore, this study used the DBN model for the structure learning of influenza spatiotemporal routes.
In particular, we estimated the DBN model by using the firstorder conditional dependencies approximation algorithm [22]. It implemented DBN learning as a twostep procedure. At the first step, it learned a DAG encoding firstorder partial dependence relationships. Then it inferred the real network structure of the DBN using the graph from the first step [18]. Once the structure learning of influenza spatiotemporal routes was completed, it could be used to infer which cities have impact on others in influenza transmission. Then the next step was to further estimate the strength of these impacts by the means of parameter learning.
The parameter learning of spatiotemporal routes by the VARMA model
As mentioned above, the parameter learning of influenza spatiotemporal routes was required to quantify how the current ILI cases in one place were impacted by the past ILI cases in other places. To this end, the multivariate time series (MTS) models were used for parameter learning [23]. Furthermore, it has been proven that the DBN is mathematically equivalent to the vectorautoregressive model (VAR) [24] (i.e., one of the most commonly used MTS model), which again shed light upon the application of the VAR model to parameter learning of influenza spatiotemporal routes. Specifically, the basic formula of VAR(p) was:
where K stood for the overall number of the involved places (i.e., K = 21 in this study since there were 21 cities in Sichuan province), ILI(i, t) was defined in Section 2.1 and ε(i, t) was defined as the residual of the fitted model for ILI cases in city i at time t. For the VAR model, it assumed that at any given time t, all the residualsε (1, t), ε (2, t), …, ε(K, t) were independent with each other. However, in the realworld situation of influenza transmission, since some important information of influenza transmission factors (e.g., population density, the effectiveness of influenza transmission in the city, festival effects and so on) may not be captured, consequently the assumption of residual independence would be violated. Therefore, this study improved the VAR model to the VARMA model as below [25]:
Compared with the VAR model, the VARMA model could compensate for the violation of residual independence assumption by adding the moving average term \( \sum \limits_{j=1}^q{\boldsymbol{B}}_j{\boldsymbol{\varepsilon}}_{tj} \) to extract the remaining dependency information between residuals. Therefore, it was plausible that the VARMA model could complete the parameter learning while reducing the error of parameter estimation as possible.
During the estimation of spatiotemporal routes, both the DBN and the VARMA models could be implemented in R 3.6.3, a free sofware environment for statistical computing and graphics. The DBN model was implemented using the {G1DBN} package and the VARMA model was built by the {MTS} package. All the packages were downloaded from the Comprehensive R Archive Network (CRAN) at http://cran.rproject.org/ and installed in advance. Additionally, the time series plots of ILI% in Sichuan and the maps of Sichuan were generated by us with R software.
Results
The ILI% in Sichuan between 2010 and 2016
The total number of outpatients in Sichuan surveillance sentinels from 2010 to 2016 was 31,898,487, and the total number of influenzalike cases in Sichuan was 784,984. As a result, the ILI% of Sichuan was 2.46% within 6 years. From 2010 to 2016, the cumulative ILI% was 2.51, 2.36, 2.38, 2.69, 2.60, 2.34 and 2.36%, respectively. The year with lowest ILI number was 2010 (84,766 visits), and highest was 2016 (137,945 visits); the lowest ILI% was in 2016 (2.36%) and the highest was in 2013 (2.69%). The 2010–2016 weekly ILI% distribution was shown in Fig. 1 and Fig. 2 below:
From Fig. 1 and Fig. 2, the weekly ILI% in Sichuan province from 2010 to 2016 was between 1.69 and 3.66%, which was consistent with the past years. From Fig. 2, it could be seen there was a slight peak in winter and spring (From the 47th to 52nd week, and from 1st to 9th week), similar to the result of the whole country [26].
Results of structure learning of possible influenza spatiotemporal routes
Combining the longest incubation period of influenza (i.e. 2 weeks [27]) and our exploratory work, we set 1week lag, 2week lag and 3week lag as lags of time order respectively, and the final influenza spatial transmission network was obtained as Fig. 3:
In the perspective of prevention for influenza, from Fig. 3, it could be observed that during 3week influenza transmission cycle, the number of suspected influenza spatiotemporal routes in Sichuan province showed a decreasing trend. Hence, it could be speculated that the period of influenza transmission cycle was longer in Western Sichuan and Chengdu Plain and shorter in Northeastern Sichuan, suggesting that we should consider extending the period for influenza prevention and control in Western Sichuan and Chengdu Plain appropriately compared with Northeastern Sichuan.
To verify the robustness of structure learning across different ages and years, we fitted DBN models within each data subset of different age groups (0–5 year age group, 5–15 year age group, 15–25 year age group, 25–60 year age group and 60+ year age group) and each year (2010, 2011, 2012, 2013, 2014, 2015 and 2016). In terms of age stratification analysis, we found that the structure learning results in each subset were almost consistent with the original result shown in Fig. 3, except for the result in 0–5 year age group and 60+ age year group. Such differences might be explained by the variance of population mobility across different age groups, because the aged and the infants were much less likely to go on a long journey than the young adults. As for subset analysis within each year, the results in 2010, 2011, 2012, 2013 and 2014 were also consistent with the structure in Fig. 3 A to some extent, but in 2015 and 2016 the results were a bit different. As the Fig. 2 implied, the ILI%s in those 2 years were lower, so it was plausible to infer that the lower sample size weakened the statistical power of structure learning in 2015 and 2016, resulting in inconsistencies in 2015 and 2016.
More importantly, it could be seen from Fig. 4 that there were four cities with singledirection arrows pointing to other cities in the province while with no city in the province pointing to themselves, i.e., Ziyang, Guangyuan, Yibin and Panzhihua. In addition, all the four cities are bordering cities in Sichuan province and essential transportation hubs connecting to other adjacent provinces. Furthermore, influenza is a typical humantohuman infectious disease, so transportation could play an important role in the transmission of influenza logically, which has already been confirmed by a large number of studies. For example, Hidenori [28] found that traffic control could delay the spread of flu during peak flu periods in a simulation study. Therefore, it was plausible to speculate that there would be potential spatiotemporal routes of influenza from bordering provinces or municipalities into Sichuan province.
In addition, to verify the rationality of this speculation, it was interesting to show the high correlation between the number of estimated influenza spatiotemporal routes and the highway transportation capacities in each subregion of Sichuan province (Table 1). Except for the correlation in numbers, the relations between influenza and transportation could again be suggested by the coincidence of the estimated influenza spatiotemporal routes and the major highways in Sichuan. Taking Chengdu Plain subregion as an example, the Ziyang→Meishan influenza spatiotemporal routes coincided with the Suizimei Expressway; the Chengdu→Ya’an influenza spatiotemporal route was to some extent in accordance with the Chengya Expressway; the Guangyuan→Mianyang← → Deyang← → Chengdu→Ya’an influenza spatiotemporal route was possibly in line with the G5 Jingkun Expressway (the Mianguang section, the Chengmian section and the Chengya section); the Nanchong ← → Mianyang ← → Deyang ← → Chengdu influenza spatiotemporal route was expected to be brought by the Chengdemiannan Expressway. Besides, Leshan→Meishan← → Ya’an influenza spatiotemporal route was supposed to be related with the Leya Expressway, which was another important highway flowing through Ya’an. More examples about the coincidence of the estimated influenza spatiotemporal routes and the major highways in other subregions in Sichuan province could be seen in the Additional file 1. All these examples indicated that highway transportation might be a key factor underlying the estimated influenza spatiotemporal routes in Sichuan province.
Results of parameter learning of possible influenza spatiotemporal routes
The estimated parameters of influenza spatiotemporal routes were summarized in Table 2, and meanwhile the specific parameter learning results using the VARMA model were in Additional file 2. It could be seen that the median of the estimated parameters lied around zero, indicating that most of the influenza spatiotemporal routes were in general not obvious. Among the four bordering cities mentioned in Section 3.2, the spatiotemporal effect of Ziyang deserved special attention. In particular, the 1week lagged spatiotemporal effect of Ziyang to Meishan was 0.35, and from Meishan to Ya’an was 0.07, implying that Ziyang played an essential role in terms of spatiotemporal effects. All these results showed that the estimated timelagged relationships of influenza cases among these cities were relatively close, suggesting that the prevention of influenza transmission among those cities should be highlighted. In addition, it could be seen from Table 2 that from a lag of 1 week to a lag of 3 weeks, the degrees of dispersion were almost stable, implying the change of spatiotemporal effects might not be obvious.
Discussion
This study proposed the concept of spatiotemporal route as well as its estimation methods to construct the influenza transmission network in a novel way. To our knowledge, this study may contribute to the infectious diseases surveillance in at least the following three ways.

(1)
This study initially proposed the concept of spatiotemporal route to better solve the problem of how to construct a potential spatial transmission network of influenza when mobility data is unavailable. On the one hand, in traditional epidemiology, there is a similar concept named the route of transmission, which mainly refers to the entire process experienced by pathogens in the external environment from the time they are discharged from infection sources to the time they invade new susceptible hosts [29]. Hence, one could judge by definition that it is impossible to answer the question of constructing a spatial transmission network by study on the route of transmission. On the other hand, in most of the latest researches about the influenza transmission network, the analytical models were based on theoretical physics or internet disciplines [30, 31], as well as the mobility information like the air transportation network. However, very few studies participated in constructing a spatial transmission network of influenza through the perspective of spatiotemporal distribution. Hence, the spatiotemporal routes not only help to clarify the parameters of interest in this study, but also provide a theoretical foundation for further researches to study the propagation and epidemic law of infectious diseases from the temporal and spatial dimensions.

(2)
This study also put forward a joint pattern based on DBN and VARMA models to estimate the spatiotemporal routes for the first time. Although neither of the two models was proposed by this study for the first time, this research has made full use of the theoretical properties and combined their advantages together. In particular, previous studies have proved the validity and robustness of DBN model when handling complicated data structure such as high dimension, high noise and nonlinearity [21], which served as a powerful guarantee to reveal the spatiotemporal correlations between the ILIs of different areas. In addition, the VARMA model has advantages in dealing with the potential confounders due to data unavailability in practice. Therefore, when the two models were combined, it was plausible that they would be well applied to infer influenza transmission network in the complicated realworld of influenza surveillance.

(3)
Another potential contribution of this study was to help making timely and reliable prediction of the spatiotemporal trend of influenza, and further determining the possible key areas for the next influenza epidemic outbreak. According to Stoto [32], a practical symptom surveillance system required continuous surveillance data (possibly multivariate data), an alert generated by the application of predictive algorithms, and a prescribed process for how to respond to the alert. To this end, it was promising that this study could help to improve the current influenza surveillance system in the following ways. Firstly, this study confirmed that our model could efficiently utilize continuous surveillance data of multiple places for influenza surveillance. Secondly, the results of this study revealed that some adjacent cities were indeed close to each other in ILI cases, which suggested that there might be some potential stable spatiotemporal routes among these cities so that the CDC could locate the key areas of the epidemics and send alarms when necessary. For example, as mentioned before, the spatiotemporal association between Ziyang to Meishan was relatively high. Therefore, if an influenza outbreak happened in Ziyang, one should pay attention to Meishan because our model reminded close relationship between them. Finally, on the basis of the previous early warning prevention and control policies, the local authorities receiving the alert could formulate specific prevention and control strategies according to the actual situation, such as identifying the key population for vaccination, enhancing inspection and quarantine as well as timely allocating health resources.
Although there were some interesting findings in this study, some limitations should also be acknowledged. First, as mentioned above, this study did not consider the possible influences affected by population density, influenza effectiveness, festival effects and other factors directly in the process of constructing the spatiotemporal routes. Although we have adopted the VARMA model which specifically dealt with residual effects as a remedy, it was definitely not as good as incorporating these factors into the analysis directly. Second, this study only analyzed the possible spreading directions of influenza between different cities from the perspective of spatiotemporal statistics. However, statistical significance cannot fully represent actual significance after all. In order to further demonstrate the issue of the spreading directions of spatiotemporal transmission, much additional work involving pathogenic detections, epidemiological investigations and so on still needs to be done. To this end, it is highly expected that this study could provide a bit inspiration and reference for future researches on surveillance of infectious diseases including but not limited to influenza.
Conclusions
This study proposed a new framework for exploring the potentially stable spatiotemporal routes among different places and measuring specific sizes of the transmission effects. It showed that the period of influenza transmission cycle was longer in the Western Sichuan and Chengdu Plain than that in Northeastern Sichuan, and there would be potential spatiotemporal routes of influenza from bordering provinces or municipalities into Sichuan province. Furthermore, this study also pointed out several estimated spatiotemporal routes with relatively high strengths of associations, which could serve as clues of hot spot areas detection for influenza surveillance. The results could be used for the detection and early warning of infectious diseases in the future.
Availability of data and materials
The data that support the findings of this study was obtained by applying it to the Sichuan CDC.
Abbreviations
 ILI:

Influenzalike illness
 DBN:

Dynamic Bayesian network
 VARMA:

Vector autoregressive moving average
 CDC:

Center for Disease Control and Prevention
 DAG:

Directed acyclic graph
 MTS:

Multivariate time series
 VAR:

Vectorautoregressive model
 CRAN:

Comprehensive R Archive Network
References
 1.
Su S, Fu X, Li G, Kerlin F, Veit M. Novel influenza D virus: epidemiology, pathology, evolution and biological characteristics. Virulence. 2017;8(8):1580–91.
 2.
Influenza (Seasonal). World health organization. 6 Nov 2018. https://www.who.int/newsroom/factsheets/detail/influenza(seasonal). Accessed 15 Feb 2019.
 3.
Sah P, AlfaroMurillo JA, Fitzpatrick MC, Neuzil KM, Meyers LA, Singer BH, et al. Future epidemiological and economic impacts of universal influenza vaccines. Proc Natl Acad Sci U S A. 2019;116(41):20786–92.
 4.
Grohskopf LA, Alyanak E, Broder KR, Walter EB, Fry AM, Jernigan DB. Prevention and control of seasonal influenza with vaccines: recommendations of the advisory committee on immunization practicesUnited States, 201920 influenza season. MMWR Recomm Rep. 2019;68(3):1–21.
 5.
Petrie JG, Ohmit SE, Truscon R, Johnson E, Braun TM, Levine MZ, et al. Modest waning of influenza vaccine efficacy and antibody titers during the 20072008 influenza season. J Infect Dis. 2016;214(8):1142–9.
 6.
Hill EM, Petrou S, de Lusignan S, Yonova I, Keeling MJ. Seasonal influenza: Modelling approaches to capture immunity propagation. PLoS Comput Biol. 2019;15(10):e1007096.
 7.
Fu X, Small M, Chen G. Propagation dynamics on complex networks. 1st ed. West Sussex: Wiley; 2014.
 8.
Pei S, Kandula S, Yang W, Shaman J. Forecasting the spatial transmission of influenza in the United States. Proc Natl Acad Sci U S A. 2018;115(11):2752–7.
 9.
Shi B, Zhan XM, Zheng JX, Qiu H, Liang D, Ye YM, et al. Identifying key bird species and geographical hotspots of avian influenza a (H7N9) virus in China. Infect Dis Poverty. 2018;7(1):97.
 10.
WHO surveillance case definitions for ILI and SARI. World Health Organization. 2014. https://www.who.int/influenza/surveillance_monitoring/ili_sari_surveillance_case_definition/en/. Accessed 30 Oct 2019.
 11.
LobatoCordero A, Quentin E, LobatoCordero G. Spatiotemporal analysis of influenza morbidity and its association with climatic and housing conditions in Ecuador. J Environ Public Health. 2019;2019:6741202.
 12.
Li X, Chan KKY, Xu B, Lu M, Xu B. Spatial, temporal and genetic dynamics characteristics of influenza B viruses in China, 19732018. Virol Sin. 2020;35(1):14–20.
 13.
Martin LJ, Dong H, Liu Q, Talbot J, Qiu W, Yasui Y. Predicting influenzalike illnessrelated emergency department visits by modelling spatiotemporal syndromic surveillance data. Epidemiol Infect. 2019;147:e312.
 14.
Xu ZW, Li ZJ, Hu WB. Global dynamic spatiotemporal pattern of seasonal influenza since 2009 influenza pandemic. Infect Dis Poverty. 2020;9(1):2.
 15.
Held L, Hofmann M, Höhle M, Schmid V. A twocomponent model for counts of infectious diseases. Biostatistics. 2006;7(3):422–37.
 16.
Alonso WJ, Viboud C, Simonsen L, Hirano EW, Daufenbach LZ, Miller MA. Seasonality of influenza in Brazil: a traveling wave from the Amazon to the subtropics. Am J Epidemiol. 2007;165(12):1434–42.
 17.
Paul M, Held L, Toschke AM. Multivariate modelling of infectious disease surveillance data. Stat Med. 2008;27(29):6250–67.
 18.
Nagarajan R, Scutari M, Lèbre S. Bayesian networks in the presence of temporal information. In: Bayesian networks in R: with applications in systems biology. 1st ed. New York: Springer Science+Business Media; 2013. p. 59–69.
 19.
Sandri M, Berchialla P, Baldi I, Gregori D, De Blasi RA. Dynamic Bayesian networks to predict sequences of organ failures in patients admitted to ICU. J Biomed Inform. 2014;48:106–13.
 20.
Kourou K, Rigas G, Exarchos KP, Papaloukas C, Fotiadis DI. Prediction of oral cancer recurrence using dynamic Bayesian networks. Annu Int Conf IEEE Eng Med Biol Soc. 2016;2016:5275–8.
 21.
Zhang T, Ma Y, Xiao X, Lin Y, Zhang X, Yin F, et al. Dynamic Bayesian network in infectious diseases surveillance: a simulation study. Sci Rep. 2019;9(1):10376.
 22.
Lèbre S. Inferring dynamic genetic networks with low order independencies. Stat Appl Genet Mol Biol. 2009;8(1):1–38.
 23.
Tsay RS. Stationary vector autoregressive time series. Vector autoregressive movingaverage time series. In: Multivariate time series analysis: with R and financial applications. 1st ed. New Jersey: Wiley; 2014. p. 27–198.
 24.
OpgenRhein R, Strimmer K. Learning causal networks from systems biology time course data: an effective model selection procedure for the vector autoregressive process. BMC Bioinformatics. 2007;8 Suppl 2:S3.
 25.
Zhang T, Yin F, Zhou T, Zhang XY, Li XS. Multivariate time series analysis on the dynamic relationship between class B notifiable diseases and gross domestic product (GDP) in China. Sci Rep. 2016;6(1):29.
 26.
Yu H, Alonso WJ, Feng L, Tan Y, Shu Y, Yang W, et al. Characterization of regional influenza seasonality patterns in China and implications for vaccination strategies: spatiotemporal modeling of surveillance data. PLoS Med. 2013;10(11):e1001552.
 27.
Zhou L, Li Q, Uyeki TM. Estimated incubation period and serial interval for humantohuman influenza A(H7N9) virus transmission. Emerg Infect Dis. 2019;25(10):1982–3.
 28.
Yasuda H, Yoshizawa N, Kimura M, Shigematsu M, Matsumoto M, Kawachi S, et al. Preparedness for the spread of influenza: prohibition of traffic, school closure, and vaccination of children in the commuter towns of Tokyo. J Urban Health. 2008;85(4):619–35.
 29.
Volovskaia ML. Epidemiology and fundamentals of infectious diseases. 1st ed. Moscow: Mir Publishers; 1990.
 30.
Rvachev LA, Longini IML Jr. A mathematical model for the global spread of influenza. Math Biosci. 1985;75:3–22.
 31.
Pei S, Shaman J. Counteracting structural errors in ensemble forecast of influenza outbreaks. Nat Commun. 2017;8(1):925.
 32.
Stoto A. Comment on: a bayesian dynamic model for influenza surveillance by Sebastianiet et al. statistics in medicine (this issue). Stat Med. 2006;25(11):1817–8.
Acknowledgements
Many thanks are expressed to deceased professor. The authors also want to thank Hongli Wan, Yang Zhong, Cheng Li, Jianxiang Tang and Wenhui Zhu for manuscript revisions.
Funding
Tao Zhang’s work was supported by Sichuan Science and Technology Program (2020YFS0015, 2020YFS0091), Health Commision of Sichuan province (20PJ092), National Natural Science Foundation of China (81602935), Chongqing Science and Technology Program (cstc2020jscxcylhX0003) and Sichuan University (2018hhf26, GSSCU2018038). Changhong Yang’s work was supported by Sichuan Science and Technology Program (2020YFS0015, 2020YFS0007, 2020YFS0001, 2020YFS0559). The funding body did not participate in the design, collection, analysis, interpretation and writing of this study.
Author information
Affiliations
Contributions
JQQ analyzed the data and was a major contributor in writing the manuscript. All authors read and approved the final manuscript, HMW and LH helped perform the analysis with constructive discussions, CHY collected and provided us with the data, and TZ contributed significantly to analysis and manuscript preparation.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The ILI surveillance was a routine surveillance activity. The analysis of ILI data was not considered as human subject research. No administrative permission was needed to assess the data.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
The association between traffic pattern and the spatiotemporal routes of Sichuan
Additional file 2.
The parameter learning results using VARMA
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Qiu, J., Wang, H., Hu, L. et al. Spatial transmission network construction of influenzalike illness using dynamic Bayesian network and vectorautoregressive moving average model. BMC Infect Dis 21, 164 (2021). https://doi.org/10.1186/s12879021057696
Received:
Accepted:
Published:
Keywords
 Influenza
 Spatial transmission network
 Dynamic Bayesian network
 Vector autoregressive moving average model
 Spatiotemporal route