AktuNotes
← Kembali
TA2 · Materi 7.2

Bühlmann and Bühlmann-Straub Models

Calculation-Intensive Bobot: 20–25% (bersama Topik 7) Klugman, Panjer & Willmot (2019), Loss Models 5th ed., Bab 16 & 17; Tse (2009), Bab 6 & 7
TA2KredibilitasBuhlmannBuhlmannStraubTeoriRisiko

📊 7.2 — Bühlmann and Bühlmann-Straub Models

Ringkasan Cepat

Topik: Model Bühlmann & Bühlmann-Straub | Bobot: ~20–25% (Topik 7) | Difficulty: Calculation-Intensive Ref: Klugman et al. (2019), Bab 16 & 17; Tse (2009), Bab 6 & 7 | Prereq: 7.1 Classical Credibility, 6.3 Bayesian Parameter Estimation


Section 0 — Pemetaan Topik

Topik TA2Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Teori Kredibilitas7.2Menghitung parameter struktural μ\mu, vv, aa; faktor kredibilitas ZZ; premi kredibilitas ZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu; ekstensi Bühlmann-Straub dengan bobot mim_i20–25% (bersama Topik 7)Calculation-Intensive7.1 Classical Credibility, 6.3 Bayesian Parameter Estimation7.3 Bayesian Credibility, 7.4 Empirical Bayesian MethodsKlugman et al. (2019), Bab 16 & 17; Tse (2009), Bab 6 & 7

Section 1 — Intuisi

Bayangkan sebuah perusahaan asuransi kendaraan bermotor memiliki ribuan nasabah, tapi setiap nasabah individu hanya memiliki beberapa tahun riwayat klaim. Data per nasabah terlalu sedikit untuk memperkirakan risiko mereka secara akurat, namun mengabaikan data individual dan memakai rata-rata portofolio penuh juga terlalu kasar — seorang pengemudi yang tidak pernah klaim selama lima tahun tentu bukan risiko yang sama dengan rata-rata populasi. Model Bühlmann menjawab pertanyaan ini: seberapa besar kepercayaan yang seharusnya kita berikan pada data individu, versus kepercayaan pada rata-rata kelompok besar?

Ide dasarnya elegan: premi terbaik untuk seorang nasabah adalah rata-rata tertimbang antara pengalaman klaimnya sendiri dan rata-rata kolektif seluruh portofolio. Bobot yang diberikan pada data individu disebut faktor kredibilitas ZZ — nilainya antara 0 dan 1. Jika nasabah punya banyak data historis yang relevan, ZZ mendekati 1 dan kita hampir sepenuhnya mempercayai data individunya. Jika datanya sangat sedikit, ZZ mendekati 0 dan kita kembali mengandalkan rata-rata kolompok. Besar kecilnya ZZ ditentukan oleh tiga besaran struktural: rata-rata keseluruhan μ\mu, variabilitas dalam satu risiko (vv, process variance), dan variabilitas antar risiko (aa, variance of hypothetical means).

Model Bühlmann-Straub memperluas ide ini untuk situasi yang lebih realistis: tiap nasabah tidak hanya memiliki jumlah periode observasi yang berbeda, tapi juga volume eksposur yang berbeda — misalnya jumlah kendaraan yang diasuransikan, atau jumlah karyawan dalam polis asuransi jiwa kelompok. Nasabah korporat yang mengasuransikan 500 kendaraan seharusnya mendapat bobot lebih besar daripada nasabah individu dengan satu kendaraan, meskipun keduanya hanya punya dua tahun riwayat klaim. Model Bühlmann-Straub menangani ini dengan parameter bobot mijm_{ij}, menghasilkan formula kredibilitas yang lebih umum dan lebih kuat dalam praktik aktuaria nyata.


Section 2 — Definisi Formal

Definisi Matematis — Premi Kredibilitas Bühlmann

Misalkan terdapat rr risiko independen, dan risiko ke-ii memiliki nn periode observasi Xi1,Xi2,,XinX_{i1}, X_{i2}, \ldots, X_{in}. Premi kredibilitas Bühlmann untuk risiko ke-ii adalah:

μ^i=ZiXˉi+(1Zi)μ\hat{\mu}_i = Z_i \bar{X}_i + (1 - Z_i)\mu

di mana faktor kredibilitas Zi=nn+kZ_i = \dfrac{n}{n + k}, dengan k=vak = \dfrac{v}{a}.

Tabel Variabel & Parameter

SimbolMaknaCatatan
rrJumlah risiko (grup/nasabah)Baris dalam data panel
nnJumlah periode observasi per risikoKolom dalam data panel (Bühlmann: sama untuk semua ii)
XijX_{ij}Observasi risiko ke-ii periode ke-jjVariabel acak; bisa klaim per unit eksposur
Xˉi\bar{X}_iRata-rata sampel risiko ke-iiXˉi=1nj=1nXij\bar{X}_i = \frac{1}{n}\sum_{j=1}^n X_{ij}
μ\muExpected value of hypothetical meansμ=E[μ(Θ)]=E[Xij]\mu = E[\mu(\Theta)] = E[X_{ij}]; rata-rata keseluruhan
vvExpected value of process variancesv=E[σ2(Θ)]=E[Var(XijΘ)]v = E[\sigma^2(\Theta)] = E[\text{Var}(X_{ij}\|\Theta)]
aaVariance of hypothetical meansa=Var(μ(Θ))=Var(E[XijΘ])a = \text{Var}(\mu(\Theta)) = \text{Var}(E[X_{ij}\|\Theta])
kkRasio struktural Bühlmannk=v/ak = v/a; semakin kecil kk, semakin besar kredibilitas
ZiZ_iFaktor kredibilitasZi=n/(n+k)(0,1)Z_i = n/(n+k) \in (0,1)
mijm_{ij}Bobot eksposur (Bühlmann-Straub)Jumlah unit/eksposur risiko ke-ii periode ke-jj
mim_iTotal bobot risiko ke-iimi=j=1nimijm_i = \sum_{j=1}^{n_i} m_{ij} (Bühlmann-Straub)

