| Literature DB >> 33536798 |
Ohmi Watanabe1, Norio Narita2, Masahito Katsuki2, Naoya Ishida1, Siqi Cai1, Hiroshi Otomo3, Kenichi Yokota3.
Abstract
PURPOSE: With the aging population in Japan, the prediction of ambulance transports is needed to save the limited medical resources. Some meteorological factors were risks of ambulance transports, but it is difficult to predict in a classically statistical way because Japan has 4 seasons. We tried to make prediction models for ambulance transports using the deep learning (DL) framework, Prediction One (Sony Network Communications Inc., Tokyo, Japan), with the meteorological and calendarial variables.Entities:
Keywords: ambulance transport; cardiopulmonary arrest; deep learning; meteorological factors; trauma
Year: 2021 PMID: 33536798 PMCID: PMC7850460 DOI: 10.2147/OAEM.S293551
Source DB: PubMed Journal: Open Access Emerg Med ISSN: 1179-1500
Figure 1(A) Map of Japan. Miyagi prefecture, where Kesennuma City Hospital (HP) locates, is in the northeastern part of Japan (colored black). (B) Map of Miyagi prefecture. The Kesennuma area, with a background population of about 100,000 people, is colored light green. The red triangle is Kesennuma City HP, where this study was performed.
Number of Ambulance Transports and Meteorological Factors
| Year | Month | Ambulance Transport (Number) | Meteorological Factors | ||||
|---|---|---|---|---|---|---|---|
| Total | CPA | Trauma | Tmean (°C) | Pa (hPa) | RH (%) | ||
| 2017 | Jan | 188 | 10 | 4 | 1.3 | 1009.5 | 64 |
| Feb | 159 | 7 | 3 | 2.3 | 1007.8 | 62 | |
| Mar | 154 | 1 | 11 | 4.0 | 1009.6 | 64 | |
| Apr | 150 | 4 | 7 | 9.9 | 1007.9 | 66 | |
| May | 168 | 10 | 7 | 15.7 | 1008.0 | 70 | |
| Jun | 164 | 5 | 10 | 17.0 | 1004.4 | 79 | |
| Jul | 151 | 4 | 15 | 23.6 | 1004.4 | 82 | |
| Aug | 153 | 2 | 7 | 21.8 | 1004.8 | 89 | |
| Sep | 178 | 5 | 12 | 19.4 | 1007.5 | 79 | |
| Oct | 134 | 4 | 6 | 13.7 | 1014.7 | 80 | |
| Nov | 161 | 9 | 9 | 8.1 | 1012.3 | 68 | |
| Dec | 172 | 5 | 14 | 2.5 | 1010.3 | 66 | |
| Total of 2017 | 1932 | 66 | 105 | 11.6 | 1008.4 | 72 | |
| 2018 | Jan | 189 | 15 | 11 | 0.9 | 1008.2 | 61 |
| Feb | 152 | 11 | 6 | 0.4 | 1011.5 | 60 | |
| Mar | 148 | 3 | 13 | 6.8 | 1012.5 | 60 | |
| Apr | 150 | 4 | 9 | 11.4 | 1009.6 | 66 | |
| May | 165 | 5 | 14 | 15.4 | 1007.0 | 71 | |
| Jun | 129 | 4 | 7 | 18.5 | 1005.3 | 81 | |
| Jul | 159 | 0 | 7 | 23.6 | 1006.7 | 87 | |
| Aug | 196 | 7 | 12 | 23.6 | 1005.1 | 83 | |
| Sep | 146 | 7 | 4 | 19.9 | 1010.1 | 84 | |
| Oct | 159 | 7 | 9 | 15.2 | 1011.7 | 77 | |
| Nov | 147 | 5 | 6 | 9.3 | 1015.7 | 71 | |
| Dec | 204 | 11 | 15 | 3.5 | 1013.2 | 66 | |
| Total of 2018 | 1944 | 79 | 113 | 12.4 | 1009.7 | 72 | |
| 2019 | Jan | 179 | 13 | 6 | 1.3 | 1010.2 | 61 |
| Feb | 130 | 7 | 3 | 2.0 | 1012.9 | 61 | |
| Mar | 151 | 9 | 6 | 5.3 | 1008.3 | 63 | |
| Apr | 148 | 7 | 10 | 8.9 | 1008.7 | 66 | |
| May | 179 | 6 | 10 | 16 | 1004.0 | 67 | |
| Jun | 158 | 5 | 10 | 17.8 | 1006.2 | 83 | |
| Jul | 168 | 8 | 13 | 21.7 | 1004.9 | 83 | |
| Aug | 201 | 5 | 7 | 24.8 | 1011.5 | 83 | |
| Sep | 185 | 2 | 16 | 21.1 | 1013.1 | 78 | |
| Oct | 190 | 11 | 11 | 15.8 | 1013.1 | 80 | |
| Nov | 200 | 9 | 18 | 8.6 | 1013.6 | 66 | |
| Dec | 183 | 13 | 9 | 4.1 | 1014.4 | 68 | |
| Total of 2019 | 2072 | 95 | 119 | 12.2 | 1010.1 | 72 | |
| Total of 3 years | 5948 | 240 | 337 | 12.4 | 1009.4 | 72 | |
Abbreviations: CPA, cardiopulmonary arrest; Pa, mean atmospheric pressure; RH, mean relative humidity; Tmean, mean ambient temperature.
Approximation by Linear, Quadratic, or Cubic Polynomial Curve Between the Number of Ambulance Transports and Meteorological Factors in the Kesennuma Area
| Variables (X) | r2 | p value | Polynomial Curve† |
|---|---|---|---|
| Tmean (°C) | 0.018 | <0.001* | Y = 0.005X2 − 0.131X + 5.872 |
| Tmax (°C) | 0.019 | <0.001* | Y = 0.005X2 − 0.152X + 5.418 |
| Tmin (°C) | 0.014 | 0.001* | Y = 0.004X2 − 0.083X + 5.479 |
| Tmax-min (°C) | 0.007 | 0.045* | Y = 0.003X3 − 0.082X2 + 0.761X + 3.331 |
| Pa (hPa) | 0.001 | 0.376 | Y = − 0.010X + 15.362 |
| Daily mean vapor pressure (hPa) | 0.007 | 0.024* | Y = 0.004 X2 − 0.111X + 5.974 |
| Mean wind speed (m/s) | 0.009 | 0.016* | Y = 0.046 X3 − 0.564X2 + 2.158X + 2.939 |
| Daily amount of the rainfall (mm) | 0.002 | 0.270 | Y = 0.0002X2 − 0.021X + 5.475 |
| Daily amount of the snowfall (mm) | <0.001 | 0.943 | Y = − 0.013X2 + 0.077X + 5.431 |
| Sunlight hours (hr) | 0.005 | 0.057 | Y = − 0.010X2 + 0.140X + 5.141 |
| RH (%) | 0.009 | 0.006* | Y = − 0.001X2 − 0.148X + 0.582 |
| THI (°C) | 0.016 | <0.001* | Y = 0.007X2 − 0.190X + 6.353 |
| Tmean (°C) | 0.021 | <0.001* | Y = 0.0004X2 - 0.019X + 0.306 |
| Tmax (°C) | 0.020 | <0.001* | Y = 0.0004X2 − 0.020X + 0.410 |
| Tmin (°C) | 0.021 | <0.001* | Y = 0.0004X2 − 0.013X + 0.276 |
| Tmax-min (°C) | 0.001 | 0.574 | Y = − 0.001X2 − 0.014X + 0.173 |
| Pa (hPa) | 0.005 | 0.019* | Y = 0.005X − 4.813 |
| Daily mean vapor pressure (hPa) | 0.015 | 0.001* | Y = 0.003X2 − 0.017X + 0.353 |
| Mean wind speed (m/s) | 0.001 | 0.255 | Y = 0.016 X + 0.176 |
| Daily amount of the rainfall (mm) | 0.001 | 0.418 | Y = − 0.001X + 0.223 |
| Daily amount of the snowfall (mm) | 0.001 | 0.982 | Y = 0.001X2 - 0.011X + 0.220 |
| Sunlight hours (hr) | 0.007 | 0.019* | Y = − 0.003X2 + 0.034X + 0.164 |
| RH (%) | 0.004 | 0.039* | Y = − 0.002X + 0.370 |
| THI (°C) | 0.019 | <0.001* | Y = 0.001X2 − 0.025X + 0.408 |
| Tmean (°C) | 0.004 | 0.147 | Y = − 0.0004X2 + 0.013X + 0.244 |
| Tmax (°C) | 0.004 | 0.095 | Y = − 0.0003X2 + 0.015X + 0.188 |
| Tmin (°C) | 0.003 | 0.250 | Y = − 0.0003X2 + 0.008X + 0.292 |
| Tmax-min (°C) | 0.002 | 0.154 | Y = 0.008 X + 0.243 |
| Pa (hPa) | <0.001 | 0.892 | Y = 0.0003X − 0.043 |
| Daily mean vapor pressure (hPa) | 0.003 | 0.164 | Y = − 0.001X2 + 0.018X + 0.204 |
| Mean wind speed (m/s) | 0.002 | 0.410 | Y = 0.003X2 − 0.017X + 0.353 |
| Daily amount of the rainfall (mm) | <0.001 | 0.764 | Y = 0.00003X2 − 0.002X + 0.312 |
| Daily amount of the snowfall (mm) | <0.001 | 0.993 | Y = − 0.001X2 +0.007X + 0.308 |
| Sunlight hours (hr) | 0.008 | 0.031* | Y = 0.001X3 − 0.018X2 + 0.092X + 0.227 |
| RH (%) | 0.003 | 0.191 | Y = − 0.0002X2 + 0.022X − 0.445 |
| THI (°C) | 0.003 | 0.201 | Y = − 0.001X2 + 0.017X + 0.204 |
Notes: *p < 0.05 for the calculated polynomial curve; †linear, quadratic, or cubic polynomial curve between the number of ambulance transports and meteorological factors were calculated by SPSS software.
Abbreviations: CPA, cardiopulmonary arrest; Pa, mean atmospheric pressure; r2, determination coefficient; RH, mean relative humidity; THI, thermo-hydrological index; Tmax, daily highest ambient temperature; Tmax-min, daily difference between the highest and lowest ambient temperature Tmean; mean ambient temperature; Tmin, daily lowest ambient temperature.
Day of the Week and Numbers of Ambulance Transports
| Day of the Week† | Number of the Days | Total Ambulance Transport | CPA | Trauma |
|---|---|---|---|---|
| Monday | 157 | 902 | 37 | 51 |
| Tuesday | 157 | 814 | 27 | 48 |
| Wednesday | 156 | 828 | 34 | 45 |
| Thursday | 156 | 832 | 35 | 42 |
| Friday | 156 | 865 | 33 | 41 |
| Saturday | 156 | 838 | 38 | 55 |
| Sunday | 156 | 869 | 36 | 55 |
| National holiday | 55 (4.6%) | 319 | 11 | 24 |
Notes: †Mann–Whitney U-test and Kruskal–Wallis test did not show a significant difference depending on the day of the week nor national holidays (all p > 0.05).
Abbreviation: CPA, cardiopulmonary arrest.
Figure 2The area under the curves (AUCs) and calculation tables of each model. (A) The model for total daily ambulance transport more than 5 (median) or not has an AUC of 0.972 (95%confident interval (95% CI) 0.960–0.984). Its sensitivity and specificity were 0.937 and 0.935. (B) The model for daily CPA transport (present or absent) has an AUC of 0.988 (95% CI 0.975–1.000). Its sensitivity and specificity were 0.958 and 0.972. (C) The model for daily trauma transport (present or absent) has an AUC of 0.947 (95% CI 0.927–0.968). Its sensitivity and specificity were 0.899 and 0.892.