Selecting loss distribution model by using statistical tests should be done to make sure it is accepted and the model selected can be applied in insurance.This study aims to select the best loss distribution model by applying statistical fitting test on seveal fitted Loss distribution model. The loss distribution model represents sample data from motor insurance claims in Malaysia. The model is obtained by using maximum likelihood estimation. The statistical tests used are Pearson: goodness-of fit statistic and Kolmogorov-Smirnov statistic. Based on results from maximum likelihood estimation and statistical tests, a single best model will be selected and the model will be applied in insurance

Ismail, Noriszura

Indonesian

UUM Repository

ANALISIS 8 (1 & 2), 87-100 (2001) PEMILIHAN MODEL KERUGIAN DENGAN UJIAN STATISTIK DAN APLIKASI MODEL DALAM INSURANS NORISZURA HJ .  ISMAIL Fakulti Sains dun Zknologi Un iversiti Kebangsaa n Malaysia ABSTRAK Pmiliban model taburan kerugian dengan ujian statistik perlu dilakukan untuk memustikan agar suatu model boleb diterima dan boleb diaplikasi dahm insurans. h j i a n  ini bertujuun untuk memilib model taburan kerugian terbaik dengun melakukan ujian statistik kesesuaian ke atas beberapa model taburan kerugian yang tehb tersuai. Model taburan kerugian ini mewakili sampel data tuntutan insurans kenderaan bermotor di Mahysia. Model ini diperoleh mehlui kaeahbpenganggaran kebolehjadian maksimum. Dua jenis ujian statistik yang akan dihkukan adulab ujian statistik kebagusan Pearson kbi kuusa-dua dun ujian statistik Kolmogorov-Smirnov. Berdusarkan basil daripada kaedab pengangman kebolebjadian maksimum dun ujian-ujian ini, suatu model tungal yang terbaik akan dipilib dun model ini akan diaplikasi dalam insurans. Kzta kunci: model taburan kerugian, ujian statistik, insurans kenderaan bermotor, kaedab kebolebjudian maksimum. ABSTRACTS Selecting loss distribution model by using statistical tests should be done to make sure it is accepted and the model selected can be applied in insurance. This study aims to select the best loss distribution model by applying statistical fitting test on sevoalfitted Loss distribution modeh. The loss distribution model represents sample data fiom motor insurance claims in Malaysia. Tbe model is obtained by using maximum likelibood estimation. The statistical tests used are Pearson: goodness-ofjt statistic and Kolmogorov-Smirnov statistic. Bused on results from maximum likelibood estimation and statistical tests, a single best model will be selected and the model will be applied in insurance. Keyw0t.h: loss distribution mydel, statistical tests, motor insurance, maxi- mum likelibood. PENGENALAN Penyesuaian taburan tuntutan insurans adalah penting dan perlu dilaksanakan oleh setiap syarikat insurans untuk memastikan agar premium yang dibayar oleh pemegang polisi setimpal dengan tuntutan yang dibayar. Secara tidak langsung, syarikat insurans juga dapat mengurangkan masalah untuk membayar tuntutan kemalangan jalanraya yang semakin meningkat. Oleh kerana terdapat banyak kaedah yang boleh diguna untuk menyesuai model berparameter dan menganggar parameter yang berkait dengannya, kebanyakan penyelidik berhadapan dengan pelbagai pilihan model yang telah tersuai. Kaedah yang sesuai harus diguna untuk mengurangkan lagi bilangan model sehingga akhirnya hanya satu model tunggal terbaik sahaja yang diperoleh. Terdapat dua pendekatan yang boleh diguna untuk memilih model yang terbaik (Klugman et al., 1998). Pendekatan pertama adalah mekanisma benar atau tidak. Untuk pendekatan ini, bagi setiap model yang dicadang, keputusan sama ada ia boleh diterima atau tidak, harus dilakukan. Pendekatan yang kedua adalah menyusun kedudukan model dan memilih model yang berada teratas dalam senarai. Walaupun pendekatan ini adalah lebih baik, pendekatan mekanisma benar atau tidak masih perlu dilakukan kerana akhirnya ianya akan menuju kepada kriteria yang perlu untuk menyusun kedudukan suatu model. Setelah model tunggal yang terbaik diperoleh berdasarkan ujian statistik dan kedudukannya yang teratas di senarai, model ini boleh diaplikasi dalam insurans. Dalam insurans, model tersuai yang diperoleh boleh diguna untuk melihat kesan inflasi, meramal kebarangkalian terjadinya suatu amaun tuntutan dengan penggunaan persentil dan menjangka tempoh terjadinya suatu amaun tuntutan tertentu. Terdapat banyak kajian yang telah dilakukan terhadap ujian statistik kesesuaian model kerugian dan aplikasi model kerugian dalam insurans. Hewitt dan Lefkowitz (1 979) telah menyesuai taburan Gamma, Loggamma, Lognormal, Gamma dan Loggamma dan Gamma dan Lognormal terhadap data kerugian dan menganggar kesan inflasi ke atas deduktibel taburan kerugian yang tersuai. Berdasarkan hasil kajian mereka, tidak semua data kerugian insurans dapat disuai. Namun, melalui pengalaman lepas, kaedah mereka boleh diterima pakai. Klugman dan Parsa (1993) pula telah menyesuai data kerugian dengan kaedah penganggaran jar& minima dan menguji kesesuaian taburan yang tersuai dengan ujian kebagusan penyesuaian. Hasil kajian mereka menunjukkan bahawa penganggaran jarak minimum dapat menepati objektif permasalahan aktuari. Selain itu, Taylor (1 994) telah memodel pengalaman kekerapan dan nisbah tuntutan insurans gadai janji dan membanding serta menguji model yang tersuai. Berdasarkan hasil kajian beliau, model yang tersuai membolehkan liabiliti tuntutan dianggar dan perkembangan nilai harta benda diunjur. Akhir sekali, Brockett et al. (1 995) telah menjana penganggaran ketumpatan kebarangkafian maksimum daripada data histogram dengan min dan median yang diketahui. Hasilnya, prosedur mereka boleh diguna ke atas data berkumpulan dan data sekunder. UJIAN STATISTIK KEBAGUSAN PEARSON Terdapat dua ujian yang boleh dilakukan untuk memastikan agar model kebarangkalian yang dicadang dapat menggambarkan kerawakan sampel setanding dengan taburan kerugian sebenar (Klugman et al., 1998). Kedua- dua ujian ini adalah ujian statistik kebagusan Pearson khi kuasa-dua dan ujian statistik Kolmogorov-Smirnov. Untuk ujian statistik kebagusan Pearson (Hogg & Klugman, 1984), biarkan F(x;8) sebagai fungsi ketumpatan taburan berparameter yang tersuai. Jika FJx) adalah fungsi ketumpatan populasi yang sebenar, kita ingin menguji hipotesis, Ho Fx( X) = F( x;O)  4 : Fx(x) f F(x;O) Sampel rawak XI, 5, ....., X, yang tidak bersandar dicerap dan dibahagi kepada r kumpulan, (co,c,], (cI,c2], ....., (cr-,,cI]. Kekerapan cerapan untuk setiap kumpulan dikira dan nilai-nilainya adalah n,, n2, ....., nr. Untuk hipotesis nol, kebarangkalian untuk kumpulan ke-j adalah, dan nilai jangkaan kekerapan untuk kumpulan ke-j adalah, bilangan ( X  x )  n F " ( 4  = E,  = np, ,  j =  1,2, ....., r dengan, n = bilangan sampel data Di bawah H,, nilai jangkaan kekerapan bagi setiap kumpulan akan bernilai np,, np,, ....., np,. Rumus untuk ujian statistik kebagusan Pearson adalah, Jika n besar, xkira2 dianggar bertaburan khi kuasa-dua dengan darjah kebebasan, v=r-k-1 (r ialah bilangan kumpulan dan k ialah bilangan anggaran parameter taburan). Oleh itu, jika x,.,: adalah kecil, F(x;8) diterima sebagai model yang munasabah. Jika xkir; 2 xa,"' dengan, H, ditolak pada aras keertian a. Untuk memastikan agar ujian ini dapat dipakai, setiap nilai Ej, j=1,2, ....., r, mestilah lebih daripada 5 (Klugman et al., 1998). Jika nilai Ej kurang daripada 5 ,  batas atas perlu dikembang atau lebih kumpulan perlu digabung untuk membentuk satu kumpulan dengan nilai Ej yang melebihi 5. Model taburan yang mempunyai nilai xLra2 yang terkecil merupakan model taburan kerugian yang terbaik. UJIAN STATISTIK KOLMOGOROV-SMIRNOV Ujian statistik kebagusan Pearson tidak sesuai diguna jika bilangan cerapan adalah kecil. Sebagai alternatif, ujian statistik Kolmogorov-Smirnov boleh diguna kerana kaedah ini menggunakan sampel data individu. Untuk kaedah ini, biarkan Fn(x) sebagai fungsi longgokan empirik. Rumusnya ialah (Khgman et a l ,  1998), Ujian statistik Kolmogorov-Smirnov, Dn, ditakrifkan sebagai perbezaan mutlak maksimum di antara fungsi longgokan tersuai, F(x; e), dengan fungsi longgokan empirik, F,(x). Rumusnya adalah, Fn(x) adalah tidak selanjar dan mempunyai lompatan pada setiap titik sampel cerapan, x,, x,, ....., xn. Jika F(x;e) adalah selanjar, nilai maksimum sisihan mesti terjadi sebaik sebelum atau selepas lompatan pada F,(x). Oleh itu, rumus untuk Dn adalah, dengan, {xi} = tertib sampel xl<x,< ..... <xn. Di bawah H,, Dn mempunyai taburan yang tidak berapa dikenali (Hogg & Klugman, 1984). Namun, jadual untuk taburan Dn boleh diperoleh dari beberapa buku statistik. Nilai kritikal, c, boleh diperoleh daripada, Nilai kritikal ini adalah baik jika n215. APLIKASI MODEL DALAM INSURANS Model taburan yang telah tersuai boleh diaplikasi dalam insurans. Antara aplikasi yang boleh dilakukan adalah: (i) melihat kesan inflasi terhadap model taburan yang tersuai; (ii) mengguna persentil untuk meramal kebarangkalian terjadinya suatu amaun tuntutan; (iii) menjangka selang masa terjadinya suatu amaun tuntutan tertentu. Melihat Kesan Inflasi Terhadap Model Taburan Yang Tersuai Kesan inflasi boleh mengakibatkan penokokan dalam kerugian (Hogg & Klugman, 1984). Model yang diperoleh harus diubah suai kepada tahap pengalaman kerugian yang semasa kerana kebanyakan model teranggar adalah berdasarkan cerapan yang dibuat beberapa tahun yang lepas. Selain itu, perunjuran perlu dilakukan untuk menggambarkan kerugian yang dijangka pada masa hadapan. / Biarkan, X = pemboleh ubah rawak kerugian berdasarkan data tuntutan Z = pemboleh ubah rawak kerugian semasa r = kadar inflasi untuk suatu tempoh atau tahun lepas X dan Z boleh dihubungkan dengan persamaan, z=x(1+ r) Untuk mengetahui kesan inflasi, fungsi longgokan yang diinginkan adalah fungsi longgokan Z iaitu (Hogg & Klugman, 1984), F,(Z) = k b ( 2  I 2) = kb(X(1f r )  2 z ) =  kb l + r  l + r  Jika 2 bertaburan selanjar, fungsi ketumpatan untuk 2 adalah, Mengguna Persentil Untuk Meramal Kebarangkalian Terjadinya Suatu Amaun Tuntutan Penganggaran persentil pula berkait dengan penganggaran titik taburan kerugian (Hogg & Klugman, 1984). Melalui penganggaran persentil, kebarangkalian bahawa kerugian melebihi suatu amaun tertentu dapat diperoleh dengan menggunakan fungsi longgokan. Jika kerugian adalah x, rumus untuk anggaran kebarangkalian ini adalah, Lb(X > x )  = 1 - F,  ( x ; 6 )  Menjangka Selang Masa Terjadinya Suatu Amaun Tuntutan Tertentu Untuk menjangka selang masa terjadinya suatu amaun tuntutan tertentu, b iarkan , N(t) = bilangan kejadian kerugian dalam selang (0, t], Oleh itu, N( 1) = N = bilangan kejadian kerugian dalam sehari. Jika N bertaburan Poisson dengan parameter h, h= bilangan kejadian kerugian dalam sehari. Jika n adalah bilangan kejadian kerugian dalam setahun (366 hari), anggaran untuk nilai jangkaan N, E(N), adalah, Untuk meramal nilai jangkaan selang masa (dalam hari) untuk kerugian sekurang-kurangnya sebanyak x, biarkan, K = bilangan kejadian kerugian dalam sehari untuk x>x W = selang masa (dalam hari) kejadian kerugian untuk x>x Oleh itu, K bertaburan Poisson dengan parameter hkb(X>x) dan W bertaburan eksponen dengan parameter hkb(x>x) (Bowers et al., 1997). Anggaran untuk nilai jangkaan selang masa untuk kerugian sekurang- kurangnya x adalah, 1 iZkb(X > x )  i(w)= A A  SKOP DAN METODOLOGI KAJIAN Data tuntutan diperoleh daripada salah sebuah syarikat insurans yang terdapat di Malaysia. Data yang diperoleh adalah sebanyak 6,979 daripada kesemua cawangan syarikat tersebut. Tempoh tuntutan adalah selama 12 bulan iaitu daripada bulan November 1999 hingga Oktober 2000. Untuk memudahkan penganggaran, purata tuntutan untuk setiap cawangan syarikat dan jenis skim insurans bagi setiap bulan telah diambil. Andaian yang diguna adalah, semua tuntutan adalah tuntutan individu dan saling tidak bersandar. Dengan cara ini, bilangan sampel data yang diperoleh dapat dikurangkan kepada sebanvak 459. Dengan kaedah kebolehjadian maksimum, anggaran untuk parameter dan negatif logaritma hngsi kebolehjadian, -in L(e), yang minimum bagi setiap taburan, diberi di dalam Jadual 1 di bawah (Noriszura & Lu Jenn Gio, 200 1). Berdasarkan kesemua 16 taburan daripada Jadual 1 , taburan Burr (a = 1.97, 8 = 8,811.74, y = 2.73), Inverse Burr (z = 0.48, 8 = 8,415.15, y =4.58), Loglogistic (y = 3.23,8=6,351.68) dan Paralogistic (a= 2.56,8= 10,149) memberikan -In L(8) yang paling minimum. Oleh itu, ujian statistik hanya akan diguna ke atas empat taburan ini. Untuk Ujian statistik kebagusan Pearson, kesemua 459 sampel data tuntutan dibahagikan kepada 10 kumpulan iaitu (0, 5001, (500, 20001, (2000y3~OOly (3500,50001, (~ooo,6500], (6500,8500], (8500, 1 10001, (1 1000, 160001, (16000, 50000], (50000,o0). Fungsi longgokan untuk empat taburan ini diberi di dalam Jadual 2 di bawah (Klugman et al., 1998). Jadual2 Fungsi Longgokan untuk Taburan Burr, Inverse Burr, Logistic dan Paralogistic Jadual 1 / Parameter Teranggar dan Nilai -In L(e) yang Minimum bagi Setiap Taburan Taburan Burr 1. Taburan Bilangan Parameter Teranggar Parameter -In LQ) yang minimum f \" 1. Pareto 3 0 = 28,636.83 = 4.05 a = 17.26 a= 1.97 e = 8.8 1 I .74 y = 2.73 = 0.48 e = 8,415.15 = 4.58 T = 21.46 0 = 234.05 0 = 625 1.68 a = 2.56 0 = 10,149 5 = 2.22 y = 3.23 8 = 3,975.47 4,398.2 1 2. Burr 3 4,38 I .44 3. Inverse Burr 3 4,375.5 1 2. Taburan Inverse Burr 4. Inverse Pareto 2 4,569.19 4,389.49 4,382.35 4,428.85 r , x  > 0 5. Loglogistic 2 6. Paralogistic 2 F ( x )  = 7. Inverse Paralogistic 2 8. Gamma Terubahsuai 3 0 = 476.97 5 = 0.70 a = 6.39 a = 3.22 0 = 10,512 a 2.202 c = 0.000000395 r =  1.64 0 = 4,383.37 = 1.23752 8 = 7,172.90 e = 2,227.84 4,403.57 Taburan Loglogistic 3. 9. Gamma 4,406.65 4,504.53 4,440.72 4,548.29 4,534.03 4,571.06 4,419.74 4,456.03 10. Inverse Gamma 1 1. Weibull 12. Inverse Weibull 13. Eksponen 14. Inverse Eksponen 15. Lognormal (3 , x  > 0 F ( x )  = 0 = 4,773.70 p = 8.71 = 0.60 8 = 14,272 = 7,172.90 4. Taburan Paralogistic 16. Inverse Gaussian U , x  > 0 1 F ( x )  = 1 -  1 + ( $  Fungsi longgokan untuk keempat-empat taburan di atas diguna untuk mendapatkan pj, Ej, j=1,2, ....., 10, dan xkira2 masing-masing. Taburan Bagi ujian statistik Kolmogorov-Smirnov pula, fungsi longgokan, F(x;8), untuk keempat-empat taburan seperti yang ditunjukkan di dalam Jadual 2 di atas dan fungsi longgokan empirik, F,(x), diguna untuk mendapatkan Dn masing-masing. D, C Kcputusan Aplikasi model taburan kerugian dalam insurans hanya akan diguna ke atas satu model tunggal terbaik. Model ini dipilih berddarkan nilai yang paling minimum bagi -In L(e), xkira2 dan Dn. HASIL KAJIAN Ujian Statistik Kebagusan Pearson Bagi ujian statistik ini, nilai jangkaan kekerapan kumpulan (0,5001 dan (50000, 00) untuk keempat-empat taburan adalah kurang daripada 5. Oleh itu, kumpulan (0,5001 digabungkan dengan kumpulan (500,20001 manakala kumpulan (1 6000,50000] digabungkan dengan kumpulan (50000, 00). Bilangan kumpulan dikurangkan daripada 10 kepada 8. Jadual3 di bawah menunjukkan nilai xkir;, x pada aras keertian 5%, dan keputusan sama ada untuk menerima atau menolak H, bagi setiap taburan. Jadual3 Nilai xKira2, xa "* dan Keputusan untuk Menerima atau Menolak H, bagi Setiap Taburan Taburan Kep u tusan x K i n  6.807 X, ,,,,'=9.4 8 8 terima H, Inverse Burr 2.712 x,,~,, =9.488 terima H, Loglogistic 22.276 x , , ~ . ~  = 1 1.070 tolak H, -  6.228 ~ , , ~ , ~ ' = l  1.070 terirna H, Paralogistic V i a n  Statistik Kolmogorov-Smirnov Untuk ujian statistik ini pula, Jadual 4 di bawah menunjukkan nilai D,,c pada aras keertian 5%, dan keputusan untuk menerima atau menolak H, bagi setiap taburan. Jadual4 Nilai Dn, C dan Keputusan untuk Menerima atau Menolak H, bagi Setiap Taburan terima H, Inverse Burr Loglogis tic 0.05068 0.06348 terima H, Paralogis tic 0.03498 0.06348 terima H, Melihat Kesan Inflasi Terhadap Model Taburan Yang Tersuai Bagi melihat kesan inflasi, taburan Inverse Burr dengan parameter ~=0.48, e=8,415.15 dan y=4.58 akan dipertimbangkan kerana taburan ini memberikan nilai -In L(e), xkir; dan Dn yang paling minimum berbanding dengan taburan lain. Fungsi ketumpatan yang diinginkan untuk mengetahui kesan inflasi adalah fungsi ketumpatan taburan Z. Oleh itu, rumusnya adalah, 8(1 +r) ,x > 0 Berdasarkan hngsi ketumpatan Z di atas, Z masih lagi bertaburan Inverse Burr, tetapi parameternya bernilai e( 1 +r), z dan 7. Jika kadar inflasi adalah 6% setahun bagi tahun 2000, taburan yang patut diguna untuk tahun 200 1 adalah taburan Inverse Burr dengan parameter 8,=8,920.065, ~,=0.48 dan yz=4.58. Mengguna Persentil Untuk Meramal Kebarangkaian Tejadinya Suatu Amaun Tuntutan K E S I M P W  Sekarang, kebarangkalian terjadinya suatu amaun tuntutan akan diramal. Jika kebarangkalian kerugian melebihi RM20,OOO bagi tahun*2OOO ingin diketahui, kebarangkaliannya adalah, ib(X> 20,000) = 1 -&(20,000) = 1 - = 0.0089797 Oleh kerana parameter T, 8 dan y dianggar dengan kaedah kebolehjadian maksimum, kebarangkalian yang dianggar di atas juga adalah penganggar kebolehjadian maksimum. Menjangka Selang Masa Terjadinya Suatu Amaun Tuntutan Tertentu Katakan nilai jangkaan selang masa (dalam hari), E(W), untuk kerugian melebihi RM20,OOO ingin dijangka. Berdasarkan skop kajian, terdapat 459 bilangan sampel data tuntutan kemalangan kenderaan dalam selang masa 366 hari. Anggaran nilai jangkaan N, E(N), adalah, A 459 A E(N)  =- = 1.254098 = A 366 Oleh itu, anggaran nilai jangkaan selang masa, E(W), untuk kerugian sekurang-kurangnya RM20,OOO ialah, = 88.80 hari 1 ib( X > 20,000) E(W) = Berdasarkan nilai -lnL(8), taburan Inverse Burr memberikan nilai yang paling minimum dan ini diikuti oleh taburan Burr, Paralogistic dan Loglogistic. Berdasarkan ujian kebagusan Pearson pula, hanya taburan Burr, Inverse Burr dan Paralogistic sahaja yang diterima sebagai model. Taburan Inverse Burr sekali lagi memberikan nilai xkir; yang paling mini- mum dan ini diikuti oleh taburan Paralogistic dan Burr. Akhir sekali, berdasarkan ujian Kolmogorov-Smirnov, keempat-empat taburan Burr, Inverse Burr, Loglogistic dan Paralogistic diterima sebagai model. Taburan Inverse Burr sekali lagi memberikan nilai Dn yang paling minimum dan ini diikuti oleh taburan Burr, Paralogistic dan Logistic. Taburan Inverse Burr ( ~ 0 . 4 8 ,  8=8,415.15, ~ 4 . 5 8 )  boleh dianggap sebagai model taburan kerugian yang terbaik untuk mewakili sampel data kerana ianya terletak paling atas dalam senarai berdasarkan penganggaraan kebolehjadian maksimum dan kedua-dua ujian yang dilakukan. RUJUKAN Bowers, N.L., Gerber, H.U., Hickman, J.C., Jones, D.A. & Nesbitt, C.J. ( 1997). Actuarial Mathematics. Schaumburg, Illinois: The Society of Actuaries. Brockett, PL., Cox, S.H., Golany, B., Phillips, F.Y. & Song, Y. (1 995). Actuarial usage of group data: an approach to incorporating secondary data. Society ofActuaries: Transations, 47, 89-1 13. Hewitt, C.C & Lefkowitz, B. (1979). Methods for fitting distributions to insurance loss data. Proc. Casualty Actuarial Society, 66( 123,  Hogg, RV., Klugman, S.A. (1984). Loss Distribution. Canada: John Wiley & Sons. Klugman, S.A., Panjer, H.H., Willmot, G.E. (1998). Loss Models: From Data to Decision. Canada: John Wiley & Sons. Klugman, S.A. & Parsa, A.R. (1 993). Minimum distance estimation of loss distributions. Proc. Casualty Actuarial Society, 80( 152), 250 269. 139- 160. Noriszura Ismail & Lu Jenn Gio (2001). Penyesuaian taburan kerugian: Kajian awal terhadap tuntutan insurans kenderaan di Malaysia. Prosiding Simposium Kebangsaan Sains Matematik ke-3 ( Belum ter b it) . Taylor, G.C. (1  994). Modelling mortgage insurance claims experience: a case study. ASTIN Bulletins, 24( l) ,  97-129. 100 ANALISIS 8 (1 & 2), 87-100 (2001) 

Pemilihan model kerugian dengan ujian statistik dan aplikasi model dalam insurans

http://repo.uum.edu.my/443/1/Noriszura_Hj._Ismail.pdf

Pemilihan model kerugian dengan ujian statistik dan aplikasi model dalam insurans

Abstract

Similar works

Full text

Available Versions

UUM Repository