隨著全球氣候暖化趨於嚴重，登革熱逐漸成為全球關切的公共
生疾病，近年臺灣登革熱疫情頻繁發生，過去相關研究著重在南部登革熱疫情的空間分析，通常以發生率來描述疫情之嚴重性。針對2015年登革熱疫情擴散全臺，本研究欲以全台為範圍，以流行風險指標進行空間分析。
使用資料為2015年台灣鄉鎮地區的登革熱資料，由於登革熱疫情在空間上布不平衡，約3成地區沒有確定病例，資料分布上也呈現不均且右偏，不適合以傳統迴歸(平均數迴歸)進行分析。因此，以分量迴歸與地理加權分量迴歸之技術進行研究，欲探究空間異質的變化情形，亦可觀察分量在空間上的變化情形，並探討與文獻中疫情嚴重等級的劃分區塊之關係。

As the global warming getting serious, the Dengue fever gradually become a global public health concern. In recent years, the epidemic of the Dengue fever frequently occurs in Taiwan. However, researchers in the past focus in the southern of Taiwan, using spatial analysis to determine the epidemic. In this paper, we will focus the epidemic of Dengue fever all over Taiwan in 2015, using epidemic risk index to process spatial analysis.
The data used for dengue fever data on Taiwan’s townships in 2015.The dengue fever epidemic is the spatial imbalance. The data distribution is also uneven and right-skewed,  there are about 30% no confirmed cases in the regions, which is not suitable for traditional regression (mean regression)
analyze. In order to explore the variation in spatial heterogeneity, and also observe the spatial variation of quantile, using quantile regression and geographically weighted quantile regression is suitable. We also discussed the relationship between the results of this thesis and the severity level of epidemic on the references.