Rumus Utama

1. Parameter Struktural — Tiga Besaran Wajib Dihitung:

μ=E[Xij],v=E[Var(XijΘ)],a=Var(E[XijΘ])\mu = E[X_{ij}], \qquad v = E[\text{Var}(X_{ij}\mid\Theta)], \qquad a = \text{Var}(E[X_{ij}\mid\Theta])

Label: Ketiganya diturunkan dari distribusi campuran. Hubungan penting: Var(Xij)=v+a\text{Var}(X_{ij}) = v + a (hukum variansi total).

2. Faktor Kredibilitas — Model Bühlmann:

Zi=nn+k,k=vaZ_i = \frac{n}{n + k}, \qquad k = \frac{v}{a}

Label: ZiZ_i sama untuk semua risiko jika jumlah periode nn sama. Z1Z \to 1 ketika nn \to \infty atau k0k \to 0.

3. Premi Kredibilitas — Model Bühlmann:

μ^i=ZiXˉi+(1Zi)μ\hat{\mu}_i = Z_i \bar{X}_i + (1 - Z_i)\mu

Label: Rata-rata tertimbang antara pengalaman individu Xˉi\bar{X}_i dan rata-rata kolektif μ\mu.

4. Faktor Kredibilitas — Model Bühlmann-Straub:

Zi=mimi+k,k=va,mi=j=1nimijZ_i = \frac{m_i}{m_i + k}, \qquad k = \frac{v}{a}, \qquad m_i = \sum_{j=1}^{n_i} m_{ij}

Label: Ekstensi Bühlmann — nn digantikan oleh total bobot eksposur mim_i. Risiko dengan eksposur lebih besar mendapat kredibilitas lebih tinggi.

5. Premi Kredibilitas — Model Bühlmann-Straub:

μ^i=ZiXˉiw+(1Zi)μ,Xˉiw=j=1nimijXijmi\hat{\mu}_i = Z_i \bar{X}_i^w + (1 - Z_i)\mu, \qquad \bar{X}_i^w = \frac{\sum_{j=1}^{n_i} m_{ij} X_{ij}}{m_i}

Label: Xˉiw\bar{X}_i^w adalah rata-rata tertimbang berbobot eksposur. Berbeda dari rata-rata aritmetika biasa.

6. Estimasi Parameter Struktural dari Data (Bühlmann):

μ^=Xˉ=1ri=1rXˉi\hat{\mu} = \bar{X} = \frac{1}{r}\sum_{i=1}^r \bar{X}_i v^=1r(n1)i=1rj=1n(XijXˉi)2\hat{v} = \frac{1}{r(n-1)} \sum_{i=1}^r \sum_{j=1}^n (X_{ij} - \bar{X}_i)^2 a^=1r1i=1r(XˉiXˉ)2v^n\hat{a} = \frac{1}{r-1}\sum_{i=1}^r (\bar{X}_i - \bar{X})^2 - \frac{\hat{v}}{n}

Label: Estimator ini adalah ANOVA-style. Perhatikan: a^\hat{a} bisa negatif secara teoritis — jika demikian, tetapkan a^=0\hat{a} = 0 (kredibilitas nol, gunakan μ\mu).

Asumsi Eksplisit

  1. Independensi kondisional: Diberikan Θi\Theta_i, observasi Xi1,,XinX_{i1}, \ldots, X_{in} saling independen dan identik terdistribusi — disebut conditional i.i.d.
  2. Independensi antar risiko: Pasangan (Θi,Xi)(\Theta_i, X_{i\cdot}) independen satu sama lain untuk iii \neq i' — tiap risiko diperlakukan sebagai unit independen.
  3. Homogenitas struktur: Semua risiko berasal dari populasi yang sama — Θi\Theta_i berdistribusi sama (i.i.d. dari distribusi prior yang sama), walaupun nilai spesifik θi\theta_i berbeda.
  4. Keberadaan momen: μ(θ)=E[XΘ=θ]\mu(\theta) = E[X\mid\Theta=\theta] dan σ2(θ)=Var(XΘ=θ)\sigma^2(\theta) = \text{Var}(X\mid\Theta=\theta) keduanya terdefinisi hingga; v=E[σ2(Θ)]>0v = E[\sigma^2(\Theta)] > 0 dan a=Var(μ(Θ))0a = \text{Var}(\mu(\Theta)) \geq 0.
  5. Bühlmann-Straub tambahan: Var(XijΘi)=σ2(Θi)/mij\text{Var}(X_{ij}\mid\Theta_i) = \sigma^2(\Theta_i)/m_{ij} — variansi berbanding terbalik dengan bobot eksposur (asumsi homoskedastisitas berbasis eksposur).

Section 3 — Jembatan Logika

