Preskoči na vsebino
Petka Zavarovanja – logotipPETKAZavarovanja
GAMLSS modeli: AIC, BIC in diagnostika za optimalno oceno škod
Nezgoda
  • Tehnični vpogledi
Nezgoda in poškodbe

GAMLSS modeli: AIC, BIC in diagnostika za optimalno oceno škod

Kot zavarovalna strokovnjakinja z 20-letnimi izkušnjami se pogosto srečujem z izzivom natančnega modeliranja škodnih dogodkov, še posebej pri kompleksnih podatkih nezgodnih zavarovanj. V tem članku se posvečam ključnim statističnim kriterijem in diagnostičnim testom, ki nam pomagajo izbrati optimalne GAMLSS modele.

Petra Guštin · 15 min branja ·

Zadnjič posodobljeno:

Kako delam z vami: pripravim natančno primerjavo pogojev, izračune kritij in pregledno specifikacijo police — vse s preverljivimi viri in številkami.

Hitra Petka — 5 točk za hitro branje

Bistvo objave v 30 sekundah.

  • GAMLSS modeli omogočajo prilagodljivo modeliranje kompleksnih distribucij škod.
  • AIC in BIC kriterija uravnotežita prileganje modela in njegovo kompleksnost.
  • Nižje vrednosti AIC in BIC so indikatorji boljšega modela.
  • Diagnostični testi residuov so ključni za preverjanje predpostavk in robustnosti modela.
  • Optimalna izbira modela temelji na kombinaciji informacijskih kriterijev in vizualne diagnostike.

Uvod v GAMLSS modele in pomen natančnega modeliranja škod

V svetu aktuarske znanosti in zavarovalništva je natančno modeliranje škodnih dogodkov temelj za določanje primernih premij, obvladovanje tveganj in načrtovanje rezervacij. Z nezgodnimi zavarovanji, kjer se srečujemo z raznolikimi in pogosto asimetričnimi distribucijami škod, postanejo tradicionalni regresijski modeli (npr. splošni linearni modeli – GLM) pogosto nezadostni. Ti namreč pogosto predpostavljajo določene oblike distribucije (npr. normalno) in konstanten razpon variance, kar v resničnem svetu ni vedno izpolnjeno. Zato se v zadnjih dveh desetletjih vse pogosteje obračamo k bolj fleksibilnim pristopom.

Prav tu pridejo do izraza splošni aditivni modeli lokacije, skale in oblike (Generalized Additive Models for Location, Scale and Shape – GAMLSS), ki jih je leta 2005 predstavil Rigby in Stasinopoulos. GAMLSS modeli nam omogočajo, da poleg povprečja (lokacije) modeliramo tudi druge parametre distribucije, kot so varianca (skala) in asimetrija ter sploščenost (oblika) v odvisnosti od prediktorjev. To pomeni, da lahko model natančneje zajame kompleksno strukturo podatkov in se prilagodi širokemu spektru distribucij, vključno s tistimi, ki so močno asimetrične ali imajo debeli rep, kar je pogosto značilno za škodne dogodke. Moja izkušnja potrjuje, da ta fleksibilnost bistveno izboljša napovedno moč in robustnost modelov, kar je ključno za zanesljive aktuarske ocene.

Razumevanje in pravilna uporaba GAMLSS modelov mi omogočata, da pri delu z zavarovalnicami in strankami ponudim poglobljeno analizo tveganj. Cilj tega članka je predstaviti ključne kriterije za evalvacijo in primerjavo različnih GAMLSS modelov, saj ni dovolj le zgraditi model, temveč tudi preveriti njegovo ustreznost in izbrati najboljšega izmed možnih. Osredotočila se bom na informacijske kriterije AIC in BIC ter na diagnostične teste residuov, ki so nepogrešljivi pri validaciji modela in preverjanju njegovih predpostavk. Te metode so neprecenljive za zagotavljanje, da so naši modeli ne le teoretično ustrezni, ampak tudi praktično uporabni in zanesljivi v realnih scenarijih.

Temeljni principi GAMLSS modelov in njihova prednost pri škodnih dogodkih

Osnovna ideja GAMLSS modelov je modeliranje vseh parametrov distribucije (npr. povprečja μ, variance σ², asimetrije ν, sploščenosti τ) kot funkcij prediktorjev, kar omogoča izjemno fleksibilnost. Za razliko od GLM, ki modelirajo le povprečje kot funkcijo prediktorjev, GAMLSS omogočajo, da so vsi parametri distribucije (ne le povprečje) odvisni od prediktorjev, in sicer preko povezovalnih funkcij (link functions). To pomeni, da lahko na primer variabilnost škodnih dogodkov niha glede na starost zavarovanca ali tip nezgode, kar je v praksi zelo pogosto. Na ta način GAMLSS modeli precej bolje zajamejo kompleksnost podatkov in morebitne heteroskedastičnosti.

