AktuNotes
← Kembali
TA2 · Materi 1.4

Tail Characteristics

Hard Bobot: 5–10% Klugman, Panjer & Willmot (2019) Loss Models 5th ed., Bab 3 (kecuali Section 3.5), 4, 5
TA2TopicModelBesarKlaimTailCharacteristicsHazardRateMeanExcessLossTeoriRisiko

📊 1.4 — Tail Characteristics

Ringkasan Cepat

Topik: Tail Characteristics | Bobot: ~5–10% | Difficulty: Hard Ref: Klugman et al. (2019) Loss Models 5th ed., Bab 3–5 | Prereq: 1.1 Moment and Probability Generating Functions, 1.2 Distribution Classes and Extreme Value

Section 0 — Pemetaan Topik

Topik TA2Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Model Besar Klaim1.4Menghitung dan menginterpretasikan momen, rasio momen, limiting tail behaviour, hazard rate, dan mean excess loss function; membandingkan dua distribusi berdasarkan berat ekor5–10%Hard1.1 Moment and Probability Generating Functions, 1.2 Distribution Classes and Extreme Value, 1.3 Techniques for Creating New Distributions3.1 Coverage Modifications on Severity and Frequency, 3.2 Loss Elimination Ratio and Inflation, 5.2 VaR and TVaRKPW (2019) Bab 3–5

Section 1 — Intuisi

Bayangkan sebuah perusahaan reasuransi sedang mengevaluasi dua portofolio asuransi properti di Indonesia — satu dari wilayah Jakarta Pusat dengan profil klaim yang relatif stabil, dan satu lagi dari kawasan industri di pesisir Cilegon yang rentan terhadap klaim besar akibat bencana. Secara rata-rata, kedua portofolio mungkin memiliki ekspektasi klaim yang serupa. Namun bagi reasuradur, yang jauh lebih penting adalah apa yang terjadi di ekor distribusi — seberapa sering dan seberapa besar klaim-klaim ekstrem yang bisa terjadi? Inilah yang dijawab oleh karakteristik ekor (tail characteristics).

Konsep pertama adalah hazard rate (atau failure rate). Fungsi ini menjawab pertanyaan: “Jika sebuah klaim sudah mencapai nilai xx, seberapa cepat peluang terjadinya klaim yang lebih besar dari xx menurun?” Distribusi dengan hazard rate yang menurun artinya makin besar klaim yang sudah terjadi, makin besar pula kemungkinan klaim akan terus bertumbuh — inilah tanda distribusi heavy-tail yang berbahaya. Sebaliknya, hazard rate yang meningkat mengindikasikan distribusi light-tail yang “aman”. Distribusi Eksponensial adalah kasus unik dengan hazard rate konstan, mencerminkan sifat memoryless-nya.

Konsep kedua yang sangat praktis di aktuaria adalah mean excess loss function (MEL), yang menjawab: “Jika sebuah klaim sudah diketahui melebihi deductible sebesar dd, berapa ekspektasi kelebihan klaim di atas dd itu?” Fungsi ini secara langsung terhubung dengan penetapan premi reasuransi excess-of-loss. Jika MEL suatu distribusi naik seiring dd, distribusi tersebut heavy-tailed — semakin besar deductible yang diterapkan, semakin besar pula ekspektasi kelebihan yang harus ditanggung reasuradur. Memahami karakteristik ekor bukan hanya soal akademis; ia adalah kompas yang digunakan aktuaris untuk menentukan berapa premi layak dikenakan untuk perlindungan terhadap risiko katastrofi.

Section 2 — Definisi Formal

Definisi Matematis — Lima Karakteristik Ekor

Misalkan XX adalah variabel acak kontinu non-negatif dengan PDF f(x)f(x), CDF F(x)F(x), dan survival function S(x)=1F(x)S(x) = 1 - F(x).

Hazard Rate (Force of Mortality / Failure Rate):

h(x)=f(x)S(x)=f(x)1F(x),x>0h(x) = \frac{f(x)}{S(x)} = \frac{f(x)}{1 - F(x)}, \quad x > 0

Mean Excess Loss Function (MEL):

e(d)=E[XdX>d]=dS(x)dxS(d),d0e(d) = E[X - d \mid X > d] = \frac{\int_d^\infty S(x)\, dx}{S(d)}, \quad d \geq 0

Limited Expected Value (LEV):

E[Xu]=0uS(x)dx=0u[1F(x)]dxE[X \wedge u] = \int_0^u S(x)\, dx = \int_0^u [1 - F(x)]\, dx
SimbolMaknaCatatan
S(x)S(x)Survival function: P(X>x)P(X > x)S(x)=1F(x)S(x) = 1 - F(x); selalu monotone turun
h(x)h(x)Hazard rate / force of mortalityh(x)0h(x) \geq 0; mengukur laju “kepunahan” distribusi
H(x)H(x)Cumulative hazard functionH(x)=0xh(t)dt=lnS(x)H(x) = \int_0^x h(t)\, dt = -\ln S(x)
e(d)e(d)Mean excess loss function di ddEkspektasi kelebihan klaim di atas threshold dd
e(0)e(0)Mean excess loss di d=0d=0e(0)=E[X]e(0) = E[X] — mean unconditional
E[Xu]E[X \wedge u]Limited Expected Value (LEV)Ekspektasi XX yang di-cap pada nilai uu
μk\mu_k'Raw moment ke-kkE[Xk]=0kxk1S(x)dxE[X^k] = \int_0^\infty k x^{k-1} S(x)\, dx
CVCVCoefficient of VariationCV=σ/μCV = \sigma / \mu; ukuran dispersi relatif
γ1\gamma_1Skewness (momen ke-3 yang distandarisasi)γ1=μ3/σ3\gamma_1 = \mu_3' / \sigma^3; positif = right-skewed