Dari Definisi ke Rumus — Mengapa Formula Tertimbang Ini Optimal

Premi kredibilitas Bühlmann bukan hanya konvensi — ia adalah estimator Best Linear Unbiased Predictor (BLUP) dari μ(Θi)\mu(\Theta_i). Artinya, di antara semua estimator berbentuk c0+jcjXijc_0 + \sum_j c_j X_{ij} (kombinasi linear dari data), premi Bühlmann meminimalkan expected squared error E[(μ^iμ(Θi))2]E[(\hat{\mu}_i - \mu(\Theta_i))^2]. Ini memberikan justifikasi matematis yang kuat: kita tidak memilih formula tertimbang karena “terlihat masuk akal” — ia memang optimal dalam kelas estimator linear.

Tiga Parameter Struktural — Cara Berpikir yang Benar
  • μ\mu adalah rata-rata “grand mean” seluruh portofolio — ini yang kita gunakan jika tidak ada data sama sekali.
  • vv (process variance) mengukur seberapa bervariasi klaim dalam satu risiko yang sama antar periode. Ini adalah keacakan inheren — bahkan risiko yang identik pun akan menghasilkan klaim berbeda tiap tahun.
  • aa (variance of hypothetical means) mengukur seberapa berbeda rata-rata risiko antar risiko. Jika semua nasabah pada dasarnya sama, a0a \approx 0 dan kredibilitas data individu sangat rendah. Jika ada heterogenitas besar, aa besar dan data individu sangat informatif.

Rasio k=v/ak = v/a adalah signal-to-noise ratio terbalik: vv adalah noise (variasi acak), aa adalah signal (perbedaan antar risiko). kk kecil → signal kuat → beri bobot besar pada data individu.

Derivasi Step-by-Step: Faktor Kredibilitas dari BLUP

Cari estimator berbentuk μ^i=α+βXˉi\hat{\mu}_i = \alpha + \beta \bar{X}_i yang meminimalkan:

E[(α+βXˉiμ(Θi))2]E\left[(\alpha + \beta \bar{X}_i - \mu(\Theta_i))^2\right]

Langkah 1 — Hitung E[Xˉi]E[\bar{X}_i] dan Var(Xˉi)\text{Var}(\bar{X}_i):

E[Xˉi]=E[E[XˉiΘi]]=E[μ(Θi)]=μE[\bar{X}_i] = E[E[\bar{X}_i \mid \Theta_i]] = E[\mu(\Theta_i)] = \mu Var(Xˉi)=E[Var(XˉiΘi)]+Var(E[XˉiΘi])=vn+a\text{Var}(\bar{X}_i) = E[\text{Var}(\bar{X}_i \mid \Theta_i)] + \text{Var}(E[\bar{X}_i \mid \Theta_i]) = \frac{v}{n} + a

Langkah 2 — Hitung Cov(Xˉi,μ(Θi))\text{Cov}(\bar{X}_i, \mu(\Theta_i)):

Cov(Xˉi,μ(Θi))=Var(μ(Θi))=a\text{Cov}(\bar{X}_i, \mu(\Theta_i)) = \text{Var}(\mu(\Theta_i)) = a

(karena Xˉi=μ(Θi)+noise\bar{X}_i = \mu(\Theta_i) + \text{noise} dan noise independen dari μ(Θi)\mu(\Theta_i))

Langkah 3 — Terapkan regresi linear optimal (proyeksi orthogonal):

β=Cov(Xˉi,μ(Θi))Var(Xˉi)=aa+v/n=anan+v=nn+v/a=nn+k=Z\beta^* = \frac{\text{Cov}(\bar{X}_i, \mu(\Theta_i))}{\text{Var}(\bar{X}_i)} = \frac{a}{a + v/n} = \frac{an}{an + v} = \frac{n}{n + v/a} = \frac{n}{n+k} = Z

Langkah 4 — Hitung intercept:

α=μβμ=(1Z)μ\alpha^* = \mu - \beta^* \mu = (1 - Z)\mu

Langkah 5 — Gabungkan:

μ^i=(1Z)μ+ZXˉi\hat{\mu}_i = (1-Z)\mu + Z\bar{X}_i

Inilah formula premi kredibilitas Bühlmann — bukan asumsi, melainkan hasil optimasi.

Tiga Larangan Fatal dalam Bühlmann/Bühlmann-Straub
  1. JANGAN tukar posisi vv dan aa saat menghitung kk: k=v/ak = v/a, bukan a/va/v. Ingat: kk besar → ZZ kecil → data individu kurang dipercaya. Jika aa besar (banyak heterogenitas), kk harus kecil → data individu lebih dipercaya. Ini konsisten dengan k=v/ak = v/a.
  2. JANGAN gunakan rata-rata aritmetika biasa Xˉi=1nXij\bar{X}_i = \frac{1}{n}\sum X_{ij} di model Bühlmann-Straub — gunakan rata-rata tertimbang eksposur Xˉiw=mijXij/mi\bar{X}_i^w = \sum m_{ij} X_{ij} / m_i.
  3. JANGAN panik jika a^<0\hat{a} < 0 dalam estimasi dari data — ini bisa terjadi secara komputasi. Dalam kasus ini, tetapkan a^=0\hat{a} = 0, yang berarti Z=0Z = 0 dan premi kredibilitas =μ^= \hat{\mu} (abaikan data individu sepenuhnya).

Section 4 — Contoh Soal

Soal A — Fundamental

