 Research article
 Open Access
 Published:
Integration of a Kalman filter in the geographically weighted regression for modeling the transmission of hand, foot and mouth disease
BMC Public Health volume 20, Article number: 479 (2020)
Abstract
Background
Hand, foot and mouth disease (HFMD) is a common infectious disease whose mechanism of transmission continues to remain a puzzle for researchers. The measurement and prediction of the HFMD incidence can be combined to improve the estimation accuracy, and provide a novel perspective to explore the spatiotemporal patterns and determinant factors of an HFMD epidemic.
Methods
In this study, we collected weekly HFMD incidence reports for a total of 138 districts in Shandong province, China, from May 2008 to March 2009. A Kalman filter was integrated with geographically weighted regression (GWR) to estimate the HFMD incidence. Spatiotemporal variation characteristics were explored and potential risk regions were identified, along with quantitatively evaluating the influence of meteorological and socioeconomic factors on the HFMD incidence.
Results
The results showed that the average error covariance of the estimated HFMD incidence by district was reduced from 0.3841 to 0.1846 compared to the measured incidence, indicating an overall improvement of over 50% in error reduction. Furthermore, three specific categories of potential risk regions of HFMD epidemics in Shandong were identified by the filter processing, with manifest filtering oscillations in the initial, local and longterm periods, respectively. Amongst meteorological and socioeconomic factors, the temperature and number of hospital beds per capita, respectively, were recognized as the dominant determinants that influence HFMD incidence variation.
Conclusions
The estimation accuracy of the HFMD incidence can be significantly improved by integrating a Kalman filter with GWR and the integration is effective for exploring spatiotemporal patterns and determinants of an HFMD epidemic. Our findings could help establish more accurate HFMD prevention and control strategies in Shandong. The present study demonstrates a novel approach to exploring spatiotemporal patterns and determinant factors of HFMD epidemics, and it can be easily extended to other regions and other infectious diseases similar to HFMD.
Background
Hand, foot and mouth disease (HFMD) is a common infectious disease caused by at least 20 enteroviruses including enterovirus 71 (EVA71) and Coxsackie virus A16 (CAV16) [1]. HFMD usually affects infants and children under five and its main symptoms include fever, mouth ulcers and blisters or vesicles on the hands, feet, and mouth. Existing vaccines are only partially effective for specific HFMD pathogens [2]. The transmission mechanism of HFMD epidemics is complicated and its spatiotemporal pattern is not yet fully understood [3]. During the last decades, HFMD has been widespread in Asian countries, such as Japan, Malaysia, and Singapore [4,5,6]. In China, the first largescale outbreaks of HFMD occurred in Linyi city, Shandong province in 2007 [7] and in Fuyang city, Anhui province in 2008 [8]. Next, in May 2008 the Ministry of Health of China listed HFMD as a statutorily notifiable infectious category C disease. China’s infectious disease automated alert and response system (CIDARS) was developed in the same year for the early detection and rapid response to the outbreaks of infectious diseases, and the system performance was satisfactory in the detection of HFM disease outbreaks, with a sensitivity of 92.7% and a specificity of 95.0% [9]. An extensive threelevel HFMD surveillance laboratory network was established in mainland China since 2008, and the surveillance data from 2008 to 2017 indicated a high incidence of HFMD occurred every 2 years and the highrisk regions were located in southern, eastern, and central China [10]. Numerous studies on HFMD epidemics were implemented in various regions, particularly in provinces with serious epidemics, such as Guangdong [11, 12], Sichuan [13, 14], Henan [15, 16], Shandong [17, 18], and others.
Previous studies have mainly focused on characteristics of the epidemic [1, 15, 19], such as spatiotemporal patterns and correlations with various risk factors. HFMD epidemics have significant temporal variations and seasonality features, which vary between regions [20,21,22,23]. HFMD epidemics were spatially dispersed across counties in mainland China in the summer and winter, while clustered in spring and autumn; they were also geographically clustered in and closely linked to regions with high levels of monthly precipitation [3, 24]. In addition, HFMD epidemics follow complicated spatiotemporal patterns and transmission mechanisms, and are associated with several types of risk factors. For example, the HFMD incidence in Singapore has been found to be affected in a nonlinear manner by the maximum temperature and rainfall, with a time lag of 1–2 weeks, and thresholds of 32 °C and 75 mm, respectively [25]. Furthermore, in Japan and Vietnam, temperature and humidity had significant effects on the HFMD incidence [21, 26]. The spatial variation of HFMD in counties across mainland China was found to be affected by a combination of climate variables, while the spatiotemporal transmission was largely driven by variations in temperature, with a 7week lag [3]. Extreme precipitation was significantly associated with childhood HFMD in Hefei, China, and the susceptible risk in urban areas was much higher than that in rural ones [27]. Highrisk areas of HFMD incidence temporally varied from northeast to southwest in Sichuan, China, and temperature and per capita gross domestic product (GDP) were the main positive driving factors [13].
Nonlinear associations have been found between the HFMD incidence and meteorological, landuse, normalized difference vegetation index (NDVI) and socioeconomic factors in Shandong, China [18]. Many other studies have also focused on exploring of HFMD spatiotemporal patterns and the associated driving factors, by considering a variety of methods [3, 11,12,13,14, 16, 18, 20, 22,23,24, 26,27,28,29,30,31,32,33,34]. However, the measurement and prediction of the HFMD incidence are usually considered separately, and rarely in an integrated fashion. The former is mainly accomplished by using case reports, while the latter requires specific quantitative models. The modeling of HFMD transmission and the corresponding analysis results could be influenced by explanatory variable selection, spatial autocorrelation, spatial stratified heterogeneity, spatiotemporal nonstationary, etc. Excluding the above factors, the unsatisfactory performance of some specific models is probably caused by both the prediction uncertainty and the measurement noise. On the other hand, the measurement and prediction could be combined recursively in the modeling of HFMD transmission. Considering both the measurement noise and the prediction uncertainty can positively improve the estimation accuracy of the HFMD incidence, and could possibly offer a fresh perspective in exploring spatiotemporal patterns and determinant factors of the epidemic. This study aims to estimate the spatiotemporal evolution of the HFMD incidence by districts using a Kalman filter integrated with geographically weighted regression (GWR), to explore the spatiotemporal variation characteristics and potential risk regions, and to quantitatively evaluate the influence of meteorological and socioeconomic factors on the HFMD variation.
Methods
Study region
Shandong is an eastern coastal province of China and is located between 34° 23′ and 38° 24′ north latitude and between 114° 48′ and 122° 42′ east longitude (Fig. 1). It extends to the Yellow Sea in the east and is bordered by the Hebei, Henan, Anhui and Jiangsu provinces from northwest to southwest. The Shandong province has a total population of approximately 100.47 million and a total land area of 157,100 km^{2}. The gross domestic product (GDP) of Shandong province was 7646.97 billion Yuan in 2018. Shandong falls in the warm temperate monsoon climate zone, with an annual average temperature and precipitation in the ranges of 11–14 °C and 550–590 mm, respectively. More than 60% of the annual rainfall in the Shandong province is registered in the summer, and high temperatures usually occur in seasons with high precipitation.
Data
From May 1st, 2008 to March 19th, 2009 (47 weeks), weekly HFMD incidence reports for a total of 138 districts in Shandong were collected from the Chinese Centre for Disease Control and Prevention. To reduce the influence of population size, weekly incidence rates were calculated to reflect the risk of the HFMD epidemic for sample locations, and the corresponding Thiessen polygons were constructed to account for spatial effects (Fig. 1). Monthly meteorological data from May 2008 to March 2009 were obtained from the China National Meteorological Information Center (http://data.cma.cn/), including the daily average, maximum, and minimum temperatures (°C), the air pressure (hPa), relative humidity (%), wind speed (m/s), precipitation (mm) and sunshine hours (h). The socioeconomic data were collected from the 2008 statistic Yearbook of Shandong province, including GDP (10,000 Yuan), ratio of the number of primary school students to the total population (%) and number of hospital beds per capita. u_{1}–u_{8} and u_{9}–u_{11} are used to denote the above eight meteorological factors and three socioeconomic factors, respectively. Spatial Kriging methods were used to calculate the weekly average meteorological factors for each sample location during the 47week study of HFMD epidemics. Both dynamic meteorological factors and static socioeconomic factors were normalized to the range of 0–1.
Geographically weighted regression model
Compared with the global multivariate regression model, local models can be more effective at describing potential local variations in relationships between dependent and independent variables. The geographically weighted regression [35, 36] is a typical local multivariate regression model extensively applied to measure the spatial relationships between variables and corresponding local variations across an entire area. Moreover, GWR model can clearly detect and interpret any nonstationary features of spatial patterns and associations, and has been widely used to estimate the epidemic risk and assess the influence of the epidemic determinants [37, 38]. The GWR model used in this study is as follows:
where y_{i} is the HFMD incidence rate at location i with coordinates u_{i} and v_{i}, α (u_{i}, v_{i}) is the corresponding intercept constant, x_{k,i} are a series of independent variables describing local variations, β_{k} (u_{i}, v_{i}) are the local regression coefficients to be estimated, which vary with location, z_{l,i} are a series of independent variables connected with the global stability, γ_{l} are the corresponding static coefficients, and ε_{i} indicates the estimation error.
To approximate the HFMD incidence rate of each sample location in Shandong province, we take the dynamic meteorological factors as the local variables x_{k} in the above GWR model, and the static socioeconomic factors as the global variables z_{l}. Therefore, every location in the study area has a set of specific coefficients to reflect the associations between the HFMD incidence rate and the global or local variables. To solve the proposed GWR model, we apply a Gaussian distancedecay function to represent the relative importance between locations and an adaptive kernel scheme to determine the bandwidth (optimal number of neighboring locations), which is calculated through an iterative optimization process according to the Akaike Information Criterion (AIC). Meanwhile, the significance of the estimated global/local coefficients was checked with pseudo t tests and the model significance was tested by variance analysis (F tests).
Kalman filter
The Kalman filter (KF) is a data fusion algorithm initially designed to solve the discretedata linear filtering problem and provides a recursive solution to estimate the state variable of a timevarying system [39, 40]. In this study, KF is used to estimate the HFMD incidence and quantitatively assess the influence of risk factors. For a specific district, we define a multivariate state space X, which includes the HFMD incidence and several static socioeconomic factors. The state space is timevarying and calculated using the following parametric formula:
where X_{t} is the state vector containing the HFMD incidence and socioeconomic factors at time t, A is the state transition matrix indicating the effects of each state variable at time t1 on the state vector at time t, U_{t} is a vector containing control variables which are dynamic meteorological and static socioeconomic factors relevant to this study, B is the control coefficient matrix indicating the effects of each control variable on the state vector, and ω_{t} is a random variable representing the process noise, which is drawn from a zeromean Gaussian distribution N(0, Q). Last, Q stands for the prediction noise variance and accounts for the prediction uncertainty compared with the real process. The prediction of the timevarying state vector could be implemented as follows:
where \( {\hat{X}}_t \) is the prediction state vector at time t and X_{t1} is the estimated (filtered) state vector at time t1. The a priori estimation error covariance of the above prediction model propagates according to the equation:
where \( {\hat{P}}_t \) is the estimation error covariance of the prediction model at time t. Furthermore, by considering the HFMD incidence Y as the most important variable in the state vector X, we define a simple linear relationship linking the measurement Y to the state vector X:
where Y_{t} is the measurement HFMD incidence at time t which is the observed incidence calculated based on the reported cases, C is the observation operator matrix, and v_{t} is a random variable representing the measurement noise which is also assumed to be drawn from a zeromean Gaussian distribution N(0, R). Similarly, R stands for the measurement noise variance and represents the measurement uncertainty.
When both the process prediction and the measurement are considered, the a priori estimation error covariance \( {\hat{P}}_t \) and the measurement noise variance R are combined to generate the Kalman gain:
where K_{t} stands for the Kalman gain at time t and is applied to compute the a posteriori estimation of the state vector at time t as the following linear combination of the a priori estimation \( {\hat{X}}_t \) and the actual measurement Y_{t}:
As a function of the state vector covariance and the measurement noise, the Kalman gain K_{t} is noticeably high if the estimation error covariance is much higher than the measurement noise and the a posteriori estimation of the state vector significantly follows the measurements. Conversely, when K_{t} is low, the filter will essentially follow the predictions. In fact, K_{t} establishes the best combination between the process prediction and the measurement in order to minimize the mean square error between the a posteriori estimation X_{t} and its true value. After the update of the state vector as described above, the a posteriori estimation error covariance can be expressed as:
where I is an identity matrix and P_{t} indicates the estimation error covariance after the prediction and the update at time t. The a priori estimations take place at each step of the recursive solution based on the last a posteriori estimations, according to Eqs. (3) and (4), the Kalman gain at each step is computed according to Eq. (6), and the a posteriori estimations which are also the a priori estimations of the next step are generated according to Eqs. (7) and (8). Beginning from the initial state, the prediction and the update appear at every single step of the KF recursive solution.
Integration of the Kalman filter with the GWR model
Weekly averages of HFMD incidences in the sample locations were collected; the corresponding spatial autocorrelation was weak, with a Moran’s I of 0.0208 (p = 0.5460, calculated in ArcGIS Pro 2.3, https://pro.arcgis.com). However, the spatial stratified heterogeneity of the HFMD incidence among counties was statistically significant, with a GeoDetector qstatistic of 0.2153 (p < 0.001) [41, 42]. Therefore, GWR model was applied to explore the global or local associations between the HFMD incidence and meteorological or socioeconomic factors. Eight meteorological factors (u_{1}–u_{8}) were applied to be local (varying) variables and three socioeconomic factors (u_{9}–u_{11}) were used as global (fixed) variables. The GWR model to estimate the spatial distribution of HFMD incidences can be described as:
where y_{i} is the incidence at location i, α_{i} is the intercept, and ε_{i} indicates the estimation error. x_{k,i} denote the local factors (u_{1}–u_{8}) and β_{k,i} are the varying coefficients of local meteorological variables at location i. z_{l,i} indicate the global factors (u_{9}–u_{11}) and γ_{l} are the static coefficients of global socioeconomic variables. For 138 monitored districts of the studied area, coefficients γ_{l} are constant, and coefficients β_{k,i} are organized into a matrix composed by 138 rows and 8 columns respectively.
The GWR analysis was accomplished in GWR 4.0.90 software (https://gwrtools.github.io/) and produced an overall coefficient of determination R^{2} of 0.2482, which was only an approximately 14% improvement compared with the global regression prediction. No significant local coefficients were found with an alpha level of 0.05 (see more in Additional file 1). Whereas specific corrections can optimize the coefficient significance of pseudo t test (e.g., a correction to avoid false positives in GWR [43]), there might be false negatives in our results. These were possibly caused by the measurement noise in the HFMD incidence, as well as the prediction noise of the GWR model. To better explore the spatiotemporal patterns and assess the determinant factors of the HFMD epidemic, we combined the Kalman filter with the GWR model (Fig. 2). The filtering allows to couple the measured and predicted incidences, and improve the incidence estimation accuracy. On the other hand, GWR model indicates the associations between incidence and determinant factors, and therefore could provide the prediction modeling of state vector varying in the Kalman filter. Notice that, during the timevarying process of Kalman filtering, the measurement and prediction would be combined together to improve the estimation accuracy recursively, and thus, the prediction effectiveness of the state space model derived by GWR is not prerequired. Furthermore, the influence sensitivity of the control variables can be evaluated during the incidence filtering process, and the corresponding determinants of HFMD incidence can be quantitively assessed.
In our proposed Kalman filter, the multivariate state vector X is composed by HFMD incidence (y) and static socioeconomic factors (u_{9}–u_{11}). Y indicates the measured incidence and the parameter C is a simple observation operator matrix that indicates the transition between the state vector X and the measured incidence Y. Notice that, the measurement Y can be a vector consist of multiple explained variables (e.g. incidence and prevalence), and therefore, Kalman filter is effective to solve multivariate multiple regression problems.
The state transition matrix A models the variation of the state vector that consists of the HFMD incidence and the static socioeconomic factors from time t1 to time t. Coefficients γ_{l} derived from Eq. (9) indicate the global relationships between the incidence (y) and static factors (u_{9}–u_{11}). With a short 1week timestep, the socioeconomic variables are constant along with the recursive solution. Thus, matrix A can be easily organized by an identity matrix and a row vector of coefficients γ_{l}.
Both dynamic meteorological and static socioeconomic factors were selected as the vector containing controls in the Kalman filter. That is to say, vector U in Eq. (3) was consisted of meteorological factors (u_{1}–u_{8}) and socioeconomic factors (u_{9}–u_{11}). The control coefficient matrix B in Eq. (3) indicates the effects of each control variable on the state vector, and therefore, coefficients β_{k,i} derived from Eq. (9) could be applied to generate the parameter B in Eq. (3).
For different districts in the study area, the global and local effects of the determinant factors on the HFMD incidence vary spatially. Therefore, as shown in Fig. 2, we integrated the GWR model into the Kalman filter, derived the spacevarying parameters A and B, and generated multiple filters for the various districts (138 filters for 138 districts). The integration has two implications: one is the variable correspondence and the other is the parameter transitivity. The HFMD incidence was the explained variable in GWR model, as well as the measurement Y in the Kalman filter. The local and global explanatory variables in the GWR model were the meteorological and socioeconomic factors, which also constitutes the control vector U of the Kalman filter. Moreover, the state vector X in the Kalman filter contains the HFMD incidence and the socioeconomic factors. For each district, the global coefficients γ_{l} and the local coefficients β_{k,i}, which indicate the associations between the incidence and determinant factors, were obtained from the GWR result. Thus, the corresponding parameter A in the Kalman filter could be constructed from the global regression coefficients in the GWR model, while the parameter B using the local regression coefficients. Different from the parameters A and C, the control coefficient matrix B is districtdependent (various Bs for districts), and the corresponding multiple filters describe the spatial variation of the HFMD incidence evolution patterns and determinant influence effects.
During the recursive filtering process, the prediction \( {\hat{X}}_t \) at time t was calculated by the state space model with the estimation X_{t1} at time t1 according to Eq. (3), and the a priori estimation error covariance at time t was calculated with the a posteriori one at time t1 according to Eq. (4). The Kalman gain K_{t} at time t could be derived from the a priori estimation error covariance and the measurement noise variance according to Eq. (6). And the a posteriori estimation X_{t} at time t was updated according to Eq. (7), which indicates the estimated HFMD incidence of a specific district at time t. The corresponding a posteriori estimation error at time t could be calculated according to Eq. (8). After the timevarying recursion, the estimations and errors of HFMD incidences of all sample locations could be implemented.
Results
Kalman filtering validation
The HFMD incidence rates of 138 monitored districts were obtained in 47 weeks (from May 1st, 2008 to March 19th, 2009). For each of the sample districts continuous weekly incidence rates were available, the week index varying from 1 to 47. The average weekly incidence by district varied with time and had a mean value of approximately 0.936 × 10^{− 4} (in a range of 0.043 × 10^{− 4}–4.851 × 10^{− 4}). Eight meteorological factors (air pressure, daily average, maximum, and minimum temperatures, precipitation, relative humidity, wind speed and sunshine hours) were selected as the local dynamic independent variables (u_{1}–u_{8}), and the global static independent variables (u_{9}–u_{11}) were the following three socioeconomic factors: GDP, ratio of primary school students and number of hospital beds per capita. Both dynamic and static variables were normalized to the range of 0–1.
To evaluate the overall efficiency of the Kalman filter for HFMD incidence assessment, weekly incidence rates and meteorological variables for the studied districts were first aggregated to weekly average values. Next, using the static socioeconomic variables, the regression coefficients were calculated with the ordinary least squares (OLS) linear regression method. Subsequently, these coefficients were applied to generate the parameters B and C within the Kalman filter model, and the initial prediction and measurement errors were assumed to be drawn from a standardized Gaussian distribution. As shown in Fig. 3a, the filtering provided an adjustment to the weekly average HFMD incidences in the 138 districts to some extent compared to the corresponding measured values, and the estimated HFMD incidences followed a similar distribution as the measurements. Figure 3b illustrates that the original measurement errors varied among districts, highvalue errors correlating to districts with highvalue measurements; the estimation errors after filtering only apparently approach zero (the blue error curve presents an approximately horizontal line around the x axis). Even in districts with highvalue measured incidence, the Kalman filter satisfactorily reduces the estimation errors. The measurement and estimation errors of the HFMD incidences in districts are mapped in Fig. 4a. The HFMD incidence errors were reduced from the range of − 3.55 × 10^{− 4}–3.64 × 10^{− 4} to − 0.21 × 10^{− 4}–0.41 × 10^{− 4}. The Kalman filter significantly reduced the incidence errors for the majority of the districts, especially for those with large measurement errors. Figure 4b illustrates the reduced error distribution after filtering: although several districts received negative error reductions, the errors that increased were small and approximately 10% of the reduced ones. Regions with large error reductions and large HFMD incidences had similar reduced error distributions and were surrounded by regions with negative error reductions (the lightyellow polygons surround the darkgreen ones). Overall, the Kalman filter plays an effective role for HFMD incidence assessment even if the filter parameters are derived from the OLS linear regression without spatial variances. The measurement error covariance was 0.5686, whereas the estimation error covariance was substantially reduced to 0.0211 after filtering.
The spatiotemporal pattern of HFMD incidence filtering
After the overall validation of the Kalman filtering for the HFMD incidence assessment, we applied this model to explore the spatiotemporal patterns of HFMD incidences for all 138 districts. The local and global coefficients of dynamic meteorological factors (u_{1}–u_{8}) and static socioeconomic factors (u_{9}–u_{11}) on the HFMD incidences of each district were separately calculated using the GWR model. The corresponding parameter B of the Kalman filter is a matrix array that includes 138 control coefficient matrices ([11 × 4]), indicating the effects of meteorological and socioeconomic factors (u_{1}–u_{11}) on the state vector for 138 districts, respectively. A total of 138 Kalman filters with spatial variations were used to assess the temporal changes of HFMD incidences in the studied districts under the determinant factors (u_{1}–u_{11}). As shown in Fig. 5a, the average errors of measured incidences started with a high initial value and varied from week 1 to week 47; the error interval of 1 standard deviation (1StdDev) around the average showed local fluctuations, which are probably related to the abnormal temporal intervals of the HFMD incidence evolution. For instance, there was a tiny error increase that appeared in the 28th week (beginning on November 6th) accompanying a substantial interval expansion; the error intervals expanded significantly even when the error mean decreased to nearly zero in weeks 46 and 47 (beginning on March 12th). Figure 5b shows that, compared to the measurements, the error means and 1StdDev intervals of measurement incidences were reduced. However, considering the abovementioned temporal anomalies, even after filtering the error means and 1StdDev intervals were still large in the first 8 weeks (beginning on May 1st). That is to say, the HFMD epidemic in Shandong probably had pronounced seasonality features, usually evolving from midMarch, increasing until late June and with a potential reversal in early November.
To explore the spatial variation of the HFMD incidence filtering, the error covariances of incidence measurements and estimations by district were analyzed as shown in Fig. 6a. The majority of districts had satisfactory reductions of error covariances after filtering and several districts received noticeable reductions even when the original error covariances were large. However, the error covariances of several districts were still significant after filtering, and the Kalman filters played a weak role in these districts (their positions are indicated by red arrows). Figure 6b illustrates the spatial variation of the reduced error covariances by district after HFMD incidence filtering. The average error covariance of measured incidences was 0.3841, whereas the average estimated incidence error covariance was reduced to 0.1846, indicating an overall improvement of over 50% error reduction. However, several districts with significant error reductions overlapped to a certain extent with districts of large estimated incidence error covariances (Fig. 6). In other words, the HFMD incidence evolutions in these districts were abnormal, deeming such areas as potential risk regions of HFMD epidemic outbreaks.
Further considerations were proposed in these specific districts, and among them (Fig. 6b), error covariances of estimated incidences were classified in natural breaks and mapped in Fig. 7a. Judging by the temporal variations of the filtered HFMD incidence errors in each district, three classes of potential risk areas were distinguished, and presented separately in Fig. 7b, c and d, respectively. The temporal measurement curves and the estimation errors in two districts of the same class were extremely similar to each other. Although the spatial aggregation feature of these abnormal districts was weak (Fig. 7a), we could still classify potential HFMD risk regions into three categories by using the Kalman filter model in association with the meteorological and socioeconomic factors. As shown in Fig. 7b, the error curves of HFMD incidence filtering greatly varied in the early period but maintained a longterm steady trend. The second type of potential risk regions is illustrated in Fig. 7c; such regions present a relatively longterm steady trend with slight variations within a few intervals. Last, the third type had significant oscillations during the longterm period and unsteady oscillations appear in unpredictable localized time intervals (Fig. 7d). Evidently, the former two types of potential HFMD risk regions raise concerns during the localized periods, especially in HFMD highincidence seasons. Although the risk regions of the latter type were probably characterized by relatively low incidences, the HFMD epidemic evolutions were unsteady in the longterm, thus more prevention and control policies (e.g. longterm epidemic surveillance) should be implemented in these specific districts. Overall, the proposed HFMD incidence filtering in Shandong showed a strong seasonal dependence and several specific potential HFMD risk regions were found without significant spatial clustering.
Influence sensitivity of determinant factors
The control coefficient matrix B of the Kalman filter was generated from the GWR results to indicate the relationships between the HFMD incidence and meteorological or socioeconomic factors. To assess the influence that each factor has on the HFMD incidence, we defined an index ζ_{j} (j = 1–11) to describe the assumed enhancement effect of determinant factors (u_{1}–u_{11}). Experiments were repeated to evaluate the influence sensitivity of each dynamic or static factor on the HFMD incidence filtering. In experiment j, ζ_{j} varied from 0 to 5 with a step of 0.5, which indicates that the enhancement effect of factor u_{j} had a step size of 50% increase, while the ζ_{i} (i ≠ j) of other factors was kept invariant. The average errors and covariances of incidence estimations by district were applied to assess the influence sensitivity of meteorological and socioeconomic factors.
Figure 8a and b demonstrate the variations of the average estimation errors and covariances of HFMD incidence filtering along with the variation of each meteorological factor. As expected, the temperature factors (u_{2}–u_{4}) played the most important roles in the relationship with HFMD incidence filtering, and the average estimation errors and covariances were both sensitive to their enhancement effects, suggesting that higher temperature variations would cause a higher HFMD variation. Air pressure (u_{1}) was a secondary determinant affecting the HFMD variation approximately 25% as strongly as the temperature factors (Table 1). The next secondary determinants were sunshine hours, relative humidity, and precipitation (u_{8}, u_{6}, u_{5}). Compared to the latter rainfall factors, the effect of sunshine hours on the HFMD incidence variation was almost twice as much (Table 1). As shown in Fig. 8a and b, the wind speed (u_{7}) played a very weak role in HFMD incidence filtering, with a relative variation of nearly zero, reflecting that the HFMD epidemic probably had little airborne contagious transmission. Figure 8c and d illustrate the influence of socioeconomic factors on the HFMD incidence filtering. The number of hospital beds per capita (u_{11}) was the dominant determinant, followed by the GDP (u_{9}), which influenced the HFMD incidence approximately 30% as strongly as the dominant factor (Table 1). The relative variation of HFMD incidence filtering with the ratio of primary school students (u_{10}) was very slight (Table 1), suggesting that the amount of susceptible population in the studied region was probably not the leading cause of the HFMD variation. Overall, the daily average, maximum, and minimum temperatures and air pressure were the dominant meteorological factors, while the number of hospital beds per capita and GDP were the dominant socioeconomic ones that influenced the HFMD incidence variation in Shandong. Concomitantly, the HFMD variation was extremely slight even at high values of the wind speed and ratio of primary school students.
Discussion
In recent years, Kalman filters have been extensively used in a variety of applications, such as land cover classification [44] or landslide susceptibility evaluation [45]. Typical applications in Earth science concentrate on remote sensing image processing [44, 46,47,48,49] and data assimilations in the fields of agriculture [50,51,52], agrology [53, 54], ecology [55], hydrology [56, 57], oceanography [58] and others. In epidemiology, Kalman filters are usually applied to the mathematical modeling of epidemic spreads for diseases such as HIV/AIDS and Ebola [59,60,61]. In the present study, a Kalman filter was used to estimate the spatiotemporal evolution of HFMD incidence in 138 districts of the Shandong province, China, by integration with a GWR model to identify the local relationships between the HFMD incidence and risk factors. The proposed integrated model showed significant improvement in the HFMD incidence estimation accuracy. The spatiotemporal variation characteristics and potential risk regions of HFMD incidence were explored, and the influence of meteorological and socioeconomic factors on the HFMD variation were assessed. The results showed that the Kalman filter was effective for the HFMD incidence assessment in Shandong and produced a reduction of error covariance from 0.5686 to 0.0211 at the provincial scale. Considering the spatial variation of Kalman filters for various districts, the error covariance was reduced from 0.3841 to 0.1846 after filtering. Furthermore, filter processing allowed to identify potential HFMD risk regions: three categories of risk regions could be distinguished, with manifest filtering oscillations in the initial, local and longterm periods, respectively. Although the detected potential risk regions did not exhibit significant spatial clustering, more attention should be paid to these districts, especially the ones in the third category, with longterm filtering oscillations.
In addition to exploring the HFMD spatiotemporal patterns, the influence sensitivity of meteorological and socioeconomic factors was determined. We found that three temperature factors were the dominant meteorological determinants of the HFMD epidemic in Shandong, although the air pressure also affected the HFMD epidemic to a certain extent; however, wind speed had no manifest effect. Intense variations of temperature or air pressure produced high variations of HFMD incidence, whereas the influence of wind speed on the epidemic incidence was negligible and unclear. The HFMD related viruses are probably sensible to temperatures and air pressure. The main transmission routes of HFMD epidemic are the intimate contacts, and the wind speed influences slightly to the epidemic spread. With an overall shortage of rainfall in Shandong province, the HFMD epidemic spread might have less sensitivity to precipitation and relative humidity in a lowvalue level. Our findings are consistent with a number of previous studies [13, 18, 25, 26, 62, 63]. The environmental temperature relates to behavioral patterns such as increased contact among young children, thereby facilitating the spread of an HFMD infection [16]. However, our results indicate that meteorological factors such as precipitation, relative humidity and sunshine hours were not strongly associated with HFMD incidence, which is partially inconsistent with some of the previous studies. For instance, precipitation was strongly correlated with HFMD incidence in Singapore [25], and the number of HFMD cases increased significantly with increasing relative humidity in Japan [26]. HFMD cases at the county level across mainland China were spatially clustered and closely linked to the amounts of monthly precipitation in the region [24]. Relative humidity and precipitation were also found as the dominant driving factors of HFMD incidence in Henan, China [16]. Moreover, compared to GDP and ratio of primary school students to the total population, the number of hospital beds per capita appeared to be more dominant in HFMD incidence in Shandong. The children behavior patterns were possibly consistent amongst the districts with various economic levels around Shandong province. The influence of the population density background was already partially reduced during the calculation The HFMD incidence. The healthcare level played a manifest role in the controls of the HFMD epidemic spread in Shandong province. This result differs from other studies as well. For instance, GDP was the primary risk factor contributing to the spatial distribution of HFMD incidence in Sichuan and Henan, China [13, 16]. Possible reasons for this discrepancy include the differences between the studied regions, different transmission mechanisms of the HFMD epidemics, seasonal variations of meteorological factors, scale effects, zoning effects and others.
This study provides a multiperspective on estimating the spread of an HFMD epidemic by combining measurement noise with prediction uncertainty and demonstrates a novel approach to exploring the spatiotemporal patterns and determinant factors of an HFMD epidemic. Nevertheless, there are several limitations to this study, described as follows. First, we generated the basic local associations between the HFMD incidence and meteorological and socioeconomic factors using a GWR model without considering an HFMD mathematical model. Also, a limited number of driving factors were selected, which could have led to an insufficient description and interpretation of the HFMD epidemic dynamic mechanism. Second, our method was trained on countylevel data from the Shandong Province of China from 2008 to 2009, and applied only for the pattern exploration and risk assessment of the HFMD epidemic. This approach could easily be extended to other regions and infectious diseases similar to HFMD, although it should be accompanied by a thorough analysis and benchmarking of the model on the new problem. Lastly, it was hypothesized that the measurement and prediction noises of the Kalman filter followed a zeromean Gaussian distribution, and the model of the state vector and control variables was linear. These assumptions might have limited the applicability of the model; appropriate improvements could include nonlinear filters and nonGaussian noise distributions such as an extended Kalman filter (EKF), an unscented Kalman filter (UKF), or a particle filter (PF).
Conclusion
This study introduces a novel perspective to explore the spatiotemporal patterns and determinant factors of an HFMD epidemic. To this purpose, a Kalman filter method integrated with the GWR model with the aim to identify the global and local relationships between HFMD incidence and dynamic meteorological and static socioeconomic factors was designed. The proposed method considers both measurement noise and prediction uncertainty, which reduces the estimation error covariance of the HFMD incidence and improves the estimation accuracy. The filter processing could help explore the spatiotemporal patterns and determinants of the HFMD epidemic. As a result, three specific categories of potential risk regions of HFMD epidemics in Shandong were identified, with temperature factors and number of hospital beds per capita as the dominant determinants of the epidemic incidence. Furthermore, our approach can be extended to other regions and other infectious diseases similar to HFMD.
Availability of data and materials
The meteorological data used in the study are publicly available from the China National Meteorological Information Center (http://data.cma.cn/). The Socioeconomic data used in the study are publicly available from the statistical Yearbook of Shandong province (http://tjj.shandong.gov.cn/). The Shandong HFMD dataset is owned by the Chinese Centre for Disease Control and Prevention (http://www.chinacdc.cn/) and not available for distribution due to the constraint in the consent. The other dataset are available from the corresponding author on reasonable request.
Abbreviations
 CDC:

Chinese Center for Disease Control and Prevention
 EKF:

Extended Kalman filter
 GDP:

Gross domestic product
 GWR:

Geographically weighted regression
 HFMD:

Hand, foot and mouth disease
 KF:

Kalman filter
 NDVI:

Normalized difference vegetation index
 OLS:

Ordinary least squares
 PF:

Particle filter
 StdDev:

Standard deviation
 UKF:

Unscented Kalman filter
References
 1.
Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, et al. Hand, foot, and mouth disease in China, 200812: an epidemiological study. Lancet Infect Dis. 2014;14:308–18.
 2.
Zhu FC, Meng FY, Li JX, Li XL, Mao QY, Tao H, et al. Efficacy, safety, and immunology of an inactivated alumadjuvant enterovirus 71 vaccine in children in China: a multicentre, randomised, doubleblind, placebocontrolled, phase 3 trial. Lancet. 2013;381:2024–32.
 3.
Wang JF, Xu CD, Tong SL, Chen HY, Yang WZ. Spatial dynamic patterns of handfootmouth disease in the People’s republic of China. Geospat Health. 2013;7:381–90.
 4.
Ang LW, Koh BKW, Chan KP, Chua LT, James L, Goh KT. Epidemiology and control of hand, foot and mouth disease in Singapore, 20012007. Ann Acad Med Singap. 2009;38:106–12.
 5.
Chua KB, Kasri AR. Hand foot and mouth disease due to enterovirus 71 in Malaysia. Virol Sin. 2011;26:221.
 6.
Hosoya M, Kawasaki Y, Sato M, Honzumi K, Hayashi A, Hiroshima T, et al. Genetic diversity of coxsackievirus A16 associated with hand, foot, and mouth disease epidemics in Japan from 1983 to 2003. J Clin Microbiol. 2007;45:112–20.
 7.
Zhang Y, Tan XJ, Wang HY, Yan DM, Zhu SL, Wang DY, et al. An outbreak of hand, foot, and mouth disease associated with subgenotype C4 of human enterovirus 71 in Shandong, China. J Clin Virol. 2009;44:262–7.
 8.
Zhang Y, Zhu Z, Yang W, Ren J, Tan X, Wang Y, et al. An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol J. 2010;7:94.
 9.
Li Z, Lai S, Zhang H, Wang L, Zhou D, Liu J, et al. Hand, foot and mouth disease in China: evaluating an automated system for the detection of outbreaks. Bull World Health Organ. 2014;92:656–63.
 10.
Ji T, Han T, Tan X, Zhu S, Yan D, Yang Q, et al. Surveillance, epidemiology, and pathogen spectrum of hand, foot, and mouth disease in mainland of China from 2008 to 2017. Biosafety Health. 2019;1:32–40.
 11.
Guo C, Yang J, Guo Y, Ou QQ, Shen SQ, Ou CQ, et al. Shortterm effects of meteorological factors on pediatric hand, foot, and mouth disease in Guangdong, China: a multicity timeseries analysis. BMC Infect Dis. 2016;16:524.
 12.
Du Z, Lawrence WR, Zhang W, Zhang D, Yu S, Hao Y. Interactions between climate factors and air pollution on daily HFMD cases: a time series study in Guangdong, China. Sci Total Environ. 2019;656:1358–64.
 13.
Liao J, Qin Z, Zuo Z, Yu S, Zhang J. Spatialtemporal mapping of hand foot and mouth disease and the longterm effects associated with climate and socioeconomic variables in Sichuan Province, China from 2009 to 2013. Sci Total Environ. 2016;563:152–9.
 14.
Song C, Shi X, Bo Y, Wang J, Wang Y, Huang D. Exploring spatiotemporal nonstationary effects of climate factors on hand, foot, and mouth disease using Bayesian spatiotemporally varying coefficients (STVC) model in Sichuan, China. Sci Total Environ. 2019;648:550–60.
 15.
Huang X, Wei H, Wu S, Du Y, Liu L, Su J, et al. Epidemiological and etiological characteristics of hand, foot, and mouth disease in Henan, China, 20082013. Sci Rep. 2015;5:8904.
 16.
Xu C, Zhang X, Xiao G. Spatiotemporal decomposition and risk determinants of hand, foot and mouth disease in Henan, China. Sci Total Environ. 2019;657:509–16.
 17.
Wang J, Hu T, Sun D, Ding S, Carr MJ, Xing W, et al. Epidemiological characteristics of hand, foot, and mouth disease in Shandong, China, 20092016. Sci Rep. 2017;7:8900.
 18.
Li L, Qiu W, Xu C, Wang J. A spatiotemporal mixed model to assess the influence of environmental and socioeconomic factors on the incidence of hand, foot and mouth disease. BMC Public Health. 2018;18:274.
 19.
Wang Y, Feng Z, Yang Y, Self S, Gao Y, Longini IM, et al. Hand, foot, and mouth disease in China patterns of spread and transmissibility. Epidemiology. 2011;22:781–92.
 20.
Liu W, Ji H, Shan J, Bao J, Sun Y, Li J, et al. Spatiotemporal dynamics of handfootmouth disease and its relationship with meteorological factors in Jiangsu Province, China. PLoS One. 2015;10:e0131311.
 21.
Phung D, Nguyen HX, Nguyen HLT, Do CM, Tran QD, Chu C. Spatiotemporal variation of handfootmouth disease in relation to socioecological factors: a multipleprovince analysis in Vietnam. Sci Total Environ. 2018;610:983–91.
 22.
Zhang X, Xu C, Xiao G. Spacetime heterogeneity of hand, foot and mouth disease in children and its potential driving factors in Henan, China. BMC Infect Dis. 2018;18:638.
 23.
Zhao J, Hu X. The complex transmission seasonality of hand, foot, and mouth disease and its driving factors. BMC Infect Dis. 2019;19:521.
 24.
Wang J, Guo YS, Christakos G, Yang WZ, Liao YL, Li ZJ, et al. Hand, foot and mouth disease: spatiotemporal transmission and climate. Int J Health Geogr. 2011;10:25.
 25.
Hii YL, Rocklov J, Ng N. Short term effects of weather on hand, foot and mouth disease. PLoS One. 2011;6:e16796.
 26.
Onozuka D, Hashizume M. The influence of temperature and humidity on the incidence of hand, foot, and mouth disease in Japan. Sci Total Environ. 2011;410:119–25.
 27.
Cheng J, Wu J, Xu Z, Zhu R, Wang X, Li K, et al. Associations between extreme precipitation and childhood hand, foot and mouth disease in urban and rural areas in Hefei, China. Sci Total Environ. 2014;497:484–90.
 28.
Xie Y, Chongsuvivatwong V, Tang Z, McNeil EB, Tan Y. Spatiotemporal clustering of hand, foot, and mouth disease at the county level in Guangxi, China. PLoS One. 2014;9:e88065.
 29.
Shi RX, Wang JF, Xu CD, Lai SJ, Yang WZ. Spatiotemporal pattern of handfootmouth disease in China: an analysis of empirical orthogonal functions. Public Health. 2014;128:367–75.
 30.
Gui J, Liu Z, Zhang T, Hua Q, Jiang Z, Chen B, et al. Epidemiological characteristics and spatialtemporal clusters of hand, foot, and mouth disease in Zhejiang Province, China, 20082012. PLoS One. 2015;10:e0139109.
 31.
Wang Y, Lai Y, Du Z, Zhang W, Feng C, Li R, et al. Spatiotemporal distribution of hand, foot, and mouth disease in Guangdong Province, China and potential predictors, 20092012. Int J Environ Res Public Health. 2019;16:1191.
 32.
Yu G, Li Y, Cai J, Yu D, Tang J, Zhai W, et al. Shortterm effects of meteorological factors and air pollution on childhood handfootmouth disease in Guilin, China. Sci Total Environ. 2019;646:460–70.
 33.
Huang R, Ning H, He T, Bian G, Hu J, Xu G. Impact of PM10 and meteorological factors on the incidence of hand, foot, and mouth disease in female children in Ningbo, China: a spatiotemporal and timeseries study. Environ Sci Pollut Res. 2019;26:17974–85.
 34.
Zhang Q, Zhou M, Yang Y, You E, Wu J, Zhang W, et al. Shortterm effects of extreme meteorological factors on childhood hand, foot, and mouth disease reinfection in Hefei, China: a distributed lag nonlinear analysis. Sci Total Environ. 2019;653:839–48.
 35.
Brunsdon C, Fotheringham AS, Charlton ME. Geographically weighted regression: a method for exploring spatial nonstationarity. Geogr Anal. 1996;28:281–98.
 36.
Fotheringham AS, Charlton ME, Brunsdon C. Geographically weighted regression: a natural evolution of the expansion method for spatial data analysis. Environ Plan A. 1998;30:1905–27.
 37.
Hong Z, Hao H, Li C, Du W, Wei L, Wang H. Exploration of potential risks of hand, foot, and mouth disease in Inner Mongolia Autonomous Region, China Using Geographically Weighted Regression Model. Sci Rep. 2018;8:17707.
 38.
Hu M, Li Z, Wang J, Jia L, Liao Y, Lai S, et al. Determinants of the incidence of hand, foot and mouth disease in china using geographically weighted regression models. PLoS One. 2012;7:e38978.
 39.
Kalman RE. A new approach to linear filtering and prediction problems. J Basic Eng. 1960;82:35–45.
 40.
Kalman RE, Bucy RS. New results in linear filtering and prediction theory. J Basic Eng. 1961;83:95–108.
 41.
Wang JF, Li XH, Christakos G, Liao YL, Zhang T, Gu X, et al. Geographical detectorsbased health risk assessment and its application in the neural tube defects study of the Heshun region, China. Int J Geogr Inf Sci. 2010;24:107–27.
 42.
Wang JF, Zhang TL, Fu BJ. A measure of spatial stratified heterogeneity. Ecol Indic. 2016;67:250–6.
 43.
da Silva AR, Fotheringham AS. The multiple testing issue in geographically weighted regression. Geogr Anal. 2016;48:233–47.
 44.
Kleynhans W, Olivier JC, Wessels KJ, van den Bergh F, Salmon BP, Steenkamp KC. Improving land cover class separation using an extended Kalman filter on MODIS NDVI timeseries data. IEEE Geosci Remote Sens Lett. 2010;7:381–5.
 45.
Gorsevski PV, Jankowski P. An optimized solution of multicriteria evaluation analysis of landslide susceptibility using fuzzy sets and Kalman filter. Comput Geosci. 2010;36:1005–20.
 46.
Samain O, Roujean JL, Geiger B. Use of a Kalman filter for the retrieval of surface BRDF coefficients with a timeevolving model based on the ECOCLIMAP land cover classification. Remote Sens Environ. 2008;112:1337–46.
 47.
Garzelli A, Nencini F. Panchromatic sharpening of remote sensing images using a multiscale Kalman filter. Pattern Recogn. 2007;40:3568–77.
 48.
Salmon BP, Kleynhans W, Olivier JC, van den Bergh F, Wessels KJ. A modified temporal criterion to metaoptimize the extended Kalman filter for land cover classification of remotely sensed time series. Int J Appl Earth Obs Geoinf. 2018;67:20–9.
 49.
Kanakaraj S, Nair MS, Kalady S. Adaptive importance sampling unscented Kalman filter based SAR image super resolution. Comput Geosci. 2019;133:104310.
 50.
Li R, Li C, Dong Y, Liu F, Wang J, Yang X, et al. Assimilation of remote sensing and crop model for LAI estimation based on ensemble Kaiman filter. Agric Sci China. 2011;10:1595–602.
 51.
Zhao Y, Chen S, Shen S. Assimilating remote sensing information with crop model using ensemble Kalman filter for improving LAI monitoring and yield estimation. Ecol Model. 2013;270:30–42.
 52.
Huang J, Sedano F, Huang Y, Ma H, Li X, Liang S, et al. Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to improve regional winter wheat yield estimation. Agric For Meteorol. 2016;216:188–202.
 53.
Huang C, Li X, Lu L, Gu J. Experiments of onedimensional soil moisture assimilation system based on ensemble Kalman filter. Remote Sens Environ. 2008;112:888–900.
 54.
Gruber A, De Lannoy G, Crow W. A Monte Carlo based adaptive Kalman filtering framework for soil moisture data assimilation. Remote Sens Environ. 2019;228:105–14.
 55.
Chen C, Huang J, Chen Q, Zhang J, Li Z, Lin Y. Assimilating multisource data into a threedimensional hydroecological dynamics model using ensemble Kalman filter. Environ Model Softw. 2019;117:188–99.
 56.
Xie X, Zhang D. Data assimilation for distributed hydrological catchment modeling via ensemble Kalman filter. Adv Water Resour. 2010;33:678–90.
 57.
Zou L, Zhan C, Xia J, Wang T, Gippel CJ. Implementation of evapotranspiration data assimilation with catchment scale distributed hydrological model via an ensemble Kalman filter. J Hydrol. 2017;549:685–702.
 58.
Shu Y, Zhu J, Wang D, Xiao X. Assimilating remote sensing and in situ observations into a coastal model of northern South China Sea using ensemble Kalman filter. Cont Shelf Res. 2011;31:S24–36.
 59.
Cazelles B, Chau NP. Using the Kalman filter and dynamic models to assess the changing HIV/AIDS epidemic. Math Biosci. 1997;140:131–54.
 60.
Cobb L, Krishnamurthy A, Mandel J, Beezley JD. Bayesian tracking of emerging epidemics using ensemble optimal statistical interpolation. Spatial Spatio Temporal Epidemiol. 2014;10:39–48.
 61.
Ndanguza D, Mbalawata IS, Haario H, Tchuenche JM. Analysis of bias in an Ebola epidemic model by extended Kalman filter approach. Math Comput Simul. 2017;142:113–29.
 62.
Chang HL, Chio CP, Su HJ, Liao CM, Lin CY, Shau WY, et al. The association between Enterovirus 71 infections and meteorological parameters in Taiwan. PLoS One. 2012;7:e46845.
 63.
Wei J, Hansen A, Liu Q, Sun Y, Weinstein P, Bi P. The effect of meteorological variables on the transmission of hand, foot and mouth disease in four major cities of Shanxi Province, China: a time series data analysis (20092013). Plos Neglect Trop Dis. 2015;9:e0003572.
Acknowledgements
Not applicable.
Funding
This study was supported by the following grants: National Natural Science Foundation of China (41531179, 41421001, 41961055); National Key R&D Program of China (2016YFC1302504); Innovation Project of LREIS (O88RA200YA); Opening Fund of Key Laboratory of Poyang Lake Wetland and Watershed Research (Jiangxi Normal University), Ministry of Education (No. PK2019001). The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. All authors had access to the data, and had final responsibility for the decision to submit for publication.
Author information
Affiliations
Contributions
BH completed the design of the study, implemented the model integration, finished the data analysis and interpretation, and wrote the article. WQ conducted the data pretreatment and analysis, conducted the GWR and GeoDetector computing. CX prepared the HFMD and meteorological or socioeconomic data and participated in the analysis of the results and in the writing. JW participated in the modeling design and in the writing and accomplished the manuscript revision. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
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.
Details of the GWR analysis.
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
Hu, B., Qiu, W., Xu, C. et al. Integration of a Kalman filter in the geographically weighted regression for modeling the transmission of hand, foot and mouth disease. BMC Public Health 20, 479 (2020). https://doi.org/10.1186/s12889020086077
Received:
Accepted:
Published:
Keywords
 Hand
 Foot and mouth disease
 Kalman filter
 Geographically weighted regression
 Spatiotemporal pattern
 Determinant factors