Association between Daily Hospital Outpatient Visits for Accidents and Daily Ambient Air Temperatures in an Industrial City

An accident is an unwanted hazard to a person. However, accidents occur. In this work, we search for correlations between daily accident rates and environmental factors. To study daily hospital outpatients who were admitted for accidents during a 5-year period, 2007–2011, we analyzed data regarding 168,366 outpatients using univariate regression models; we also used multivariable regression models to account for confounding factors. Our analysis indicates that the number of male outpatients admitted for accidents was approximately 1.31 to 1.47 times the number of female outpatients (P < 0.0001). Of the 12 parameters (regarding air pollution and meteorology) considered, only daily temperature exhibited consistent and significant correlations with the daily number of hospital outpatient visits for accidents throughout the 5-year analysis period. The univariate regression models indicate that older people (greater than 66 years old) had the fewest accidents per 1-degree increase in temperature, followed by young people (0–15 years old). Middle-aged people (16–65 years old) were the group of outpatients that were more prone to accidents, with an increase in accident rates of 0.8–1.2 accidents per degree increase in temperature. The multivariable regression models also reveal that the temperature variation was the dominant factor in determining the daily number of outpatient visits for accidents. Our further multivariable model analysis of temperature with respect to air pollution variables show that, through the increases in emissions and concentrations of CO, photochemical O3 production and NO2 loss in the ambient air, increases in vehicular emissions are associated with increases in temperatures. As such, increases in hospital visits for accidents are related to vehicular emissions and usage. This finding is consistent with clinical experience which shows about 60% to 80% of accidents are related to traffic, followed by accidents occurred in work place.


Introduction
Safety is a top priority when a person participates in any form of daily activity. We are fortunate if we are safe every day. However, accidents do occur. Is there a background trend that reveals a hidden fact about all accidents that occur? Ambient temperature (heat) has been demonstrated to be related to aggression, with an increase in temperature leading to an increased occurrence of human violence [1]. In professional baseball, a positive and significant relationship between temperature and the number of batters hit by pitches per game has been observed; this result suggests that higher temperatures cause pitchers to become more aggressive when pitching to batters [2]. A direct linear increase in horn honking with increasing temperature has also been observed [3]. However, a study of accident risk in a Swedish town concluded that high temperature (and rain) do not increase the risk of accidents for low-speed buses [4].
In this work, we used the daily number of hospital outpatient visits for accidents and daily temperature data to study the correlation between the occurrence of accidents and temperature in an industrial area in northern Taiwan during the period of 2007-2011. The use of daily hospital medical records ensures better data quality over a long period [5]. This long studying period enabled us to obtain and test the statistical significance of our results. In the regression model analysis presented below, we also included and analyzed correlations between the number of hospital visits for accidents and other environmental factors (meteorological factors in addition to temperature and air pollution factors, which have been demonstrated to have health effects).

