Laohasiriwong W, Puttanapong N and Luenam A. A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.12688/f1000research.12128.1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
1Faculty of Public Health and Research and Training Center for Enhancing Quality of Life for Working Age People, Khon Kaen University, Khon Kaen, Thailand 2Faculty of Economics, Thammasat University, Khlong Luang, Pathum Thani, Thailand 3Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
Background: The Centers for Disease Control and Prevention reported that deaths from chronic respiratory diseases (CRDs) in Thailand increased by almost 13% in 2010, along with an increased burden related to the disease. Evaluating the geographical heterogeneity of CRDs is important for surveillance. Previous studies have indicated that socioeconomic status has an effect on disease, and that this can be measured with variables such as night-time lights (NTLs) and industrial density (ID). However, there is no understanding of how NTLs and ID correlate with CRDs. We compared spatial heterogeneity obtained by using local cluster detection methods for CRDs and by correlating NTLs and ID with CRDs.
Methods: We applied the spatial scan statistic in SaTScan, as well as local indices of spatial association (LISA), Getis and Ord’s local Gi*(d) statistic, and Pearson correlation. In our analysis, data were collected on gender, age, household income, education, family size, occupation, region, residential area, housing construction materials, cooking fuels, smoking status and previously diagnosed CRDs by a physician from the National Socioeconomic Survey, which is a cross-sectional study conducted by the National Statistical Office of Thailand in 2010.
Results: According to our findings, the spatial scan statistic, LISA, and the local Gi*(d) statistic revealed similar results for areas with the highest clustering of CRDs. However, the hotspots for the spatial scan statistic covered a wider area than LISA and the local Gi*(d) statistic. In addition, there were persistent hotspots in Bangkok and the perimeter provinces. NTLs and ID have a positive correlation with CRDs.
Conclusions: This study demonstrates that all the statistical methods used could detect spatial heterogeneity of CRDs. NTLs and ID can serve as new parameters for determining disease hotspots by representing the population and industrial boom that typically contributes to epidemics.
Corresponding author:
Wongsa Laohasiriwong
Competing interests:
No competing interests were disclosed.
Grant information:
This research was financially supported by the Research and Training Center for Enhancing Quality of Life for Working Age People, Khon Kaen University, Thailand.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Chronic respiratory diseases (CRDs) are a public health issue worldwide1. According to the World Health Organization, CRDs are the world’s leading cause of death. In Thailand, the Centers for Disease Control and Prevention reported that deaths from CRDs increased by almost 13% in 2010, and there was an increased burden resulting from these diseases2.
Spatial clustering methods are important for statistical consideration, to develop models for prediction of disease risk positions. Disease risk positions are areas located near one another in statistically significant clusters. In this study, we indicate whether the geographical aggregation of CRDs is explained by chance or is statistically significant. It is likely that geographical areas share similar disease risk factors, since they are in a similar environment. Local spatial clustering methods are useful to identify the characteristics of clusters in terms of location, size, and prevalence of the disease4. These methods differ from standard methods3,4 used to identify the presence of spatial clustering in a whole study area. Local clustering methods in spatial epidemiology that are commonly applied to understanding spatial clustering (heterogeneity) are the Getis and Ord’s local Gi*(d) statistic, spatial scan statistics, and local indices of spatial association (LISA)5–8.
At present, there are no studies comparing the performance of these methods in epidemiological studies involving CRDs. In this study, we used local clustering methods to detect the spatial heterogeneity of CRDs and compared the performance of these local spatial clustering methods, including the local Gi*(d) statistic, the spatial scan statistic, and LISA. Spatial scan statistics performed in SaTScan circular version (version 9.4.2) were considered the best choice for detecting small, compact clusters9. A local version of Moran’s I was one of our selective cluster detection methods, due to its great performance on outlier detection, despite not being a good method for detecting large clusters3,10. The local Gi*(d) statistic was used to indicate locations surrounded by a cluster of high or low values. It shows areas where lower than average values tend to be found near each other, or areas where higher than average values tend to be found near each other, including the value at the location in which the spatial autocorrelation is being measured.
In a socioeconomic context, previous literature has recommended that night-time lights (NTLs) can substitute for other variables, such as urbanization, density, and economic growth. Data on NTLs have been used to indicate the actual scope of urban incorporations and appraise population scale in urban areas11–15 and population density16. These data have been employed to monitor the speed of urbanization17,18, as well as electricity used19. Moreover, NTLs have been utilized as a substitute for Gross Domestic Product (GDP), which works as an indirect indicator of administrative areas11,20–24. Remotely sensed NTL data can also be applied in epidemiological studies25,26. A 2011 study entitled “Explaining seasonal fluctuations of measles in Niger using night-time lights imagery” in Niger, by Bharti et al., found the concentration of the population affected by outbreaks of measles, and recommended that NTLs should be applied to public health studies27. Some studies in Israel have found that the incidence of breast and lung cancer in females is associated with NTL intensity in Israel28.
In Thailand, NTLs are used to estimate Gross Provincial Products (GPPs), because the luminosity of NTLs is normally dependent on the amount of economic activity in each area. The results showed that NTLs and GPP growth are significantly correlated and can represent the relationship between economic values and spatial inequality29. Industrial density (ID) may also be related to disease. However, previous studies have not investigated how NTLs and ID correlate with CRDs.
The government of Thailand has utilized various measures to attempt to prevent, diagnose, and treat CRDs. However, spatial clustering of diseases, particularly at a national level, has rarely been applied to examine characteristics of CRD clusters in terms of their location, size, and magnitude in Thailand. In addition, there has been limited use of socioeconomic indicators, such as NTLs and ID, in the epidemiology of CRDs.
In this study, commonly used methods for spatial clustering of CRDs were compared, including the use of the spatial scan statistic, LISA, and the local Gi*(d) statistic. We also analysed whether there was a correlation between NTLs and ID with CRDs. Our findings could help decision makers implement more relevant disease control policies using the performance of these local cluster detection methods, as well as the correlation of CRDs with NTLs and ID. The availability of a suitable strategy for spatial clustering of CRDs may allow for a better allocation of health care resources.
Methods
Study area
Thailand occupies an area of 514,000 square kilometres, which consists of an area of 511,770 square kilometres of land and 2,230 square kilometres of water. Thailand shares a border with Myanmar, Cambodia, Laos and Malaysia. In 2010, there were 76 provinces, 878 districts (Amphoe), 7,225 sub-districts (Tambon), and 74,965 villages (Mooban) in Thailand.
Ethical approval
The Ethical Committee of Khon Kaen University deemed this study to be exempt from requiring ethical approval (reference no. HE 582315).
Data collection and processing
Spatial data. The geographical coordinates of administrative areas of Thailand were retrieved from DIVA-GIS online (http://www.diva-gis.org/gdata), which is publicly available. This database of provincial maps of Thailand was processed with Quantum GIS for further GIS based analysis. For GIS based analysis, we applied zonal statistics to calculate mean NTLs and ID. The spatial distribution of CRD prevalence was obtained from the province-level polygon map, which contains information regarding the latitudes and longitudes of each province. The calculated mean NTLs, ID and CRD prevalence were further visualized with province-level layers on the polygon map and labelled with the administrative code of each province.
CRD, NTL and ID data. This study used data from the National Socioeconomics Survey (NSS), a cross-sectional study conducted by the National Statistical Office (NSO), Thailand, in 2010. Using stratified two-stage sampling, the survey selected a nationally representative sample to respond to a structured questionnaire. The questionnaires collected information on gender, age, household income, education, family size, occupation, region, residential area, housing construction materials, cooking fuels, smoking status and previous diagnosis of CRDs. A total of 17,040 individuals who met the inclusion criteria (aged between 18–59 years and diagnosed with having CRDs by a physician) were included in this analysis. The NSO administrative board officially allowed the research team to use the data (reference no.050601/1441).
The NTLs of Thailand used in this study were based on global stable lights imagery from 2010, acquired from the Operational Linescan System (OLS) sensor on-board satellite F18 under the Defense Meteorological Satellite Program (DMSP). All NTL data are publicly available (https://www.ngdc.noaa.gov/eog/dmsp/downloadV4composites.html). The ID data were calculated from GPPs, which represent the official statistics on provincial economic activity published by the National Economic and Social Development Board (NESDB). GPP data are also publicly available (http://www.nesdb.go.th/nesdb_en/more_news.php?cid=156&filename=index).
Local spatial pattern detection methods
The Open GIS software Quantum GIS (version 2.8.5) was used to conduct exploratory spatial data analysis30. Spatial autocorrelation analysis was conducted using GeoDa (version 1.6.6) and the local Gi*(d) statistic31. Stata version 10.0 (StataCorp, CollegeStation, TX) was used to calculate CRD prevalence. SaTScan (version 9.0.1) was used to determine the presence of statistically significant CRD spatial clusters and identify their approximate locations.
Local indices of spatial association. The local Moran’s I investigates the local level of spatial autocorrelation of provinces with a high and low prevalence of CRDs32. Computation of LISA assesses the local version of Moran’s I for each location to determine the variation in spatial autocorrelation over the study area. Its significance is evaluated in five categories: High-High, Low-Low, Low-High, High-Low, and Not Significant. Spatial autocorrelation occurs when a high prevalence of CRDs correlates with a high prevalence in neighbouring areas (also known as hotspots), or when a low prevalence of CRDs correlates with a low prevalence in neighbouring areas (also known as cold spots)33. This study set the spatial weight matrix of 3 k-Nearest Neighbours around each province, meaning that the clustering relation was identified with three neighbouring provinces. The statistical significance level was 0.05, and the simulation used 999 permutations34 to evaluate the sensitivity of the results.
Getis and Ord’s local Gi*(d) statistic. The local Gi*(d) statistic is used to test the statistical significance of local clusters and to determine the spatial dependence and relative magnitude between observations. In this study, the local Gi*(d) statistic was used to test statistically significant CRD prevalence and local autocorrelation, as well as to determine the dependence of neighbouring observation35 and local clustering. We used k-Nearest Neighbours, and contiguity was set as three provinces for a polygon contiguity spatial weight matrix, which was created depending on three provinces that share common boundaries and vertices. A 0.05 significance level and 999 permutations were used to identify significant clusters of local autocorrelation. A high value for the local Gi*(d) statistic indicates a hotspot, and a low value indicates a cold spot36,37.
Spatial scan statistic. The spatial scan statistic, calculated in SaTScan38, was applied to analyse the geographic distribution of CRD prevalence in Thailand and determine whether there are any spatial geographical clusters of CRD prevalence as for high or low CRDs. We used the retrospective purely spatial Poisson model to identify spatially high or low clusters of CRDs, since it fit the assumption that the number of events in an area was Poisson distributed according to a known underlying population at risk. Purely spatial analysis, which ignores the time dimension of the cases, was performed to detect CRD clusters in the study areas. The purely spatial Poisson model imposed circular windows on the map and let the circle move over the area. Each circle moved until reaching a maximum population at risk. The radius of each circle increased continuously from to a maximum radius. Therefore, the window never encompassed more than 50% of the total population at risk, i.e., the maximum spatial cluster size, in this investigation. The likelihood function was maximized over all window positions and sizes, and the one with the maximum likelihood constituted the most likely cluster. For each window, we used a Monte Carlo simulation to test the null hypothesis that no significant clusters were created, supposing that the relative risk (RR) of CRDs was no different within the window compared to outside the window. The maximum number of random Monte Carlo replications was defined as 999. The likelihood function was maximized in the overall window locations and sizes39. The window with the maximum likelihood value of the spatial window was set to 50% of the population at risk.
Results
Pearson correlation assessed the relationship between NTLs and ID with CRDs in all provinces. There was a moderate positive correlation of 0.34 (95% CI: 0.12-0.52; p-value 0.002) and 0.36 (95% CI: 0.15-0.54; p-value 0.001) with CRDs, which indicated that NTLs and ID had similar relationships to CRDs (Table 1 and Figure 1).
Table 1. Pearson correlation of NTLs and ID with CRDs, and 95% confidence interval (CI) for each factor in 2010.
Factors
CRDs
95% CI
p-value
NTLs
0.34
0.12-0.52
0.002
ID
0.36
0.15-0.54
0.001
Figure 1. Night-time lights (NTLs) and industrial density (ID) maps for Thailand, 2010.
(A) NTLs and (B) ID showed high values, representative of high population density.
The Univariate Moran’s I scatter of annual CRD prevalence among provinces in Thailand in 2010 showed a slightly positive spatial autocorrelation, as the Moran's I was 0.072, with statistical significance at 0.05. The Moran’s I indicated clustering patterns. There were hotspots in Bangkok and the perimeter provinces, namely, Samut Prakan, Nonthaburi, Pathum Thani and Nakhon Pathom. Cold spots and low prevalence areas, were located in the Nan and Phitsanulok provinces. Two high/low clusters were located in Trang, Songkhla, and part of the Phatthalung provinces, indicating that the distribution of the CRD in affected provinces was spatially autocorrelated (low clustering), even though the overall tendencies were not clear. This could be due to the nature of local Moran's I indices; the statistic cannot accommodate extreme values or outliers in the data set, such as the very low CRD prevalence in some provinces that were also included in the model. Therefore, its sensitivity in cluster detection may be reduced. Using results from other local spatial analysis methods could help explain the autocorrelations.
Considering NTLs and ID with CRDs in terms of bivariate Moran’s I scatter for annualized prevalence in 2010, the results showed a positive spatial autocorrelation, with Moran's I values of 0.345 and 0.374, respectively. These values were statistically significant, with a significance level of 0.05. Based on these results, hotspots were located in Bangkok and perimeter provinces, consisting of Samut Prakan, Nonthaburi, Pathum Thani and Nakhon Pathom. The results showed a high, significant correlation between the NTLs with CRDs, suggesting that brighter lights, together with high population densities of exclusively wealthier people, can be used to define the spatial extent of CRD development. Therefore, ID, which reflects urbanization, was correlated with CRDs (Figure 3 and Figure 4).
For the local Gi*(d) statistics, hotspots were located in Bangkok and perimeter provinces, consisting of Samut Prakan, Nonthaburi, Pathum Thani and Nakhon Pathom (Figure 5).
Figure 2.
(A) LISA (Univariate) cluster map, and (B) Moran scatter plot matrix of CRDs (p < 0.05) in 2010. There were 3 hotspots located in Bangkok and the perimeter provinces.
Figure 3.
(A) LISA (Bivariate; Night-time light with CRDs) cluster map, and (B) Moran scatter plot matrix of CRDs (p < 0.05) in 2010. We analysed NTLs with CRDs in terms of bivariate Moran’s I scatter for annualized prevalence. There were 5 hotspots located in Bangkok and the perimeter provinces.
Figure 4.
(A) LISA (bivariate; industrial density with CRDs) cluster map, and (B) Moran’s I scatter plot matrix of CRDs (p < 0.05) in 2010. We analysed ID with CRDs in terms of bivariate Moran’s I scatter for annualized prevalence. There were 5 hotspots located in Bangkok and the perimeter provinces.
Figure 5. Local Gi*(d) Cluster map of CRDs in 2010.
There were 5 hotspots located in Bangkok and the perimeter provinces.
The spatial scan statistics revealed two spatial clusters. The primary cluster was located in 23 provinces. There were 157 cases compared to 213 expected cases. Thus, the ratio between observed and expected cases was 0.74. The p-value was smaller than 0.05, which indicated that the cluster was highly significant. The relative risk for the population inside the cluster compared to the population outside the cluster was 0.65, indicating that the risk of CRDs within this area was higher than locations outside it. This number is the estimated risk within the cluster, divided by the estimated risk outside the cluster. The analysis also revealed that a secondary cluster, which is the most likely cluster, was located in 32 provinces. The number of observed cases was 348, compared to 297 expected cases. Thus, the ratio between observed and expected cases was 1.17. According to these calculations, the relative risk inside the cluster is 1.37. The p-value was smaller than 0.05, indicating that the cluster was highly significant. The relative risk for the population inside the cluster compared to the population outside the cluster was 1.37, indicating that the risk of CRDs in this area was higher than locations outside it (Figure 6).
Figure 6. Spatial scan statistic calculated in SaTScan, to detect CRD prevalence in 2010.
The spatial scan statistic, calculated in SaTScan, was used to detect spatial high or low CRD prevalence. We used the retrospective, purely spatial Poisson model. Based on these results, the spatial scan statistics revealed two spatial clusters.
In summary, as seen in Figure 1–Figure 6, spatial scan statistics, LISA, and local Gi*(d) statistics revealed similar results for areas with the highest clustering of CRDs. Nevertheless, the clusters for spatial scan statistics covered a wider area of the Bangkok, Central, Northeast and South clusters, while the clusters for LISA and the local Gi*(d) statistic were wider and more dispersed in Bangkok and the perimeter provinces.
Discussion
Comparison of the local spatial pattern detection methods
In this study, NTL and ID data were utilized to assess growth in provinces based on geography, economic growth, urbanization, or health metrics, as well as concentration of NTLs and ID. Areas with a high CRD prevalence also experienced high light growth and a high concentration of ID, and therefore NTLs and ID were probably induced CRDs. In other words, the results indicated that CRDs are affected by population density and industrial concentration. This finding was consistent with a previous study that stated that it is possible to use NTLs in public health studies26. The findings of this study suggest that NTLs and ID can substitute for some variables, such as urbanization, density, economic growth, and industrial concentration. NTLs and ID serve as a new tool for specifying disease hotspots by representing the population and industrial booms typically contributing to epidemics. In urban areas with migration, NTLs can indicate where populations are clustering, through the captured expansion and the increase brightness of lighted areas. The technique indicates fluctuations in population density, which affect epidemic risk and can replace current methods of outbreak surveillance26,40.
LISA, spatial scan statistics, and local Gi*(d) statistics are common methods applied to investigate local clusters of diseases. These techniques help identify clusters and assess their statistical significance41. In this study, all of these methods provided comparable results in the detection of geographic areas of CRDs in both high and low- rate clusters. However, there were some inconsistencies. The identification of clusters by spatial scan statistics was more localized, compared to LISA and local Gi*(d) statistics, possibly due to their dissimilar role in defining clusters. LISA detected the neighbouring provinces with rates significantly correlated to each other, and areas with similar values were thus determined as a cluster32. In this study, LISA explored the correlation between the value of a certain area and the average value of neighbouring areas. Spatial scan statistic methods were used, and larger clusters were detected; the calculation of the maximum likelihood ratio of cases was made together with the underlying population in the area to investigate the cluster of provinces considered to have higher or lower rates. Using spatial scan statistics methods, larger clusters were detected, and the calculation of maximum likelihood ratio of cases was made simultaneously with the underlying population in the area to investigate the cluster of provinces and consider them for higher or lower rates7,42. Meanwhile, the local Gi*(d) statistic and LISA reported that the two methods shared similar results in terms of rising clusters (hotspots). Hanson and Wieczorek7 compared LISA to spatial scan statistics by exploring alcohol-related mortality in New York, USA (which identified some differences between two methods). It was concluded that the spatial scan statistic is the more sensitive method. A multi-method approach is adopted for cluster detection, since different methods tend to identify different characteristics in clusters, i.e., LISA identifies the core of the cluster, and the scan statistic identifies its extent. Furthermore, Jacquez and Greiling43 used LISA to analyse spatial clustering for diagnosis of breast, lung and colorectal cancers in Long Island, USA, and to identify significant spatial patterns for all of these diseases. Their analysis confirms the clustering of breast cancer mortality previously revealed by Kulldorff’s spatial scan statistic44, but they found that the two methods identified slightly different clusters. As a result, they recommended that a combination of statistics be used when studying local clustering to assure that different aspects of spatial patterns are fully identified and that the results from the suite of analyses are logically consistent.
Due to its own principle for cluster detection, each method selected different geographic areas. The intersection of the selections indicated that these three methods resulted in different features in the same clusters, suggesting that various methods should be applied to identify clusters of CRDs at the provincial level in Thailand. These spatial cluster detection methods can be utilized collectively rather than individually. Furthermore, because each method has its own advantages and disadvantages, no single method is deemed the “gold standard” for cluster analysis45.
The strengths of this study included the methods used to analyse CRD cases, including spatial scan statistics, LISA, and local Gi*(d) statistics. All of these have been thoroughly investigated and compared to determine the appropriate methods for evaluating spatial clustering and cluster detection. NTLs and ID were included to investigate a correlation with CRDs. This study contained a large sample size that represents a nationally representative sample. Therefore, the results could represent the Thai population in general.
Abbreviations: CRDs: Chronic respiratory diseases, NSO: National Statistical Office, NSS: National Socioeconomics Survey, LISA: Local Indices of Spatial Association, ID: Industrial density, NTLs: Night time lights.
Dataset 1.Pearson correlation with CRDs dataset.
Dataset 2.NTLs and ID dataset.
Dataset 3.LISA (Univariate) cluster CRDs dataset.
Dataset 4.LISA (Bivariate; Night-Time Light with CRDs) cluster dataset.
Dataset 5.LISA (Bivariate; Industrial Density with CRDs) cluster dataset.
Dataset 6.Local Gi*(d) cluster CRDs dataset.
Dataset 7.Spatial scan statistics in SaTScan analysis of CRDs dataset.
Conclusions
This study utilized a geographic information system and spatial analyses, which can be applied to several epidemiological studies to investigate and clarify the spatial heterogeneity of CRDs in highly affected provinces in Thailand. LISA, spatial scan statistics, and local Gi*(d) statistics were calculated, with the aim of revealing the spatial characteristics of the CRDs. Spatial scan statistics performed much better in outlier detection in terms of power, compared to LISA and local Gi*(d) statistics. Based on our simulation, there was a large relative risk difference (two spatial clusters) and significant spatial heterogeneity. The findings suggested heterogeneity in the spatial pattern of CRDs in the study. Purely spatial retrospective analysis revealed persistence of CRD clusters in some geographical locations of the provinces on an annual basis. CRDs were affected by the concentration of NTLs and ID. We believe that NTLs and ID serve as a new tool for specifying disease hotspots by representing the population and industrial booms that typically contribute to epidemics. The results of the study provide helpful information on the common epidemiological situation of CRDs. Demonstrating the existence of CRDs hotspots in different provinces may enable the Ministry of Public Health or provincial health officers to launch remedial measures in the affected areas and formulate strategies for more effective handling of CRDs.
Data availability
Data used in this study were from the NSS; permission to use these data can be requested from the NSO. This study got approval from the NSS (reference no.050601/1441) to use the data on gender, age, household income, education, family size, occupation, region, residential area, housing construction materials, cooking fuels, smoking status and previously diagnosed CRDs by a physician. NTL data and ID data are publicly available.
This research was financially supported by the Research and Training Center for Enhancing Quality of Life for Working Age People, Khon Kaen University, Thailand.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgements
The authors are grateful to all contributors to this research, and especially to the NSO for the data.
Faculty Opinions recommended
References
1.
World Health Organization: Global surveillance, prevention and control of chronic respiratory diseases: a comprehensive approach. GARD_Manual. 2007; World Health Organization: Geneva, Switzerland, 2007. Reference Source
2.
Centers for disease control and prevention Prevention: Global Health-Thailand. 2011, [cited 2017 Jan 26]. Reference Source
3.
Waller LA, Gotway CA:
Applied spatial statistics for public health data. John Wiley & Sons; 2004. Publisher Full Text
4.
Pfeiffer DU, Robinson TP, Stevenson M, et al.:
Spatial Analysis in Epidemiology. Oxford University Press. 2008; 142. Publisher Full Text
6.
Kulldorff M, Tango T, Park PJ:
Power comparisons for disease clustering tests.
Comput Stat Data Anal.
2003; 42(4): 665–84. Publisher Full Text
7.
Hanson CE, Wieczorek WF:
Alcohol mortality: a comparison of spatial clustering methods.
Soc Sci Med.
2002; 55(5): 791–802. PubMed Abstract
| Publisher Full Text
8.
Fotheringham AS, Zhan FB:
A comparison of three exploratory methods for cluster detection in spatial point patterns.
Geogr Anal.
1996; 28(3): 200–18. Publisher Full Text
9.
Kulldorff M, Huang L, Pickle L, et al.:
An elliptic spatial scan statistic.
Stat Med.
2006; 25(22): 3929–43. PubMed Abstract
| Publisher Full Text
10.
Kulldorff M:
Tests of Spatial Randomness Adjusted for an Inhomogeneity.
J Am Stat Assoc.
2006; 101(475): 1289–305. Publisher Full Text
11.
Elvidge CD, Baugh KE, Kihn EA, et al.:
Mapping city lights with nighttime data from the DMSP Operational Linescan System.
Photogramm Eng Remote Sensing.
1997; 63(6): 727–34. Reference Source
12.
Elvidge CD, Baugh KE, Kihn EA, et al.:
Relation between satellite observed visible-near infrared emissions, population, economic activity and electric power consumption.
Int J Remote Sens.
1997; 18(6): 1373–9. Publisher Full Text
13.
Imhoff ML, Lawrence WT, Stutzer DC, et al.:
A technique for using composite DMSP/OLS “city lights” satellite data to map urban area.
Remote Sens Environ.
1997; 61(3): 361–70. Publisher Full Text
14.
Elvidge CD, Baugh KE, Dietz JB, et al.:
Radiance calibration of DMSP-OLS low-light imaging data of human settlements.
Remote Sens Environ.
1999; 68(1): 77–88. Publisher Full Text
15.
Small C, Elvidge CD, Balk D, et al.:
Spatial scaling of stable night lights.
Remote Sens Environ.
2011; 115(2): 269–80. Publisher Full Text
16.
Zhuo L, Ichinose T, Zheng J, et al.:
Modelling the population density of China at the pixel level based on DMSP/OLS non‐radiance‐calibrated night‐time light images.
Int J Remote Sens.
2009; 30(4): 1003–18. Publisher Full Text
17.
Zhang Q, Seto KC:
Mapping urbanization dynamics at regional and global scales using multi-temporal DMSP/OLS nighttime light data.
Remote Sens Environ.
2011; 115(9): 2320–9. Publisher Full Text
18.
Elvidge CD, Baugh KE, Anderson SJ, et al.:
The Night Light Development Index (NLDI): a spatially explicit measure of human development from satellite data.
Soc Geogr.
2012; 7(1): 23–35. Publisher Full Text
19.
Townsend AC, Bruce DA:
The use of night-time lights satellite imagery as a measure of Australia's regional electricity consumption and population distribution.
Int J Remote Sens.
2010; 31(16): 4459–80. Publisher Full Text
20.
Florida R, Gulden T, Mellander C:
The rise of the mega-region.
Econ Soc.
2008; 1(3): 459–76. Publisher Full Text
22.
Doll CH, Muller JP, Elvidge CD:
Night-time imagery as a tool for global mapping of socioeconomic parameters and greenhouse gas emissions.
J Human Environ.
2000; 29(3): 157–62. Publisher Full Text
23.
Sutton PC, Elvidge CD, Ghosh T:
Estimation of gross domestic product at sub-national scales using nighttime satellite imagery.
Int J Ecol Econ Stat.
2007; 8(S07): 5–21.
Reference Source
24.
Florida R, Mellander C, Gulden T:
Global metropolis: assessing economic activity in urban centers based on nighttime satellite images.
Prof Geogr.
2012; 64(2): 178–87. Publisher Full Text
25.
Li D, Zhao X, Li X:
Remote sensing of human beings–a perspective from nighttime light.
Geo Spat Inf Sci.
2016; 19(1): 69–79. Publisher Full Text
27.
Bharti N, Tatem AJ, Ferrari MJ, et al.:
Explaining seasonal fluctuations of measles in Niger using nighttime lights imagery.
Science.
2011; 334(6061): 1424–7. PubMed Abstract
| Publisher Full Text
| Free Full Text
28.
Kloog I, Haim A, Stevens RG, et al.:
Light at night co-distributes with incident breast but not lung cancer in the female population of Israel.
Chronobiol Int.
2008; 25(1): 65–81. PubMed Abstract
| Publisher Full Text
29.
Chaiwat T:
Night Lights, Economic Growth, and Spatial Inequality of Thailand. Puey Ungphakorn Institute for Economic Research. 2016.
Reference Source
30.
Steiniger S, Hunter AJ:
The 2012 free and open source GIS software map–A guide to facilitate research, development, and adoption.
Comput Environ Urban Syst.
2013; 39: 136–50. Publisher Full Text
31.
Anselin L, Syabri I, Kho Y:
GeoDa: an introduction to spatial data analysis.
Geogr Anal.
2006; 38(1): 5–22. Publisher Full Text
32.
Anselin L:
Local Indicators of Spatial Association—LISA.
Geogr Anal.
1995; 27(2): 93–115. Publisher Full Text
33.
Anselin L, Syabri I, Kho Y:
GeoDa: An introduction to spatial data analysis. 2004. Reference Source
34.
Anselin L:
An introduction to EDA with GeoDa. Spatial Analysis Laboratory (SAL) Department of Agricultural and Consumer Economics, University of Illinois, Urbana-Champaign, IL. 2003. Reference Source
35.
Hinman SE, Blackburn JK, Curtis A:
Spatial and temporal structure of typhoid outbreaks in Washington, D.C., 1906–1909: evaluating local clustering with the Gi* statistic.
Int J Health Geogr.
2006; 5(1): 13. PubMed Abstract
| Publisher Full Text
| Free Full Text
36.
Getis A, Ord JK:
The analysis of spatial association by use of distance statistics.
Geogr Anal.
1992; 24(3): 189–206. Publisher Full Text
37.
Getis A, Ord JK:
Local spatial statistics: an overview. Spatial analysis: modelling in a GIS environment. 1996; 374. Reference Source
38.
Kulldorff M:
A spatial scan statistic. In: Communications in Statistics - Theory and Methods. 1997; 26(6): 1481–96. Publisher Full Text
39.
Kulldorff M:
SaTScanTM User Guide for version 8. 0. 2009.
40.
Princeton University: Nighttime images help track disease from the sky.
ScienceDaily.
2011, [cited 2017 Jan 30]. Reference Source
41.
Anselin L:
Review of Cluster Analysis Software. The North American Association of Central Cancer Registries, Inc. 2004. Reference Source
42.
Nødtvedt A, Guitian J, Egenvall A, et al.:
The spatial distribution of atopic dermatitis cases in a population of insured Swedish dogs.
Prev Vet Med.
2007; 78(3–4): 210–22. PubMed Abstract
| Publisher Full Text
44.
Kulldorff M, Feuer EJ, Miller BA, et al.:
Breast cancer clusters in the northeast United States: a geographic analysis.
Am J Epidemiol.
1997; 146(2): 161–70. PubMed Abstract
| Publisher Full Text
45.
Sasson C, Cudnik MT, Nassel A, et al.:
Identifying High-risk Geographic Areas for Cardiac Arrest Using Three Methods for Cluster Analysis.
Acad Emerg Med.
2012; 19(2): 139–46. PubMed Abstract
| Publisher Full Text
46.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 1 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
47.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 2 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
48.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 3 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
49.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 4 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
50.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 5 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
51.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 6 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
52.
Laohasiriwong W, Puttanapong N, Luenam A:
Dataset 7 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand.
F1000Research.
2017. Data Source
1
Faculty of Public Health and Research and Training Center for Enhancing Quality of Life for Working Age People, Khon Kaen University, Khon Kaen, Thailand 2
Faculty of Economics, Thammasat University, Khlong Luang, Pathum Thani, Thailand 3
Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
This research was financially supported by the Research and Training Center for Enhancing Quality of Life for Working Age People, Khon Kaen University, Thailand.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Laohasiriwong W, Puttanapong N and Luenam A. A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.12688/f1000research.12128.1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.
Share
Open Peer Review
Current Reviewer Status:
?
Key to Reviewer Statuses
VIEWHIDE
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations
A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Martinez EZ. Reviewer Report For: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.13126.r29646)
The authors present interesting data about the possible associations between socioeconomic status (night-time lights and industrial density) and chronic respiratory diseases. The manuscript is well written, the methods and the results sections are adequate, and should be of great interest
... Continue reading
The authors present interesting data about the possible associations between socioeconomic status (night-time lights and industrial density) and chronic respiratory diseases. The manuscript is well written, the methods and the results sections are adequate, and should be of great interest to the readers. Overall, it is an important study, but their limitations are not mentioned. For example, ecological studies do not permit cause-and-effect conclusions or evidence that risk-modifying factors are related to risk of the disease, since they do not link individual exposure histories to individual outcome events. Thus, I believe the authors can provide a paragraph on the limitations of the study.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.
Martinez EZ. Reviewer Report For: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.13126.r29646)
Amornrat Luenam, Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
06 Mar 2018
Author Response
There are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second
...
Continue readingThere are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second limitation is caused by the characteristic of data and analytical method that do not permit the cause-and-effect conclusion, particularly the connection between individual exposure histories and individual outcome events.
There are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second limitation is caused by the characteristic of data and analytical method that do not permit the cause-and-effect conclusion, particularly the connection between individual exposure histories and individual outcome events.
Competing Interests:No competing interests were disclosed.Close
Amornrat Luenam, Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
06 Mar 2018
Author Response
There are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second
...
Continue readingThere are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second limitation is caused by the characteristic of data and analytical method that do not permit the cause-and-effect conclusion, particularly the connection between individual exposure histories and individual outcome events.
There are two aspects of limitations. The first is based on the availability of data. Therefore, there are limited resolutions in both spatial and time dimensions of data. The second limitation is caused by the characteristic of data and analytical method that do not permit the cause-and-effect conclusion, particularly the connection between individual exposure histories and individual outcome events.
Competing Interests:No competing interests were disclosed.Close
Tipayamongkholgul M. Reviewer Report For: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.13126.r28274)
In the methods, the author used SatScan and used 50% of population at risk to be a radius at the province level. Using such technique at the provincial level (large geographic area) may lead to inaccurate cluster
... Continue reading
In the methods, the author used SatScan and used 50% of population at risk to be a radius at the province level. Using such technique at the provincial level (large geographic area) may lead to inaccurate cluster as you can see in the result that one cluster cover the whole region.
In the result, Moran’ I present very low spatial correlation (0.07), although it presented statistically significant. The correlational value cannot accept the spatial correlation that exists, and therefore all following spatial data analysis may not be necessary.
LISA and Local G were used for different purpose, I don’t think comparison between this two is needed.
The number of CRD may relate with the accessibility to healthcare service in Bangkok metropolitan and some surrounding area. The author may explain whether this point may affect the study result.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Partly
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Partly
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
Competing Interests: No competing interests were disclosed.
I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.
Tipayamongkholgul M. Reviewer Report For: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand [version 1; peer review: 2 approved with reservations]. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.13126.r28274)
Amornrat Luenam, Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
06 Mar 2018
Author Response
Explanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study
...
Continue readingExplanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study to locate potential circular cluster areas based on the recommendation of Kulldorff M. Kulldorff M (2016) suggested that 50 percent is the default maximum, a smaller maximum may be requested. However, the recommended value is 50 percent (Kulldorff M. 2018, Kulldorff M, Nagarwalla.1995) which reported that a window sized up to 50% of the population at risk can generally reduce negative clusters. Whereas 90%, is more appropriately interpreted as a lower disease rate in the 10% of the area outside the ‘cluster’ rather than as an excess disease rate covering almost the whole study region, therefore, it did not appropriate with our study. References: Abdurrob A., Kulldorff M. SaTScan Tutoria # 3 Advanced Options. 2016. Available: https://www.satscan.org/tutorials/advancedoptions/SaTScanTutorialAdvancedOptions.pdf Kulldorff M, 2018. SaTScanTM User Guide for version 9.5 Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Statistics in medicine. 1995;14(8):799–810.
Explanation to question 2 and 3
We conducted Global Moran I, G* test and LISA based on contents and suggestions introduced by Anselin (1995). Specifically, all three tests were undertaken because of the following reasons as clarified in Anselin (1995).
(1) Global Moran I has a limitation. The statistical insignificance of Global Moran I might not indicate that the local association does not exist (Getis and Ord 1992, p. 201). Therefore, the methods of G* and LISA have been developed in order to overcome this incomplete capability of Global Moran I, and both are capable of identifying the local associations. (2) Based on statistical inference, G* can identify only two cases of localized associations, which are hotspots of high-high and cold spots of low-low. However, LISA has extended this capability, and it can indicate four categories of hotspots, which are high-high, high-low, low-high and low-low. Specifically, LISA can identify both cases of positive spatial correlation (e.g. high-high and low-low) and those of negative spatial association (e.g. high-low and low-high) Based on these clarifications, this paper included results obtained from Global Moran I, G* test, and LISA in order to generate the robust results, statically examining spatial associations for both aspects - the global vs. local and univariate vs. bivariate.
References: Anselin, L. Local Indicators of Spatial Association—LISA. Geographical Analysis. 1995; 27: 93–115. Getis, A., and K. Ord. ‘The Analysis of Spatial Association by Use of Distance Statistics.” Geographical Analysis. 1992; 24:189-206.
Explanation to question 4 We study the CRDs morbidity which is mainly linked to the epidemic of tobacco exposure and indoor and outdoor air pollution in Asian countries (Tan WC, Ng TP, 2008), whereas the mortality may relate with the accessibility to healthcare service in Bangkok metropolitan and some surrounding area. References: Tan WC, Ng TP. COPD in Asia: where East meets West. Chest. 2008; 133(2):517-527.
Explanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study to locate potential circular cluster areas based on the recommendation of Kulldorff M. Kulldorff M (2016) suggested that 50 percent is the default maximum, a smaller maximum may be requested. However, the recommended value is 50 percent (Kulldorff M. 2018, Kulldorff M, Nagarwalla.1995) which reported that a window sized up to 50% of the population at risk can generally reduce negative clusters. Whereas 90%, is more appropriately interpreted as a lower disease rate in the 10% of the area outside the ‘cluster’ rather than as an excess disease rate covering almost the whole study region, therefore, it did not appropriate with our study. References: Abdurrob A., Kulldorff M. SaTScan Tutoria # 3 Advanced Options. 2016. Available: https://www.satscan.org/tutorials/advancedoptions/SaTScanTutorialAdvancedOptions.pdf Kulldorff M, 2018. SaTScanTM User Guide for version 9.5 Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Statistics in medicine. 1995;14(8):799–810.
Explanation to question 2 and 3
We conducted Global Moran I, G* test and LISA based on contents and suggestions introduced by Anselin (1995). Specifically, all three tests were undertaken because of the following reasons as clarified in Anselin (1995).
(1) Global Moran I has a limitation. The statistical insignificance of Global Moran I might not indicate that the local association does not exist (Getis and Ord 1992, p. 201). Therefore, the methods of G* and LISA have been developed in order to overcome this incomplete capability of Global Moran I, and both are capable of identifying the local associations. (2) Based on statistical inference, G* can identify only two cases of localized associations, which are hotspots of high-high and cold spots of low-low. However, LISA has extended this capability, and it can indicate four categories of hotspots, which are high-high, high-low, low-high and low-low. Specifically, LISA can identify both cases of positive spatial correlation (e.g. high-high and low-low) and those of negative spatial association (e.g. high-low and low-high) Based on these clarifications, this paper included results obtained from Global Moran I, G* test, and LISA in order to generate the robust results, statically examining spatial associations for both aspects - the global vs. local and univariate vs. bivariate.
References: Anselin, L. Local Indicators of Spatial Association—LISA. Geographical Analysis. 1995; 27: 93–115. Getis, A., and K. Ord. ‘The Analysis of Spatial Association by Use of Distance Statistics.” Geographical Analysis. 1992; 24:189-206.
Explanation to question 4 We study the CRDs morbidity which is mainly linked to the epidemic of tobacco exposure and indoor and outdoor air pollution in Asian countries (Tan WC, Ng TP, 2008), whereas the mortality may relate with the accessibility to healthcare service in Bangkok metropolitan and some surrounding area. References: Tan WC, Ng TP. COPD in Asia: where East meets West. Chest. 2008; 133(2):517-527.
Competing Interests:No competing interests were disclosed.Close
Amornrat Luenam, Faculty of Public Health, Khon Kaen University, Khon Kaen, Thailand
06 Mar 2018
Author Response
Explanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study
...
Continue readingExplanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study to locate potential circular cluster areas based on the recommendation of Kulldorff M. Kulldorff M (2016) suggested that 50 percent is the default maximum, a smaller maximum may be requested. However, the recommended value is 50 percent (Kulldorff M. 2018, Kulldorff M, Nagarwalla.1995) which reported that a window sized up to 50% of the population at risk can generally reduce negative clusters. Whereas 90%, is more appropriately interpreted as a lower disease rate in the 10% of the area outside the ‘cluster’ rather than as an excess disease rate covering almost the whole study region, therefore, it did not appropriate with our study. References: Abdurrob A., Kulldorff M. SaTScan Tutoria # 3 Advanced Options. 2016. Available: https://www.satscan.org/tutorials/advancedoptions/SaTScanTutorialAdvancedOptions.pdf Kulldorff M, 2018. SaTScanTM User Guide for version 9.5 Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Statistics in medicine. 1995;14(8):799–810.
Explanation to question 2 and 3
We conducted Global Moran I, G* test and LISA based on contents and suggestions introduced by Anselin (1995). Specifically, all three tests were undertaken because of the following reasons as clarified in Anselin (1995).
(1) Global Moran I has a limitation. The statistical insignificance of Global Moran I might not indicate that the local association does not exist (Getis and Ord 1992, p. 201). Therefore, the methods of G* and LISA have been developed in order to overcome this incomplete capability of Global Moran I, and both are capable of identifying the local associations. (2) Based on statistical inference, G* can identify only two cases of localized associations, which are hotspots of high-high and cold spots of low-low. However, LISA has extended this capability, and it can indicate four categories of hotspots, which are high-high, high-low, low-high and low-low. Specifically, LISA can identify both cases of positive spatial correlation (e.g. high-high and low-low) and those of negative spatial association (e.g. high-low and low-high) Based on these clarifications, this paper included results obtained from Global Moran I, G* test, and LISA in order to generate the robust results, statically examining spatial associations for both aspects - the global vs. local and univariate vs. bivariate.
References: Anselin, L. Local Indicators of Spatial Association—LISA. Geographical Analysis. 1995; 27: 93–115. Getis, A., and K. Ord. ‘The Analysis of Spatial Association by Use of Distance Statistics.” Geographical Analysis. 1992; 24:189-206.
Explanation to question 4 We study the CRDs morbidity which is mainly linked to the epidemic of tobacco exposure and indoor and outdoor air pollution in Asian countries (Tan WC, Ng TP, 2008), whereas the mortality may relate with the accessibility to healthcare service in Bangkok metropolitan and some surrounding area. References: Tan WC, Ng TP. COPD in Asia: where East meets West. Chest. 2008; 133(2):517-527.
Explanation to question 1. We used the SaTScan, using the Poisson probability model and constrained to clusters no larger than 50% of the population at risk and 50% of this study to locate potential circular cluster areas based on the recommendation of Kulldorff M. Kulldorff M (2016) suggested that 50 percent is the default maximum, a smaller maximum may be requested. However, the recommended value is 50 percent (Kulldorff M. 2018, Kulldorff M, Nagarwalla.1995) which reported that a window sized up to 50% of the population at risk can generally reduce negative clusters. Whereas 90%, is more appropriately interpreted as a lower disease rate in the 10% of the area outside the ‘cluster’ rather than as an excess disease rate covering almost the whole study region, therefore, it did not appropriate with our study. References: Abdurrob A., Kulldorff M. SaTScan Tutoria # 3 Advanced Options. 2016. Available: https://www.satscan.org/tutorials/advancedoptions/SaTScanTutorialAdvancedOptions.pdf Kulldorff M, 2018. SaTScanTM User Guide for version 9.5 Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Statistics in medicine. 1995;14(8):799–810.
Explanation to question 2 and 3
We conducted Global Moran I, G* test and LISA based on contents and suggestions introduced by Anselin (1995). Specifically, all three tests were undertaken because of the following reasons as clarified in Anselin (1995).
(1) Global Moran I has a limitation. The statistical insignificance of Global Moran I might not indicate that the local association does not exist (Getis and Ord 1992, p. 201). Therefore, the methods of G* and LISA have been developed in order to overcome this incomplete capability of Global Moran I, and both are capable of identifying the local associations. (2) Based on statistical inference, G* can identify only two cases of localized associations, which are hotspots of high-high and cold spots of low-low. However, LISA has extended this capability, and it can indicate four categories of hotspots, which are high-high, high-low, low-high and low-low. Specifically, LISA can identify both cases of positive spatial correlation (e.g. high-high and low-low) and those of negative spatial association (e.g. high-low and low-high) Based on these clarifications, this paper included results obtained from Global Moran I, G* test, and LISA in order to generate the robust results, statically examining spatial associations for both aspects - the global vs. local and univariate vs. bivariate.
References: Anselin, L. Local Indicators of Spatial Association—LISA. Geographical Analysis. 1995; 27: 93–115. Getis, A., and K. Ord. ‘The Analysis of Spatial Association by Use of Distance Statistics.” Geographical Analysis. 1992; 24:189-206.
Explanation to question 4 We study the CRDs morbidity which is mainly linked to the epidemic of tobacco exposure and indoor and outdoor air pollution in Asian countries (Tan WC, Ng TP, 2008), whereas the mortality may relate with the accessibility to healthcare service in Bangkok metropolitan and some surrounding area. References: Tan WC, Ng TP. COPD in Asia: where East meets West. Chest. 2008; 133(2):517-527.
Competing Interests:No competing interests were disclosed.Close
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations -
A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 1 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178700)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 2 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178701)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 3 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178702)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 4 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178703)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 5 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178704)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 6 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178705)
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
How to fix it
Save downloaded CSV file
Open spreadsheet program (e.g. Excel)
Click the ‘Data’ tab at the top
Click the ‘From text’ icon (top left)
Browse for downloaded CSV file, click ‘Import’
Ensure ‘Delimited’ radio button is selected, click ‘Next’
Check one of the appropriate delimiter checkboxes (you can visualize the formatting by looking at the data preview below these options)
Laohasiriwong W, Puttanapong N and Luenam A. Dataset 7 in: A comparison of spatial heterogeneity with local cluster detection methods for chronic respiratory diseases in Thailand. F1000Research 2017, 6:1819 (https://doi.org/10.5256/f1000research.12128.d178706)
Adjust parameters to alter display
View on desktop for interactive features
Includes Interactive Elements
View on desktop for interactive features
Competing Interests Policy
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Examples of 'Non-Financial Competing Interests'
Within the past 4 years, you have held joint grants, published or collaborated with any of the authors of the selected paper.
You have a close personal relationship (e.g. parent, spouse, sibling, or domestic partner) with any of the authors.
You are a close professional associate of any of the authors (e.g. scientific mentor, recent student).
You work at the same institute as any of the authors.
You hope/expect to benefit (e.g. favour or employment) as a result of your submission.
You are an Editor for the journal in which the article is published.
Examples of 'Financial Competing Interests'
You expect to receive, or in the past 4 years have received, any of the following from any commercial organisation that may gain financially from your submission: a salary, fees, funding, reimbursements.
You expect to receive, or in the past 4 years have received, shared grant support or other funding with any of the authors.
You hold, or are currently applying for, any patents or significant stocks/shares relating to the subject matter of the paper you are commenting on.
Stay Updated
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Comments on this article Comments (0)