Soal: Sebuah portofolio asuransi memiliki tiga kelompok risiko. Data klaim per periode disajikan berikut:

RisikoTahun 1Tahun 2Tahun 3
A468
B10128
C756

Diketahui parameter struktural: μ=7\mu = 7, v=4v = 4, a=4a = 4. Hitung faktor kredibilitas dan premi kredibilitas untuk setiap risiko menggunakan model Bühlmann.

Solusi Soal A

Pendekatan: Hitung kk, lalu ZZ, lalu premi kredibilitas untuk masing-masing risiko.

1. Identifikasi Variabel

  • r=3r = 3 risiko, n=3n = 3 periode (sama untuk semua risiko)
  • μ=7\mu = 7, v=4v = 4, a=4a = 4
  • Rata-rata sampel: XˉA=(4+6+8)/3=6\bar{X}_A = (4+6+8)/3 = 6, XˉB=(10+12+8)/3=10\bar{X}_B = (10+12+8)/3 = 10, XˉC=(7+5+6)/3=6\bar{X}_C = (7+5+6)/3 = 6

2. Identifikasi Model Model Bühlmann — jumlah periode sama untuk semua risiko (n=3n = 3), tanpa bobot eksposur berbeda.

3. Setup Persamaan

k=va=44=1,Z=nn+k=33+1=0.75k = \frac{v}{a} = \frac{4}{4} = 1, \qquad Z = \frac{n}{n+k} = \frac{3}{3+1} = 0.75

4. Eksekusi Aljabar

Premi kredibilitas untuk masing-masing risiko:

μ^A=0.75×6+0.25×7=4.5+1.75=6.25\hat{\mu}_A = 0.75 \times 6 + 0.25 \times 7 = 4.5 + 1.75 = 6.25 μ^B=0.75×10+0.25×7=7.5+1.75=9.25\hat{\mu}_B = 0.75 \times 10 + 0.25 \times 7 = 7.5 + 1.75 = 9.25 μ^C=0.75×6+0.25×7=4.5+1.75=6.25\hat{\mu}_C = 0.75 \times 6 + 0.25 \times 7 = 4.5 + 1.75 = 6.25

5. Verification Setiap premi kredibilitas berada di antara Xˉi\bar{X}_i dan μ=7\mu = 7. Risiko A dan C keduanya memiliki Xˉ=6<7\bar{X} = 6 < 7, sehingga premi mereka ditarik ke atas mendekati 7. Risiko B dengan Xˉ=10>7\bar{X} = 10 > 7 ditarik ke bawah mendekati 7. Rata-rata tertimbang dari ketiga premi: (6.25+9.25+6.25)/3=7.257(6.25 + 9.25 + 6.25)/3 = 7.25 \neq 7 — ini wajar karena rata-rata sampel tidak persis 7.

Hasil: Z=0.75Z = 0.75; μ^A=6.25\hat{\mu}_A = 6.25, μ^B=9.25\hat{\mu}_B = 9.25, μ^C=6.25\hat{\mu}_C = 6.25.

Exam Tips — Soal A

Target waktu: 3 menit. Common trap: Menggunakan k=a/vk = a/v (terbalik). Shortcut: Hafal k=v/ak = v/a; jika v=av = a, maka k=1k = 1 dan Z=n/(n+1)Z = n/(n+1) — sangat mudah dihitung.


Soal B — Exam-Typical

Soal: Sebuah portofolio memiliki r=4r = 4 risiko, masing-masing diamati selama n=5n = 5 tahun. Data klaim tahunan:

RisikoT1T2T3T4T5Xˉi\bar{X}_i
1243513.0
28697108.0
3546555.0
4354444.0

Estimasi parameter struktural μ^\hat{\mu}, v^\hat{v}, a^\hat{a} dari data, lalu hitung premi kredibilitas untuk setiap risiko.

Solusi Soal B

Pendekatan: Estimasi tiga parameter struktural menggunakan rumus ANOVA, lalu terapkan formula premi Bühlmann.

1. Identifikasi Variabel

  • r=4r = 4, n=5n = 5
  • Xˉ1=3.0\bar{X}_1 = 3.0, Xˉ2=8.0\bar{X}_2 = 8.0, Xˉ3=5.0\bar{X}_3 = 5.0, Xˉ4=4.0\bar{X}_4 = 4.0
  • Xˉ=(3+8+5+4)/4=20/4=5.0\bar{X} = (3+8+5+4)/4 = 20/4 = 5.0

2. Identifikasi Model Model Bühlmann dengan estimasi parameter dari data. Gunakan rumus ANOVA tiga-langkah.

3. Setup Persamaan — Tiga Estimator Struktural

μ^=Xˉ=5.0\hat{\mu} = \bar{X} = 5.0 v^=1r(n1)i=1rj=1n(XijXˉi)2\hat{v} = \frac{1}{r(n-1)}\sum_{i=1}^r \sum_{j=1}^n (X_{ij} - \bar{X}_i)^2 a^=1r1i=1r(XˉiXˉ)2v^n\hat{a} = \frac{1}{r-1}\sum_{i=1}^r (\bar{X}_i - \bar{X})^2 - \frac{\hat{v}}{n}

4. Eksekusi Aljabar

Hitung within-group sum of squares untuk v^\hat{v}:

Risiko 1: (23)2+(43)2+(33)2+(53)2+(13)2=1+1+0+4+4=10(2-3)^2+(4-3)^2+(3-3)^2+(5-3)^2+(1-3)^2 = 1+1+0+4+4 = 10