Ethics Statement
The statistics regarding the daily number of hospital visits used in this work are openly published each year by the hospital [6][7][8][9]. Only the frequencies of visits for diseases were used in the study. No details of patient's personal information were involved in this work. Consequently, neither ethics committee approval nor written consent was required.
Outpatient data from weekends (Saturdays and Sundays), national holidays (e.g., the Chinese New Year), and natural hazards (e.g., closure of public offices due to typhoons) were excluded from the analysis to eliminate artificial effects of zero outpatient visits for accidents during these periods. Because hospital emergency rooms are open continuously, hospital visits for accidents outside the normal working hours are admitted through emergency rooms. Consequently, we also conducted regression model analyses of all emergency visits for accidents. eng/). The highest population densities are in the 4 districts (Tao-Yuan, 416,403 people, 11,965 people km −2 ; Chung-Li, 380,829 people, 4978 people km −2 ; Ping-Jen, 212,073 people, 4441 people km −2 ; Bade, 135,885 people, 4031 people km −2 ) located between the Tao-Yuan and Chung-Li ambient air monitoring stations (Fig 1). These 4 districts contain more than 1 million people (1,145,190), and account for 55% of entire population in Tao-Yuan City. Most people live in areas with flat terrain, which are also characterized by high densities of industrial emissions (Fig 1). Hence, except for exceptional industrial emissions (such as industrial volatile organic compounds, hydrogen sulfide, etc), monitoring data from Chung-Li station is representative of Tao-Yuan City. For example, the R 2 values for temperatures between Chung-Li and Tao  Spatial domain of Tao-Yuan City (marked by red colored curve), and terrain heights (shaded colors, in the unit of meters) distributed over northern Taiwan. Geographical locations of the Landseed Hospital (red corss); Chung-Li (a blue cross close to the red corss), Tao-Yuan (a blue cross on the right), and other ambient air monitoring stations (blue circles); and stationary air pollution sources (black crosses).
The hourly ambient data included the concentrations of 7 air pollutants and 5 meteorological parameters. The 7 air pollutants concentrations included those of particulate matter with particle diameters of less than 10 μm (PM 10 ), particulate matter with particle diameters of less than 2.5 μm (PM 2.5 ), ozone (O 3 ), carbon monoxide (CO), nitrogen oxide (NO), nitrogen dioxide (NO 2 ), and sulfur dioxide (SO 2 ). The 5 meteorological parameters were temperature, precipitation, relative humidity, wind speed, and wind direction [10].
In a manner similar to the analysis of daily outpatient visits, data from the weekend and national holidays were excluded. This ensures that the time-series data regarding outpatient data, air pollution, and meteorological data covered consistent time periods. For the regression model analysis of hospital emergency visits, no public holidays were excluded because emergency rooms were open at all times. All data were analyzed for the period of 2007-2011. Air pollution and meteorological data used in this work can be openly accessed on the website of the Envrionmental Protection Administrtion (Taiwan) at http://www.epa.gov.tw.

Test of Normal Distribution for Data
Data Distribution. The time-series data regarding hospital outpatient visits for accidents, air pollutants, and meteorological parameters were tested for normality. The procedures used to test the data regarding air pollutants, meteorological factors, and outpatient visits for normality are similar [10]. Test of normality for air pollution and temperature data were presented in a previous work [10]. The following are the steps that were used to test the data regarding outpatient visits for accidents for normality. The daily outpatient visit data for each age group were binned into 15 bins. For the outpatients with ages of 0-15 years and greater than 66 years old, the bin interval was 3, and the bin ranges were from 0 to 45 visits; for the 16-65 years-old outpatients, the bin interval was 12, and the bin ranges were from 0 to 180 visits. We then computed histograms of days (y-axis) versus daily outpatient visits binned into these 15 bins (xaxis). With these histograms determined, we computed standard deviations and mean values for each age group in each year of the period of 2007-2011.
Theoretical Normal Distribution. The mean value and standard deviation for an age group of outpatients in a year determined from the above method was then used to compute the theoretical normal distribution for the selected 15 bin intervals.
Test of Data Normality. Finally, with the data distribution and theoretical normal distribution determined, we could test the data for normality using a method such as the chi-squared test [11]. We used the chi-squared test for determining the goodness of fit to determine whether there was any difference between the observed (real-world) data and the expected (theoretical) data. The binned observed data were from the medical, meteorological, and air pollution data. The expected (theoretical) data were computed using the theoretical normal distribution with the mean and standard values determined from observations [11].

