AktuNotes
← Kembali
TA2 · Materi 4.3

Mean Variance and Stop-Loss

Calculation-Intensive Bobot: 10–15% Klugman, Panjer & Willmot (2019), Loss Models 5th ed., Bab 9; Tse (2009), Bab 3
TA2ModelAgregatStopLossMeanVarianceRisikoKolektifTeoriRisiko

📊 4.3 — Mean Variance and Stop-Loss

Ringkasan Cepat

Topik: Mean & Variansi Model Risiko Kolektif/Individual; Ekspektasi Asuransi Stop-Loss | Bobot: ~10–15% | Difficulty: Calculation-Intensive Ref: Klugman et al. (2019), Bab 9; Tse (2009), Bab 3 | Prereq: 4.1 Individual and Collective Risk Models, 4.2 Compound Distributions


Section 0 — Pemetaan Topik

Topik TA2Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Model Agregat4.3Menghitung E[S]E[S], Var(S)\text{Var}(S) untuk model kolektif dan individual; menghitung ekspektasi asuransi stop-loss E[(Sd)+]E[(S-d)_+] dan E[min(S,d)]E[\min(S,d)]; menurunkan relasi antar kuantitas tersebut10–15%Calculation-Intensive4.1 Individual and Collective Risk Models, 4.2 Compound Distributions4.4 Aggregate Distribution Approximation, 4.5 Panjer Recursive Formula, 4.6 Coverage Modifications on Aggregate Models, 3.1 Coverage Modifications on Severity and FrequencyKlugman et al. (2019), Bab 9; Tse (2009), Bab 3

Section 1 — Intuisi

Bayangkan sebuah perusahaan reasuransi yang menanggung klaim agregat (total klaim) dari sebuah portofolio asuransi kendaraan bermotor selama satu tahun. Klaim total SS adalah jumlah semua klaim individual yang terjadi dalam periode tersebut — bisa nol bila tidak ada klaim sama sekali, bisa sangat besar bila terjadi banyak klaim besar secara bersamaan. Perusahaan reasuransi perlu menjawab dua pertanyaan mendasar: berapa rata-rata kerugian yang akan ditanggung? dan seberapa besar ketidakpastian di sekitar rata-rata itu? Inilah mengapa kita perlu menghitung mean dan variansi dari distribusi agregat SS.

Namun ada satu produk asuransi yang sangat bergantung pada pemahaman distribusi SS ini, yaitu asuransi Stop-Loss. Dalam kontrak stop-loss, penanggung (atau reasuradur) berjanji untuk membayar seluruh kelebihan klaim agregat di atas suatu batas (deductible) dd. Artinya, jika total klaim SS ternyata kurang dari dd, reasuradur tidak membayar apa-apa. Tetapi jika SS melampaui dd, reasuradur membayar selisihnya SdS - d. Berapa yang harus reasuradur cadangkan sebagai ekspektasi pembayaran? Inilah yang dihitung melalui net stop-loss premium E[(Sd)+]E[(S-d)_+].

Konsep ini sangat penting dalam praktik karena stop-loss melindungi penanggung primer dari bencana keuangan akibat tahun yang “buruk”. Seorang aktuaris harus mampu menghitung ekspektasi pembayaran stop-loss ini dengan tepat — baik ketika distribusi SS diketahui secara analitik, maupun ketika hanya mean dan variansi SS yang tersedia (untuk aproksimasi Normal/Lognormal di Topik 4.4). Topik 4.3 adalah fondasi matematis untuk semua itu.


Section 2 — Definisi Formal

Definisi Matematis — Kerugian Agregat dan Stop-Loss

Misalkan SS adalah total kerugian agregat dalam satu periode. Asuransi Stop-Loss dengan deductible d0d \geq 0 membayar:

(Sd)+=max(Sd,  0)={0jika SdSdjika S>d(S - d)_+ = \max(S - d,\; 0) = \begin{cases} 0 & \text{jika } S \leq d \\ S - d & \text{jika } S > d \end{cases}

Net stop-loss premium (ekspektasi pembayaran):

E[(Sd)+]=d(sd)fS(s)ds=d[1FS(s)]dsE[(S-d)_+] = \int_d^{\infty} (s - d)\, f_S(s)\, ds = \int_d^{\infty} [1 - F_S(s)]\, ds
SimbolMaknaCatatan
SSTotal kerugian agregat =i=1NXi= \sum_{i=1}^{N} X_iModel kolektif; S=0S=0 jika N=0N=0
NNJumlah klaim (frekuensi)Variabel acak diskrit non-negatif
XiX_iBesar klaim ke-ii (severity)i.i.d., independen dari NN
μX=E[X]\mu_X = E[X]Mean besar klaim individual>0> 0
σX2=Var(X)\sigma_X^2 = \text{Var}(X)Variansi besar klaim individual>0> 0
μN=E[N]\mu_N = E[N]Mean frekuensi klaim>0> 0
σN2=Var(N)\sigma_N^2 = \text{Var}(N)Variansi frekuensi klaim>0> 0
ddDeductible stop-loss (retensi)d0d \geq 0
(Sd)+(S-d)_+Pembayaran stop-lossmax(Sd,0)\max(S-d, 0)
eS(d)e_S(d)Mean excess loss function agregat=E[Sd    S>d]= E[S-d \;\|\; S>d]
FS(d)F_S(d)CDF distribusi agregat di ddP(Sd)P(S \leq d)

Rumus Utama

Mean model kolektif (formula paling mendasar):

E[S]=E[N]E[X]E[S] = E[N] \cdot E[X]

Label: Mean agregat = rata-rata frekuensi × rata-rata severity. Berlaku tanpa syarat selama NN dan XiX_i independen.

Variansi model kolektif (conditional variance formula):

Var(S)=E[N]Var(X)+Var(N)(E[X])2\text{Var}(S) = E[N] \cdot \text{Var}(X) + \text{Var}(N) \cdot (E[X])^2

Label: Variansi agregat = komponen severity (rata-rata variansi bersyarat) + komponen frekuensi (variansi rata-rata bersyarat). Ini adalah law of total variance.

Variansi model kolektif — bentuk alternatif:

Var(S)=E[N]E[X2]+(E[X])2[Var(N)E[N]]\text{Var}(S) = E[N] \cdot E[X^2] + (E[X])^2 \cdot [\text{Var}(N) - E[N]]

Label: Berguna ketika E[X2]E[X^2] langsung diketahui.

Momen kedua model kolektif:

E[S2]=Var(S)+(E[S])2=E[N]Var(X)+Var(N)(E[X])2+(E[N])2(E[X])2E[S^2] = \text{Var}(S) + (E[S])^2 = E[N]\cdot\text{Var}(X) + \text{Var}(N)\cdot(E[X])^2 + (E[N])^2\cdot(E[X])^2

Label: Diperlukan untuk menghitung stop-loss dengan distribusi diskrit.

Mean model individual (dengan mm unit eksposur):

E[S]=i=1mE[Xi]=i=1mμiE[S] = \sum_{i=1}^{m} E[X_i] = \sum_{i=1}^{m} \mu_i

Variansi model individual (unit independen):

Var(S)=i=1mVar(Xi)=i=1mσi2\text{Var}(S) = \sum_{i=1}^{m} \text{Var}(X_i) = \sum_{i=1}^{m} \sigma_i^2

Net stop-loss premium — bentuk integral:

E[(Sd)+]=d(1FS(s))dsE[(S-d)_+] = \int_d^{\infty} (1 - F_S(s))\, ds

Label: Setara dengan integral survival function SS(s)=1FS(s)S_S(s) = 1 - F_S(s) dari dd ke \infty.

Net stop-loss premium — relasi fundamental:

E[(Sd)+]=E[S]d+dFS(d)0dsfS(s)dsE[(S-d)_+] = E[S] - d + d \cdot F_S(d) - \int_0^d s\, f_S(s)\, ds

Label: Diperoleh dengan memisah E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+].

