分層廣義線性模型在準備金評估中的建模研究

閆春 李延星 陳祥輝 邱藝偉

摘 要 考慮到賠付流量三角形數(shù)據(jù)同一事故年反復(fù)觀測的縱向特征以及數(shù)據(jù)結(jié)構(gòu)的層次性，建立了分層廣義線性模型.與通常的隨機模型相比，分層廣義線性模型不但可以選擇條件反應(yīng)變量的分布而且風(fēng)險參數(shù)分布范圍也更加廣泛.利用h-似然函數(shù)估計分層廣義線性模型的模型參數(shù)，降低了計算量.為使模型具有可比性，評估模型的預(yù)測精度，推導(dǎo)了模型預(yù)測誤差的估計式.為充分利用已知賠付信息，將賠付額和賠付次數(shù)兩種賠付信息納入未決賠款準備金評估模型，建立了兩階段分層廣義線性模型.在線性預(yù)測量中考慮了各種固定效應(yīng)和隨機效應(yīng)以及模型結(jié)構(gòu)的散布參數(shù)，改進了線性預(yù)估量結(jié)構(gòu).研究表明：分層廣義線性模型對于數(shù)據(jù)的各種分布及形式都具有很好的適應(yīng)性，更加符合保險實務(wù)現(xiàn)實的賠付規(guī)律.

關(guān)鍵詞 保險數(shù)學(xué)；分層廣義線性模型；h-似然函數(shù)；預(yù)測誤差

中圖分類號 F840.65 文獻標識碼 A

Abstract By considering the longitudinal characteristics of repeated measurements over time of loss for a given accident year in the loss runoff triangles and regarding the loss runoff triangles as hierarchical data， this paper established a hierarchical generalized linear model. Compared with the usual stochastic model， the hierarchical generalized linear model can choose the distribution of conditional response variables and the distribution of risk parameters is more extensive. Using the h-likelihood function to estimate the model parameters of the hierarchical generalized linear model， the calculation amount was reduced. In order to make the model to be comparable， the prediction accuracy of the model was evaluated， and the estimation formula of the model prediction error was derived. In order to make full use of the known payment information， a two-stage hierarchical generalized linear model was established to consider two payment information -the number of compensation and the amount of compensation. The linear prediction considers the distribution parameters of fixed effects， random effects and the distribution parameters of the model structure. The research shows that the hierarchical generalized linear model has good adaptability to all kinds of distribution and form of the data， and it is more consistent with the real payment rule in the insurance practice.

Key words insurance mathematics； hierarchical generalized linear model； h-likelihood function； prediction error

1 引 言

未決賠款準備金的評估方法是非壽險精算的重要組成部分，其評估模型的理論與實踐正在經(jīng)歷著巨大的變革.為維護保險行業(yè)的健康發(fā)展，保監(jiān)會規(guī)定了保險實務(wù)中常用的幾種評估未決賠款準備金的確定性方法.但賠款額本質(zhì)上是一個隨機變量，確定性評估方法只能基于歷史數(shù)據(jù)對未決賠款準備金進行點估計，無法得到估計誤差及波動程度.在全球經(jīng)濟市場環(huán)境不穩(wěn)定、巨額索賠風(fēng)險頻繁發(fā)生的情況下，引入可靠性強、能夠度量準備金波動性以及模擬準備金預(yù)測分布的隨機模型顯得尤為重要.為了實施新的會計準則，保監(jiān)會財務(wù)會計部在2010年1月舉辦了關(guān)于非壽險未決賠款責(zé)任準備金評估方法的系列培訓(xùn)，新出臺的保險合同相關(guān)會計處理規(guī)定鼓勵財險公司采用隨機性模型與方法（如Mack模型、GLM模型、Bootstrap方法等） 評估未決賠款準備金.

但是隨機性模型并非完美無瑕，仍有不少亟需完善的問題.即使經(jīng)典的廣義線性模型（Generalized Linear Models，GLM）仍存在不少缺陷，如對于在某些水平上數(shù)據(jù)量很少的分類解釋變量GLM參數(shù)估計的標準誤差會很大.對于相互獨立且具有層次性結(jié)構(gòu)的數(shù)據(jù)，若直接應(yīng)用GLM會產(chǎn)生過多的待估參數(shù).此外，GLM線性預(yù)估量僅僅包含固定效應(yīng)，忽略了不同事故年賠付數(shù)據(jù)之間的異質(zhì)性以及相同事故年不同進展年間賠付數(shù)據(jù)的相關(guān)性.為了克服了GLM的缺陷，有些學(xué)者引入了廣義線性混合模型（Generalized Linear Mixed Models， GLMM），該模型在線性預(yù)測部分引入了隨機效應(yīng)，不僅考慮了不同事故年間的互異性，還反映了同一事故年不同進展年賠付數(shù)據(jù)由于共同縱向特征而產(chǎn)生的組內(nèi)數(shù)據(jù)的相關(guān)性.近年來GLMM在精算中的應(yīng)用越來越受到國內(nèi)外學(xué)者的重視，如Kelvin（2003）[1]等對索賠頻率分別使用GLM和GLMM建模，通過對比分析指出了GLMM的優(yōu)勢.Antonio和 Beirlant（ 2008）[2]基于GLMM 的半?yún)?shù)回歸模型結(jié)合貝葉斯方法模擬了索賠準備金的預(yù)測分布.Klinker（2011） [3]進一步研究了GLMM在Buhlmann-Straub信度模型理論上的應(yīng)用.謝遠濤和楊娟（2014）[4]基于操作時間來重新設(shè)計流量三角形，并利用GLMM評估準備金.劉新紅和孟生旺（2014）[5]對增量已決賠款建立GAMLSS模型，并將此模型應(yīng)用于一組具有明顯異方差的車險數(shù)據(jù)，證明擬合效果優(yōu)于均值回歸模型.