GAMLSS modeli mi omogočajo, da delam z zelo širokim naborom distribucij – ne le z eksponentno družino (kot pri GLM), temveč z več kot 100 različnimi parametriziranimi distribucijami. To vključuje distribucije z debelimi repi, asimetrične distribucije in tiste, ki so primerne za diskretne ali zvezne podatke. Pri modeliranju škodnih dogodkov, kjer so zneski pogosto asimetrično porazdeljeni in imajo precejšnje ekstreme, je to ključnega pomena. Na primer, za zneske škod se pogosto uporabljajo distribucije, kot so Gamma, inverzna Gaussova, Weibull ali celo bolj kompleksne, kot so zmešane distribucije ali distribucije z ničelnim inflacijo. Možnost izbire najustreznejše distribucije za dane podatke bistveno poveča natančnost modela in zmanjša pristranskost ocene.

Poleg tega GAMLSS modeli podpirajo aditivne prediktorje, kar pomeni, da lahko vključimo gladke funkcije (splines) za modeliranje nelinearnih odnosov med prediktorji in odzivno spremenljivko. To je še posebej koristno, ko ne poznamo točne oblike odnosa, na primer med starostjo zavarovanca in pogostostjo škod. V moji praksi sem opazila, da lahko to zmanjša tveganje napačne specifikacije modela in s tem izboljša njegovo napovedno moč. Kombinacija teh lastnosti – fleksibilnost distribucij, modeliranje vseh parametrov distribucije in uporaba gladkih funkcij – postavlja GAMLSS modele v ospredje pri modeliranju kompleksnih aktuarskih podatkov.

Informacijski kriteriji: AIC (Akaike Information Criterion)

Ko imamo na voljo več kandidatnih GAMLSS modelov, se soočimo z vprašanjem, kateri je 'najboljši'. Pri tem nam pomagajo informacijski kriteriji, ki poskušajo uravnotežiti prileganje modela podatkom (likelihood) in njegovo kompleksnost (število parametrov). Eden izmed najbolj razširjenih kriterijev je Akaike Information Criterion (AIC). Moja izkušnja kaže, da je AIC pogosto prva izbira pri primerjavi modelov, saj je relativno enostaven za izračun in interpretacijo.

Formula za AIC je naslednja: AIC = -2 * log(L) + 2 * k, kjer je L vrednost funkcije verjetnosti (likelihood) modela, k pa število parametrov v modelu. Nižja vrednost AIC indikator boljšega modela. Pomembno je razumeti, da AIC ne meri absolutne 'dobrote' modela, temveč omogoča primerjavo med modeli. Model z najnižjo vrednostjo AIC je tisti, ki je po Akaikejevem mnenju najboljši kompromis med prileganjem in kompleksnostjo. Vendar moramo biti previdni; AIC ima tendenco izbirati bolj kompleksne modele, še posebej pri velikih vzorcih, saj je 'kazen' za število parametrov razmeroma majhna.

Pri uporabi AIC je nujno, da primerjamo le modele, ki so vgnezdni ali vsaj modelirani na istih podatkih. Čeprav se AIC široko uporablja, sem v praksi spoznala, da ni edini in vedno optimalni kriterij. Zato ga vedno dopolnjujem z drugimi kriteriji in diagnostiko, da dobim celovito sliko o ustreznosti modela. Prav tako sem ugotovila, da pri manjših vzorcih obstaja popravljena verzija, AICc (AIC corrected), ki je bolj primerna za takšne situacije, saj bolje obravnava manjšo velikost vzorca in zmanjšuje tendenco k prekomerni kompleksnosti modela.

Informacijski kriteriji: BIC (Bayesian Information Criterion)

Drug pomemben informacijski kriterij je Bayesian Information Criterion (BIC), znan tudi kot Schwarz Information Criterion (SIC). BIC je, podobno kot AIC, namenjen primerjavi modelov in izbiri najboljšega. Vendar pa BIC uvaja strožjo 'kazen' za število parametrov kot AIC, kar pomeni, da ima tendenco izbirati bolj parsimoniozne (enostavnejše) modele, še posebej pri velikih vzorcih podatkov. Moje izkušnje kažejo, da je BIC še posebej uporaben, ko želimo resnično zmanjšati tveganje overfittinga in preferiramo robustnejši, a morda nekoliko manj natančno prilegajoč model, ki se bolje generalizira na nove podatke.

Formula za BIC je: BIC = -2 * log(L) + k * log(n), kjer je L vrednost funkcije verjetnosti modela, k število parametrov in n število opazovanj (velikost vzorca). Vidimo, da se 'kazen' za število parametrov (k * log(n)) povečuje z velikostjo vzorca (n) hitreje kot pri AIC (2 * k). Zato bo BIC pogosteje izbral model z manj parametri, kar lahko privede do bolj robustnega in interpretibilnega modela, ki se manj nagiba k prekomernemu prileganju šumu v podatkih. Tako kot pri AIC, tudi pri BIC velja, da nižja vrednost označuje boljši model.