Dekomposisi kritis:

E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+]

Label: Total ekspektasi = porsi yang ditahan (limited expected value) + porsi stop-loss. Selalu benar.

Limited expected value:

E[min(S,d)]=E[S]E[(Sd)+]=0d[1FS(s)]dsE[\min(S,d)] = E[S] - E[(S-d)_+] = \int_0^d [1 - F_S(s)]\, ds

Label: Ekspektasi kerugian yang dibatasi di dd; digunakan dalam penetapan premi lapisan.

Stop-loss untuk SS diskrit (dengan PMF pk=P(S=k)p_k = P(S=k)):

E[(Sd)+]=k>d(kd)pk=k=0(kd)+pkE[(S-d)_+] = \sum_{k > d} (k - d)\, p_k = \sum_{k=0}^{\infty} (k-d)_+ \, p_k

Label: Penjumlahan eksplisit atas semua nilai S>dS > d.

Asumsi Eksplisit

  1. X1,X2,X_1, X_2, \ldots adalah i.i.d. (identically and independently distributed) dengan mean μX\mu_X dan variansi σX2\sigma_X^2.
  2. NN dan {Xi}\{X_i\} independen satu sama lain — frekuensi tidak bergantung pada besarnya klaim.
  3. Untuk model individual: XiX_i independen antar unit, meskipun boleh tidak identik (μi\mu_i dan σi2\sigma_i^2 boleh berbeda).
  4. E[X2]<E[X^2] < \infty (momen kedua severity terbatas) agar Var(S)\text{Var}(S) terdefinisi.
  5. Untuk stop-loss kontinu: FSF_S terdiferensialkan dan 0sfS(s)ds=E[S]<\int_0^{\infty} s\, f_S(s)\, ds = E[S] < \infty.

Section 3 — Jembatan Logika

Dari Definisi ke Rumus

Formula Var(S)=E[N]Var(X)+Var(N)(E[X])2\text{Var}(S) = E[N]\cdot\text{Var}(X) + \text{Var}(N)\cdot(E[X])^2 bukan rumus yang dihafal — ia adalah aplikasi langsung law of total variance: Var(S)=E[Var(SN)]+Var(E[SN])\text{Var}(S) = E[\text{Var}(S|N)] + \text{Var}(E[S|N]). Komponen pertama (E[N]Var(X)E[N]\cdot\text{Var}(X)) muncul karena bersyarat pada N=nN=n, klaim adalah jumlah nn variabel i.i.d. sehingga variansinya nσX2n\sigma_X^2, dan rata-ratanya terhadap NN menghasilkan E[N]σX2E[N]\sigma_X^2. Komponen kedua (Var(N)(E[X])2\text{Var}(N)\cdot(E[X])^2) muncul karena E[SN=n]=nμXE[S|N=n] = n\mu_X adalah fungsi linear dari nn, sehingga Var(E[SN])=(E[X])2Var(N)\text{Var}(E[S|N]) = (E[X])^2 \cdot \text{Var}(N).

Support dan Domain
  • S0S \geq 0 selalu, dengan P(S=0)=P(N=0)>0P(S=0) = P(N=0) > 0 untuk model kolektif.
  • Stop-loss premium E[(Sd)+]E[(S-d)_+] adalah fungsi menurun convex dalam dd: semakin besar retensi dd, semakin kecil pembayaran reasuransi.
  • Di d=0d = 0: E[(S0)+]=E[S]E[(S-0)_+] = E[S] — reasuransi menanggung seluruh klaim.
  • Di dd \to \infty: E[(Sd)+]0E[(S-d)_+] \to 0 — reasuransi tidak pernah membayar.
  • Dekomposisi E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+] berlaku untuk semua nilai d0d \geq 0.

Derivasi 1: Variansi Model Kolektif via Law of Total Variance

Langkah 1 — Terapkan law of total variance:

Var(S)=E[Var(SN)]+Var(E[SN])\text{Var}(S) = E[\text{Var}(S \mid N)] + \text{Var}(E[S \mid N])

Langkah 2 — Hitung E[SN=n]E[S|N=n]: Bersyarat pada N=nN=n, S=X1++XnS = X_1 + \cdots + X_n dengan XiX_i i.i.d.