目錄
第一章 緒論. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 
1.1 研究背景與動機 . . . . . . . . . . . . . . . . . . . . . . . . . . 1 
1.2 研究目的 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 
1.3 研究架構 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 
第二章 文獻探討. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 
2.1 登革熱概說 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 
2.1.1 登革熱病毒與病媒 . . . . . . . . . . . . . . . . . . . . 6 
2.1.2 登革熱發病症狀 . . . . . . . . . . . . . . . . . . . . . . 6 
2.1.3 登革熱流行 . . . . . . . . . . . . . . . . . . . . . . . . 7 
2.2 危險因子 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 
2.2.1 環境因子 . . . . . . . . . . . . . . . . . . . . . . . . . 8 
2.2.2 社會因子 . . . . . . . . . . . . . . . . . . . . . . . . . 8 
2.3 登革熱分布狀況 . . . . . . . . . . . . . . . . . . . . . . . . . . 9 
2.3.1 疾病地圖 . . . . . . . . . . . . . . . . . . . . . . . . . 9 
2.3.2 空間相關 . . . . . . . . . . . . . . . . . . . . . . . . . 9 
2.3.3 空間迴歸 . . . . . . . . . . . . . . . . . . . . . . . . . 12 
第三章 資料介紹. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 
3.1 資料來源 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 
3.2 反應變數 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 
3.3 解釋變數 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 
3.3.1 環境因子 . . . . . . . . . . . . . . . . . . . . . . . . . 17 
3.3.2 社會因子 . . . . . . . . . . . . . . . . . . . . . . . . . 20 
3.4 篩選結果 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 
第四章 研究方法. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 
4.1 地理加權迴歸 . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4.2 分量迴歸 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 
4.2.1 估計 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 
4.3 地理加權分量迴歸 . . . . . . . . . . . . . . . . . . . . . . . . 30 
4.3.1 估計 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 
4.4 區域型Moran’s I(LISA) . . . . . . . . . . . . . . . . . . . . . . 33 
第五章 研究結果. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 
5.1 敘述性探索分析 . . . . . . . . . . . . . . . . . . . . . . . . . . 34 
5.1.1 反應變數 . . . . . . . . . . . . . . . . . . . . . . . . . 34 
5.1.2 解釋變數 . . . . . . . . . . . . . . . . . . . . . . . . . 41 
5.2 分量迴歸 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 
5.2.1 變數確認 . . . . . . . . . . . . . . . . . . . . . . . . . 47 
5.2.2 檢測分量 . . . . . . . . . . . . . . . . . . . . . . . . . 48 
5.2.3 結果分析 . . . . . . . . . . . . . . . . . . . . . . . . . 49
5.2.4 Moran’s Index . . . . . . . . . . . . . . . . . . . . . . . 50 
5.3 空間分析 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 
5.3.1 空間異質性 . . . . . . . . . . . . . . . . . . . . . . . . 51 
5.3.2 空間自相關與局部共線性的檢驗 . . . . . . . . . . . . 55 
5.3.3 變數顯著性 . . . . . . . . . . . . . . . . . . . . . . . . 57 
5.3.4 模型預測 . . . . . . . . . . . . . . . . . . . . . . . . . 82 
5.3.5 對應指標 . . . . . . . . . . . . . . . . . . . . . . . . . 84 
第六章 結論. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 
6.1 總結與討論 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 
6.2 未來研究 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
圖目錄
圖 2.1 LISA聚類地圖範例 . . . . . . . . . . . . . . . . . . . . . . 11
圖 3.1 相等發生率之不同時間特性 . . . . . . . . . . . . . . . . . 15
圖 5.1 反應變數分布 . . . . . . . . . . . . . . . . . . . . . . . . . 37 
圖 5.2 反應變數地圖分布 . . . . . . . . . . . . . . . . . . . . . . . 38 
圖 5.3 LISA聚類圖 . . . . . . . . . . . . . . . . . . . . . . . . . . 39 
圖 5.4 LISA兩等級分級 . . . . . . . . . . . . . . . . . . . . . . . . 39 
圖 5.5 疫情分類圖 . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 
圖 5.6 解釋變數地圖分布 . . . . . . . . . . . . . . . . . . . . . . . 42 
圖 5.7 解釋變數地圖分布 . . . . . . . . . . . . . . . . . . . . . . . 43 
圖 5.8 解釋變數地圖分布 . . . . . . . . . . . . . . . . . . . . . . . 44 
圖 5.9 分量0.5敏感性族群顯著分布圖 . . . . . . . . . . . . . . . . 59 
圖 5.10 分量0.9敏感性族群顯著分布圖 . . . . . . . . . . . . . . . 60 
圖 5.11 分量0.5弱勢族群顯著分布圖 . . . . . . . . . . . . . . . . 61 
圖 5.12 分量0.9弱勢族群顯著分布圖 . . . . . . . . . . . . . . . . 62 
圖 5.13 分量0.5都市化顯著分布圖 . . . . . . . . . . . . . . . . . . 63 
圖 5.14 分量0.9都市化顯著分布圖 . . . . . . . . . . . . . . . . . . 64 
圖 5.15 分量0.5綠覆率顯著分布圖 . . . . . . . . . . . . . . . . . . 66 
圖 5.16 分量0.9綠覆率顯著分布圖 . . . . . . . . . . . . . . . . . . 67 
圖 5.17 分量0.5布氏級數顯著分布圖 . . . . . . . . . . . . . . . . 68 
圖 5.18 分量0.9布氏級數顯著分布圖 . . . . . . . . . . . . . . . . 69 
圖 5.19 分量0.5週平均溫度大於18度之週數顯著分布圖 . . . . . . 70 
圖 5.20 分量0.9週平均溫度大於18度之週數顯著分布圖 . . . . . . 71 
圖 5.21 分量0.5 1-5月最長濕潤日數顯著分布圖 . . . . . . . . . . 73 
圖 5.22 分量0.9 1-5月最長濕潤日數顯著分布圖 . . . . . . . . . . 74
圖 5.23 分量0.5 1-5月最長乾燥日數顯著分布圖 . . . . . . . . . . 75 
圖 5.24 分量0.9 1-5月最長乾燥日數顯著分布圖 . . . . . . . . . . 76 
圖 5.25 分量0.5年平均降雨量顯著分布圖 . . . . . . . . . . . . . . 77 
圖 5.26 分量0.9年平均降雨量顯著分布圖 . . . . . . . . . . . . . . 78 
圖 5.27 分量0.75綠覆率顯著分布圖 . . . . . . . . . . . . . . . . . 80 
圖 5.28 分量0.95綠覆率顯著分布圖 . . . . . . . . . . . . . . . . . 81
表目錄
表 3.1 資料一覽 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 
表 3.2 布氏指數與布氏級數轉換 . . . . . . . . . . . . . . . . . . . 17 
表 3.3 環境因子定義與計算 . . . . . . . . . . . . . . . . . . . . . 18 
表 3.4 社會因子定義與計算 . . . . . . . . . . . . . . . . . . . . . 21 
表 3.5 特徵比例 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 
表 3.6 轉軸後的成分矩陣 . . . . . . . . . . . . . . . . . . . . . . . 23 
表 3.7 變異數膨脹因子 . . . . . . . . . . . . . . . . . . . . . . . . 25
表 4.1 分類基準與描述 . . . . . . . . . . . . . . . . . . . . . . . . 33
表 5.1 各變數基礎統計量 . . . . . . . . . . . . . . . . . . . . . . . 34 
表 5.2 各分位數下反應變數值 . . . . . . . . . . . . . . . . . . . . 36 
表 5.3 反應變數之分量迴歸係數 . . . . . . . . . . . . . . . . . . . 49 
表 5.4 全域Moran’s Index . . . . . . . . . . . . . . . . . . . . . . 50 
表 5.5 流行頻率在不同分量下各變數的空間異質性 . . . . . . . . 52 
表 5.6 流行延時在不同分量下各變數的空間異質性 . . . . . . . . 53 
表 5.7 傳播強度在不同分量下各變數的空間異質性 . . . . . . . . 54 
表 5.8 全域與局部Moran’s Index . . . . . . . . . . . . . . . . . . . 56 
表 5.9 局部變異數膨脹 . . . . . . . . . . . . . . . . . . . . . . . . 57 
表 5.10 殘差比較 . . . . . . . . . . . . . . . . . . . . . . . . . . . 83 
表 5.11 疫情分類之基礎統計量 . . . . . . . . . . . . . . . . . . . 84 
表 5.12 疫情分類之基礎統計量 . . . . . . . . . . . . . . . . . . . 85
表 1 環境因子溫度面向相關係數 . . . . . . . . . . . . . . . . . . 89 
表 2 環境因子雨量面向相關係數 . . . . . . . . . . . . . . . . . . 90 
表 3 社會因子之相關係數矩陣 . . . . . . . . . . . . . . . . . . . . 91
表 4 流行頻率之變數確認 . . . . . . . . . . . . . . . . . . . . . . 92 
表 5 流行延時之變數確認 . . . . . . . . . . . . . . . . . . . . . . 93 
表 6 傳染強度之變數確認 . . . . . . . . . . . . . . . . . . . . . . 94 
表 7 流行頻率在不同分量下的迴歸係數之檢定 . . . . . . . . . . . 95 
表 8 流行延時在不同分量下的迴歸係數之檢定 . . . . . . . . . . . 96 
表 9 傳播強度在不同分量下的迴歸係數之檢定 . . . . . . . . . . . 97