V praksi se pogosto zgodi, da AIC in BIC izbereta različne modele. Takrat je pomembno, da razmislimo o kontekstu in ciljih modeliranja. Če je naš cilj predvsem napovedovanje in nas ne skrbi prevelika kompleksnost, bi se morda odločili za model z nižjim AIC. Če pa želimo razumeti osnovne procese, identificirati ključne prediktorje in imeti bolj interpretibilen model, je model z nižjim BIC pogosto boljša izbira. Medsebojna primerjava obeh kriterijev mi pomaga pri bolj premišljeni odločitvi, saj mi omogoča vpogled v kompromis med prileganjem in kompleksnostjo. Včasih je treba razmisliti tudi o kombinaciji teh pristopov, na primer z uporabo AIC za izbiro družine modelov in nato BIC za natančnejšo selekcijo znotraj te družine.

Primerjava AIC in BIC: Kateri kriterij izbrati in zakaj?

Odločitev med AIC in BIC ni vedno enostavna in je odvisna od specifičnih ciljev analize. Kot sem že omenila, AIC nagrajuje bolj kompleksne modele in je uporaben, ko je cilj maksimalna napovedna natančnost, tudi če to pomeni nekoliko manj interpretibilen model. BIC, po drugi strani, bolj kaznuje kompleksnost in je preferiran, ko iščemo bolj parsimoniozen, stabilen in interpretibilen model, ki se bo dobro obnesel tudi na novih, nevidnih podatkih. To je še posebej pomembno v zavarovalništvu, kjer je robustnost modelov ključna za dolgoročno profitabilnost in obvladovanje tveganj.

Matematično gledano, razlika v kaznovanju parametrov je v členu: 2 * k (za AIC) proti k * log(n) (za BIC). Za n > 7 bo log(n) > 2, kar pomeni, da bo BIC vedno bolj kaznoval kompleksnost kot AIC, ko je število opazovanj večje od 7. V aktuarski praksi imamo običajno veliko število opazovanj, zato je ta razlika pomembna. To pogosto pomeni, da bo BIC izbral enostavnejši model. Moja izkušnja je, da je pri velikih naborih podatkov, značilnih za zavarovalništvo, BIC pogosto boljša izbira za zmanjšanje tveganja prekomernega prileganja in izbiro modela, ki se bolje generalizira.

Kljub temu pa nobeden od kriterijev ni absolutno 'pravilen' ali 'napačen'. Če se odločiva za model, ki ga izbere AIC, je smiselno dodatno preveriti njegovo stabilnost in robustnost, na primer z validacijo na ločenem testnem naboru podatkov. Če pa izberemo model, ki ga predlaga BIC, se prepričamo, da ta še vedno zajema dovolj informacij iz podatkov in da ni preveč poenostavljen. Pogosto uporabljam oba kriterija in iščem ravnovesje. Model, ki ima relativno nizko vrednost obeh kriterijev in hkrati izkazuje dobre diagnostične lastnosti residuov, je pogosto najboljša izbira. Včasih je smiselno razmisliti tudi o drugih kriterijih, kot je AICC (popravljen AIC za majhne vzorce), vendar se v zavarovalništvu, kjer so vzorci običajno veliki, ta potreba pojavlja redkeje.

Diagnostični testi residuov: preverjanje predpostavk modela

Izbira modela na podlagi informacijskih kriterijev je pomemben korak, vendar ni dovolj. Ključno je tudi preveriti, ali model resnično ustreza podatkom in ali so izpolnjene njegove predpostavke. Pri GAMLSS modelih to pomeni diagnostiko residuov. Residuumi so razlike med dejanskimi in napovedanimi vrednostmi, in če je model dober, naj bi bili ti resuduumi naključni, neodvisni in ustrezali določeni distribuciji (običajno standardizirani normalni distribuciji, odvisno od transformacije). Kot strokovnjakinja vem, da je ta faza kritična za preprečevanje zmotnih zaključkov.

Ena izmed ključnih metod diagnostike so kvantilni-kvantilni (Q-Q) grafi standardiziranih residuov. Ti grafi primerjajo kvantile opazovanih residuov s kvantili teoretične distribucije (npr. standardizirane normalne distribucije). Če se točke na grafu približno nahajajo vzdolž referenčne diagonale (y=x), to kaže na dobro prileganje modela in skladnost residuov s predpostavljeno distribucijo. Odstopanja od diagonale nakazujejo na težave, kot so neustrezna distribucija, prisotnost odstopajočih vrednosti ali napačno specifikacijo modela. Vedno preverjam tudi graf razpršenosti residuov glede na napovedane vrednosti, da bi identificirala morebitno heteroskedastičnost ali nelinearnosti, ki jih model ni zajel. Idealen scenarij je naključen razpršenost brez vidnih vzorcev.