然而GLMM的主要缺點是隨機效應(yīng)必須服從正態(tài)分布，且隨機效應(yīng)項的引入使得邊際似然函數(shù)的計算更加復(fù)雜，往往需要高階微積分運算，算法實現(xiàn)具有局限性.針對GLMM的缺陷，學(xué)者們進一步將模型推廣為隨機效應(yīng)服從指數(shù)類分布族的分層廣義線性模型（Hierarchical Generalized Linear Models，HGLM）.該模型的基本思想在于模型的某些參數(shù)本身需要建模，分層模型中固定效應(yīng)為通過樣本數(shù)據(jù)直接估計的模型參數(shù)，隨機效應(yīng)則是通過模型的超參數(shù)間接估計的模型參數(shù).GLM等隨機評估方法大多沒有體現(xiàn)賠付數(shù)據(jù)隨事故年反復(fù)觀測的縱向特征，而HGLM作為處理縱向數(shù)據(jù)的一種自然方式，將賠付流量三角形視作分層數(shù)據(jù)，以每個事故年相應(yīng)的賠付數(shù)據(jù)為一個“目標”來評估索賠準備金.近年來有關(guān)廣義線性分層模型的研究不斷涌現(xiàn)，1996年Y. Lee和J.A.Nelder[6]首次提出了HGLM的概念，并于2006[7]年對HGLM進一步推廣給出了雙廣義線性分層模型.

隨后Y. Lee和M. Noh于2012[8]年又給出了雙廣義線性分層模型的隨機效應(yīng)方差.Payne R W（2014）[9]回顧了廣義非線性模型及算法，依據(jù)HGLM相關(guān)理論確定了分層廣義非線性模型.HGLM通過h-似然函數(shù)來估計模型參數(shù)，它是反應(yīng)變量和風(fēng)險參數(shù)聯(lián)合密度函數(shù)的對數(shù).目前國內(nèi)關(guān)于分層模型的研究剛剛起步，其中張連增、段白鴿[9]（2013）依據(jù)損失進展模型LDF模型、Cape Cod模型以及Clark（2003） [10] 提出的兩種體現(xiàn)損失進展過程的非線性增長曲線——Loglogistic增長曲線和Weibull增長曲線建立了非線性分層模型，后來又將其與貝葉斯理論相融合建立了貝葉斯非線性分層模型（2013）[11]，并做了分層模型在非壽險精算學(xué)中應(yīng)用的研究評述（2013）[12]，指出了HGLM評估準備金的研究前景.除此之外，非壽險領(lǐng)域中尚未見更一般意義下的HGLM理論及應(yīng)用研究.

HGLM不但可以處理具有相關(guān)性的縱向數(shù)據(jù)、空間聚類數(shù)據(jù)、甚至更寬泛的聚類數(shù)據(jù)，為不滿足獨立性假設(shè)且具有層次結(jié)構(gòu)的非壽險賠付數(shù)據(jù)提供了處理“大規(guī)模分類”問題的一種自然方式，而且HGLM能提供統(tǒng)一的貝葉斯信度建?？蚣埽N種優(yōu)勢注定了HGLM技術(shù)成為當(dāng)前國際非壽險精算理論研究的熱點領(lǐng)域.隨著國內(nèi)財險公司對GLM在非壽險精算中的應(yīng)用，相信經(jīng)過一段時間的實踐探索，國內(nèi)財險公司必然會迫切需要對更合適準備金評估和非壽險定價的HGLM應(yīng)用研究.鑒于此，建立了分層廣義線性模型，利用h-似然函數(shù)估計分層廣義線性模型的模型參數(shù)，降低了計算量.通過推導(dǎo)模型預(yù)測誤差的估計式來評估模型的預(yù)測精度.利用賠付額和賠付次數(shù)兩種賠付信息建立了未決賠款準備金評估的兩階段分層廣義線性模型，并改進了模型的線性預(yù)估量使之更加符合真實的賠付規(guī)律.提出的未決賠款準備金評估的HGLM，理論上有助于完善未決賠款準備金評估的建模理論和方法，改進模型穩(wěn)定性、適用范圍及預(yù)測能力，豐富未決賠款準備金評估領(lǐng)域內(nèi)的研究內(nèi)容.在實際應(yīng)用價值方面，提出的這一改進模型為精算人員提供構(gòu)建評估模型的新思路與方法，無論對保險公司還是保險監(jiān)管部門都將提供理論支持和實務(wù)參考.