Linear Regression Model Analysis
The processed time-series data regarding the daily outpatient visits were combined with each of the air pollutant-related and meteorological parameters for univariate regression model analysis to calculate correlation coefficients R [11]. Because there were 5 years of data , and 3 age groups of outpatients for each year, a total of 15 univariate regression model calculations (3 age groups per year x 5 years of data) were conducted for each air pollutionrelated and meteorological parameter.
For each pair of the outpatient-air pollution time-series data, a total of 15 correlation coefficients R were determined via 15 univariate regression model calculations. For each pair of the outpatient-meteorology time-series data, a total of 15 correlation coefficients R were determined via another 15 univariate regression model calculations. A total of 7 air pollutants required 105 univariate regression model calculations; and a total of 5 meteorological factors required 75 univariate regression model calculations. Hence, a total of 180 univariate regression model calculations were performed, and 180 correlations coefficients R were produced. These correlation coefficients R provided a good quantitative measure for analysis and comparison. The significances of each correlation coefficient R were measured using a P-value, which was obtained using the chi-squared test [11]. The univariate regression model was written in the Fortran language [12]. The model calculations were performed on the computers running an open source Linux operating system, Fedora Core [13]. The univariate regression model can also be run on the Windows operation system using the open source software package Cygwin [14]. The Cygwin package is a marvelous software package that enables development of Linux code in a Windows environment on a notebook computer. A detailed description of the algorithms and example code that were used in the univariate regression model calculations presented in this work can be found in reference 11.