Poleg vizualne diagnostike uporabljamo tudi formalne statistične teste. Eden takšnih je Shapiro-Wilk test, ki testira normalnost distribucije residuov. Ničelna hipoteza tega testa je, da so resuduumi normalno porazdeljeni. Če je p-vrednost testa nizka (npr. < 0.05), zavrnemo ničelno hipotezo in sklepamo, da resuduumi niso normalno porazdeljeni. Vendar je pomembno biti previden; pri velikih vzorcih je test Shapiro-Wilk zelo občutljiv in pogosto zavrne normalnost, tudi če so odstopanja minimalna. Zato ga vedno interpretiram v povezavi z vizualno diagnostiko in splošnim razumevanjem podatkov. Drugi testi, kot so test homoskedastičnosti (npr. Breusch-Pagan test) ali testi avtokorelacije (npr. Durbin-Watson test), so prav tako relevantni, odvisno od specifičnih predpostavk modela in narave podatkov. Vsak od teh testov mi pomaga zgraditi celovito sliko o ustreznosti in robustnosti modela.

Napredna diagnostika: Z-score residuumi in robusnost modela

Pri GAMLSS modelih je diagnostika standardiziranih residuov še bolj niansirana. Ker GAMLSS modeli omogočajo modeliranje vseh parametrov distribucije, lahko ustvarimo t.i. 'z-score' residuume, ki so enakovredni standardiziranim residuumom v GLM. Ti z-score residuumi so izračunani tako, da so idealno neodvisni in sledijo standardni normalni distribuciji (N(0,1)). Preverjanje te predpostavke je ključnega pomena. Če z-score residuumi odstopajo od standardne normalne distribucije, to kaže na morebitno neustrezno izbiro distribucije za odzivno spremenljivko ali napačno specifikacijo modela parametrov distribucije (npr. skale ali oblike).

Za grafično analizo teh z-score residuov pogosto uporabljam t.i. 'worm plots' (Weight Residuals for Model Checking). Ti grafi so razširitev Q-Q grafov in so še posebej uporabni za vizualno preverjanje prileganja z-score residuov standardni normalni distribuciji. Worm plot razdeli residuume v več intervalov (običajno glede na vrednosti prediktorjev) in za vsak interval prikaže de-trendirani Q-Q graf. Če je model dober, naj bi vse krivulje na worm plotu ostale znotraj dveh standardnih odklonov od ravne črte (ki predstavlja idealno prileganje). Odstopanja kažejo na lokalne težave s prileganjem modela v določenih segmentih podatkov, kar mi pomaga natančno določiti, kje je model pomanjkljiv.

Poleg vizualne diagnostike, robustnost modela preverjamo tudi z analizo vplivnih točk in odstopajočih vrednosti. V GAMLSS kontekstu lahko posamezne opazovalne enote močno vplivajo na ocene parametrov. Identifikacija in obravnava teh točk je pomembna za zmanjšanje pristranskosti in povečanje stabilnosti modela. Uporabljamo različne metrike, kot so Cookova razdalja ali DFFITS, za oceno vpliva posameznih opazovanj. V moji praksi sem ugotovila, da je kombinacija vizualnih in kvantitativnih metod ključna za celovito in zanesljivo evalvacijo GAMLSS modelov, saj nam omogoča, da vidimo tako 'veliko sliko' kot tudi specifične probleme znotraj modela. To mi omogoča, da z zaupanjem predstavim rezultate in priporočila strankam.

Interpretacija rezultatov in izbira optimalnega modela

Po izračunu informacijskih kriterijev in izvedbi diagnostičnih testov se soočim z interpretacijo rezultatov. Idealno bi bilo, da bi imel en model najnižji AIC in BIC, ter da bi njegovi residuumi izkazovali dobro prileganje na Q-Q grafih in da bi bili statistični testi za residuume ne-zavrnljivi. Vendar se v praksi to pogosto ne zgodi. Lahko se zgodi, da AIC izbere en model, BIC pa drugega. Takrat je potrebna premišljena presoja in razumevanje specifičnih potreb analize.

Kot sem že poudarila, je pomembno, da poleg numeričnih kriterijev upoštevamo tudi strokovno znanje. V zavarovalništvu se to nanaša na aktuarsko intuicijo in razumevanje narave tveganj. Ali je model, ki ga je izbral AIC (morda bolj kompleksen), smiseln z aktuarskega stališča? Ali parametri, ki jih ocenjuje, odražajo pričakovane odnose med dejavniki tveganja in škodami? Če model s sicer nižjim AIC ali BIC prinaša neintuitivne rezultate, je to rdeča zastavica, ki mi pravi, da moram dodatno raziskati in preveriti specifikacijo modela ali podatke.