Rumus Utama

Hubungan fundamental: hazard rate ↔ survival function:

S(x)=exp ⁣(0xh(t)dt)=eH(x)S(x) = \exp\!\left(-\int_0^x h(t)\, dt\right) = e^{-H(x)}

Label: Jika diketahui h(x)h(x), distribusi sepenuhnya ditentukan — ini adalah representasi alternatif yang ekuivalen.

PDF dari hazard rate:

f(x)=h(x)S(x)=h(x)eH(x)f(x) = h(x) \cdot S(x) = h(x)\, e^{-H(x)}

Label: Berguna untuk mengidentifikasi distribusi dari bentuk hazard rate yang diberikan.

Momen via survival function (alternatif integrasi langsung):

E[Xk]=0kxk1S(x)dxE[X^k] = \int_0^\infty k x^{k-1} S(x)\, dx

Label: Seringkali lebih mudah dihitung daripada 0xkf(x)dx\int_0^\infty x^k f(x)\, dx, khususnya untuk distribusi heavy-tail.

Mean excess loss function — bentuk integral:

e(d)=dS(x)dxS(d)=E[X]E[Xd]S(d)e(d) = \frac{\int_d^\infty S(x)\, dx}{S(d)} = \frac{E[X] - E[X \wedge d]}{S(d)}

Label: Hubungan MEL dengan LEV sangat penting — soal sering memberikan E[Xd]E[X \wedge d] dan meminta e(d)e(d).

Hubungan MEL dengan mean unconditional:

e(0)=E[X]=0S(x)dxe(0) = E[X] = \int_0^\infty S(x)\, dx

Label: Kasus khusus d=0d=0: MEL di threshold nol adalah mean distribusi.

Rasio momen (moment ratio) ke-kk:

E[Xk](E[X])k\frac{E[X^k]}{(E[X])^k}

Label: Digunakan untuk membandingkan berat ekor dua distribusi dengan mean yang berbeda; rasio lebih besar = ekor lebih berat.

Hazard rate distribusi-distribusi penting:

DistribusiHazard Rate h(x)h(x)Sifat Ekor
Eksponensial(θ)(\theta)1/θ1/\theta (konstan)Memoryless; benchmark light-tail
Weibull(τ,θ)(\tau, \theta)(τ/θ)(x/θ)τ1(\tau/\theta)(x/\theta)^{\tau-1}τ>1\tau > 1: naik (light-tail); τ<1\tau < 1: turun (heavy-tail)
Gamma(α,θ)(\alpha, \theta)Meningkat ke 1/θ1/\thetaLight-tail untuk α>1\alpha > 1; mendekati Eksponensial
Pareto(α,θ)(\alpha, \theta)α/(x+θ)\alpha/(x + \theta)Selalu menurun — heavy-tail kanonik
Lognormal(μ,σ2)(\mu, \sigma^2)Berbentuk unimodal, lalu turunHeavy-tail (ekor lebih berat dari Eksponensial)

MEL distribusi-distribusi penting:

DistribusiMEL e(d)e(d)Sifat MEL
Eksponensial(θ)(\theta)θ\theta (konstan)Memoryless: MEL tidak bergantung dd
Pareto(α,θ)(\alpha, \theta)(d+θ)/(α1)(d + \theta)/(\alpha - 1)Naik linear — bukti heavy-tail
Gamma(α,θ)(\alpha, \theta)Menurun ke θ\thetaLight-tail: semakin besar dd, MEL mendekati θ\theta

Asumsi Eksplisit

  1. XX adalah variabel acak kontinu non-negatif dengan F(0)=0F(0) = 0.
  2. S(d)>0S(d) > 0 untuk semua dd dalam domain yang relevan (agar MEL terdefinisi).
  3. MEL e(d)e(d) terdefinisi hanya jika E[X]<E[X] < \infty — distribusi yang tidak memiliki momen pertama tidak memiliki MEL yang bermakna.
  4. Untuk limiting tail behaviour, perbandingan dilakukan pada xx \to \infty — rasio S1(x)/S2(x)S_1(x)/S_2(x) atau perbandingan h1(x)h_1(x) vs h2(x)h_2(x).
  5. Distribusi dikatakan heavy-tailed relatif satu sama lain — “berat” dan “ringan” adalah konsep komparatif, bukan absolut.

Section 3 — Jembatan Logika

Dari Definisi ke Rumus — Tiga Cara Membaca Ekor

Ada tiga “lensa” yang saling melengkapi untuk mengkarakterisasi ekor distribusi:

Lensa 1 — Hazard rate h(x)h(x): Melihat laju “peluruhan” distribusi secara lokal di titik xx. Analoginya: seberapa cepat intensitas hujan melemah ketika curah hujan sudah mencapai tingkat xx? Hazard rate menurun berarti distribusi “semakin lambat melemah” — tanda heavy-tail.

Lensa 2 — MEL e(d)e(d): Melihat ekspektasi kelebihan, kondisional pada sudah melampaui dd. Analoginya: “klaim sudah melebihi Rp500 juta — berapa lagi yang harus kita cadangkan?” MEL yang naik berarti semakin besar threshold, semakin besar ekspektasi tambahan — juga tanda heavy-tail.