Multivariable Model Analysis
Multivariable regression models, which were based on the principles described in reference [15], were developed [10] and used to model the variation in the number of outpatient visits with respect to air pollutants and meteorological parameters for each year of the time period of 2007-2011. There were 7 air pollutant concentrations and 5 meteorological parameters (temperature, precipitation, wind direction, wind speed, and relative humidity) considered in the multivariable model analysis. The 7 air pollutants were PM 10 , PM 2.5 , O 3 , CO, NO, NO 2 , and SO 2 . Table 1 presents a list of the multivariable regression models used in this work. Given the input time-series data of daily outpatient visits, air pollution variables, and meteorological variables, one can determine a correlation coefficient β i of a variable i for each multivariable model presented in Table 1 using a linear algebra method [11]. Here i = 1 indicates air pollution variable PM 10 , i = 2 indicates air pollution variable PM 2.5 , and so forth. The first 7 variables are related to air pollution factors, and the rest of the 5 variables are related to meteorological factors.
As for the univariate regression model described previously, the multivariable regression models were written in the Fortran programming language on a Linux platform running the open source operating system Fedora Core and can also be run using a Windows platform running the open source Cygwin package on a notebook computer. The Fedora Core operating  system, Fortran compiler, and Cygwin package can be freely downloaded from open sources on the Web. Reference [15] provides plenty of useful examples and data, which are ideal for testing the multivariable regression model code developed in Fortran language [10]. Once the accuracy of these codes is verified, these multivariable regression models will be powerful tools for big data analysis, which requires efficient handling of the heterogeneous data sets that can easily exceed the capacity of commercial software packages. In this work, there were 3 age groups of outpatients for a period of 5 years and a list of 12 variables for each age group of outpatients in each year to be modeled. Hence, a total of 180 (3 age groups x 5 years x 12 variables) multivariable regression model calculations were conducted to produce the analysis results presented in this work. These results include the coefficients β i for each variable i, F values, R 2 , and P-values to measure the statistical significance of the multivariable regression model calculations. Table 2 presents a list of annual outpatient visits for accidents during 2007-2011. The frequency of visits varied between 32,189 and 34,770. These data consistently demonstrated that there were more male accident outpatients (ranging between 18,271 and 20,049) than female outpatients (ranging between 13,610 and 15,038) each year. The number of male accident outpatients was approximately 1.31 to 1.47 times the number for females. More than 90% of patients that visited Landseed Hospital were from Tao-Yuan County (now called Tao-Yuan City) [6][7][8][9]. The population in the Tao-Yuan County is approximately equally divided between males and females, i.e., the male-to-females ranged between 1.03 to 1.01 during the study [ Table 3]. Even accounting for the fact that accident attendance rates depend also on the distance from an accident and emergency department [16], the gender distribution of outpatient visits for accidents differed significantly from the population distribution (P < 0.0001; Oneway ANOVA test, http://in-silo.net/tools/statistics/anova). One factor that could possibly cause more males than females to have accidents is that males tend to engage in more risky behaviors and get hurt more often than females [17]. Alcohol abuse was reported as a primary cause of accidents [18], and males may exhibit more frequent alcohol abuse than females. Moreover, there are more male than female drivers in the high-traffic urban areas that are typical of Tao-Yuan City, and more curbside parking can increase the risk of injury of child pedestrians by motor vehicles [19]. A study regarding injuries and the risk of disability in teenagers and young adults in England, Scotland, and Wales also observed that more males (62%) than females (26%) reported at least one accident that resulted in hospital treatment since age 16 [20]. The seasonal trends in temperature, with low temperatures occurred during the winter season (December, January, and February), and high temperatures occurred during the summer season (June, July, and August), are readily identified from the plots. The temperatures varied between 12°C and 38°C. In addition to the pronounced seasonal variations in temperatures, day-to-day fluctuations in temperatures were also evident in the data.

Outpatient Visits During 2007-2011
In contrast with the sharp seasonal temperature variations, no distinctive seasonal variations in the number of outpatients for accidents were determined from the outpatient data. The number daily outpatient visits for accidents exhibited large day-to-day variations. Fewer accidents occurred in older people (greater than 65 years old) than the young people (0-15 years old and 16-65 years old). These results are consistent with the results of a previous work [21]. indicate that the outpatient data and the respective theoretical normal distribution were significantly different [11]. Our calculations yielded large vales of PROB, thereby indicating that outpatients for accidents were distributed nearly normally.  The linear equations indicated that the basic number of outpatients per day was approximately 14 to 17 per day when the temperature approached 0°C. For the outpatients with ages of greater than 66 years old, the univariate regression models also indicated positive correlations. The linear equations indicated that the basic number of outpatients was 5 to 10 per day as the temperature approached 0°C, and the correlation coefficients R varied between 0.03 and 0.35. Statistically significant correlations (i.e., P < 0.1) occur when correlation coefficients R are large (i.e., R > 0.14). The larger the values for R, the smaller the P-values and the higher the significance of the statistical correlation and the univariate regression model analysis.

Univariate Regression Model Analysis
Statistically These x coefficients in the linear equations indicated that the older people (greater than 66 years old) had the fewest accident per 1 degree increase in temperature, followed by the young people (0-15 years old). The middle-aged people (16-65 years old) were the group of outpatients that were most prone to accidents, with an increase of 1 accident per degree increase in temperature. Fig 5 presents the regression model analyses of the number of hospital emergency room visits for accidents versus temperature during the study period for the three age groups of patients. Positive correlations between the numbers of emergency room visits and temperature were consistently obtained. These results are consistent with the positive correlations between the number of outpatients for accidents and temperature (Fig 4). We note that the coefficients R for the correlations between the number of emergency room visits and temperature were typically greater than the coefficients for outpatient visits for the 0-15 and 16-65 years-old patient age groups. However, for the patients aged greater than 66 years old, the correlation coefficients R for outpatients were greater than the coefficients for emergency room visits. Because hospital emergency rooms are open continuously, the consistently greater correlation coefficients R between the number of hospital visits for accidents and temperature compared with outpatient visits indicate significant correlations for the 0-15 and 16-65 years-old patient age groups.
In order to access time-lag effects on the correlation results shown here, Table 4 compares the correlation coefficients R for the daily hospital emergency room visits that are 0-day, 1-day, and 2-day lagging the daily temperature. For all the 45 scenarios analyzed (3 aged groups of outpatients X 5 years of data X 3 scenarios (0-day, 1-day, and 2-day delays)), 91% of scenarios show smaller correlation coefficients R for the 1-day and 2-day delay when compared with the no delay scenarios. Only 9% of the scenarios show equal correlations coefficients R when either the 1-day or 2-day delay are compared with the 0-day delay (no delay) scenarios. No scenario is found to have higher R from the 1-day and the 2-day delayed visits compared with the no delay visit.
These results indicate the immediate effect of temperatures on the hospital emergency room visits. Interestingly, similar R occurs between the 1-day delay and 0-day delays for the 16-65 years old outpatients at some years (2008, 2009), but not for the young (0-15 years old) and the old (>65 years old) outpatients. This could reflect that, when the needs to immediately visit hospital emergency rooms occur, a more delayed character for the 16-65 years old outpatients compared with other two groups of outpatients. This figure clearly shows that temperature was the most distinctive and persistent factor that determined the number of outpatient visits for accidents compared with other air pollution and meteorological factors. For the 16-65 years-old outpatients, the correlation coefficients R were consistently positive, with R values greater than 0.22 throughout the entire period of 2007-2011. The P-values for this group of outpatients during this 5-year period were very small, less than 0.01 (Fig 6b), thereby indicating statistically significant correlations.
For the other groups of outpatients, the coefficients for the correlations with temperature were also consistently positive. The larger the correlation coefficients R, the smaller the P-value (P less than 0.10), thus indicating the significance of the correlation.
For 12 of the parameters considered, Fig 6 reveals that only the temperature exhibited a persistent and significant correlation with the number of daily hospital outpatient visits for accidents throughout the 5-year analysis period. This is a key result of this work.

Multivariable Regression Model Analysis: Significance of Temperature Factor
To further examine whether there are other potential confounding factors that could affect the correlation between temperatures and daily outpatient visits for accidents, we conducted a multivariable analysis for the 7 air pollution-related variables and 5 meteorological variables considered in this work. Table 5 presents a list of coefficients β, F values, R 2 , and P values for the 12 multivariable regression models listed in Table 1. For multivariable regression models 1 through 7, the statistical significances of the model analyses were low, as indicated by the large P values. In model 8, when the temperature effect was considered in the model, the P-value decreased to 0.001, thus indicating statistical significance of the regression model, and R 2 was equal to 0.11. Further inclusion of other meteorological factors, as demonstrated by models 9 through 12, slightly increased the R 2 values to 0.15 while maintaining statistical significance, i.e., P values of less than 0.001. Correlations between Daily Accident Rates and Daily Ambient Air Temperatures Based on the changes of R 2 and P values, the multivariable regression models revealed that the temperature variable was the dominant factor that associated air pollution-related and meteorological variables with the daily number of outpatient visits for accidents.
The multivariable regression models for 2008 (Table 6) indicated a significant change in R 2 values from 0.04 to 0.16, and the P value decreased from 0.161 to less than 0.001 when the temperature factor was included in the model; for 2009 (Table 7), the R 2 value changed from 0.08 to 0.10, and the P value decreased from 0.005 to 0.001; for 2010 (Table 8), the R 2 value increased from 0.02 to 0.09, and the P value decreased from 0.600 to 0.003 when the temperature was included in the model.
For 2011 (Table 9), the R 2 values increased from 0.13 to 0.19, whereas the P value remain less than 0.001 in model 8, i.e., when the temperature variable was included in the model. Possible confounding factors occurred when CO was included, i.e., in model 4 (R 2 increased from 0.01 to 0.04, and P decreased from 0.382 to 0.067), and when NO 2 was included, i.e., in model 6 (R 2 increased from 0.05 to 0.13, and P decreased from 0.020 to less than 0.001) in the year 2011 (Table 9).
For these 5 years of analysis with the multivariable regression models, the R 2 values varied between 0.10 and 0.19, with P values less than 0.017. These results indicated that the multivariable regression model analyses were statistically significant. The coefficients β for the Model β 1 (PM 10 ) β 7 (SO 2 ) β 2 (PM 2.5 ) β 8 (TEMP) β 3 (O 3 ) β 9 (RAIN) β 4 (CO) β 10 (WD) β 5 (NO) β 11 (WS) β 6 (NO 2 ) β 12 (RH) F R 2 P temperature variable vary between 0.7 and 1.4, consistent with the results from the univariate regression models. As similar to Table 4 for the univariate model analysis, we run the multivariable models through the 180 scenarios (1 aged group of outpatient X 5 years of data X 12 multivariable models X 3 scenarios (0-day, 1-day, and 2-day delays)) to investigate the sensitivity of the time-lag effects on the results shown above. Table 10 shows a list of R 2 from 180 scenarios. Examinations of these results show no increases in R 2 when the 1-day and 2-days delays in the hospital emergency room visits are taken into account. Also, the temperature variable (included in models 8 and after) remains to be the dominating factor in associating with the hospital emergency room visits. Table 9 shows the occurrence of confounding factors CO and NO 2 in 2011. In order to investigate the possible causes of these two confounding factors, we run multivariable models for temperature with respect to variables from air pollutants for the period 2007-2011. Table 11 shows results of β for each air pollutant when the following multivariable model is used: T = β 1 (PM 10 )+β 2 (PM 2.5 )+β 3 (O 3 )+β 4 (CO)+β 5 (NO)+β 6 (NO 2 )+β 7 (SO 2 ). Table 12 shows evolution of R 2 when each of the air pollutant variables is successively included in the multivariable models. The β coefficients (Table 11) show that temperature is positively associated with Model β 1 (PM 10 ) β 7 (SO 2 ) β 2 (PM 2.5 ) β 8 (TEMP) β 3 (O 3 ) β 9 (RAIN) β 4 (CO) β 10 (WD) β 5 (NO) β 11 (WS) β 6 (NO 2 ) β 12 (RH) F R 2 P O 3 and CO, negatively associated with NO and NO 2 , and the associations between temperature and other pollutants (PM 10 , PM 2.5 , and SO 2 ) are weak. As demonstrated in Table 12, when the multivariable model for temperature contains only PM 10 and PM 2.5 two variables, the R 2 values are low (less than 0.04) and statistically not significant. As O 3 been included in the model, R 2 values rise to 0.10 and 0.11, indicating that O 3 levels are positively associated with temperature. This is consistent with higher temperatures often occur during the day-time and cloud-free conditions, resulting in atmospheric conditions that are conducive for photochemical O 3 production.
The inclusion of CO variable further enhances R 2 values to 0.19 to 0.31, indicating that high temperatures are associated with high CO levels. Since vehicular emissions are major sources of CO in urban areas in Taiwan [22], this positive association indicate enhanced vehicular emissions and vehicular usages as temperatures become higher.
The next significant increases in R 2 values occur when NO 2 variable is included in the model. As NO 2 is also a main pollutant from vehicular emissions, high R 2 values demonstrate again the importance of vehicular sources in the study area. The negative coefficient for NO 2 variable, while positive coefficients for O 3 and CO clearly indicate effective photochemical production of O 3 from NO 2 while emissions of CO are sustained in the urban environment. Model β 1 (PM 10 ) β 7 (SO 2 ) β 2 (PM 2.5 ) β 8 (TEMP) β 3 (O 3 )β 9 (RAIN) β 4 (CO) β 10 (WD) β 5 (NO) β 11 (WS) β 6 (NO 2 ) β 12 (RH) F R 2 P Detailed causes of each accident are not available. However, the multivariable models predict that increase in vehicular usage and emissions are strongly associated with increase in temperatures, which in turn are strongly associated with increase in the number of hospital outpatient visits for accidents. Hence, it is likely that the causes of accidents are related to the increases in vehicular usages and emissions during the time periods with high temperatures. More vehicles and higher driving speeds can result in increase in CO and NO 2 emissions, increase in O 3 levels [23], more accidents, and more hospital usage. This analysis is consistent with the hospital clinical experience, which showed that most (about 60%-80%) of the patients for accidents were related to traffic accidents, followed by accidents occurred in work places. Table 12 shows exceptional high R 2 values of 0.50 in for model when NO 2 is included as a model variable compared with other years. As a result, NO 2 variable also exhibits distinctive association with the hospital outpatients for accidents in 2011 than other years.

Conclusions
An accident is an unwanted hazard for any person. However, accidents do occur. It is important to analyze accident data accumulated from people seeking medical care to understand the history of accidents that have occurred such that we can develop methods to reduce the risk of accidents in the future. In this work, we searched for correlations between daily accident rates and Model β 1 (PM 10 ) β 7 (SO 2 ) β 2 (PM 2.5 ) β 8 (TEMP) β 3 (O 3 ) β 9 (RAIN) β 4 (CO) β 10 (WD) β 5 (NO) β 11 (WS) β 6 (NO 2 ) β 12 (RH) F R 2 P Our analysis demonstrated that the number of male outpatients admitted for accidents was approximately 1.31 to 1.47 times the number of females; the difference was statistically significant because the population is approximately equally distributed between males and females (P < 0.0001). For all 12 (air pollution and meteorology) of the parameters considered, only the daily temperature exhibited consistent and significant correlations with the daily number of hospital outpatient visits for accidents throughout the 5-year analysis period.
The univariate regression models indicated that older people (greater than 66 years old) had fewer accidents per 1-degree increase in temperature, followed by young people (0-15 years old). Middle-aged people (16-65 years old) were the group of outpatients that were most prone to accidents, with an increase of 0.8-1.2 accidents per degree increase in temperature.
The multivariable regression models also revealed that the temperature variable was the dominant factor in associating air pollution and meteorological variables with the daily number of outpatient visits for accidents. The results of the multivariable regression model analyses were statistically significant. The accident rates of 0.7-1.4 additional accidents per 1-degree increase in temperature were consistent with the results of the univariate regression models. Model β 1 (PM 10 ) β 7 (SO 2 ) β 2 (PM 2.5 ) β 8 (TEMP) β 3 (O 3 ) β 9 (RAIN) β 4 (CO) β 10 (WD) β 5 (NO) β 11 (WS) β 6 (NO 2 ) β 12 (RH) F R 2 P Our further multivariable model analysis of temperature with respect to air pollution variables show that, through the increases in emissions and concentrations of CO, photochemical O3 production and NO2 loss in the ambient air, increases in vehicular emissions are associated with increases in temperatures. As such, increases in hospital visits for accidents are related to vehicular emissions and usage. This finding is consistent with clinical experience which shows about 60% to 80% of accidents are related to traffic, followed by accidents occurred in work place.
The accident data analyzed in this work were selected according to the ICD-9 codes 800-949. Follow-up works to analyze the detailed nature of each accident as documented according to the ICD-9 codes to further explore how the nature of accidents on changes in meteorology and climate factors are ongoing.
Supporting Information S1 File. Daily meteorology, air pollution, and hospital visit data used in this work. The daily meteorological data (temperature, relative humidity, wind direction, wind speed, rainfall), air pollution data (PM 10 , PM 2.5 , O 3 , CO, NO, NO 2 , SO 2 ), and hospital visits for accidents (outpatient visits, emergency room visits) for each year of the 2007-2011 study period are included in the uploaded data file. There is also a README file, containing the detailed file structures for the data set, and the description for each data file. Interested readers can contact Kuo-Ying Wang at kuoying@mail.atm.ncu.edu.tw to access data as described above. Alternatively, interested readers can use following methods to obtain data (meteorology and air pollution) from Taiwan Environmental Protection Administration (EPA): 1) Click on http://taqm.epa.gov.tw/ taqm/en/default.aspx; 2) Look on a list of items from the left column shown from the pop-up window. Search and click on an item called "Data Service", which will lead to "Air Quality". Then click on it to see "Hourly Value". Then click on "Hourly Value" to download hourly data for all ambient air monitoring sites in Taiwan