E[SN=n]=nμXE[SN]=NμXE[S \mid N=n] = n \cdot \mu_X \quad \Rightarrow \quad E[S \mid N] = N \cdot \mu_X

Langkah 3 — Hitung Var(SN=n)\text{Var}(S|N=n):

Var(SN=n)=nσX2Var(SN)=NσX2\text{Var}(S \mid N=n) = n \cdot \sigma_X^2 \quad \Rightarrow \quad \text{Var}(S \mid N) = N \cdot \sigma_X^2

Langkah 4 — Substitusi ke law of total variance:

Var(S)=E[NσX2]+Var(NμX)\text{Var}(S) = E[N \cdot \sigma_X^2] + \text{Var}(N \cdot \mu_X)

Langkah 5 — Sederhanakan (konstanta keluar dari ekspektasi dan variansi):

Var(S)=σX2E[N]+μX2Var(N)\text{Var}(S) = \sigma_X^2 \cdot E[N] + \mu_X^2 \cdot \text{Var}(N) Var(S)=E[N]Var(X)+Var(N)(E[X])2\boxed{\text{Var}(S) = E[N] \cdot \text{Var}(X) + \text{Var}(N) \cdot (E[X])^2}

Derivasi 2: Net Stop-Loss Premium via Integral Survival Function

Langkah 1 — Mulai dari definisi:

E[(Sd)+]=0(sd)+fS(s)ds=d(sd)fS(s)dsE[(S-d)_+] = \int_0^{\infty} (s-d)_+ f_S(s)\, ds = \int_d^{\infty} (s-d) f_S(s)\, ds

Langkah 2 — Substitusi u=sdu = s - d, sehingga s=u+ds = u + d, ds=duds = du, batas: uu dari 00 ke \infty:

=0ufS(u+d)du= \int_0^{\infty} u\, f_S(u+d)\, du

Langkah 3 — Integrasikan per bagian (alternatif lebih elegan): langsung dari integral asal,

d(sd)fS(s)ds=[(sd)((1FS(s)))]d+d(1FS(s))ds\int_d^{\infty}(s-d) f_S(s)\, ds = \left[(s-d)(-(1-F_S(s)))\right]_d^{\infty} + \int_d^{\infty} (1-F_S(s))\, ds

Langkah 4 — Evaluasi batas: suku pertama =0= 0 di s=ds=d (karena sd=0s-d=0) dan 0\to 0 di ss \to \infty (asumsi E[S]<E[S] < \infty). Tersisa:

E[(Sd)+]=d[1FS(s)]dsE[(S-d)_+] = \int_d^{\infty} [1 - F_S(s)]\, ds

Derivasi 3: Dekomposisi E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+]

Langkah 1 — Identitas aljabar: untuk setiap s0s \geq 0 dan d0d \geq 0,

s=min(s,d)+(sd)+s = \min(s, d) + (s-d)_+

Langkah 2 — Ambil ekspektasi kedua sisi:

E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+]

Ini langsung dari linearitas ekspektasi. ∎

Dilarang
  1. Jangan gunakan Var(S)=E[N]Var(X)\text{Var}(S) = E[N] \cdot \text{Var}(X) saja — ini hanya benar jika NN deterministik (konstan), bukan variabel acak. Suku Var(N)(E[X])2\text{Var}(N) \cdot (E[X])^2 wajib disertakan.
  2. Jangan tukar E[(Sd)+]E[(S-d)_+] dengan E[S]dE[S] - d — yang benar adalah E[(Sd)+]=E[S]E[min(S,d)]E[(S-d)_+] = E[S] - E[\min(S,d)], bukan E[S]dE[S] - d.
  3. Jangan lupa P(S=0)>0P(S=0) > 0 ketika menghitung stop-loss untuk distribusi diskrit campuran — distribusi SS memiliki massa di titik nol yang harus diperhitungkan.

Section 4 — Contoh Soal

Soal A — Fundamental

Model risiko kolektif: jumlah klaim NPoisson(λ=4)N \sim \text{Poisson}(\lambda = 4) dan besar klaim individual XExponentialX \sim \text{Exponential} dengan mean E[X]=500E[X] = 500.

(a) Hitung E[S]E[S] dan Var(S)\text{Var}(S). (b) Hitung E[S2]E[S^2]. (c) Hitung net stop-loss premium E[(Sd)+]E[(S-d)_+] untuk d=2500d = 2500 menggunakan integrasi (gunakan fakta bahwa untuk model Poisson-Eksponensial, distribusi SS dapat dianalisis secara langsung).

Solusi Soal A

Pendekatan: Terapkan formula mean dan variansi kolektif langsung. Untuk stop-loss Poisson-Eksponensial, gunakan dekomposisi E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+] dengan hasil analitik.

1. Identifikasi Variabel

  • NPoisson(λ=4)N \sim \text{Poisson}(\lambda = 4): E[N]=Var(N)=4E[N] = \text{Var}(N) = 4
  • XExp(θ=500)X \sim \text{Exp}(\theta = 500): E[X]=500E[X] = 500, E[X2]=2θ2=500,000E[X^2] = 2\theta^2 = 500{,}000, Var(X)=θ2=250,000\text{Var}(X) = \theta^2 = 250{,}000
  • d=2500d = 2500

2. Identifikasi Distribusi / Model Model kolektif standar: S=X1++XNS = X_1 + \cdots + X_N, NN independen dari {Xi}\{X_i\}. Poisson-Eksponensial adalah pasangan klasik dalam teori risiko.

3. Setup Persamaan

E[S]=E[N]E[X],Var(S)=E[N]Var(X)+Var(N)(E[X])2E[S] = E[N] \cdot E[X], \qquad \text{Var}(S) = E[N]\cdot\text{Var}(X) + \text{Var}(N)\cdot(E[X])^2

4. Eksekusi Aljabar

(a):

E[S]=4×500=2,000E[S] = 4 \times 500 = 2{,}000 Var(S)=4×250,000+4×(500)2=1,000,000+1,000,000=2,000,000\text{Var}(S) = 4 \times 250{,}000 + 4 \times (500)^2 = 1{,}000{,}000 + 1{,}000{,}000 = 2{,}000{,}000