Risiko 2: (88)2+(68)2+(98)2+(78)2+(108)2=0+4+1+1+4=10(8-8)^2+(6-8)^2+(9-8)^2+(7-8)^2+(10-8)^2 = 0+4+1+1+4 = 10

Risiko 3: (55)2+(45)2+(65)2+(55)2+(55)2=0+1+1+0+0=2(5-5)^2+(4-5)^2+(6-5)^2+(5-5)^2+(5-5)^2 = 0+1+1+0+0 = 2

Risiko 4: (34)2+(54)2+(44)2+(44)2+(44)2=1+1+0+0+0=2(3-4)^2+(5-4)^2+(4-4)^2+(4-4)^2+(4-4)^2 = 1+1+0+0+0 = 2

Total within-SS =10+10+2+2=24= 10+10+2+2 = 24

v^=244×4=2416=1.5\hat{v} = \frac{24}{4 \times 4} = \frac{24}{16} = 1.5

Hitung between-group sum of squares untuk a^\hat{a}:

i=14(XˉiXˉ)2=(35)2+(85)2+(55)2+(45)2=4+9+0+1=14\sum_{i=1}^4 (\bar{X}_i - \bar{X})^2 = (3-5)^2 + (8-5)^2 + (5-5)^2 + (4-5)^2 = 4+9+0+1 = 14 a^=14411.55=1430.3=4.6670.3=4.367\hat{a} = \frac{14}{4-1} - \frac{1.5}{5} = \frac{14}{3} - 0.3 = 4.667 - 0.3 = 4.367

Hitung kk dan ZZ:

k=v^a^=1.54.367=0.3436,Z=55+0.3436=55.3436=0.9357k = \frac{\hat{v}}{\hat{a}} = \frac{1.5}{4.367} = 0.3436, \qquad Z = \frac{5}{5 + 0.3436} = \frac{5}{5.3436} = 0.9357

Premi kredibilitas:

μ^1=0.9357×3.0+0.0643×5.0=2.807+0.321=3.13\hat{\mu}_1 = 0.9357 \times 3.0 + 0.0643 \times 5.0 = 2.807 + 0.321 = 3.13 μ^2=0.9357×8.0+0.0643×5.0=7.486+0.321=7.81\hat{\mu}_2 = 0.9357 \times 8.0 + 0.0643 \times 5.0 = 7.486 + 0.321 = 7.81 μ^3=0.9357×5.0+0.0643×5.0=4.679+0.321=5.00\hat{\mu}_3 = 0.9357 \times 5.0 + 0.0643 \times 5.0 = 4.679 + 0.321 = 5.00 μ^4=0.9357×4.0+0.0643×5.0=3.743+0.321=4.06\hat{\mu}_4 = 0.9357 \times 4.0 + 0.0643 \times 5.0 = 3.743 + 0.321 = 4.06

5. Verification Z=0.936Z = 0.936 tinggi — masuk akal karena n=5n=5 cukup besar dan aa jauh lebih besar dari vv (banyak heterogenitas antar risiko). Risiko 3 memiliki Xˉ3=μ^\bar{X}_3 = \hat{\mu}, sehingga premi kredibilitas-nya tepat sama dengan grand mean.

Hasil: μ^=5\hat{\mu}=5, v^=1.5\hat{v}=1.5, a^=4.367\hat{a}=4.367, Z=0.936Z=0.936; premi: 3.13, 7.81, 5.00, 4.06.

Exam Tips — Soal B

Target waktu: 5–6 menit. Common trap: (1) Lupa mengurangi v^/n\hat{v}/n saat menghitung a^\hat{a} — ini koreksi bias yang sering terlewat. (2) Membagi between-SS dengan rr alih-alih r1r-1. Shortcut: Kerjakan dalam urutan ketat: within-SS → v^\hat{v} → between-SS → a^\hat{a}kkZZ → premi. Jangan lewati langkah.


Soal C — Challenging

Soal: Model Bühlmann-Straub. Tiga risiko diamati selama 3 tahun dengan eksposur (jumlah polis) berbeda:

Risiko(mi1,Xi1)(m_{i1}, X_{i1})(mi2,Xi2)(m_{i2}, X_{i2})(mi3,Xi3)(m_{i3}, X_{i3})
1(100, 0.04)(150, 0.06)(200, 0.05)
2(200, 0.08)(200, 0.07)(400, 0.09)
3(50, 0.03)(100, 0.05)(50, 0.04)

XijX_{ij} adalah klaim per polis per tahun. Diketahui v=0.0005v = 0.0005, a=0.0004a = 0.0004. Hitung premi kredibilitas Bühlmann-Straub untuk setiap risiko.

Solusi Soal C

Pendekatan: Hitung total eksposur mim_i dan rata-rata tertimbang Xˉiw\bar{X}_i^w untuk setiap risiko, lalu gunakan formula Bühlmann-Straub.

1. Identifikasi Variabel

  • v=0.0005v = 0.0005, a=0.0004a = 0.0004, k=v/a=0.0005/0.0004=1.25k = v/a = 0.0005/0.0004 = 1.25
  • Eksposur dan klaim per periode diberikan untuk r=3r = 3 risiko, ni=3n_i = 3 periode

2. Identifikasi Model Bühlmann-Straub — eksposur berbeda antar periode dan antar risiko. Gunakan rata-rata tertimbang eksposur.