Lensa 3 — Limiting tail behaviour: Membandingkan survival function dua distribusi saat xx \to \infty. Jika S1(x)/S2(x)S_1(x)/S_2(x) \to \infty, distribusi 1 jauh lebih heavy-tailed. Eksponensial memiliki ekor seperti ex/θe^{-x/\theta}, sedangkan Pareto seperti xαx^{-\alpha} — untuk xx besar, xαx^{-\alpha} jauh lebih besar dari ex/θe^{-x/\theta}, sehingga Pareto selalu heavy-tailed relatif terhadap Eksponensial.

Support dan Domain
  • Hazard rate h(x)h(x) terdefinisi hanya di mana S(x)>0S(x) > 0. Untuk distribusi dengan support (0,ω)(0, \omega) terbatas (seperti Uniform), h(x)h(x) \to \infty saat xωx \to \omega^-.
  • MEL e(d)e(d) terdefinisi hanya jika S(d)>0S(d) > 0 dan E[X]<E[X] < \infty. Distribusi Pareto dengan α1\alpha \leq 1 tidak memiliki E[X]E[X] yang finite — MEL tidak terdefinisi!
  • Untuk limiting tail comparison: selalu bandingkan survival function (bukan PDF) saat xx \to \infty — perbandingan PDF bisa menyesatkan karena PDF bisa kecil di ekor tapi integral-nya (yaitu probabilitas ekor) masih besar.

Derivasi: Rumus MEL via Integrasi by Parts

Ini adalah derivasi paling sering diuji — menunjukkan mengapa e(d)=dS(x)dxS(d)e(d) = \frac{\int_d^\infty S(x)\, dx}{S(d)}.

Langkah 1 — Definisi ekspektasi bersyarat:

e(d)=E[XdX>d]=E[(Xd)1X>d]P(X>d)=d(xd)f(x)dxS(d)e(d) = E[X - d \mid X > d] = \frac{E[(X-d)\,\mathbf{1}_{X>d}]}{P(X > d)} = \frac{\int_d^\infty (x-d)\, f(x)\, dx}{S(d)}

Langkah 2 — Fokus pada pembilang; substitusi u=xdu = x - d, du=dxdu = dx:

d(xd)f(x)dx=0uf(u+d)du\int_d^\infty (x-d)\, f(x)\, dx = \int_0^\infty u\, f(u+d)\, du

Langkah 3 — Integrasi by parts pada pembilang asli: uparts=xdu_{\text{parts}} = x - d, dv=f(x)dxdv = f(x)\,dx:

d(xd)f(x)dx=[(xd)(S(x))]d+dS(x)dx\int_d^\infty (x-d)\, f(x)\, dx = \Big[(x-d)(-S(x))\Big]_d^\infty + \int_d^\infty S(x)\, dx

Langkah 4 — Evaluasi batas: jika momen pertama finite, maka (xd)S(x)0(x-d)S(x) \to 0 saat xx \to \infty (distribusi tidak memiliki massa di \infty), dan di x=dx = d: (dd)S(d)=0(d-d)S(d) = 0. Sehingga suku pertama gugur:

d(xd)f(x)dx=dS(x)dx\int_d^\infty (x-d)\, f(x)\, dx = \int_d^\infty S(x)\, dx

Langkah 5 — Bagi dengan S(d)S(d):

e(d)=dS(x)dxS(d)e(d) = \frac{\int_d^\infty S(x)\, dx}{S(d)}

∎ Inilah bentuk yang paling berguna secara komputasi.

Derivasi: MEL Pareto dan Sifat Heavy-Tail-nya

Misalkan XPareto(α,θ)X \sim \text{Pareto}(\alpha, \theta) dengan S(x)=(θθ+x)αS(x) = \left(\frac{\theta}{\theta + x}\right)^\alpha.

Langkah 1 — Hitung dS(x)dx\int_d^\infty S(x)\, dx:

d(θθ+x)αdx\int_d^\infty \left(\frac{\theta}{\theta+x}\right)^\alpha dx

Substitusi t=θ+xt = \theta + x, dt=dxdt = dx, batas t:θ+dt: \theta+d \to \infty:

=θαθ+dtαdt=θαtα+1α+1θ+d=θα(θ+d)1αα1(α>1)= \theta^\alpha \int_{\theta+d}^\infty t^{-\alpha}\, dt = \theta^\alpha \cdot \frac{t^{-\alpha+1}}{-\alpha+1}\Bigg|_{\theta+d}^\infty = \frac{\theta^\alpha (\theta+d)^{1-\alpha}}{\alpha - 1} \quad (\alpha > 1)

Langkah 2 — Hitung S(d)S(d):

S(d)=(θθ+d)αS(d) = \left(\frac{\theta}{\theta+d}\right)^\alpha

Langkah 3 — Bagi:

e(d)=θα(θ+d)1α/(α1)θα(θ+d)α=(θ+d)1α(α1)(θ+d)α=θ+dα1e(d) = \frac{\theta^\alpha (\theta+d)^{1-\alpha}/(\alpha-1)}{\theta^\alpha (\theta+d)^{-\alpha}} = \frac{(\theta+d)^{1-\alpha}}{(\alpha-1)(\theta+d)^{-\alpha}} = \frac{\theta + d}{\alpha - 1}

Kesimpulan: e(d)=d+θα1e(d) = \frac{d + \theta}{\alpha - 1} adalah fungsi linear naik dalam dd — ini adalah bukti formal bahwa Pareto heavy-tailed. Semakin besar deductible, semakin besar ekspektasi kelebihan yang harus dibayar reasuradur.