Končna izbira optimalnega GAMLSS modela je torej iterativen proces, ki vključuje: 1) postavitev več kandidatnih modelov, 2) izračun AIC in BIC, 3) izvedbo obsežne diagnostike residuov (Q-Q grafi, worm plots, statistični testi), 4) preverjanje robustnosti modela (npr. vplivne točke), in 5) aktuarsko interpretacijo rezultatov. Cilj ni le najti model, ki se matematično najbolje prilega podatkom, ampak tudi takšnega, ki je smiseln, robusten in nam omogoča zanesljive napovedi in razumevanje procesov, kar je temelj za strateške odločitve v zavarovalništvu. Le tako lahko dosežemo tisto stopnjo natančnosti, ki jo od nas pričakujejo regulativni okviri, kot so Zavarovalniška regulativa (npr. Solvency II direktiva, ki je v Sloveniji implementirana preko ZZavar-1) in dobra aktuarska praksa.

Kaj je krito in kaj ni krito pri nezgodnem zavarovanju (relevantno za modeliranje škod)

Pri modeliranju škodnih dogodkov, še posebej z GAMLSS modeli, je ključno razumevanje, kaj sploh spada v obseg 'škode', ki jo obravnavamo. Splošni pogoji nezgodnih zavarovanj, ki so določeni s strani posamezne zavarovalnice (v skladu z Zakonom o zavarovalništvu – ZZavar-1 in Zakonom o zavarovalnih pogodbah – ZZavPog), jasno določajo, kateri dogodki so kriti in kateri ne. Moja vloga je, da stranke seznanim s temi podrobnostmi, saj je od tega odvisna tudi pravilna kategorizacija in zbiranje podatkov za modeliranje.

**Kaj je običajno krito:**

* **Smrt kot posledica nezgode:** Zavarovalnica izplača dogovorjeno zavarovalno vsoto dedičem. Pri tem je ključna definicija nezgode, kot nenadnega, od zavarovančeve volje neodvisnega dogodka, ki je povzročil smrt. Pogosto se modelira frekvenca in verjetnost takšnih dogodkov.

* **Trajna invalidnost kot posledica nezgode:** Zavarovalnica izplača določen odstotek zavarovalne vsote glede na stopnjo invalidnosti, ki je določena po lestvici, specifični za posamezno zavarovalnico. Ta element je pogosto predmet kompleksnega modeliranja, saj vključuje oceno stopnje invalidnosti in zneska odškodnine.

* **Dnevna odškodnina za čas zdravljenja:** Izplačilo določenega zneska za vsak dan zdravljenja zaradi nezgode (bolnišnično zdravljenje, zdravljenje na domu). Tukaj modeliramo trajanje in višino dnevnih odškodnin.

* **Stroški zdravljenja in reševanja:** Kritje stroškov, ki niso pokriti iz obveznega zdravstvenega zavarovanja (ZZVZZ), kot so prevoz, specialistični pregledi, rehabilitacija. Modeliramo variabilnost teh stroškov.

* **Bolezenska odškodnina (če je vključena):** Nekatere police vključujejo tudi kritje za določene bolezni, ki so posledica nesreče ali so s pogodbo izrecno določene. To je običajno razširitev standardnega nezgodnega zavarovanja.

**Kaj običajno ni krito:**

* **Bolezni, ki niso posledica nezgode:** Razen če je vključeno dodatno kritje za specifične bolezni, standardno nezgodno zavarovanje ne krije bolezni. Na primer, srčni infarkt ali možganska kap, če nista neposredna posledica nezgode, nista krita. Modeliranje GAMLSS je v tem kontekstu osredotočeno na dejanske nezgode, ne na bolezni. Zakon o pokojninskem in invalidskem zavarovanju (ZPIZ-2) določa splošno invalidnost, ki je lahko tudi posledica bolezni, vendar se to področje loči od nezgodnega zavarovanja.

* **Nezgode, ki so posledica določenih aktivnosti:** Sem spadajo določeni ekstremni športi, vojne aktivnosti, namerne samopoškodbe ali nezgode pod vplivom alkohola/drog. Podatki o takšnih dogodkih so običajno izključeni iz učnih naborov za GAMLSS modele ali pa so obravnavani kot posebni primeri.

* **Nezgode zaradi malomarnosti ali hude malomarnosti:** V nekaterih primerih, če je bila nezgoda povzročena z izjemno malomarnostjo, zavarovalnica lahko zmanjša ali zavrne izplačilo.

* **Škoda na premoženju:** Nezgodno zavarovanje je namenjeno kritju poškodb osebe, ne njenega premoženja. Za premoženje obstajajo druga zavarovanja.

Pri GAMLSS modeliranju moramo vedno jasno definirati obseg podatkov in pogojev kritja, da se izognemo napačnim napovedim in napačni interpretaciji rezultatov. Upoštevanje teh določb mi omogoča, da zgradim robustnejše modele, ki dejansko odražajo realno zavarovalno tveganje.