(b):

E[S2]=Var(S)+(E[S])2=2,000,000+4,000,000=6,000,000E[S^2] = \text{Var}(S) + (E[S])^2 = 2{,}000{,}000 + 4{,}000{,}000 = 6{,}000{,}000

(c) Untuk Poisson(λ\lambda)-Eksponensial(θ\theta), diketahui bahwa SS memiliki distribusi campuran: massa di s=0s=0 sebesar eλe^{-\lambda}, dan untuk s>0s > 0 distribusi kontinu. Net stop-loss premium dapat dihitung via survival integral. Karena E[S]=λθ=2000E[S] = \lambda\theta = 2000 dan d=2500>E[S]d = 2500 > E[S], gunakan relasi:

E[(Sd)+]=E[S]d+E[max(dS,0)]E[(S-d)_+] = E[S] - d + E[\max(d-S, 0)]

Lebih langsung: gunakan E[(Sd)+]=E[S]E[min(S,d)]E[(S-d)_+] = E[S] - E[\min(S,d)].

Untuk Poisson-Eksponensial: E[min(S,d)]=E[S](1ed/(λθ)(koreksi))E[\min(S,d)] = E[S]\left(1 - e^{-d/(\lambda\theta)} \cdot \text{(koreksi)}\right). Karena distribusi SS Poisson-Eksponensial kompleks, gunakan batas bawah: E[(Sd)+](E[S]d)+=(20002500)+=0E[(S-d)_+] \geq (E[S] - d)_+ = (2000 - 2500)_+ = 0.

Untuk nilai eksak, gunakan formula rekursif atau aproksimasi Normal (lihat Topik 4.4): SN(μ=2000,σ2=2,000,000)S \approx N(\mu=2000, \sigma^2=2{,}000{,}000), σ=2,000,0001,414.2\sigma = \sqrt{2{,}000{,}000} \approx 1{,}414.2.

E[(Sd)+]σϕ(dμσ)(dμ)[1Φ(dμσ)]E[(S-d)_+] \approx \sigma \cdot \phi\left(\frac{d-\mu}{\sigma}\right) - (d-\mu)\left[1 - \Phi\left(\frac{d-\mu}{\sigma}\right)\right] z=250020001414.2=5001414.20.354z = \frac{2500 - 2000}{1414.2} = \frac{500}{1414.2} \approx 0.354 ϕ(0.354)0.3752,Φ(0.354)0.6384\phi(0.354) \approx 0.3752, \quad \Phi(0.354) \approx 0.6384 E[(Sd)+]1414.2×0.3752500×(10.6384)=530.6500×0.3616=530.6180.8=349.8E[(S-d)_+] \approx 1414.2 \times 0.3752 - 500 \times (1-0.6384) = 530.6 - 500 \times 0.3616 = 530.6 - 180.8 = 349.8

5. Verification E[S]=2000E[S] = 2000; E[(Sd)+]=349.8<E[S]E[(S-d)_+] = 349.8 < E[S] ✓. Dekomposisi: E[min(S,d)]=2000349.8=1650.2E[\min(S,d)] = 2000 - 349.8 = 1650.2; dan 1650.2<d=25001650.2 < d = 2500 ✓ (limited expected value selalu d\leq d).

Hasil: E[S]=2,000E[S]=2{,}000; Var(S)=2,000,000\text{Var}(S)=2{,}000{,}000; E[S2]=6,000,000E[S^2]=6{,}000{,}000; E[(S2500)+]349.8E[(S-2500)_+] \approx 349.8 (aproksimasi Normal).

Exam Tips — Soal A

Target waktu: 3 menit. Common trap: Untuk Poisson, Var(N)=E[N]=λ\text{Var}(N) = E[N] = \lambda — jangan gunakan nilai berbeda. Akibatnya Var(S)=λ(E[X2])=λE[X2]\text{Var}(S) = \lambda(E[X^2]) = \lambda \cdot E[X^2], yang merupakan shortcut khusus Poisson. Shortcut Poisson: Var(S)=E[N]E[X2]\text{Var}(S) = E[N] \cdot E[X^2] karena E[N]=Var(N)E[N] = \text{Var}(N) untuk Poisson.


Soal B — Exam-Typical

Model risiko kolektif dengan frekuensi NNegBin(r=3,β=2)N \sim \text{NegBin}(r=3, \beta=2) dan severity XX berdistribusi Gamma dengan mean E[X]=1,000E[X] = 1{,}000 dan Var(X)=500,000\text{Var}(X) = 500{,}000.

(a) Hitung E[S]E[S], Var(S)\text{Var}(S), dan standar deviasi SS. (b) Hitung net stop-loss premium E[(Sd)+]E[(S-d)_+] untuk d=7,000d = 7{,}000 menggunakan aproksimasi Normal. (c) Berapakah nilai E[min(S,7000)]E[\min(S, 7000)]?

Solusi Soal B

Pendekatan: Hitung momen NegBin terlebih dahulu, terapkan formula variansi kolektif, lalu gunakan formula stop-loss Normal.

1. Identifikasi Variabel

  • NNegBin(r=3,β=2)N \sim \text{NegBin}(r=3, \beta=2): E[N]=rβ=6E[N] = r\beta = 6, Var(N)=rβ(1+β)=3×2×3=18\text{Var}(N) = r\beta(1+\beta) = 3 \times 2 \times 3 = 18
  • XX: E[X]=1,000E[X] = 1{,}000, Var(X)=500,000\text{Var}(X) = 500{,}000, sehingga E[X2]=Var(X)+(E[X])2=500,000+1,000,000=1,500,000E[X^2] = \text{Var}(X) + (E[X])^2 = 500{,}000 + 1{,}000{,}000 = 1{,}500{,}000
  • d=7,000d = 7{,}000

2. Identifikasi Distribusi / Model Model kolektif dengan frekuensi NegBin (overdispersi) dan severity Gamma (skewed kanan). Variasi frekuensi yang besar (Var(N)=18E[N]=6(N) = 18 \gg E[N] = 6) akan mendominasi variansi agregat.