Dilarang
  1. Jangan menghitung e(d)e(d) dengan formula Pareto untuk α1\alpha \leq 1 — MEL tidak terdefinisi karena E[X]=E[X] = \infty untuk α1\alpha \leq 1. Ini bukan hanya kesalahan teknis; secara aktuaria artinya distribusi tersebut tidak dapat di-insurance-kan dengan cara konvensional.
  2. Jangan membandingkan berat ekor hanya dari nilai PDF — PDF distribusi heavy-tail bisa sangat kecil di titik tertentu namun memiliki probabilitas ekor yang jauh lebih besar daripada distribusi light-tail. Selalu gunakan survival function S(x)S(x), hazard rate h(x)h(x), atau MEL e(d)e(d) untuk perbandingan ekor.
  3. Jangan gunakan formula e(d)=E[X]E[Xd]S(d)e(d) = \frac{E[X] - E[X \wedge d]}{S(d)} tanpa menghitung E[Xd]E[X \wedge d] terlebih dahulu — LEV E[Xd]E[X \wedge d] memerlukan perhitungan tersendiri; jangan substitusi E[Xd]=dF(d)E[X \wedge d] = d \cdot F(d) (ini salah kecuali XX deterministik).

Section 4 — Contoh Soal

Soal A — Fundamental

Variabel acak XEksponensial(θ=2000)X \sim \text{Eksponensial}(\theta = 2000) merepresentasikan besar klaim asuransi kendaraan (dalam ribuan rupiah). Hitung: (a) hazard rate h(x)h(x), (b) MEL e(d)e(d) untuk d=3000d = 3000, dan (c) interpretasikan hasil secara aktuaria.

Solusi Soal A

Pendekatan: Gunakan formula langsung hazard rate dan MEL untuk distribusi Eksponensial. Sifat memoryless akan terkonfirmasi.

1. Identifikasi Variabel

  • XEksponensial(θ=2000)X \sim \text{Eksponensial}(\theta = 2000)
  • f(x)=12000ex/2000f(x) = \frac{1}{2000} e^{-x/2000}
  • S(x)=ex/2000S(x) = e^{-x/2000}
  • d=3000d = 3000

2. Identifikasi Distribusi / Model Eksponensial adalah distribusi memoryless benchmark — hazard rate konstan dan MEL konstan. Distribusi ini sering digunakan sebagai distribusi klaim dasar sebelum dilakukan pemilihan model yang lebih kompleks.

3. Setup Persamaan

h(x)=f(x)S(x),e(d)=dS(x)dxS(d)h(x) = \frac{f(x)}{S(x)}, \qquad e(d) = \frac{\int_d^\infty S(x)\,dx}{S(d)}

4. Eksekusi Aljabar

(a) Hazard rate:

h(x)=12000ex/2000ex/2000=12000h(x) = \frac{\frac{1}{2000}e^{-x/2000}}{e^{-x/2000}} = \frac{1}{2000}

Konstan untuk semua x>0x > 0 — tidak bergantung pada xx.

(b) MEL dengan formula integral:

3000ex/2000dx=[2000ex/2000]3000=2000e3000/2000=2000e1.5\int_{3000}^\infty e^{-x/2000}\, dx = \Big[-2000\, e^{-x/2000}\Big]_{3000}^\infty = 2000\, e^{-3000/2000} = 2000\, e^{-1.5} S(3000)=e3000/2000=e1.5S(3000) = e^{-3000/2000} = e^{-1.5} e(3000)=2000e1.5e1.5=2000e(3000) = \frac{2000\, e^{-1.5}}{e^{-1.5}} = 2000

5. Verification e(d)=θ=2000e(d) = \theta = 2000 untuk semua dd — ini membuktikan sifat memoryless: apapun threshold deductible yang diterapkan, ekspektasi kelebihan klaim selalu sama dengan mean distribusi. Juga: e(0)=E[X]=2000e(0) = E[X] = 2000 ✓.

Hasil: h(x)=1/2000h(x) = 1/2000 (konstan), e(3000)=2000e(3000) = 2000. Interpretasi aktuaria: meskipun deductible Rp3 juta diterapkan, reasuradur masih menghadapi ekspektasi klaim tambahan Rp2 juta — sama persis seperti tanpa deductible. Ini mencerminkan sifat memoryless Eksponensial.

Exam Tips — Soal A

Target waktu: 2–3 menit. Common trap: Beberapa kandidat mencoba menghitung E[XdX>d]E[X - d \mid X > d] langsung dari definisi bersyarat — ini benar tapi lebih panjang. Gunakan formula e(d)=dS(x)dx/S(d)e(d) = \int_d^\infty S(x)\,dx / S(d) yang lebih efisien. Shortcut: Untuk Eksponensial, langsung tulis h(x)=1/θh(x) = 1/\theta dan e(d)=θe(d) = \theta tanpa perhitungan — ini adalah fakta yang wajib dihafal.


Soal B — Exam-Typical

Distribusi klaim XPareto(α=3,θ=6000)X \sim \text{Pareto}(\alpha = 3, \theta = 6000). Sebuah polis memiliki deductible biasa d=2000d = 2000. Hitung: (a) MEL e(2000)e(2000), (b) E[X2000]E[X \wedge 2000] (LEV), dan (c) verifikasi hubungan e(d)=E[X]E[Xd]S(d)e(d) = \frac{E[X] - E[X \wedge d]}{S(d)}.

Solusi Soal B

Pendekatan: Gunakan formula MEL Pareto untuk (a), hitung LEV dari definisi untuk (b), lalu verifikasi konsistensi keduanya.