Praktični primer: Modeliranje stroškov nezgodnih škod

Dovolite mi, da ponazorim uporabo GAMLSS modelov in evalvacijskih kriterijev s poenostavljenim praktičnim primerom iz moje prakse, seveda brez kakršnihkoli konkretnih številk ali imen. Recimo, da sem imela nalogo modelirati višino stroškov zdravljenja po nezgodi (v EUR) za zavarovance določene zavarovalnice. Podatki so pokazali močno asimetrično distribucijo z dolgim desnim repom, kar je značilno za zdravstvene stroške – večina primerov ima nizke stroške, medtem ko je nekaj primerov z zelo visokimi stroški. Tradicionalni GLM z Gaussovo distribucijo (ki predpostavlja simetrijo) bi bil v tem primeru neustrezen.

Odločila sem se, da preizkusim več GAMLSS modelov z različnimi distribucijami za stroške (npr. Gamma, inverzna Gaussova, log-normalna) in različnimi strukturami prediktorjev (npr. starost zavarovanca, spol, tip nezgode, kraj nezgode – nekateri kot gladke funkcije, nekateri kot kategorijske spremenljivke), vplivajočimi na parametre lokacije (μ) in skale (σ). Po izgradnji petih kandidatnih modelov, sem izračunala AIC in BIC za vsakega. Model A (Gamma distribucija, starost kot gladka funkcija za μ in σ) je imel AIC=15678 in BIC=15752. Model B (inverzna Gaussova distribucija, starost kot gladka funkcija samo za μ) je imel AIC=15690 in BIC=15745. Model C (log-normalna distribucija, starost in tip nezgode za μ) pa AIC=15685 in BIC=15720.

Na podlagi informacijskih kriterijev se je Model C z log-normalno distribucijo in enostavnejšo strukturo za σ (konstantna) izkazal za najboljšega po BIC, medtem ko je Model A imel najnižji AIC. To je bil klasičen primer, kjer sta se AIC in BIC razšla. Pri pregledu diagnostike residuov sem opazila, da so z-score residuumi Modela C na worm plotu izkazovali boljše prileganje standardni normalni distribuciji v vseh segmentih kot Model A, ki je kazal rahla odstopanja v segmentih z višjimi napovedanimi vrednostmi. Tudi Shapiro-Wilk test je bil manj verjeten za zavrnitev normalnosti residuov pri Modelu C. Čeprav je imel Model A nekoliko nižji AIC, sem se zaradi boljše diagnostike residuov, bolj parsimonioznega modela po BIC in bolj robustnega prileganja na kritičnih delih distribucije odločila za Model C. Ta odločitev mi je omogočila zgraditi zanesljiv model za oceno višine stroškov zdravljenja po nezgodi, ki je bil nato uporabljen za optimizacijo premij in rezervacij. To je resnično primer, kako se aktuarska znanost in praktične odločitve prepletajo.

Regulativni okviri in vpliv na modeliranje škod

V zavarovalništvu nisem omejena le na izbiro najboljšega statističnega modela, temveč moram delovati tudi znotraj strogih regulativnih okvirov. V Sloveniji so to predvsem Zakon o zavarovalništvu (ZZavar-1), Zakon o zavarovalnih pogodbah (ZZavPog) in na širšem evropskem nivoju direktiva Solvency II. Ti predpisi ne določajo neposredno, katere statistične modele naj uporabljamo, določajo pa standarde glede zanesljivosti, natančnosti in transparentnosti aktuarskih izračunov, kar posredno vpliva na izbiro in evalvacijo GAMLSS modelov.

ZZavar-1, na primer, določa zahteve glede kapitalne ustreznosti zavarovalnic, izračuna tehničnih rezervacij in upravljanja s tveganji. Pri izračunu tehničnih rezervacij je ključno zanesljivo ocenjevanje bodočih škodnih dogodkov. Model, ki ga uporabljamo, mora biti dovolj robusten in natančen, da zagotavlja ustrezno višino rezervacij in s tem solventnost zavarovalnice. Nezadostna evalvacija modela, ki bi privedla do podcenjenih rezervacij, bi lahko imela resne finančne in regulativne posledice. Zato je izbira optimalnega GAMLSS modela, podkrepljena z AIC, BIC in diagnostičnimi testi, ključna za skladnost z regulativo.

Agencija za zavarovalni nadzor (AZN) kot regulator v Sloveniji redno preverja aktuarske metodologije in modele zavarovalnic. Uporaba naprednih tehnik, kot so GAMLSS modeli, in transparentna predstavitev evalvacijskih kriterijev (AIC, BIC, diagnostika residuov) prispevata k večji verodostojnosti in skladnosti z regulatornimi pričakovanji. Kot strokovnjakinja v zavarovalništvu, se zavedam, da ni dovolj le razviti statistično eleganten model, ampak mora biti ta model tudi regulativno sprejemljiv in ga moram biti sposobna prepričljivo predstaviti nadzornim organom. To pomeni, da morajo biti vsi koraki, od izbire distribucije do evalvacije, skrbno dokumentirani in utemeljeni.