參考文獻
[1] Almeida M. C. D., Caiaffa W. T., Assuncao R. M. and Proietti F. A.,(2007) Spatial vulnerability to dengue in a Brazilian urban area during a 7-year surveillance, Journal of Urban Health-Bulletin of the New York Academy of Medicine: 84(3) 334-345.
[2] Brunsdon, C., Fotheringham, S., Charlton, M., (1998) Geographically weighted regression. Journal of the Royal Statistical Society: Series D (The Statistician), 47(3), 431-443.
[3] Chen, V. Y. J., Wu, P. C., Yang, T. C., Su, H. J., Examining nonstationary effects of social determinants on cardiovascular mortality after cold surges in Taiwan. Science of the Total Environment, 408(9), 2042 2049.
[4] Chen, V.Y.J., Deng, W.S., Yang, T.C. and Matthews, S.A. (2012). Geographically Weighted Quantile Regression (GWQR): An Application to U.S. Mortality Data. Geographically Analysis, 44(2): 134-150.
[5] David Wheeler, Michael Tiefelsdorf (2005)Multicollinearity and Correlation Among Local Regression Coefficients in Geographically Weighted Regression.J Geograph Syst 7: 161–187
[6] David Wheeler, Diagnostic tools and a remedial method for collinearity in geographically weighted regression. Environment and Planning A 2007, volume 39, pages 2464-2481
[7] DJ Gubler et al. (1945) Dengue and Dengue Hemorragic Fever
[8] Fotheringham, A. S., Brunsdon, C. and Charlton, M. (2002). Geographically weighted regression: the analysis of spatially varying relationships.
[9] H. Wu, L. Gao, Z. Zhang.,(2014) Analysis of crash data using quantile regression for counts. J. Transp. Eng., Vol.140, No.4. Apr.2014.
[10] Jansen CC, Beebe NW., (2010) The dengue vector Aedes aegypti: what comes next. Microbes and Infection, 12: 272-279.
[11] Jerry M. Spiegel, Mariano Bonet, AnaMaria Ibarra, Nino Pagliccia, Veronic Ouellette, Annalee Yassi. (2007) Social and environmental determinants of Aedes aegypti infestation in Central Havana: results of a case–control study nested in an integrated dengue surveillance programme in Cuba. Tropical Medicine and International Health Volume 12, Issue 4.
[12] Jos´e A. F Machado,J. M. C. Santos Silva (2005) Quantiles for Counts, Journal of the American Statistical Association, 100:472, 1226-1237.
[13] Koenker R., Bassett(1978).Regression Quantile. Econometrica, 46,33-50.
[14] Koenker R., Pin N.G., Stephen Portnoy. Quantile smoothing splines. Biometrika, Volume 81, Issue 4, 1 December 1994, Pages 673–680.
[15] Koenker R. (2005) Quantile Regression. Cambridge University Press, Cambridge.
[16] Lifson AR (1996) Mosquitoes, models, and dengue. Lancet. 1996 May 4; 347(9010):1201-2.
[17] Lin C.H., Wen T.H. (2011, Jul). Using Geographically Weighted Regression (GWR) to Explore Spatial Varying Relationships of Immature Mosquito and Human Densities with the Incidence of Dengue. International Journal of Environmental Research and Public Health , 8(7): 2798
2815.
[18] Luc Anselin (1995) Local Indicators of Spatial Association—LISA. Geographical Analysis Volume 27, Issue 2.
[19] Ma, L., and Pohlman, L. (2008). Return forecasts and optimal portfolioconstruction: A quantile regression approach.
[20] Normile (2013)D. Surprising new dengue virus throws a spanner in disease control efforts. Science. 342 (6157): 415.
[21] P. A. P. Moran (1950) Notes on Continuous Stochastic Phenomena. Biometrika Vol. 37, No. 1/2, pp. 17-23
[22] Pei-Chih Wu, Jinn-Guey Lay, How-Ran Guoc, Chuan-Yao Lin, Shih Chun Lung, Huey-Jen Suc. (2009) Higher temperature and urbanization affect the spatial patterns of dengue fever transmission in subtropical Taiwan. Science of The Total Environment Volume 407, Issue 7,Pages 2224-2233.
[23] Rosa-Freitas, M. G., Tsouris, P., Sibajev, A., de Souza Weimann, E.T., Marques, A. U., Ferreira, R. L., Luitgards-Moura, J. F., (2003) Exploratory temporal and spatial distribution analysis of dengue notifications in Boa Vista, Roraima, Brazilian Amazon, 1999-2001., Dengue Bulletin, 27, 63-80.
[24] Patz, J. A. , Public health risk assessment linked to climatic and ecological change, Human and Ecological Risk Assessment: An International Journal (2001), 7(5), 1317-1327.
[25] S Promprou, M Jaroensutasinee, K Jaroensutasinee (2005) Climatic Factors Affecting Dengue Haemorrhagic Fever Incidence in Southern Thailand.
[26] Variyam, J., Blaylock, J., Smallwood, D., (2002) Characterizing the Distribution of Macronutrient Intake among U.S. Adults: A Quantile Regression Approach, American Journal of Agricultural Economics. 84:454-466.
[27] Wen, T.H., Lin, N.H., Lin, C.H., King, C.C. Su, M.D. (2006). Spatial mapping of temporal risk characteristics to improve environmental health risk identification: A case study of a dengue epidemic in Taiwan. Science of the Total Environment, 367 (2-3), 631-640.
[28] Wu, H. (2011). A framework for developing road risk indices using quantile regression based crash prediction model. Ph.D. dissertation, Dept. of Civil, Architectural, and Environmental Engineering, Univ. of Texas at Austin.
[29] 王一安 (2016) 地理加權卜瓦松迴歸於台灣各鄉鎮區登革熱資料之分 析, 淡江大學統計學系碩士班學位論文2016: 1-62.
[30] 王麗容 , (1994) 社會變遷中的親職教育需求、觀念與策略 , 國立台灣 大學社會學刊 (1994) 23 : 191-216.
[31] 吳樹寰 (2013)臺灣地區性氣候因子及聖嬰/反聖嬰事件對登革熱發生 之時間影響趨勢, 中國醫藥大學公共衛生學系碩士班學位論文
[32] 林政宏, (2007) 台灣地區登革熱擴散之空間分析 , 臺灣大學地理環境資 源學研究所學位論文
[33] 高醫 (1988)登革熱的症狀與診斷.醫訊月刊, 第二十二卷(第十一期)。
[34] 亞拉妮 (2013) 降雨對於登革熱發病率的潛在衝擊：印尼與台灣的案例研究
[35] 徐筱瑜、賴淑寬、郭俊賢、吳智文、劉定萍, (2013) 2003-2011 年臺灣 登革出血熱個案流行病學分析. 疫情報導 (2013),29(21) : 319-328.
[36] 莊家彰,管中閔 (2005), 台灣與美國股市價量關係的分量迴歸分析, 經濟 論文, 33:4, 379-404
[37] 陳薇琳 (2012) 臺灣登革熱疫情之時空群集研究。國立陽明大學生物醫 學資訊研究所碩士班學位論文
[38] 陳郁雰 (2014) 台灣鄉鎮市區自殺死亡率之地理加權迴歸分析 , 淡江大 學統計學系碩士班學位論文
[39] 陳慈忻, (2017) 登革熱疫情特徵的時空脆弱度因素:多層次模型的分析 , 國立台灣大學地理環境資源學系碩士班學位論文
[40] 陳錦生、黃正中, 登革熱病媒蚊之生態學研究 I.發育零點與生命表之 比較,玉山生物學報 (1988) 7: 95-114.
[41] 黃基森 (1991) 臺灣地區斑蚊生態及其與登革熱流行之關係，中華昆蟲特刊
[42] 黃基森、徐爾烈 , 高雄市登革熱病媒蚊分布及孳生環境調查 , 中華昆 蟲 (1994),14: 233-244.
[43] 許惠美, (1993) 台灣南部的登革熱 , 台南市衛生局疾病管制課
[44] 連日清、王正雄、周欽賢, (1988) 醫學昆蟲學 , 南山堂出版社.
[45] 溫在弘、金傳春(協同)、謝心怡(2006)建立時間-空間擴散動態模型評 估登革熱流行的空間危險區域與防疫策略
[46] 溫在弘、陳慈忻 (2016)，氣候變遷下的登革熱擴散風險與調適策略， 《氣候變遷下的國家發展藍圖》，台北：國立臺灣大學。
[47] 衛生福利部疾病管制署 , 登革熱病媒蚊生態及習性介紹 (2012) : 1-25.
[48] 衛生福利部疾病管制署 , 傳染病防治工作手冊 (2016).
[49] 賴進貴、葉高華、陳汶軍, (2005) 地理資訊系統應用於登革熱疫情防 治之檢討與建議 , 環境與世界 11 : 65-81.
[50] 龍紀萱, (2011) 原住民長期照護服務模式之探討 ,社區發展季刊 136 :1-14.
[51] 顏吟真 (2012) 地理加權自相關分量迴歸 , 淡江大學統計學系碩士班學 位論文: 1-62.
[52] 內政資料開放平台：https://data.moi.gov.tw/MoiOD/default/Index.aspx
[53] 政府開放平台：https://data.gov.tw/
[54] 財政部財政資訊中心:https://www.fia.gov.tw
[55] 衛生福利部：https://www.cdc.gov.tw/rwd
[56] 社會經濟資料庫：https://segis.moi.gov.tw/STAT/Web/Portal/STAT
PortalHome.aspx
[57] Luc Anselin (2004) Geoda操作手冊