1. Identifikasi Variabel

  • XPareto(α=3,θ=6000)X \sim \text{Pareto}(\alpha = 3, \theta = 6000); α>1\alpha > 1 sehingga E[X]E[X] finite ✓
  • S(x)=(60006000+x)3S(x) = \left(\frac{6000}{6000+x}\right)^3
  • E[X]=θα1=60002=3000E[X] = \frac{\theta}{\alpha - 1} = \frac{6000}{2} = 3000
  • d=2000d = 2000

2. Identifikasi Distribusi / Model Pareto adalah distribusi heavy-tail kanonik. MEL-nya naik linear — makin besar deductible, makin besar ekspektasi kelebihan. Ini berbanding terbalik dengan Eksponensial (MEL konstan).

3. Setup Persamaan

e(d)=d+θα1,E[Xu]=θα1[1(θθ+u)α1]e(d) = \frac{d + \theta}{\alpha - 1}, \qquad E[X \wedge u] = \frac{\theta}{\alpha-1}\left[1 - \left(\frac{\theta}{\theta+u}\right)^{\alpha-1}\right] e(d)=E[X]E[Xd]S(d)e(d) = \frac{E[X] - E[X \wedge d]}{S(d)}

4. Eksekusi Aljabar

(a) MEL Pareto:

e(2000)=2000+600031=80002=4000e(2000) = \frac{2000 + 6000}{3 - 1} = \frac{8000}{2} = 4000

(b) LEV Pareto pada u=2000u = 2000:

E[X2000]=60002[1(60008000)2]=3000[1(34)2]=3000[1916]=3000×716=1312.5E[X \wedge 2000] = \frac{6000}{2}\left[1 - \left(\frac{6000}{8000}\right)^{2}\right] = 3000\left[1 - \left(\frac{3}{4}\right)^2\right] = 3000\left[1 - \frac{9}{16}\right] = 3000 \times \frac{7}{16} = 1312.5

(c) Verifikasi hubungan:

S(2000)=(60008000)3=(34)3=2764S(2000) = \left(\frac{6000}{8000}\right)^3 = \left(\frac{3}{4}\right)^3 = \frac{27}{64} E[X]E[X2000]S(2000)=30001312.527/64=1687.527/64=1687.5×6427=10800027=4000\frac{E[X] - E[X \wedge 2000]}{S(2000)} = \frac{3000 - 1312.5}{27/64} = \frac{1687.5}{\,27/64\,} = 1687.5 \times \frac{64}{27} = \frac{108000}{27} = 4000 \checkmark

5. Verification e(2000)=4000>E[X]=3000e(2000) = 4000 > E[X] = 3000: MEL melebihi mean unconditional — ini adalah bukti heavy-tail. Untuk Eksponensial, e(d)=E[X]e(d) = E[X] selalu; untuk Pareto, e(d)>E[X]e(d) > E[X] segera saat d>0d > 0.

Hasil: e(2000)=4000e(2000) = 4000, E[X2000]=1312.5E[X \wedge 2000] = 1312.5, keduanya konsisten via hubungan LEV–MEL.

Exam Tips — Soal B

Target waktu: 4–5 menit. Common trap: Formula LEV Pareto sering disalahtulis — perhatikan bahwa eksponen di LEV adalah α1\alpha - 1, bukan α\alpha. Shortcut: Hafal formula LEV Pareto: E[Xu]=θα1 ⁣[1(θθ+u)α1]E[X \wedge u] = \frac{\theta}{\alpha-1}\!\left[1 - \left(\frac{\theta}{\theta+u}\right)^{\alpha-1}\right]. Ini muncul sangat sering di soal modifikasi coverage (Topik 3).


Soal C — Challenging

Dua model sedang dibandingkan untuk distribusi klaim asuransi kebakaran: Model A adalah XAPareto(α=4,θ=9000)X_A \sim \text{Pareto}(\alpha = 4, \theta = 9000) dan Model B adalah XBEksponensial(θ=3000)X_B \sim \text{Eksponensial}(\theta = 3000). Keduanya memiliki mean yang sama. Bandingkan: (a) hazard rate kedua model untuk xx besar, (b) MEL e(d)e(d) untuk d=0,3000,9000d = 0, 3000, 9000, dan (c) tentukan model mana yang lebih konservatif untuk penetapan premi reasuransi excess-of-loss dengan attachment point d=9000d = 9000.

Solusi Soal C

Pendekatan: Bandingkan hazard rate dan MEL secara sistematis; gunakan perilaku ekor untuk menarik kesimpulan aktuaria.

1. Identifikasi Variabel

  • Model A: XAPareto(4,9000)X_A \sim \text{Pareto}(4, 9000); E[XA]=9000/(41)=3000E[X_A] = 9000/(4-1) = 3000
  • Model B: XBEksponensial(3000)X_B \sim \text{Eksponensial}(3000); E[XB]=3000E[X_B] = 3000
  • Kedua model memiliki mean =3000= 3000 ✓ — perbandingan fair
  • Attachment point reasuransi: d=9000d = 9000

2. Identifikasi Distribusi / Model Pareto: heavy-tail, hazard rate menurun. Eksponensial: light-tail (memoryless), hazard rate konstan. Dengan mean sama, perbedaan sepenuhnya terletak pada perilaku ekor. Ini adalah soal perbandingan ekor klasik.

3. Setup Persamaan

Hazard rate:

hA(x)=αθ+x=49000+x,hB(x)=13000h_A(x) = \frac{\alpha}{\theta + x} = \frac{4}{9000 + x}, \qquad h_B(x) = \frac{1}{3000}

MEL:

eA(d)=d+90003,eB(d)=3000e_A(d) = \frac{d + 9000}{3}, \qquad e_B(d) = 3000

4. Eksekusi Aljabar

(a) Perbandingan hazard rate:

Titik persilangan: hA(x)=hB(x)h_A(x) = h_B(x) saat 49000+x=13000\frac{4}{9000+x} = \frac{1}{3000}, yaitu x=3000x = 3000.

  • Untuk x<3000x < 3000: hA(x)>hB(x)h_A(x) > h_B(x) — Pareto memiliki hazard rate lebih tinggi (lebih “aman”) di ekor kiri
  • Untuk x>3000x > 3000: hA(x)<hB(x)h_A(x) < h_B(x) — Pareto memiliki hazard rate lebih rendah (ekor kanan lebih berat)

Untuk xx besar: hA(x)4x0h_A(x) \approx \frac{4}{x} \to 0, sedangkan hB(x)=1/3000h_B(x) = 1/3000 tetap konstan. Pareto semakin lambat “punah” — bukti heavy-tail.

(b) Tabel MEL:

ddeA(d)e_A(d)eB(d)e_B(d)Selisih eAeBe_A - e_B
009000/3=30009000/3 = 30003000300000
3000300012000/3=400012000/3 = 400030003000+1000+1000
9000900018000/3=600018000/3 = 600030003000+3000+3000

(c) Premi reasuransi excess-of-loss di d=9000d = 9000:

Premi reasuransi proporsional dengan e(d)×S(d)e(d) \times S(d) (ekspektasi klaim yang ditanggung reasuradur):

SA(9000)=(900018000)4=(12)4=116=0.0625S_A(9000) = \left(\frac{9000}{18000}\right)^4 = \left(\frac{1}{2}\right)^4 = \frac{1}{16} = 0.0625 SB(9000)=e9000/3000=e30.0498S_B(9000) = e^{-9000/3000} = e^{-3} \approx 0.0498

Ekspektasi klaim reasuransi =e(d)×S(d)= e(d) \times S(d):

Model A: 6000×0.0625=375\text{Model A: } 6000 \times 0.0625 = 375 Model B: 3000×0.0498149.4\text{Model B: } 3000 \times 0.0498 \approx 149.4

5. Verification Model A menghasilkan ekspektasi klaim reasuransi 2.5×\approx 2.5\times lebih besar dari Model B meskipun mean sama. Ini mengkonfirmasi bahwa Pareto jauh lebih berbahaya dari perspektif reasuransi. Fakta bahwa eA(d)e_A(d) naik sementara eB(d)e_B(d) konstan adalah tanda definitif perbedaan berat ekor.

Hasil: Model A (Pareto) secara signifikan lebih heavy-tailed — hazard rate menurun ke nol sementara Eksponensial konstan. MEL Pareto di d=9000d = 9000 adalah 6000 vs 3000 untuk Eksponensial. Model A menghasilkan ekspektasi klaim reasuransi 2.5× lebih besar — Model A lebih konservatif untuk penetapan premi reasuransi.

Exam Tips — Soal C

Target waktu: 5–6 menit. Common trap: Menyimpulkan bahwa model dengan mean sama pasti “setara” untuk tujuan reasuransi — ini salah! Distribusi dengan mean sama bisa sangat berbeda perilaku ekornya. Shortcut: Untuk membandingkan dua distribusi: (1) cek apakah MEL naik atau turun/konstan, dan (2) bandingkan hazard rate di xx besar. Kedua tes ini sudah cukup untuk menentukan distribusi mana yang heavy-tailed.

Section 5 — Verifikasi & Sanity Check

Cross-check MEL dengan Mean

Selalu verifikasi bahwa e(0)=E[X]e(0) = E[X]:

e(0)=0S(x)dxS(0)=E[X]1=E[X]e(0) = \frac{\int_0^\infty S(x)\, dx}{S(0)} = \frac{E[X]}{1} = E[X]

Jika e(0)E[X]e(0) \neq E[X] dari formula yang digunakan, ada kesalahan dalam derivasi MEL. Ini adalah sanity check paling cepat dan paling wajib dilakukan.

Cross-check Hazard Rate via Survival Function

Jika hazard rate h(x)h(x) telah dihitung, verifikasi via:

S(x)=exp ⁣(0xh(t)dt)S(x) = \exp\!\left(-\int_0^x h(t)\, dt\right)

Misal: jika h(x)=ch(x) = c (konstan), maka S(x)=ecxS(x) = e^{-cx} — ini adalah Eksponensial dengan mean 1/c1/c. Jika h(x)=τxτ1/θτh(x) = \tau x^{\tau-1}/\theta^\tau, maka S(x)=e(x/θ)τS(x) = e^{-(x/\theta)^\tau} — ini adalah Weibull. Jika hasil rekonstruksi S(x)S(x) tidak cocok dengan S(x)S(x) yang diberikan, ada kesalahan dalam h(x)h(x).

Metode Alternatif

Untuk menghitung E[Xk]E[X^k] distribusi heavy-tail, formula alternatif via survival function seringkali lebih mudah:

E[Xk]=0kxk1S(x)dxE[X^k] = \int_0^\infty k x^{k-1} S(x)\, dx

Khususnya untuk Pareto(α,θ)(\alpha, \theta) dengan α>k\alpha > k:

E[Xk]=k!θkj=1k(αj)=θkΓ(αk)Γ(k+1)Γ(α)E[X^k] = \frac{k!\, \theta^k}{\prod_{j=1}^{k}(\alpha - j)} = \frac{\theta^k \Gamma(\alpha - k) \Gamma(k+1)}{\Gamma(\alpha)}

