近年來隨著科技快速發展、醫療技術進步，人們面臨非預期的生存風險，長壽風險， 近來台灣的壽險轉型推動「保障型保單」，包含了高齡化商品、年金險等，死亡率的精確估計對於保險公司來說越發重要。
本研究使用資料為年齡介於60至99歲的女性，時間為1940年至2017年的多國資料，共澳大利亞、加拿大、丹麥、芬蘭、法國、冰島、義大利、荷蘭、挪威、匈牙利、西班牙、瑞士、瑞典、美國、英國，總計十五國的死亡率資料，以估計美國女性高齡死亡率為目標，建立兩階段因子關聯結構模型。
本研究第一階段，按年齡從60歲至99歲共分40個子集，每個子集建立縮減SUR模型，以AIC選取SUR模型最佳落後項，並加入機器學習(Boosting)估計值作為外生變數，最後以個別t檢定篩選掉不顯著的變數，以建立各特定年齡的邊際分配。
接續將縮減SUR模型的估計殘差標準化，得到標準化估計殘差，以此建立各年齡間，關聯結構下的殘差分配，參考Chen, MacMinn, and Sun (2017)，先以單因子模型縮減參數，再根據模型建立潛在變數，模擬樣本，最後建立出各年齡殘差間的聯合機率分配，在此也參考Oh, Patton (2013)所使用的GAS結構，讓因子關聯結構能有效地動態調整。
經過上述步驟模型建立完成後，本研究於樣本內、樣本外跟原模型進行比較，若以AIC衡量模型表現，發現本研究模型表現最佳，並且能有效改善原模型殘差間存在正相關性的問題。
本研究推論，因為在第一階段本研究使用了多國資料，亦採用演算法基礎、非線性的機器學習方式協助估計，理應比原模型捕捉更多的資訊，讓各年齡間的殘差相關性有效降低。但仍然可以發現，年齡在85歲以上，甚至是95至99歲的超高年齡族群間，仍存在明顯相關，如何處理超高齡死亡率間的相依性，是未來研究的一大方向。

In this study, we used the multiple countries female mortality data from 1940 to 2017, and aged 60 to 99. With the goal of estimating the mortality rate of the US females at an advanced age, we established a two-stage factor copula model in this study.
In the first stage of this study, divided data into 40 subsets according to age from 60 to 99 years old, and a reduced SUR model was established for each subset. And the Boosting estimated value was added as the exogenous variables. Also, individual t-tests are used to filter out insignificant variables to establish the marginal distribution of each specific age.
Then standardize the estimated residuals to establish the distribution of residuals under the correlation structure between ages. Refer to Chen, MacMinn, and Sun (2017), We estimated a single factor model and establish the latent variables according to the model. Finally, we established the joint probability distribution among the residuals. We used the GAS structure referred by Oh, Patton (2013), so that the factor structure can be effectively adjusted dynamically.
Through the above steps, we compared our model with the original one in and out of the sample. We measured performances of the model by AIC, it’s found that the model in this study performs the best and can effectively improve the positive correlation between the residuals of the original.
Because, in the first stage, we used data from multiple countries and also used algorithm-based and non-linear machine learning methods to assist in the estimation, it should capture more information than the original model. However, it can still be found that there is still a significant correlation between the ages of 85 and over, even between the ages of 95 and 99. How to deal with the interdependence of the mortality rate of the very elderly is a major direction of future research.

目錄
摘要 i
Abstract ii
誌謝 iii
目錄 iv
圖目錄 v
表目錄 vi
附錄 vi
第一章、緒論 1
1.1 研究背景與動機 1
1.2 研究設計 2
第二章、文獻回顧 3
2.1 死亡率模型 3
2.2 Lee-Carter 模型 3
2.3 因子關聯結構 (Factor Copula) 4
2.4 Generalized Autoregressive Score Model (GAS Model) 7
第三章、研究方法 8
3.1 研究設計 8
3.2 研究方法 9
3.3 模型設計 14
第四章、實證研究 16
4.1 死亡率資料介紹 16
4.2 資料檢定 18
4.3 機器學習(Boosting)估計 20
4.4 縮減SUR模型(Reduced SUR Model) 21
4.5 因子關聯結構(Factor Copula) 25
4.6 模型表現比較 27
第五章、總結 31
5.1 結論與探討 31
5.2 研究限制與後續工作 31
參考文獻 32
附錄 34

1. H. Chen, R. D. MacMinn, and T. Sun, “Mortality Dependence and Longevity Bond Pricing:A Dynamic Factor Copula Mortality Model with the GAS Structure,” The Journal of Riskand Insurance, vol. 84, no. S1R, pp. 393–415, 2017.
2. D. H. Oh and A. J. Patton, “Time-Varying Systemic Risk:Evidence From a Dynamic CopulaModel of CDS Spreads,” Journal of Business and Economic Statisticss, vol. 36, no. 2,pp. 181–195, 2018.
3. D. Creal, S. J. Koopman, and A. Lucas, “Generalized Autoregressive Score Models WithApplications,” Journal of Applied Econometrics, vol. 28, pp. 777–795, 2013.
4. R. D. Lee and L. R. Carter, “Modeling and forecasting U.S. mortality. ,” Journal of the American Statistical Association, vol. 87, no. 419, pp. 659–671, 1992.
5. H. Chen, R. D. MacMinn, and T. Sun, “Multi-population mortality models:A factor copula approach,” Insurance: Mathematics and Economics, no. 63, pp. 135–146, 2015.
6. J. Patton, “On the Out-of-Sample Importance of Skewness and Asymmetric Dependence for Asset Allocation,” Journal of Financial Econometrics, vol. 1, no. 2, pp. 130–168,2004.
7. D. H. Oh and A. J. Patton, “Modelling Dependence in High Dimensions with Factor Copulas,”Journal of Business and Economic Statistics, vol. 35, no. 1, pp. 139–154, 2017.
8. T. Bollerslev, “Gernalized autoregressive conditonal heteroskedasticity ,” Journal of Econometics, no. 31, pp. 307–327, 1986.
9. Philippe Deprez, Pavel V. Shevchenko, Mario V. Wüthrich,“Machine Learning Techniques for Mortality Modeling”, European Actuarial Journal ,volume 7, pp.337–352(2017)
10. Ronald Richman, Mario V. Wuthrichy,“A Neural Network Extension of the Lee-Carter Model to Multiple Populations”, Annals of Actuarial Science, June 2019
11. Li, Hong & Lu, Yang, 2017. "Coherent Forecasting Of Mortality Rates: A Sparse Vector-Autoregression Approach," ASTIN Bulletin, Cambridge University Press, vol. 47(2), pages 563-600, May.
12. Feng, Lingbing & Shi, Yanlin & Chang, Le, 2021. "Forecasting mortality with a hyperbolic spatial temporal VAR model," International Journal of Forecasting, Elsevier, vol. 37(1), pages 255-273.
13. Guibert, Quentin & Lopez, Olivier & Piette, Pierrick, 2019. "Forecasting mortality rate improvements with a high-dimensional VAR," Insurance: Mathematics and Economics, Elsevier, vol. 88(C), pages 255-272.
14. Hong Li & Yanlin Shi, 2021. "Mortality Forecasting with an Age-Coherent Sparse VAR Model," Risks, MDPI, Open Access Journal, vol. 9(2), pages 1-19, February.