Spatiotemporal Heterogeneity Detection
Between September 1, 2012, and August 31, 2013, a total of 8,014 cases in 53 weeks of bacillary dysentery in 103 counties were reported in Shandong Province. The highest number of cases occurred in summer (June to August), with a monthly incidence of 1.22 per 10,000 people. The lowest number of cases appeared in winter (December to February), with a monthly incidence of 0.35 per 10,000 people.
Geographically, the relative risks (RRs) differed dramatically, with the q statistic value of 0.51, indicating that there exists apparent spatial heterogeneity. Figure 3 presents the spatial RRs of bacillary dysentery by county level from 2012 to 2013. The spatial RRs of counties in eastern and northern Shandong mostly located in coast areas were higher, denoting that these counties have relatively higher bacillary dysentery risk. Conversely, counties in western and southern Shandong mainly belong to inland regions, presented relatively lower disease risk.
Additionally, the overall temporal trend presents an increase (Fig. 4), and meanwhile there remains seasonality, in which the highest disease risk occurred in summer (June to August) and the lowest disease risk occurred in winter (December to February), indicating that the risks of bacillary dysentery have obvious temporal heterogeneity, demonstrated by the q statistic value of 0.51.
Among the 103 counties in Shandong, 13 (12.62%) and 12 (11.65%) counties were perceived as hot and cold spots, respectively. Another 78 (75.73%) counties were considered to be neither hot nor cold spots. Figure 5 presents that hotspot areas were mainly distributed in the eastern and northern coast areas.
Risk Factor Analysis
The bacillary dysentery risk is obviously related to seasonal changes (Fig. 2 and 4), indicating that meteorological factors play a dominant role in the temporal evolution of the disease. In this study, the GeoDetector model was introduced to quantify individual meteorological factors and their interactive impacts. The result show that, the top factors with the higher determinant power are average temperature, precipitation and wind speed in three different regions.
Specifically, in middle mountainous areas, average temperature presents the most dramatically relationship with bacillary dysentery, with a q value of 0.66, it was the dominant factor explaining the temporal variation of the bacillary dysentery incidence (Table 1).
The other selected potential meteorological risk factors also have a non-negligible impact, such as, wind speed and precipitation had significant association with a higher extent of deviations, with q values of 0.48 and 0.22, respectively (Table 1).
The results of interactive impacts in GeoDetector denote that combined effect between randomly two meteorological factors also play an important role in the transmission of bacillary dysentery. Taking an example, in central regions, the determinant power of average temperature and relative humidity is 0.82, the determinant power of average temperature and sun hour is 0.82, and the determinant power of average temperature and precipitation is 0.81 (Table 1). Among these meteorological factors, comparing with their independent influence, all shows “bivariate enhance” effect.
Meanwhile, in western hilly areas, there also was a strong relationship between bacillary dysentery and average temperature, with the value of q being 0.47 (Table 2). Precipitation and wind speed have a similar determinant power with a higher extent of deviations, and q values were 0.32 and 0.32, respectively (Table 2).
The results of interaction in GeoDetector indicate that interactive impact between randomly two meteorological factors also play an important role in the transmission of bacillary dysentery. For example, in western regions, the determinant power of average temperature and wind speed is 0.64, the determinant power of wind speed and sun hour is 0.62, and the determinant power of sun hour and precipitation is 0.58 (Table 2). Comparing with their independent influence among these meteorological factors, all shows “bivariate enhance” effect.
Additionally, in eastern coastal areas, wind speed presents the most significant impact on the bacillary dysentery, with the value of q was 0.28. And then average temperature and sun hour also show obvious association with bacillary dysentery, with the q value was 0.25 and 0.22, respectively (Table 3).
And the results of interaction in GeoDetector present that interactive effect between randomly two meteorological factors play an important role in the transmission of bacillary dysentery. For example, in eastern regions, the determinant power of sun hour and wind speed is 0.71, the determinant power of wind speed and precipitation is 0.69, and the determinant power of sun hour and average temperature is 0.59 (Table 3). Comparing with their independent influence among these meteorological factors, also shows “bivariate enhance” effect. These results indicated that a hot and moist environment is more likely to promote the transmission of bacillary dysentery.