3. Setup Persamaan

mi=j=1nimij,Xˉiw=jmijXijmi,Zi=mimi+km_i = \sum_{j=1}^{n_i} m_{ij}, \quad \bar{X}_i^w = \frac{\sum_j m_{ij} X_{ij}}{m_i}, \quad Z_i = \frac{m_i}{m_i + k}

4. Eksekusi Aljabar

Risiko 1:

m1=100+150+200=450m_1 = 100 + 150 + 200 = 450 Xˉ1w=100(0.04)+150(0.06)+200(0.05)450=4+9+10450=23450=0.05111\bar{X}_1^w = \frac{100(0.04) + 150(0.06) + 200(0.05)}{450} = \frac{4 + 9 + 10}{450} = \frac{23}{450} = 0.05111 Z1=450450+1.25=450451.25=0.99723Z_1 = \frac{450}{450 + 1.25} = \frac{450}{451.25} = 0.99723

Risiko 2:

m2=200+200+400=800m_2 = 200 + 200 + 400 = 800 Xˉ2w=200(0.08)+200(0.07)+400(0.09)800=16+14+36800=66800=0.08250\bar{X}_2^w = \frac{200(0.08) + 200(0.07) + 400(0.09)}{800} = \frac{16 + 14 + 36}{800} = \frac{66}{800} = 0.08250 Z2=800800+1.25=800801.25=0.99844Z_2 = \frac{800}{800 + 1.25} = \frac{800}{801.25} = 0.99844

Risiko 3:

m3=50+100+50=200m_3 = 50 + 100 + 50 = 200 Xˉ3w=50(0.03)+100(0.05)+50(0.04)200=1.5+5+2200=8.5200=0.04250\bar{X}_3^w = \frac{50(0.03) + 100(0.05) + 50(0.04)}{200} = \frac{1.5 + 5 + 2}{200} = \frac{8.5}{200} = 0.04250 Z3=200200+1.25=200201.25=0.99379Z_3 = \frac{200}{200 + 1.25} = \frac{200}{201.25} = 0.99379

Estimasi μ\mu (grand mean tertimbang):

μ^=imiXˉiwimi=450(0.05111)+800(0.08250)+200(0.04250)1450\hat{\mu} = \frac{\sum_i m_i \bar{X}_i^w}{\sum_i m_i} = \frac{450(0.05111) + 800(0.08250) + 200(0.04250)}{1450} =23+66+8.51450=97.51450=0.06724= \frac{23 + 66 + 8.5}{1450} = \frac{97.5}{1450} = 0.06724

Premi kredibilitas:

μ^1=0.99723×0.05111+0.00277×0.06724=0.05097+0.00019=0.05116\hat{\mu}_1 = 0.99723 \times 0.05111 + 0.00277 \times 0.06724 = 0.05097 + 0.00019 = 0.05116 μ^2=0.99844×0.08250+0.00156×0.06724=0.08237+0.00010=0.08247\hat{\mu}_2 = 0.99844 \times 0.08250 + 0.00156 \times 0.06724 = 0.08237 + 0.00010 = 0.08247 μ^3=0.99379×0.04250+0.00621×0.06724=0.04223+0.00042=0.04265\hat{\mu}_3 = 0.99379 \times 0.04250 + 0.00621 \times 0.06724 = 0.04223 + 0.00042 = 0.04265

5. Verification ZZ sangat tinggi untuk semua risiko (>0.993) karena mik=1.25m_i \gg k = 1.25 — volume eksposur besar membuat data individu sangat dominan. Premi risiko 2 paling tinggi (konsisten dengan Xˉ2w=0.0825\bar{X}_2^w = 0.0825, tertinggi). Semua premi berada di antara Xˉiw\bar{X}_i^w dan μ^=0.0672\hat{\mu} = 0.0672 — terverifikasi.

Hasil: μ^1=0.0512\hat{\mu}_1 = 0.0512, μ^2=0.0825\hat{\mu}_2 = 0.0825, μ^3=0.0427\hat{\mu}_3 = 0.0427 (klaim per polis per tahun).

Exam Tips — Soal C

Target waktu: 6–7 menit. Common trap: (1) Menggunakan rata-rata aritmetika biasa Xˉi=(Xi1+Xi2+Xi3)/3\bar{X}_i = (X_{i1}+X_{i2}+X_{i3})/3 alih-alih rata-rata tertimbang eksposur — ini kesalahan paling kritis di Bühlmann-Straub. (2) Lupa menghitung grand mean μ^\hat{\mu} secara tertimbang. Shortcut: Hitung total klaim absolut jmijXij\sum_j m_{ij} X_{ij} terlebih dahulu, baru bagi dengan mim_i — lebih aman daripada langsung menggunakan desimal kecil.


Section 5 — Verifikasi & Sanity Check

Cross-check Faktor Kredibilitas — Monoton dalam nn (atau mim_i)

Faktor kredibilitas ZZ harus:

Z=nn+k(0,1),Z ketika n,Z1 ketika nZ = \frac{n}{n+k} \in (0, 1), \quad Z \nearrow \text{ ketika } n \nearrow, \quad Z \to 1 \text{ ketika } n \to \infty

Jika dua risiko memiliki kk yang sama tapi periode berbeda, risiko dengan nn lebih besar pasti memiliki ZZ lebih besar. Jika hasil menunjukkan sebaliknya, ada kesalahan.