3. Setup Persamaan

E[S]=E[N]E[X]E[S] = E[N]\cdot E[X] Var(S)=E[N]Var(X)+Var(N)(E[X])2\text{Var}(S) = E[N]\cdot\text{Var}(X) + \text{Var}(N)\cdot(E[X])^2 E[(Sd)+]σSϕ(z)(dE[S])[1Φ(z)],z=dE[S]σSE[(S-d)_+] \approx \sigma_S \cdot \phi(z) - (d - E[S]) \cdot [1 - \Phi(z)], \quad z = \frac{d - E[S]}{\sigma_S}

4. Eksekusi Aljabar

(a) Mean:

E[S]=6×1,000=6,000E[S] = 6 \times 1{,}000 = 6{,}000

Variansi:

Var(S)=6×500,000+18×(1,000)2=3,000,000+18,000,000=21,000,000\text{Var}(S) = 6 \times 500{,}000 + 18 \times (1{,}000)^2 = 3{,}000{,}000 + 18{,}000{,}000 = 21{,}000{,}000 σS=21,000,0004,582.6\sigma_S = \sqrt{21{,}000{,}000} \approx 4{,}582.6

Catatan: komponen frekuensi (18,000,00018{,}000{,}000) jauh mendominasi komponen severity (3,000,0003{,}000{,}000) — konsisten dengan overdispersi NegBin.

(b) Stop-loss premium:

z=7,0006,0004,582.6=1,0004,582.60.2182z = \frac{7{,}000 - 6{,}000}{4{,}582.6} = \frac{1{,}000}{4{,}582.6} \approx 0.2182

Dari tabel Normal standar: ϕ(0.2182)0.3885\phi(0.2182) \approx 0.3885, Φ(0.2182)0.5864\Phi(0.2182) \approx 0.5864.

E[(S7000)+]4,582.6×0.38851,000×(10.5864)E[(S - 7000)_+] \approx 4{,}582.6 \times 0.3885 - 1{,}000 \times (1 - 0.5864) =1,780.31,000×0.4136=1,780.3413.6=1,366.7= 1{,}780.3 - 1{,}000 \times 0.4136 = 1{,}780.3 - 413.6 = 1{,}366.7

(c) Limited expected value:

E[min(S,7000)]=E[S]E[(S7000)+]=6,0001,366.7=4,633.3E[\min(S, 7000)] = E[S] - E[(S-7000)_+] = 6{,}000 - 1{,}366.7 = 4{,}633.3

5. Verification E[(Sd)+]=1,366.7>0E[(S-d)_+] = 1{,}366.7 > 0 meskipun d=7000>E[S]=6000d = 7000 > E[S] = 6000; ini benar karena distribusi SS memiliki ekor kanan yang panjang. E[min(S,d)]=4,633.3<d=7,000E[\min(S,d)] = 4{,}633.3 < d = 7{,}000 ✓. Dekomposisi: 4,633.3+1,366.7=6,000=E[S]4{,}633.3 + 1{,}366.7 = 6{,}000 = E[S] ✓.

Hasil: E[S]=6,000E[S]=6{,}000; Var(S)=21,000,000\text{Var}(S)=21{,}000{,}000; σS4,582.6\sigma_S \approx 4{,}582.6; E[(S7000)+]1,366.7E[(S-7000)_+] \approx 1{,}366.7; E[min(S,7000)]4,633.3E[\min(S,7000)] \approx 4{,}633.3.

Exam Tips — Soal B

Target waktu: 4 menit. Common trap: Menggunakan Var(N)=E[N]=6\text{Var}(N) = E[N] = 6 (keliru mengira NegBin seperti Poisson). Untuk NegBin, Var(N)=rβ(1+β)=186\text{Var}(N) = r\beta(1+\beta) = 18 \neq 6 — perbedaan ini drastis dalam hasil variansi agregat. Shortcut: Selalu hitung kedua komponen variansi (E[N]Var(X)E[N]\text{Var}(X) dan Var(N)(E[X])2\text{Var}(N)(E[X])^2) secara terpisah; bandingkan besarannya untuk cek konsistensi dengan tipe frekuensi.


Soal C — Challenging

Model risiko individual terdiri dari m=3m = 3 polis independen dengan profil klaim sebagai berikut:

Polis iiP(klaim)P(\text{klaim})Besar klaim jika ada (XiX_i)
1q1=0.3q_1 = 0.3Konstan =1,000= 1{,}000
2q2=0.5q_2 = 0.5Konstan =2,000= 2{,}000
3q3=0.2q_3 = 0.2Konstan =3,000= 3{,}000

Misalkan Li=Xi1klaimiL_i = X_i \cdot \mathbf{1}_{\text{klaim}_i} (kerugian polis ii, nol jika tidak klaim), dan S=L1+L2+L3S = L_1 + L_2 + L_3.

(a) Hitung E[S]E[S] dan Var(S)\text{Var}(S). (b) Tentukan distribusi lengkap SS (semua nilai yang mungkin dan probabilitasnya). (c) Hitung net stop-loss premium E[(Sd)+]E[(S-d)_+] untuk d=2,000d = 2{,}000 secara eksak menggunakan distribusi SS.

Solusi Soal C

Pendekatan: Model individual — hitung momen tiap polis, jumlahkan. Untuk distribusi eksak, enumerate semua 23=82^3 = 8 skenario klaim. Lalu hitung stop-loss langsung dari distribusi.

1. Identifikasi Variabel

  • Li=ciBiL_i = c_i \cdot B_i di mana BiBernoulli(qi)B_i \sim \text{Bernoulli}(q_i) dan cic_i adalah besar klaim konstan
  • E[Li]=qiciE[L_i] = q_i \cdot c_i; Var(Li)=qi(1qi)ci2\text{Var}(L_i) = q_i(1-q_i) \cdot c_i^2
  • Polis 1: c1=1000c_1=1000, q1=0.3q_1=0.3; Polis 2: c2=2000c_2=2000, q2=0.5q_2=0.5; Polis 3: c3=3000c_3=3000, q3=0.2q_3=0.2