5 未決賠款準備金評估的

兩階段分層廣義線性模型

有了模型的基本結(jié)構(gòu)、參數(shù)估計方法以及預(yù)測誤差，就可以建立未決賠款準備金的評估模型.通常的未決賠款準備金評估模型，往往直接對各單元的累計賠付總額的分布進行假設(shè)，并對累計賠付額進行參數(shù)估計，忽略了索賠次數(shù)和案均賠款等重要數(shù)據(jù)信息.因此，試圖建立一個兩階段分層廣義線性模型，該模型分別對索賠次數(shù)與案均賠款構(gòu)造分層廣義線性模型，更符合非壽險未決賠款準備金評估的實際需求.此外，改進了模型的線性預(yù)估量，使之不但符合隨進展年的推延先按冪升再按指數(shù)降的賠付規(guī)律，而且減少了待估計參數(shù)得的個數(shù).

6 結(jié) 論

在分層廣義線性模型的全新視角下建立了未決賠款準備金的評估模型，鑒于流量三角形層次性和相關(guān)性的數(shù)據(jù)結(jié)構(gòu)，在線性預(yù)測量中考慮了固定效應(yīng)和隨機效應(yīng)，并且拓展了風(fēng)險參數(shù)的分布范圍.利用h-似然函數(shù)估計分層廣義線性模型的模型參數(shù)，減少了估計過程的計算量.利用參數(shù)誤差與過程誤差的和作為均方誤差來評估模型的預(yù)測精度，使得模型的精確程度具有可比性.此外，將賠付額和賠付次數(shù)兩種賠付信息納入評估模型，建立了兩階段分層廣義線性模型，在線性預(yù)測量中考慮了各種固定效應(yīng)和隨機效應(yīng)以及模型結(jié)構(gòu)的散布參數(shù)，改進了線性預(yù)估量結(jié)構(gòu)，并充分利用了已知賠付信息，這使得分層廣義線性模型能夠更好的適應(yīng)各類數(shù)據(jù)的結(jié)構(gòu)與形式，具有更高的靈活性，在保險實務(wù)中可以處理各種賠付規(guī)律的案件，較之傳統(tǒng)的廣義線性模型具有更高的實用性.

參考文獻

[1] KELVIN K W， YAU Karen C H，YIP H K YUEN. Modelling repeated insurance claim frequency data using the generalized linear mixed model [J]. Journal of Applied Statistics， 2003， 30（8）： 857-865.

[2] ANTONIO K， BEIRLANT J. Issues in claims reserving and credibility： a semiparametric approach with mixed models[J]. Journal of Risk and Insurance，2008，75（ 3） ： 643-676.

[3] KLINKER F. Generalized linear mixed models for ratemaking： a means of introducing credibility into a generalized linear model setting [J]. Casualty Actuarial Society E-Forum， 2011（2）：1-25.

[4] 謝遠濤，楊娟. 基于操作時間和廣義線性混合模型的準備金評估技術(shù)研究[J]. 保險研究，2014（3）：54-62.

[5] 劉新紅，孟生旺. 基于藤Copula的GAMLSS模型與非壽險準備金評估[J]. 經(jīng)濟數(shù)學(xué)，2014（4）：68-74.

[6] LEE Y， NELDER J A. Hierarchical generalized linear models[J].Journal of the Royal Statistical Society，Series B，1996，58 （ 4 ） ：619-678.

[7] GELMAN A， HILL J. Data analysis using regression and multilevel hierarchical models[M].Cambridge： Cambridge University Press， 2007.

[8] LEE Y， NOH M. Modelling random effect variance with double hierarchical generalized linear models[J]. Statistical Modelling， 2012， 12（6）：487-502.

[9] 段白鴿，張連增. 索賠準備金評估的非線性分層增長曲線模型研究[J]. 財經(jīng)理論與實踐，2013（4）：23-29.

[10]CLARK D R. LDF Curve-fitting and stochastic reserving： a maximum likelihood approach [J]. CAS Forum， 2003（3）：41-92.

[11]段白鴿，張連增. 索賠準備金評估的貝葉斯非線性分層模型[J]. 山西財經(jīng)大學(xué)學(xué)報，2013（10）：20-31.

[12]段白鴿，張連增. 分層模型在非壽險精算學(xué)中的應(yīng)用研究評述[J]. 統(tǒng)計研究，2013（05）：98-105.

[13]BOOTH J G， HOBERT J P. Standard errors of prediction in generalized linear mixed models[J]. Journal of the American Statistical Association， 2012， 93（441）：262-272.

[14]LEE Y， HA I D.Orthodox BLUP versus h-likelihood methods for inferences about random effects in Tweedie mixed models. Statistics and Computing 2010， 20（3）： 295-303.