Cross-check Premi Kredibilitas — Selalu di Antara Xˉi\bar{X}_i dan μ\mu

Untuk setiap risiko ii:

min(Xˉi,μ)μ^imax(Xˉi,μ)\min(\bar{X}_i, \mu) \leq \hat{\mu}_i \leq \max(\bar{X}_i, \mu)

Karena μ^i=ZXˉi+(1Z)μ\hat{\mu}_i = Z\bar{X}_i + (1-Z)\mu dengan Z(0,1)Z \in (0,1), premi kredibilitas selalu berada di antara dua ujung. Ini dapat dicek secara visual tanpa kalkulator.

Cross-check a^\hat{a} — Relasi dengan Variansi Total

Berlaku: Var(Xij)=v+a\text{Var}(X_{ij}) = v + a. Setelah menghitung v^\hat{v} dan a^\hat{a}, verifikasi bahwa estimasi variansi total dari data (jika bisa dihitung) konsisten dengan v^+a^\hat{v} + \hat{a}. Juga, jika a^<0\hat{a} < 0 secara komputasi → tetapkan a^=0\hat{a} = 0 dan Z=0Z = 0.

Metode Alternatif

Hubungan dengan pendekatan Bayesian: Dalam kasus khusus distribusi konjugat (misalnya Poisson-Gamma), premi kredibilitas Bühlmann identik dengan Bayes estimator posterior mean. Lihat 7.3 Bayesian Credibility untuk bukti formal. Ini berarti dalam kasus konjugat, kedua pendekatan memberikan jawaban yang sama — verifikasi silang yang berguna.


Section 6 — Visualisasi Mental

Visualisasi Bühlmann — Tarik-menarik antara Data dan Prior:

                  μ = 7 (grand mean / prior)

RISIKO B          │←────── tarik ke μ ───────→│ B̄ = 10
(data tinggi)     │         (1-Z)=0.25         │
                  │                             │
                  ├─────────────────────────────┤
                  │  ← premi B = 9.25 →         │
                  │         Z=0.75              │

RISIKO A          │←────── tarik ke μ ───────→│ Ā = 6
(data rendah)     │         (1-Z)=0.25         │
                  │                             │
                  ├─────────────────────────────┤
                  │  → premi A = 6.25 ←         │
                         Z=0.75

Visualisasi Bühlmann-Straub — Kredibilitas Berbasis Volume:

Risiko 2 (m=800): ████████████████████████ Z = 0.998  ← hampir 100% percaya data
Risiko 1 (m=450): █████████████           Z = 0.997
Risiko 3 (m=200): ██████                  Z = 0.994

         ↑ semakin besar volume eksposur, semakin tinggi Z

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Formula
Panjang “tarikan” ke μ\muBobot prior (1Z)=k/(n+k)(1-Z) = k/(n+k)
Panjang “tarikan” ke Xˉi\bar{X}_iBobot data Z=n/(n+k)Z = n/(n+k)
Lebar bar eksposur (BS)mim_i — total bobot risiko
kk kecil → bar panjangk=v/ak=v/a kecil berarti signal kuat, ZZ besar

Section 7 — Jebakan Umum

Kesalahan Parametrisasi — kk Terbalik

Salah: k=a/vk = a/v sehingga kk besar ketika heterogenitas antar risiko besar → menghasilkan ZZ kecil padahal seharusnya besar. Benar: k=v/ak = v/a. Mnemonic: “variasi dalam risiko di atas, antar risiko di bawah” — k=vdalamaantark = \frac{v_{\text{dalam}}}{a_{\text{antar}}}.

Cara verifikasi: ketika aa \to \infty (risiko sangat heterogen), k0k \to 0, Z1Z \to 1 (percayai penuh data individu). Ini intuitif dan hanya konsisten jika k=v/ak = v/a.

Kesalahan Konseptual — Empat Miskonsepsi Kritis
  1. ZZ adalah probabilitas.” — Salah. ZZ adalah bobot optimasi, bukan probabilitas. Ia bisa sangat dekat 0 atau 1 tetapi tidak pernah persis 0 atau 1 (kecuali dalam batas).
  2. “Model Bühlmann membutuhkan asumsi distribusi tertentu.” — Salah. Bühlmann hanya memerlukan momen orde dua (mean dan variansi) — distribusi tidak perlu dispesifikasikan. Ini kekuatan besar dibanding MLE.
  3. “Bühlmann-Straub sama dengan Bühlmann hanya dengan nn yang berbeda.” — Kurang tepat. Perbedaannya bukan hanya jumlah periode, tapi volume eksposur per periode yang mempengaruhi variansi kondisional.
  4. a^\hat{a} negatif adalah kesalahan hitung.” — Tidak selalu. Secara matematis a^\hat{a} bisa negatif dari formula estimasi; ini berarti data tidak mendukung adanya heterogenitas antar risiko, sehingga a^=0\hat{a} = 0 dipaksakan.
Kesalahan Interpretasi Soal — Bühlmann vs Bühlmann-Straub
  • Jika soal hanya memberikan nilai klaim tanpa bobot eksposur → gunakan Bühlmann biasa.
  • Jika soal memberikan pasangan (mij,Xij)(m_{ij}, X_{ij}) atau menyebut “jumlah polis”, “jumlah karyawan”, “volume eksposur” → gunakan Bühlmann-Straub.
  • Jika soal menyebut XijX_{ij} adalah “klaim per unit eksposur” dan eksposur berbeda antar periode → wajib Bühlmann-Straub.
  • Jika soal menyebut “estimate the parameters” dari data → gunakan rumus ANOVA estimasi; jangan asumsikan vv dan aa diberikan.