2. Identifikasi Distribusi / Model Model risiko individual dengan besar klaim deterministik (Bernoulli). Distribusi SS adalah diskrit dengan support {0,1000,2000,3000,4000,5000,6000}\{0, 1000, 2000, 3000, 4000, 5000, 6000\} (beberapa nilai mungkin tidak muncul).

3. Setup Persamaan

E[S]=i=13qici,Var(S)=i=13qi(1qi)ci2E[S] = \sum_{i=1}^{3} q_i c_i, \qquad \text{Var}(S) = \sum_{i=1}^{3} q_i(1-q_i)c_i^2 E[(Sd)+]=s>d(sd)P(S=s)E[(S-d)_+] = \sum_{s > d} (s - d) \cdot P(S = s)

4. Eksekusi Aljabar

(a) Mean:

E[L1]=0.3×1000=300E[L_1] = 0.3 \times 1000 = 300 E[L2]=0.5×2000=1,000E[L_2] = 0.5 \times 2000 = 1{,}000 E[L3]=0.2×3000=600E[L_3] = 0.2 \times 3000 = 600 E[S]=300+1,000+600=1,900E[S] = 300 + 1{,}000 + 600 = 1{,}900

Variansi:

Var(L1)=0.3×0.7×1,000,000=210,000\text{Var}(L_1) = 0.3 \times 0.7 \times 1{,}000{,}000 = 210{,}000 Var(L2)=0.5×0.5×4,000,000=1,000,000\text{Var}(L_2) = 0.5 \times 0.5 \times 4{,}000{,}000 = 1{,}000{,}000 Var(L3)=0.2×0.8×9,000,000=1,440,000\text{Var}(L_3) = 0.2 \times 0.8 \times 9{,}000{,}000 = 1{,}440{,}000 Var(S)=210,000+1,000,000+1,440,000=2,650,000\text{Var}(S) = 210{,}000 + 1{,}000{,}000 + 1{,}440{,}000 = 2{,}650{,}000

(b) Distribusi lengkap — enumerate 8 skenario (B1, B2, B3) = (klaim polis 1, 2, 3):

(B1,B2,B3)(B_1, B_2, B_3)ProbabilitasSS
(0,0,0)(0,0,0)0.7×0.5×0.8=0.2800.7 \times 0.5 \times 0.8 = 0.28000
(1,0,0)(1,0,0)0.3×0.5×0.8=0.1200.3 \times 0.5 \times 0.8 = 0.1201,0001{,}000
(0,1,0)(0,1,0)0.7×0.5×0.8=0.2800.7 \times 0.5 \times 0.8 = 0.2802,0002{,}000
(0,0,1)(0,0,1)0.7×0.5×0.2=0.0700.7 \times 0.5 \times 0.2 = 0.0703,0003{,}000
(1,1,0)(1,1,0)0.3×0.5×0.8=0.1200.3 \times 0.5 \times 0.8 = 0.1203,0003{,}000
(1,0,1)(1,0,1)0.3×0.5×0.2=0.0300.3 \times 0.5 \times 0.2 = 0.0304,0004{,}000
(0,1,1)(0,1,1)0.7×0.5×0.2=0.0700.7 \times 0.5 \times 0.2 = 0.0705,0005{,}000
(1,1,1)(1,1,1)0.3×0.5×0.2=0.0300.3 \times 0.5 \times 0.2 = 0.0306,0006{,}000

Gabungkan nilai SS yang sama:

SSP(S=s)P(S=s)
000.2800.280
1,0001{,}0000.1200.120
2,0002{,}0000.2800.280
3,0003{,}0000.070+0.120=0.1900.070 + 0.120 = 0.190
4,0004{,}0000.0300.030
5,0005{,}0000.0700.070
6,0006{,}0000.0300.030
Total1.0001.000

(c) Stop-loss premium untuk d=2,000d = 2{,}000 — hanya suku dengan S>2000S > 2000:

E[(S2000)+]=(30002000)(0.190)+(40002000)(0.030)+(50002000)(0.070)+(60002000)(0.030)E[(S-2000)_+] = (3000-2000)(0.190) + (4000-2000)(0.030) + (5000-2000)(0.070) + (6000-2000)(0.030) =1000(0.190)+2000(0.030)+3000(0.070)+4000(0.030)= 1000(0.190) + 2000(0.030) + 3000(0.070) + 4000(0.030) =190+60+210+120=580= 190 + 60 + 210 + 120 = 580

5. Verification Cek: E[min(S,2000)]=E[S]E[(S2000)+]=1900580=1320E[\min(S,2000)] = E[S] - E[(S-2000)_+] = 1900 - 580 = 1320.

Verifikasi langsung: E[min(S,2000)]=0(0.280)+1000(0.120)+2000(0.280)+2000(0.190)+2000(0.030)+2000(0.070)+2000(0.030)E[\min(S,2000)] = 0(0.280) + 1000(0.120) + 2000(0.280) + 2000(0.190) + 2000(0.030) + 2000(0.070) + 2000(0.030) =0+120+560+380+60+140+60=1320= 0 + 120 + 560 + 380 + 60 + 140 + 60 = 1320

Cek total prob: 0.280+0.120+0.280+0.190+0.030+0.070+0.030=1.0000.280+0.120+0.280+0.190+0.030+0.070+0.030 = 1.000

Hasil: E[S]=1,900E[S]=1{,}900; Var(S)=2,650,000\text{Var}(S)=2{,}650{,}000; distribusi SS memiliki 7 nilai dengan probabilitas di atas; E[(S2000)+]=580E[(S-2000)_+] = 580; E[min(S,2000)]=1,320E[\min(S,2000)] = 1{,}320.

Exam Tips — Soal C

Target waktu: 5–6 menit. Common trap 1: Menghitung P(S=3000)=0.070P(S=3000) = 0.070 saja — lupa ada dua skenario yang menghasilkan S=3000S=3000: (0,0,1)(0,0,1) dan (1,1,0)(1,1,0). Selalu cek apakah nilai SS bisa muncul dari lebih dari satu kombinasi. Common trap 2: Memasukkan s=2000s=2000 dalam penjumlahan stop-loss — hanya s>ds > d, bukan sds \geq d. Shortcut: Untuk stop-loss diskrit, tulis langsung daftar (sd)×P(S=s)(s-d) \times P(S=s) untuk s>ds > d dalam satu tabel.