Prihodnji trendi: Umetna inteligenca in strojno učenje v kombinaciji z GAMLSS

Razvoj na področju umetne inteligence (AI) in strojnega učenja (ML) prinaša nove možnosti tudi v aktuarski znanosti in modeliranju škod. Čeprav so GAMLSS modeli že sami po sebi zelo fleksibilni, se pojavljajo hibridni pristopi, ki združujejo moč GAMLSS z nekaterimi tehnikami ML, predvsem za izboljšanje napovedne moči in obvladovanje izjemno kompleksnih interakcij v podatkih. Menim, da je to smer, v katero se bo aktuarska znanost v prihodnosti vse bolj razvijala.

Ena izmed možnih smeri je uporaba ML algoritmov (npr. gradient boosting trees ali nevronske mreže) za identifikacijo in transformacijo pomembnih prediktorjev, ki se nato vključijo v GAMLSS model. Na ta način lahko ML pomaga pri 'inženiringu značilnosti' (feature engineering) in s tem izboljša prileganje in interpretibilnost GAMLSS modela. Drugi pristopi vključujejo uporabo ML za modeliranje posameznih parametrov distribucije znotraj GAMLSS okvira, kar bi omogočilo še večjo fleksibilnost pri zajemanju kompleksnih, nelinearnih odnosov, ki so jih morda tradicionalne gladke funkcije spregledale.

Vendar pa je pri integraciji ML v GAMLSS ključnega pomena ohranitev interpretibilnosti in robustnosti modela. Medtem ko so ML modeli pogosto izjemno natančni v napovedovanju, so lahko manj transparentni ('black-box' modeli), kar je v reguliranem okolju zavarovalništva problematično. GAMLSS modeli nam omogočajo, da ohranimo statistično rigoroznost in interpretacijo parametrov, hkrati pa izkoristimo nekatere prednosti ML. Evalvacija teh hibridnih modelov bo zahtevala še bolj sofisticirane kombinacije informacijskih kriterijev (prilagojenih za kompleksnost ML komponent) in naprednih diagnostičnih testov, da se zagotovi njihova ustreznost in zanesljivost. Kot vedno, bo tudi v prihodnosti, cilj doseči ravnovesje med naprednostjo, natančnostjo in transparentnostjo.

Zaključek: Pomen celovite evalvacije modelov

V tem obsežnem pregledu sem poskušala poudariti ključni pomen celovite in sistematične evalvacije GAMLSS modelov pri modeliranju škodnih dogodkov, še posebej v kontekstu nezgodnih zavarovanj. Uporaba informacijskih kriterijev, kot sta AIC in BIC, nam omogoča kvantificirano primerjavo različnih modelov in nas usmerja k izbiri modela, ki optimalno uravnoteži prileganje podatkom in kompleksnost. Vendar, kot sem večkrat poudarila, ti kriteriji sami po sebi niso dovolj.

Kritična vloga diagnostičnih testov residuov, vključno s Q-Q grafi, worm plots in statističnimi testi, je nepogrešljiva za preverjanje predpostavk modela, identifikacijo morebitnih neskladij in zagotavljanje robustnosti. Le kombinacija numeričnih kriterijev in podrobne vizualne ter statistične diagnostike mi omogoča, da z gotovostjo izberem in priporočim model, ki je tako statistično ustrezen kot tudi praktično uporaben in zanesljiv za sprejemanje poslovnih odločitev.

V moji 20-letni karieri sem spoznala, da je aktuarska znanost nenehno učenje in prilagajanje novim izzivom in metodologijam. GAMLSS modeli predstavljajo pomemben korak naprej pri modeliranju kompleksnih distribucij škod, vendar je njihova dejanska vrednost odvisna od rigorozne evalvacije. Zato bom še naprej sledila najnovejšim trendom in metodam, da svojim strankam zagotovim najboljše in najzanesljivejše zavarovalne rešitve. Če imate vprašanja ali potrebujete poglobljeno analizo vašega zavarovalnega portfelja, me prosim kontaktirajte – z veseljem vam bom pomagala.

Primer iz prakse

Modeliranje trajne invalidnosti: Primer zavarovanca Marka