Red Flags — Keyword yang Wajib Memicu Prosedur Khusus
  • “Jumlah polis/karyawan/unit berbeda” → Bühlmann-Straub; hitung mi=mijm_i = \sum m_{ij} dan Xˉiw\bar{X}_i^w.
  • “Estimate vv and aa from the data” → Gunakan rumus ANOVA tiga langkah; ingat koreksi bias v^/n-\hat{v}/n pada a^\hat{a}.
  • "a^<0\hat{a} < 0" → Tetapkan a^=0\hat{a} = 0, premi kredibilitas =μ^= \hat{\mu} untuk semua risiko.
  • “Hubungan dengan Bayesian” → Dalam kasus konjugat, Bühlmann \equiv Bayes estimator; lihat 7.3 Bayesian Credibility.
  • "nn \to \infty"Z1Z \to 1, premi kredibilitas Xˉi\to \bar{X}_i (mendekati MLE individual).

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. Tiga parameter struktural wajib:
μ=E[Xij],v=E[Var(XijΘ)],a=Var(E[XijΘ])\mu = E[X_{ij}], \quad v = E[\text{Var}(X_{ij}\mid\Theta)], \quad a = \text{Var}(E[X_{ij}\mid\Theta])
  1. Faktor kredibilitas Bühlmann:
k=va,Z=nn+kk = \frac{v}{a}, \qquad Z = \frac{n}{n+k}
  1. Premi kredibilitas (berlaku untuk keduanya):
μ^i=ZiXˉi+(1Zi)μ\hat{\mu}_i = Z_i \bar{X}_i + (1-Z_i)\mu
  1. Bühlmann-Straub — ganti nn dengan mim_i dan Xˉi\bar{X}_i dengan Xˉiw\bar{X}_i^w:
mi=jmij,Xˉiw=jmijXijmi,Zi=mimi+km_i = \sum_j m_{ij}, \quad \bar{X}_i^w = \frac{\sum_j m_{ij} X_{ij}}{m_i}, \quad Z_i = \frac{m_i}{m_i + k}
  1. Estimasi a^\hat{a} dari data — jangan lupa koreksi bias:
a^=1r1i(XˉiXˉ)2v^n\hat{a} = \frac{1}{r-1}\sum_i (\bar{X}_i - \bar{X})^2 - \frac{\hat{v}}{n}

Kapan Digunakan

  • Soal menyebut “credibility premium”, “Bühlmann”, atau “Bühlmann-Straub”
  • Soal memberikan data panel multi-risiko multi-periode dan meminta premi
  • Soal memberikan vv, aa, dan μ\mu secara eksplisit → langsung hitung kk, ZZ, premi
  • Soal meminta estimasi vv dan aa dari data → gunakan ANOVA estimator
  • Soal memberikan eksposur berbeda per periode/risiko → Bühlmann-Straub

Kapan TIDAK Boleh Digunakan

  • Hanya satu risiko tanpa pembanding → tidak ada aa yang bisa diestimasi
  • Soal meminta estimasi MLE atau method of moments tanpa konteks kredibilitas → gunakan 6.1 Parameter Estimation Methods
  • Soal meminta distribusi posterior penuh → gunakan 7.3 Bayesian Credibility
  • Data tidak tersusun dalam struktur panel (multi-risiko, multi-periode) → model kredibilitas tidak applicable

Quick Decision Tree

graph TD
    A["Ada data multi-risiko multi-periode?"] -->|"Ya"| B["Apakah eksposur per periode berbeda?"]
    A -->|"Tidak"| Z["Gunakan metode estimasi lain\n(MLE, MoM)"]
    B -->|"Ya — ada bobot m_ij"| C["Model Bühlmann-Straub<br>Hitung m_i dan X-bar-w_i"]
    B -->|"Tidak — semua n sama"| D["Model Bühlmann biasa<br>Hitung X-bar_i aritmetika"]
    C --> E["v dan a diberikan?"]
    D --> E
    E -->|"Ya"| F["Langsung hitung k = v/a<br>lalu Z dan premi"]
    E -->|"Tidak — estimasi dari data"| G["Hitung within-SS -> v-hat<br>Hitung between-SS -> a-hat<br>Koreksi: a-hat = ... - v-hat/n"]
    G --> H["a-hat negatif?"]
    H -->|"Ya"| I["Tetapkan a-hat = 0<br>Z = 0, premi = mu"]
    H -->|"Tidak"| F
    F --> J["Premi = Z * X-bar_i + (1-Z) * mu"]

Follow-up Options
  1. “Berikan contoh soal estimasi vv dan aa dari data panel lengkap Bühlmann-Straub”
  2. “Jelaskan hubungan 7.2 Bühlmann and Bühlmann-Straub Models dengan 7.3 Bayesian Credibility secara matematis”
  3. “Buat flashcard 1-halaman untuk formula Bühlmann vs Bühlmann-Straub side-by-side”

📖 Ref: Klugman, Panjer & Willmot (2019), Loss Models 5th ed., Bab 16 & 17; Tse (2009), Bab 6 & 7 | 🗓️ 2026-04-17 | #TA2 #Kredibilitas #Buhlmann