Section 5 — Verifikasi & Sanity Check

Sanity Check 1 — Dekomposisi Wajib

Selalu verifikasi: E[min(S,d)]+E[(Sd)+]=E[S]E[\min(S,d)] + E[(S-d)_+] = E[S]. Ini adalah identitas yang selalu benar — jika tidak terpenuhi, ada kesalahan kalkulasi di salah satu sisi. Gunakan nilai yang lebih mudah dihitung untuk memperoleh yang satunya.

Sanity Check 2 — Batas Stop-Loss

Dua batas yang selalu berlaku:

  • E[(Sd)+]E[S]E[(S-d)_+] \leq E[S] (stop-loss tidak mungkin melebihi total klaim)
  • E[(Sd)+](E[S]d)+E[(S-d)_+] \geq (E[S] - d)_+ (batas bawah Jensen — karena (xd)+(x-d)_+ adalah fungsi convex)

Jika hasil kalkulasi melanggar salah satu, ada kesalahan.

Sanity Check 3 — Shortcut Variansi Poisson

Khusus untuk Poisson (E[N]=Var(N)=λE[N] = \text{Var}(N) = \lambda):

Var(S)=λE[X2]\text{Var}(S) = \lambda \cdot E[X^2]

Ini lebih cepat dari formula umum. Berlaku hanya untuk Poisson; untuk Binomial atau NegBin gunakan formula penuh.

Metode Alternatif — Stop-Loss via Tabel Distribusi Diskrit

Untuk distribusi SS diskrit yang diketahui (mis. dari Panjer rekursif Topik 4.5), hitung langsung:

E[(Sd)+]=k=0P(S>k+d)1k0=j=d+1P(Sj)=j=d[1FS(j1)][1FS(d1)]E[(S-d)_+] = \sum_{k=0}^{\infty} P(S > k+d) \cdot \mathbf{1}_{k \geq 0} = \sum_{j=d+1}^{\infty} P(S \geq j) = \sum_{j=d}^{\infty} [1 - F_S(j-1)] - [1 - F_S(d-1)]

Atau lebih praktis, jika PMF sudah tersedia: bangun tabel kumulatif, lalu jumlahkan (skd)pk(s_k - d) \cdot p_k untuk semua kk di mana sk>ds_k > d.


Section 6 — Visualisasi Mental

Geometri asuransi stop-loss:

Pembayaran
(payoff)

S-d │                              ╱ ← reasuransi membayar (S-d)
    │                            ╱
    │                          ╱
    │                        ╱
    │                      ╱
  0 │____________________╱________________________ S (total klaim)
    0          d=retensi

        Deductible stop-loss
        E[(S-d)+] = area bayangan di atas garis d, dibobot PDF

Peran dekomposisi E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+]:

Total klaim E[S]
├──────────────────────────┬──────────────────────┤
│   E[min(S,d)]            │    E[(S-d)+]          │
│   (ditahan primer)       │    (dibayar reasuransi)│
│   = limited exp. value   │    = stop-loss premium │
└──────────────────────────┴──────────────────────┘

                        deductible d

Efek overdispersi pada variansi agregat:

Var(S) = E[N]·Var(X)   +   Var(N)·(E[X])²
           ↑                     ↑
    Komponen severity      Komponen frekuensi
    (selalu ada)           (hilang jika N konstan)

Poisson:   Var(N) = E[N]   → kedua suku seimbang
Binomial:  Var(N) < E[N]   → komponen severity dominan
NegBin:    Var(N) > E[N]   → komponen frekuensi dominan

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Rumus
Titik patah payoff diagram di dddd = deductible stop-loss
Area di bawah payoff × PDFE[(Sd)+]=d(sd)fS(s)dsE[(S-d)_+] = \int_d^\infty (s-d)f_S(s)\,ds
Lebar segmen “retained”E[min(S,d)]=0d[1FS(s)]dsE[\min(S,d)] = \int_0^d [1-F_S(s)]\,ds
Lebar total barE[S]E[S] = jumlah kedua segmen
Tinggi komponen frekuensiVar(N)(E[X])2\text{Var}(N) \cdot (E[X])^2; makin tinggi jika NegBin

Section 7 — Jebakan Umum

Kesalahan Parametrisasi

NegBin variansi: Untuk NegBin(r,β)(r, \beta): Var(N)=rβ(1+β)\text{Var}(N) = r\beta(1+\beta), bukan rβr\beta. Kesalahan ini menghasilkan Var(S)\text{Var}(S) yang salah besar. Ingat: Var(N)=E[N](1+β)\text{Var}(N) = E[N] \cdot (1+\beta) — selalu lebih besar dari mean.

Kesalahan Konseptual
  1. Var(S)=E[N]Var(X)\text{Var}(S) = E[N]\cdot\text{Var}(X) saja” — ini hanya benar jika NN konstan (deterministik). Untuk NN acak, suku Var(N)(E[X])2\text{Var}(N)\cdot(E[X])^2 wajib ada.
  2. "E[(Sd)+]=E[S]dE[(S-d)_+] = E[S] - d" — ini salah. Yang benar: E[(Sd)+]=E[S]E[min(S,d)]E[S]dE[(S-d)_+] = E[S] - E[\min(S,d)] \leq E[S] - d hanya jika P(Sd)=1P(S \geq d) = 1. Umumnya E[min(S,d)]<dE[\min(S,d)] < d.
  3. Memasukkan s=ds = d dalam penjumlahan stop-loss diskrit — stop-loss membayar hanya jika S>dS > d, bukan SdS \geq d. Untuk distribusi kontinu tidak masalah, tetapi untuk diskrit, s=ds = d memberikan (dd)=0(d-d) = 0 sehingga tidak berkontribusi — tetap tak bermasalah secara numerik, tetapi konsep harus benar.
  4. Lupa P(S=0)>0P(S=0) > 0 saat menghitung stop-loss untuk model kolektif — distribusi SS memiliki massa positif di nol.