Untuk k=1k=1: E[X]=θ/(α1)E[X] = \theta/(\alpha-1). Untuk k=2k=2: E[X2]=2θ2/((α1)(α2))E[X^2] = 2\theta^2/((\alpha-1)(\alpha-2)).

Section 6 — Visualisasi Mental

Diagram perbandingan hazard rate empat distribusi:

h(x)

  |  Weibull (τ>1) — naik
  |       /
  |      /
  |-----/------------- Eksponensial — konstan (1/θ)
  |    /  \
  |   /    \
  |  /      Gamma (α>1) — naik menuju 1/θ
  | /         \_______________
  |/                          ──── Pareto — turun ke 0
  +----------------------------------→ x
  0           x*

          titik silang Pareto & Eksponensial

Diagram perbandingan MEL tiga distribusi:

e(d)

  |                      / Pareto — naik linear
  |                    /
  |                  /
  |----------------/-------------- Eksponensial — konstan
  |               /
  |              /   Gamma — turun menuju θ
  |             /  \___________
  |            /
  +----------------------------------→ d
  0         E[X]

         e(0) = E[X] untuk semua distribusi

Interpretasi visual:

  • Titik awal semua kurva MEL: e(0)=E[X]e(0) = E[X] — semuanya mulai dari mean
  • Pareto: kemiringan positif — setiap kenaikan threshold menambah ekspektasi kelebihan
  • Eksponensial: garis horizontal — memoryless property secara visual
  • Gamma: kurva menurun menuju θ\theta — distribusi “terkendali” untuk klaim besar

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Rumus
Kurva hazard rate menurunh(x)<0h'(x) < 0 ↔ distribusi heavy-tailed
Kurva MEL naike(d)>0e'(d) > 0 ↔ distribusi heavy-tailed
Titik awal MEL di E[X]E[X]e(0)=E[X]=0S(x)dxe(0) = E[X] = \int_0^\infty S(x)\,dx
Kemiringan MEL Pareto =1/(α1)= 1/(\alpha-1)dddePareto(d)=1α1\frac{d}{dd}e_{\text{Pareto}}(d) = \frac{1}{\alpha-1}
Area bawah survival function0S(x)dx=E[X]\int_0^\infty S(x)\,dx = E[X]

Section 7 — Jebakan Umum

Kesalahan Parametrisasi
  • LEV Pareto: Formula E[Xu]=θα1 ⁣[1(θθ+u)α1]E[X \wedge u] = \frac{\theta}{\alpha-1}\!\left[1 - \left(\frac{\theta}{\theta+u}\right)^{\alpha-1}\right] berlaku untuk parameterisasi KPW. Beberapa sumber menggunakan Pareto dengan parameterisasi berbeda (misal Pareto Tipe II / Lomax), menghasilkan formula berbeda. Selalu verifikasi dengan E[X]=E[X]=θ/(α1)E[X \wedge \infty] = E[X] = \theta/(\alpha-1).
  • Weibull hazard rate: h(x)=τθ(xθ)τ1h(x) = \frac{\tau}{\theta}\left(\frac{x}{\theta}\right)^{\tau-1} dalam parameterisasi KPW. Parameterisasi lain menggunakan λ=θτ\lambda = \theta^{-\tau} — hasilnya berbeda. Cek: h(1)=τ/θh(1) = \tau/\theta vs h(1)=λτh(1) = \lambda\tau.
Kesalahan Konseptual
  1. Heavy-tail ≠ large mean: Distribusi bisa memiliki mean kecil tapi ekor sangat berat (atau sebaliknya). Sifat heavy-tail sepenuhnya ditentukan oleh bentuk ekor (hazard rate, MEL), bukan magnitude mean.
  2. MEL undefined untuk α1\alpha \leq 1: Untuk Pareto(α1,θ)(\alpha \leq 1, \theta), E[X]=E[X] = \infty, sehingga e(d)e(d) tidak terdefinisi. Jangan serta-merta menghitung e(d)=(d+θ)/(α1)e(d) = (d+\theta)/(\alpha-1) tanpa verifikasi α>1\alpha > 1.
  3. Hazard rate menurun ≠ distribusi tidak berguna: Meskipun hazard rate menurun mengindikasikan heavy-tail (yang lebih “berisiko”), distribusi ini tetap digunakan luas dalam aktuaria — justru karena lebih realistis untuk data klaim ekstrem.
  4. Salah arah implikasi MEL: MEL naik → heavy-tail (benar). Tetapi bukan: heavy-tail → MEL selalu ada (salah — perlu E[X]<E[X] < \infty).
Kesalahan Interpretasi Soal
  • “Tentukan apakah distribusi A lebih heavy-tailed dari B” → jangan hanya bandingkan mean atau variansi. Gunakan (1) perbandingan h(x)h(x) untuk xx besar, atau (2) apakah MEL naik/turun/konstan, atau (3) rasio SA(x)/SB(x)S_A(x)/S_B(x) saat xx \to \infty.
  • “Hitung MEL pada dd → pastikan hitung S(d)S(d) secara eksplisit, bukan asumsi S(d)=1S(d) = 1. Ini kesalahan umum ketika dd besar.
  • “Distribusi dengan hazard rate konstan” → ini hanya Eksponensial (di antara distribusi kontinu positif umum). Jangan asumsikan distribusi lain bisa memiliki hazard rate konstan.