Brez ustreznega zavarovanja
Marko, 45-letni zavarovanec, utrpi hudo nezgodo na delovnem mestu, ki mu pusti trajne posledice. Zavarovalnica, ki uporablja preveč poenostavljen GLM model za oceno invalidnosti, ki ne upošteva vseh demografskih dejavnikov in specifik nezgode, podceni verjetnost določene stopnje invalidnosti in s tem potrebnih rezervacij. Ob izplačilu ugotovijo, da so pričakovani stroški in obveznosti bistveno višji od modeliranih, kar povzroči nepredvidene obremenitve za zavarovalnico in potencialno nezadostno kritje za Marka. Model je zaradi neustrezne distribucije in predpostavk precenil število manjših škod in podcenil ekstreme, kar je bila finančna obremenitev za zavarovalnico in nezadovoljstvo pri stranki.
Z ustreznim zavarovanjem
Zavarovalnica je, na podlagi podrobne GAMLSS analize, ki je upoštevala starost, spol, vrsto poklica in mehanizem nezgode, zmožna natančneje modelirati distribucijo stopenj trajne invalidnosti. Uporabljen je bil GAMLSS model z Beta distribucijo (za delež invalidnosti) in log-normalno distribucijo (za znesek odškodnine), kjer so bili parametri μ in σ odvisni od prediktorjev. Evalvacija z AIC, BIC in podrobna diagnostika residuov (worm plots) je potrdila robustnost in ustreznost modela. Zavarovalnica je tako lahko natančno določila pričakovane obveznosti in rezervacije za primer Markove invalidnosti, kar je zagotovilo pravočasno in ustrezno izplačilo, hkrati pa optimalno upravljanje s tveganji za zavarovalnico.

Primer je ilustrativen in povzet po tipičnih situacijah iz prakse. Kritja, izključitve in postopki se med zavarovalnicami razlikujejo.

Pogosta vprašanja

Zakaj so GAMLSS modeli boljši od GLM za škodne dogodke?
GAMLSS modeli so bolj fleksibilni, saj omogočajo modeliranje vseh parametrov distribucije (lokacije, skale, oblike) kot funkcij prediktorjev, ne le povprečja. To je ključno za asimetrične in kompleksne distribucije, značilne za škodne dogodke, kjer GLM pogosto ne zadoščajo.
Kaj pomenita AIC in BIC in kateri je boljši?
AIC (Akaike Information Criterion) in BIC (Bayesian Information Criterion) sta informacijska kriterija za primerjavo statističnih modelov. Nižja vrednost pomeni boljši model. AIC nagrajuje kompleksnost, BIC pa bolj favorizira enostavnejše, robustnejše modele, še posebej pri velikih vzorcih. Izbira je odvisna od ciljev analize.
Zakaj so diagnostični testi residuov pomembni?
Diagnostični testi residuov (npr. Q-Q grafi, worm plots, Shapiro-Wilk test) so ključni za preverjanje predpostavk modela in identifikacijo morebitnih neskladij med modelom in podatki. Zagotavljajo, da je model robusten, zanesljiv in da so njegovi zaključki verodostojni. Brez njih, se lahko zmotno zanašamo na neustrezen model.
Ali lahko GAMLSS modele uporabimo tudi za druge vrste zavarovanj?
Seveda! GAMLSS modeli so izjemno vsestranski in se lahko uporabljajo tudi za modeliranje škod v premoženjskih zavarovanjih, življenjskih zavarovanjih (npr. stopnje umrljivosti), zdravstvenih zavarovanjih in drugih področjih, kjer so podatki kompleksni in ne ustrezajo preprostim distribucijam. Njihova fleksibilnost je velika prednost.
Kako zakonodaja vpliva na izbiro in evalvacijo modelov?
Zakonodaja (npr. ZZavar-1, Solvency II) ne predpisuje specifičnih modelov, ampak zahteva zanesljive in natančne aktuarske izračune za določanje rezervacij in kapitalne ustreznosti. To pomeni, da morajo biti izbrani modeli ustrezno evalvirani, dokumentirani in utemeljeni, kar poudarja pomen AIC, BIC in diagnostike residuov.

Viri in reference

  • Uradni list RS – Zakon o zavarovalništvu (ZZavar-1)
  • Uradni list RS – Zakon o zavarovalnih pogodbah (ZZavPog)
  • Agencija za zavarovalni nadzor (AZN) – smernice in poročila
  • Rigby, R. A., & Stasinopoulos, D. M. (2005). Generalized additive models for location, scale and shape. Applied Statistics, 54(3), 507-554.
  • Akaike, H. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control, 19(6), 716-723.
  • Schwarz, G. (1978). Estimating the dimension of a model. The Annals of Statistics, 6(2), 461-464.

Nadaljujte branje o tej temi

Povezave so izbrane samodejno glede na steber zaščite in ključne besede te objave.

Priporočeno branje

Petkin letni pregled polic

Rezervirajte vaš redni letni 15-minutni pregledni pogovor

Če vas je ta tema pritegnila, jo pri pregledu obravnavava konkretno na vaših policah — v okviru področja Nezgoda. Pogovor je informativne narave, brez ponudbe in brez obveznosti.

Vsebina objave je splošna informacija in ne osebno svetovanje; veljajo pogoji posamezne police.