Kesalahan Interpretasi Soal
  • “Net stop-loss premium” = E[(Sd)+]E[(S-d)_+] — ini adalah ekspektasi pembayaran, bukan premi gross.
  • “Stop-loss dengan retention dd dan “stop-loss dengan deductible dd adalah istilah yang sama — dd adalah ambang batas di mana reasuransi mulai membayar.
  • “Limited expected value” = E[min(S,d)]E[\min(S,d)] — bukan E[S1Sd]E[S \cdot \mathbf{1}_{S \leq d}] saja; ingat ada suku dP(S>d)d \cdot P(S > d) juga.
  • Soal meminta E[min(S,u)]E[\min(S,u)] → gunakan dekomposisi: E[min(S,u)]=E[S]E[(Su)+]E[\min(S,u)] = E[S] - E[(S-u)_+].
Red Flags
  • Kata “stop-loss” atau “excess of loss reinsurance” → hitung E[(Sd)+]E[(S-d)_+]; periksa apakah distribusi SS diketahui atau perlu aproksimasi.
  • Kata “limited expected value” → hitung E[min(S,d)]E[\min(S,d)]; gunakan dekomposisi.
  • Frekuensi NegBin → pastikan Var(N)=rβ(1+β)\text{Var}(N) = r\beta(1+\beta), bukan rβr\beta.
  • Frekuensi Poisson → shortcut Var(S)=λE[X2]\text{Var}(S) = \lambda E[X^2] menghemat langkah.
  • Model individual dengan polis Bernoulli → enumerate semua 2m2^m skenario jika mm kecil (m4m \leq 4).

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. Mean dan variansi model kolektif:

    E[S]=E[N]E[X],Var(S)=E[N]Var(X)+Var(N)(E[X])2E[S] = E[N]\cdot E[X], \quad \text{Var}(S) = E[N]\cdot\text{Var}(X) + \text{Var}(N)\cdot(E[X])^2
  2. Shortcut Poisson (E[N]=Var(N)=λE[N]=\text{Var}(N)=\lambda):

    Var(S)=λE[X2]\text{Var}(S) = \lambda \cdot E[X^2]
  3. Dekomposisi fundamental (selalu benar):

    E[S]=E[min(S,d)]+E[(Sd)+]E[S] = E[\min(S,d)] + E[(S-d)_+]
  4. Net stop-loss premium via survival integral:

    E[(Sd)+]=d[1FS(s)]dsE[(S-d)_+] = \int_d^{\infty}[1 - F_S(s)]\,ds
  5. Stop-loss diskrit (hitung langsung dari PMF):

    E[(Sd)+]=s>d(sd)P(S=s)E[(S-d)_+] = \sum_{s > d}(s - d)\cdot P(S = s)

Kapan Digunakan

  • Soal meminta E[S]E[S] dan/atau Var(S)\text{Var}(S) dari model kolektif atau individual.
  • Soal menyebutkan stop-loss, excess of loss, atau reinsurance retention → hitung E[(Sd)+]E[(S-d)_+].
  • Soal meminta limited expected value E[min(S,d)]E[\min(S,d)] → gunakan dekomposisi.
  • Sebagai input untuk aproksimasi Normal/Lognormal di 4.4 Aggregate Distribution Approximation.
  • Sebagai komponen dalam soal 4.6 Coverage Modifications on Aggregate Models yang melibatkan deductible agregat.

Kapan TIDAK Boleh Digunakan

  • Formula model kolektif tidak berlaku untuk model individual dengan polis berbeda-beda — gunakan μi\sum \mu_i dan σi2\sum \sigma_i^2 untuk model individual.
  • Shortcut Poisson Var(S)=λE[X2]\text{Var}(S) = \lambda E[X^2] hanya untuk frekuensi Poisson — jangan gunakan untuk Binomial atau NegBin.
  • Formula stop-loss kontinu d[1FS]ds\int_d^\infty [1-F_S]\,ds tidak langsung berlaku untuk distribusi SS yang campuran (ada massa di nol) — perlu diperlakukan dengan hati-hati.

Quick Decision Tree

graph TD
    A["Soal model agregat S"] --> B{"Jenis model?"}
    B --> |"Kolektif: S = X1+...+XN"| C{"Tipe frekuensi N?"}
    B --> |"Individual: S = L1+...+Lm"| D["E[S] = sum mu_i<br>Var(S) = sum sigma_i^2"]
    C --> |"Poisson (lambda)"| E["E[S] = lambda*E[X]<br>Var(S) = lambda*E[X^2]<br>(shortcut)"]
    C --> |"NegBin atau lainnya"| F["E[S] = E[N]*E[X]<br>Var(S) = E[N]*Var(X)<br>+ Var(N)*(E[X])^2"]
    E --> G{"Ada stop-loss?"}
    F --> G
    D --> G
    G --> |"Ya, deductible d"| H{"Distribusi S diketahui?"}
    G --> |"Tidak"| I["Selesai: laporkan<br>E[S] dan Var(S)"]
    H --> |"Ya (diskrit/analitik)"| J["Hitung eksak:<br>sum (s-d)*P(S=s) untuk s > d"]
    H --> |"Tidak, gunakan aproksimasi"| K["Gunakan Normal/Lognormal<br>lihat Topik 4.4"]
    J --> L["Cek: E[min(S,d)] + E[(S-d)+] = E[S]"]
    K --> L

Follow-up Options
  1. “Berikan contoh soal stop-loss dengan distribusi SS yang dihitung via Panjer rekursif”
  2. “Jelaskan hubungan 4.3 Mean Variance and Stop-Loss dengan 4.4 Aggregate Distribution Approximation — kapan aproksimasi Normal vs Lognormal lebih baik?”
  3. “Buat flashcard 1-halaman untuk topik ini”

📖 Ref: Klugman, Panjer & Willmot (2019), Loss Models 5th ed., Bab 9; Tse (2009), Bab 3 | 🗓️ 2026-04-17 | #TA2 #ModelAgregat #StopLoss