Red Flags
  • Soal menyebut “excess-of-loss reinsurance” atau “deductible” → MEL e(d)e(d) hampir pasti digunakan; lihat juga 3.1 Coverage Modifications on Severity and Frequency
  • Soal meminta “membandingkan dua distribusi” tanpa spesifikasi ukuran → kemungkinan besar meminta perbandingan berat ekor via hazard rate atau MEL
  • Soal menyebut α1\alpha \leq 1 untuk Pareto → peringatan: E[X]E[X] tidak finite, MEL tidak terdefinisi, banyak operasi standar gagal
  • Soal memberikan h(x)h(x) dan meminta distribusi → rekonstruksi S(x)=eh(t)dtS(x) = e^{-\int h(t)\,dt}, kemudian f(x)=h(x)S(x)f(x) = h(x) S(x); kenali bentuk distribusi dari hasil S(x)S(x)
  • Soal menyebut “limiting tail behaviour” atau “ekor distribusi untuk xx besar” → perbandingan S(x)S(x) saat xx \to \infty; untuk Pareto S(x)xαS(x) \sim x^{-\alpha} (power decay), untuk Eksponensial S(x)ex/θS(x) \sim e^{-x/\theta} (exponential decay) — power decay selalu lebih lambat

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. Hazard rate dan survival function saling menentukan:
h(x)=f(x)S(x),S(x)=exp ⁣(0xh(t)dt)h(x) = \frac{f(x)}{S(x)}, \qquad S(x) = \exp\!\left(-\int_0^x h(t)\,dt\right)
  1. MEL via survival function (bentuk paling berguna):
e(d)=dS(x)dxS(d)=E[X]E[Xd]S(d)e(d) = \frac{\int_d^\infty S(x)\,dx}{S(d)} = \frac{E[X] - E[X \wedge d]}{S(d)}
  1. Titik awal MEL selalu mean:
e(0)=E[X]e(0) = E[X]
  1. Dua tanda definitif heavy-tail:
h(x) menurun (atau h(x)0)    heavy-tailh(x) \text{ menurun (atau } h(x) \to 0) \iff \text{heavy-tail} e(d) naik seiring d    heavy-taile(d) \text{ naik seiring } d \iff \text{heavy-tail}
  1. Formula hafalan kritis:
eEksponensial(d)=θ(konstan),ePareto(d)=d+θα1(naik linear)e_{\text{Eksponensial}}(d) = \theta \quad \text{(konstan)}, \qquad e_{\text{Pareto}}(d) = \frac{d + \theta}{\alpha - 1} \quad \text{(naik linear)}

Kapan Digunakan

  • Soal meminta “karakterisasi ekor”, “perbandingan distribusi”, atau “sifat distribusi” → gunakan hazard rate dan/atau MEL
  • Soal menyebut “excess-of-loss” atau “deductible” dan meminta ekspektasi kelebihan klaim → langsung gunakan e(d)e(d)
  • Soal memberikan h(x)h(x) dan meminta PDF/CDF → rekonstruksi via S(x)=ehS(x) = e^{-\int h}
  • Soal meminta momen via survival function sebagai metode alternatif → gunakan E[Xk]=0kxk1S(x)dxE[X^k] = \int_0^\infty k x^{k-1} S(x)\,dx

Kapan TIDAK Boleh Digunakan

  • MEL e(d)e(d) tidak terdefinisi jika E[X]=E[X] = \infty — selalu cek α>1\alpha > 1 untuk Pareto, dan α>0\alpha > 0 untuk Gamma
  • Hazard rate tidak dapat dibandingkan langsung antar distribusi tanpa memperhatikan domain — pastikan keduanya dievaluasi pada xx yang sama
  • Jangan gunakan “mean sama = ekor sama” sebagai argumen — mean hanya menangkap sentral distribusi, bukan perilaku ekornya

Quick Decision Tree

graph TD
    A["Soal tentang karakteristik ekor distribusi?"] -->|"Hitung h(x)"| B["h(x) = f(x) / S(x)<br>Rekonstruksi S(x) = exp(-integral h)"]
    A -->|"Hitung MEL e(d)"| C["Cek E[X] finite?"]
    A -->|"Bandingkan dua distribusi"| D["Bandingkan h(x) untuk x besar<br>ATAU bandingkan e(d): naik vs konstan"]
    C -->|"Ya (misal alpha > 1 untuk Pareto)"| E["e(d) = integral S / S(d)<br>ATAU = (E[X] - E[X^d]) / S(d)"]
    C -->|"Tidak (alpha <= 1)"| F["MEL tidak terdefinisi!<br>E[X] = infinity"]
    D -->|"h(x) turun ke 0"| G["Heavy-tail (Pareto, Lognormal)"]
    D -->|"h(x) konstan"| H["Eksponensial — memoryless"]
    D -->|"h(x) naik"| I["Light-tail (Gamma alpha>1, Weibull tau>1)"]
    E -->|"e(d) naik"| G
    E -->|"e(d) konstan"| H
    E -->|"e(d) turun"| I

Follow-up Options
  1. “Berikan contoh soal variasi MEL untuk distribusi Gamma dengan perbandingan Eksponensial”
  2. “Jelaskan hubungan 1.4 Tail Characteristics dengan 3.2 Loss Elimination Ratio and Inflation
  3. “Buat flashcard 1-halaman: hazard rate dan MEL untuk lima distribusi penting TA2”

📖 Ref: Klugman, Panjer & Willmot (2019) Loss Models 5th ed., Bab 3–5 | 🗓️ 2026-04-16 | #TA2 #TailCharacteristics #HazardRate #MeanExcessLoss #ModelBesarKlaim