AktuNotes
← Kembali
TA2 · Materi 7.3

Bayesian Credibility

Hard Bobot: 20–25% (bersama Topik 7) Klugman et al. (2019) Bab 16, 17, 18; Tse (2009) Bab 6, 7, 8, 9 (kecuali 9.4)
TA2KredibilitasBayesianTeoriRisiko

📊 7.3 — Bayesian Credibility

Ringkasan Cepat

Topik: Bayesian Credibility | Bobot: ~20–25% (Topik 7 keseluruhan) | Difficulty: Hard Ref: Klugman et al. (2019) Bab 16–18; Tse (2009) Bab 6–9 | Prereq: 7.2 Bühlmann and Bühlmann-Straub Models, 6.3 Bayesian Parameter Estimation


Section 0 — Pemetaan Topik

Topik TA2Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Teori Kredibilitas7.3Hitung estimator Bayes; posterior mean; hubungan Bayes–Bühlmann20–25% (Topik 7)Hard7.2 Bühlmann and Bühlmann-Straub Models, 6.3 Bayesian Parameter Estimation7.1 Classical Credibility, 7.4 Empirical Bayesian MethodsKlugman et al. (2019) Bab 16, 17, 18; Tse (2009) Bab 6–9

Section 1 — Intuisi

Bayangkan seorang manajer risiko di perusahaan asuransi kendaraan yang harus menetapkan premi untuk seorang nasabah baru. Di satu sisi, ia punya data pengalaman klaim nasabah tersebut selama beberapa tahun. Di sisi lain, ia juga punya keyakinan awal (prior belief) tentang profil risiko nasabah berdasarkan karakteristik populasinya — usia, jenis kendaraan, riwayat SIM. Masalahnya: mana yang harus ia percaya lebih? Data individu yang mungkin masih sedikit, atau keyakinan awal yang lebih stabil tapi mungkin tidak tepat untuk nasabah ini?

Inilah inti dari Bayesian Credibility: sebuah kerangka matematis untuk mengombinasikan informasi prior (keyakinan sebelum data ada) dengan informasi dari data observasi. Pendekatan Bayesian tidak membuang salah satu — ia menemukan bobot optimal antara keduanya secara formal. Hasilnya adalah estimator posterior yang merupakan prediksi terbaik kita tentang kerugian masa depan nasabah tersebut.

Yang membuat pendekatan ini elegan adalah hubungannya dengan model Bühlmann yang telah kita pelajari. Dalam kondisi tertentu — khususnya ketika distribusi prior dan likelihood membentuk conjugate pair — estimator Bayesian menghasilkan formula premi yang persis sama dengan formula kredibilitas Bühlmann ZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu. Artinya, Bühlmann bukan hanya sebuah pendekatan aproksimasi yang “kebetulan berguna” — ia adalah manifestasi alami dari logika Bayesian yang lengkap dalam keluarga distribusi konjugat.


Section 2 — Definisi Formal

Definisi Matematis — Estimator Bayesian

Diberikan parameter risiko Θ\Theta dengan distribusi prior π(θ)\pi(\theta) dan data klaim X=(X1,X2,,Xn)\mathbf{X} = (X_1, X_2, \ldots, X_n), estimator Bayesian dari kerugian masa depan Xn+1X_{n+1} adalah:

μ^Bayes=E[Xn+1X]=E[μ(Θ)X]\hat{\mu}^{Bayes} = E[X_{n+1} \mid \mathbf{X}] = E[\mu(\Theta) \mid \mathbf{X}]

yaitu mean posterior dari μ(Θ)=E[XΘ]\mu(\Theta) = E[X \mid \Theta] setelah mengamati data X\mathbf{X}.


Tabel Variabel & Parameter

SimbolMaknaCatatan
Θ\ThetaParameter risiko (variabel acak)Mewakili profil risiko tak terobservasi
π(θ)\pi(\theta)Distribusi prior atas Θ\ThetaKeyakinan sebelum melihat data
f(xθ)f(x \mid \theta)Fungsi likelihood: distribusi XX given Θ=θ\Theta = \thetaDistribusi klaim untuk profil risiko tetap
π(θx)\pi(\theta \mid \mathbf{x})Distribusi posteriorPrior diperbarui setelah melihat data
μ(θ)\mu(\theta)E[XΘ=θ]E[X \mid \Theta = \theta]Mean hipotetis untuk profil risiko θ\theta
μ\muE[μ(Θ)]=E[X]E[\mu(\Theta)] = E[X]Mean kolektif (grand mean)
v(θ)v(\theta)Var(XΘ=θ)\text{Var}(X \mid \Theta = \theta)Expected process variance untuk θ\theta
vvE[v(Θ)]E[v(\Theta)]Expected value of process variance (EVPV)
aaVar(μ(Θ))\text{Var}(\mu(\Theta))Variance of hypothetical means (VHM)
kkv/av/aRasio Bühlmann
ZZn/(n+k)n/(n+k)Faktor kredibilitas Bühlmann

Rumus Utama

Teorema Bayes untuk distribusi posterior:

π(θx)=f(xθ)π(θ)f(x)f(xθ)π(θ)\pi(\theta \mid \mathbf{x}) = \frac{f(\mathbf{x} \mid \theta)\, \pi(\theta)}{f(\mathbf{x})} \propto f(\mathbf{x} \mid \theta)\, \pi(\theta)

Label: Posterior ∝ Likelihood × Prior


Estimator Bayesian (mean posterior):

μ^Bayes=E[μ(Θ)X=x]=μ(θ)π(θx)dθ\hat{\mu}^{Bayes} = E[\mu(\Theta) \mid \mathbf{X} = \mathbf{x}] = \int \mu(\theta)\, \pi(\theta \mid \mathbf{x})\, d\theta

Label: Rata-rata tertimbang μ(θ)\mu(\theta) menggunakan distribusi posterior sebagai bobot


Hubungan Bühlmann–Bayes (kasus conjugate):

E[μ(Θ)X]=ZXˉ+(1Z)μE[\mu(\Theta) \mid \mathbf{X}] = Z \bar{X} + (1-Z)\mu Z=nn+k,k=vaZ = \frac{n}{n + k}, \quad k = \frac{v}{a}

Label: Ketika prior dan likelihood membentuk pasangan conjugate, estimator Bayesian persis sama dengan premi Bühlmann


Pasangan Conjugate penting dalam aktuaria:

Poisson(λ) likelihood+Gamma(α,β) priorGamma posterior\text{Poisson}(\lambda) \text{ likelihood} + \text{Gamma}(\alpha, \beta) \text{ prior} \Rightarrow \text{Gamma posterior} Normal(μ,σ2) likelihood+Normal(μ0,σ02) priorNormal posterior\text{Normal}(\mu, \sigma^2) \text{ likelihood} + \text{Normal}(\mu_0, \sigma_0^2) \text{ prior} \Rightarrow \text{Normal posterior} Binomial(m,q) likelihood+Beta(α,β) priorBeta posterior\text{Binomial}(m, q) \text{ likelihood} + \text{Beta}(\alpha, \beta) \text{ prior} \Rightarrow \text{Beta posterior}

Label: Conjugate pairs memungkinkan komputasi posterior secara analitik


Asumsi Eksplisit

  1. Kondisional independen: Diberikan Θ=θ\Theta = \theta, variabel X1,X2,,XnX_1, X_2, \ldots, X_n saling independen dan identically distributed (i.i.d.) dengan distribusi f(xθ)f(x \mid \theta).
  2. Prior tertentu: Distribusi prior π(θ)\pi(\theta) telah ditentukan sebelum data diobservasi dan merepresentasikan heterogenitas populasi.
  3. Squared-error loss: Estimator Bayesian optimal di bawah fungsi loss kuadratik adalah mean posterior.
  4. Stasioneritas: Proses yang menghasilkan klaim diasumsikan stasioner — Θ\Theta tidak berubah antar periode untuk individu yang sama.
  5. Kasus conjugate (untuk hubungan Bühlmann): Distribusi prior dan likelihood membentuk pasangan konjugat sehingga posterior berada dalam keluarga yang sama dengan prior.

Section 3 — Jembatan Logika

Dari Distribusi Prior ke Estimator Posterior

Inti Bayesian Credibility adalah pertanyaan: setelah melihat data x\mathbf{x}, apa keyakinan terbaik kita tentang μ(Θ)\mu(\Theta)? Jawabannya adalah mean dari distribusi posterior π(θx)\pi(\theta \mid \mathbf{x}). Distribusi posterior ini lahir dari perkalian prior π(θ)\pi(\theta) dan likelihood f(xθ)f(\mathbf{x} \mid \theta) — Bayes rule. Momen dari posterior inilah yang menjadi estimator optimal kita di bawah squared-error loss.

Mengapa Posterior Mean Optimal?

Di bawah squared-error loss L(μ^,μ(θ))=(μ^μ(θ))2L(\hat{\mu}, \mu(\theta)) = (\hat{\mu} - \mu(\theta))^2, estimator yang meminimalkan expected loss adalah:

μ^=argminμ^E[(μ^μ(Θ))2X]=E[μ(Θ)X]\hat{\mu}^* = \arg\min_{\hat{\mu}} E\left[(\hat{\mu} - \mu(\Theta))^2 \mid \mathbf{X}\right] = E[\mu(\Theta) \mid \mathbf{X}]

Ini adalah hasil standar dari decision theory: estimator Bayes optimal adalah conditional expectation dari kuantitas yang ingin diestimasikan.


Derivasi: Kasus Poisson–Gamma (Paling Umum di Exam)

Setup: Frekuensi klaim XiΘ=θPoisson(θ)X_i \mid \Theta = \theta \sim \text{Poisson}(\theta), dan prior ΘGamma(α,β)\Theta \sim \text{Gamma}(\alpha, \beta) (parametrisasi: mean =αβ= \alpha\beta).

Langkah 1 — Likelihood:

f(xθ)=i=1neθθxixi!enθθxif(\mathbf{x} \mid \theta) = \prod_{i=1}^{n} \frac{e^{-\theta}\theta^{x_i}}{x_i!} \propto e^{-n\theta} \theta^{\sum x_i}

Langkah 2 — Prior:

π(θ)θα1eθ/β\pi(\theta) \propto \theta^{\alpha - 1} e^{-\theta/\beta}

Langkah 3 — Posterior (∝ likelihood × prior):

π(θx)enθθxiθα1eθ/β=θ(α+xi)1exp ⁣(θ ⁣(n+1β))\pi(\theta \mid \mathbf{x}) \propto e^{-n\theta} \theta^{\sum x_i} \cdot \theta^{\alpha - 1} e^{-\theta/\beta} = \theta^{(\alpha + \sum x_i) - 1} \exp\!\left(-\theta\!\left(n + \frac{1}{\beta}\right)\right)

Ini berbentuk Gamma dengan parameter baru:

α=α+i=1nxi,β=βnβ+1\alpha^* = \alpha + \sum_{i=1}^n x_i, \qquad \beta^* = \frac{\beta}{n\beta + 1}

Langkah 4 — Mean posterior:

E[ΘX]=αβ=(α+xi)βnβ+1E[\Theta \mid \mathbf{X}] = \alpha^* \beta^* = \frac{(\alpha + \sum x_i)\beta}{n\beta + 1}

Langkah 5 — Reformulasi sebagai kredibilitas Bühlmann:

E[ΘX]=nβnβ+1Xˉ+1nβ+1αβE[\Theta \mid \mathbf{X}] = \frac{n\beta}{n\beta + 1} \cdot \bar{X} + \frac{1}{n\beta + 1} \cdot \alpha\beta =ZXˉ+(1Z)μ= Z\bar{X} + (1-Z)\mu

di mana Z=nβnβ+1=nn+1/βZ = \frac{n\beta}{n\beta + 1} = \frac{n}{n + 1/\beta}, sehingga k=1/βk = 1/\beta.

Untuk Poisson–Gamma: v=E[Var(XΘ)]=E[Θ]=αβv = E[\text{Var}(X\|\Theta)] = E[\Theta] = \alpha\beta dan a=Var(Θ)=αβ2a = \text{Var}(\Theta) = \alpha\beta^2, sehingga k=v/a=1/βk = v/a = 1/\beta. Hasil ini konsisten sempurna dengan formula Bühlmann.

Dilarang
  1. Jangan menyamakan prior mean μ\mu dengan sample mean Xˉ\bar{X} — keduanya berbeda secara konseptual dan numerik.
  2. Jangan lupa update kedua parameter posterior (bukan hanya α\alpha^*) ketika menghitung mean posterior dalam kasus Gamma.
  3. Jangan menggunakan hubungan Bühlmann–Bayes di luar kasus conjugate tanpa klarifikasi — hubungan tersebut hanya eksak untuk pasangan conjugate; secara umum Bühlmann adalah aproksimasi linear terbaik dari estimator Bayesian.

Section 4 — Contoh Soal

Soal A — Fundamental

Sebuah portofolio asuransi memiliki frekuensi klaim XΘ=θPoisson(θ)X \mid \Theta = \theta \sim \text{Poisson}(\theta). Distribusi prior atas Θ\Theta adalah Gamma(α=3,β=2)\text{Gamma}(\alpha = 3, \beta = 2) (sehingga mean prior =6= 6). Seorang tertanggung menghasilkan n=4n = 4 klaim dalam satu tahun. Hitung estimator Bayesian untuk Θ\Theta.

Solusi Soal A

Pendekatan: Gunakan update Poisson–Gamma. Hitung parameter posterior lalu ambil mean-nya.

1. Identifikasi Variabel

  • α=3\alpha = 3, β=2\beta = 2 → mean prior =6= 6
  • n=1n = 1 tahun dengan total klaim xi=4\sum x_i = 4

2. Identifikasi Distribusi / Model Pasangan Poisson–Gamma: posterior adalah Gamma dengan parameter terupdate.

3. Setup Persamaan

α=α+xi,β=βnβ+1\alpha^* = \alpha + \sum x_i, \qquad \beta^* = \frac{\beta}{n\beta + 1}

4. Eksekusi Aljabar

α=3+4=7\alpha^* = 3 + 4 = 7 β=212+1=23\beta^* = \frac{2}{1 \cdot 2 + 1} = \frac{2}{3} E[ΘX=4]=αβ=7×23=1434,667E[\Theta \mid X = 4] = \alpha^* \cdot \beta^* = 7 \times \frac{2}{3} = \frac{14}{3} \approx 4{,}667

5. Verification Mean posterior 4,6674{,}667 berada di antara mean prior 66 dan data observasi 44 — masuk akal. Dengan n=1n = 1 yang kecil, posterior lebih dekat ke prior.

Hasil: Estimator Bayesian =14/34,667= 14/3 \approx 4{,}667 klaim per tahun.

Exam Tips — Soal A

Target waktu: 2 menit. Common trap: Mengira nn adalah jumlah klaim, padahal nn adalah jumlah periode. Di sini n=1n = 1 periode dengan xi=4\sum x_i = 4 klaim. Shortcut: Hafal rumus β=β/(nβ+1)\beta^* = \beta/(n\beta + 1) langsung.


Soal B — Exam-Typical

Frekuensi klaim tahunan XiΘ=θPoisson(θ)X_i \mid \Theta = \theta \sim \text{Poisson}(\theta). Prior: ΘGamma(α=4,β=0,5)\Theta \sim \text{Gamma}(\alpha = 4, \beta = 0{,}5). Sebuah tertanggung diamati selama n=3n = 3 tahun dengan klaim masing-masing 2,5,32, 5, 3.

(a) Hitung distribusi posterior ΘX\Theta \mid \mathbf{X}. (b) Hitung estimator Bayesian dan nyatakan sebagai kombinasi kredibilitas ZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu. (c) Tentukan faktor kredibilitas ZZ dan bandingkan dengan formula Bühlmann.

Solusi Soal B

Pendekatan: Update posterior Poisson–Gamma, kemudian reformulasi sebagai premi Bühlmann.

1. Identifikasi Variabel

  • Prior: Γ(α=4,β=0,5)\Gamma(\alpha = 4,\, \beta = 0{,}5) → mean prior μ=αβ=2\mu = \alpha\beta = 2
  • Data: n=3n = 3, Xˉ=(2+5+3)/3=10/3\bar{X} = (2+5+3)/3 = 10/3, xi=10\sum x_i = 10

2. Identifikasi Distribusi / Model Pasangan Poisson–Gamma, posterior tetap Gamma.

3. Setup Persamaan

α=α+xi=4+10=14\alpha^* = \alpha + \sum x_i = 4 + 10 = 14 β=βnβ+1=0,53(0,5)+1=0,52,5=0,2\beta^* = \frac{\beta}{n\beta + 1} = \frac{0{,}5}{3(0{,}5) + 1} = \frac{0{,}5}{2{,}5} = 0{,}2

4. Eksekusi Aljabar

(a) Posterior: ΘXGamma(14,  0,2)\Theta \mid \mathbf{X} \sim \text{Gamma}(14,\; 0{,}2)

(b) Mean posterior:

E[ΘX]=14×0,2=2,8E[\Theta \mid \mathbf{X}] = 14 \times 0{,}2 = 2{,}8

Formula kredibilitas:

Z=nβnβ+1=3×0,53×0,5+1=1,52,5=0,6Z = \frac{n\beta}{n\beta + 1} = \frac{3 \times 0{,}5}{3 \times 0{,}5 + 1} = \frac{1{,}5}{2{,}5} = 0{,}6 ZXˉ+(1Z)μ=0,6×103+0,4×2=2+0,8=2,8Z\bar{X} + (1-Z)\mu = 0{,}6 \times \frac{10}{3} + 0{,}4 \times 2 = 2 + 0{,}8 = 2{,}8 \checkmark

(c) Faktor Bühlmann: Untuk Poisson–Gamma, k=1/β=2k = 1/\beta = 2:

ZBu¨hlmann=nn+k=33+2=0,6Z^{Bühlmann} = \frac{n}{n+k} = \frac{3}{3+2} = 0{,}6

Identik dengan ZZ dari Bayes. Hubungan terbukti.

5. Verification Posterior mean 2,82{,}8 berada di antara μ=2\mu = 2 dan Xˉ=3,33\bar{X} = 3{,}33. Dengan Z=0,6Z = 0{,}6, data lebih mendominasi — wajar karena n=3n=3 sudah cukup untuk menggeser dari prior.

Hasil: μ^Bayes=2,8\hat{\mu}^{Bayes} = 2{,}8; Z=0,6Z = 0{,}6; formula Bühlmann dan Bayes memberikan hasil identik.

Exam Tips — Soal B

Target waktu: 4 menit. Common trap: Salah menghitung μ=αβ\mu = \alpha\beta (mean prior) vs Xˉ\bar{X} (sample mean) saat menyusun formula ZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu. Shortcut: k=1/βk = 1/\beta langsung dari parameter prior untuk kasus Poisson–Gamma.


Soal C — Challenging

Klaim asuransi jiwa XiQ=qBinomial(m=5,q)X_i \mid Q = q \sim \text{Binomial}(m = 5, q), di mana qq adalah peluang klaim per eksposur. Prior: QBeta(α=2,β=8)Q \sim \text{Beta}(\alpha = 2, \beta = 8).

Sebuah kelompok diamati selama n=4n = 4 tahun dengan jumlah klaim: 1,3,2,21, 3, 2, 2 (dari m=5m = 5 eksposur per tahun).

(a) Hitung distribusi posterior QXQ \mid \mathbf{X}. (b) Hitung estimator Bayesian E[QX]E[Q \mid \mathbf{X}]. (c) Nyatakan sebagai kombinasi kredibilitas. Berapa ZZ dan kk? (d) Hitung v=E[v(Q)]v = E[v(Q)] dan a=Var(μ(Q))a = \text{Var}(\mu(Q)), verifikasi k=v/ak = v/a.

Solusi Soal C

Pendekatan: Pasangan Binomial–Beta. Update posterior Beta, hitung mean, reformulasi sebagai kredibilitas.

1. Identifikasi Variabel

  • Prior: QBeta(α=2,β=8)Q \sim \text{Beta}(\alpha = 2, \beta = 8) → mean =α/(α+β)=2/10=0,2= \alpha/(\alpha+\beta) = 2/10 = 0{,}2
  • m=5m = 5, n=4n = 4 tahun, klaim: 1,3,2,21, 3, 2, 2xi=8\sum x_i = 8, Xˉ=2\bar{X} = 2
  • μ(q)=E[XQ=q]=mq=5q\mu(q) = E[X \mid Q = q] = mq = 5q

2. Identifikasi Distribusi / Model Pasangan Binomial–Beta: posterior adalah Beta dengan parameter terupdate.

3. Setup Persamaan

Likelihood untuk xi\sum x_i klaim dari nmnm eksposur:

f(xq)qxi(1q)nmxif(\mathbf{x} \mid q) \propto q^{\sum x_i}(1-q)^{nm - \sum x_i}

Prior: π(q)qα1(1q)β1\pi(q) \propto q^{\alpha-1}(1-q)^{\beta-1}

Posterior: π(qx)qα+xi1(1q)β+nmxi1\pi(q \mid \mathbf{x}) \propto q^{\alpha + \sum x_i - 1}(1-q)^{\beta + nm - \sum x_i - 1}

4. Eksekusi Aljabar

(a) Posterior:

α=α+xi=2+8=10\alpha^* = \alpha + \sum x_i = 2 + 8 = 10 β=β+nmxi=8+(4)(5)8=8+12=20\beta^* = \beta + nm - \sum x_i = 8 + (4)(5) - 8 = 8 + 12 = 20 QXBeta(10,20)Q \mid \mathbf{X} \sim \text{Beta}(10, 20)

(b) Mean posterior:

E[QX]=αα+β=1030=130,333E[Q \mid \mathbf{X}] = \frac{\alpha^*}{\alpha^* + \beta^*} = \frac{10}{30} = \frac{1}{3} \approx 0{,}333

(c) Reformulasi kredibilitas: Di sini kita bekerja pada skala X=mQX = mQ. Mean prior μ=mαα+β=5×0,2=1\mu = m \cdot \frac{\alpha}{\alpha+\beta} = 5 \times 0{,}2 = 1. Sample mean Xˉ=2\bar{X} = 2.

E[μ(Q)X]=mE[QX]=5×13=53E[\mu(Q) \mid \mathbf{X}] = m \cdot E[Q \mid \mathbf{X}] = 5 \times \frac{1}{3} = \frac{5}{3}

Menyatakan sebagai ZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu:

53=Z2+(1Z)1=1+Z    Z=23\frac{5}{3} = Z \cdot 2 + (1-Z) \cdot 1 = 1 + Z \implies Z = \frac{2}{3}

Karena Z=n/(n+k)=4/(4+k)=2/3Z = n/(n+k) = 4/(4+k) = 2/3, maka k=2k = 2.

(d) Verifikasi k=v/ak = v/a:

μ(q)=mq=5q\mu(q) = mq = 5qv(q)=Var(XQ=q)=mq(1q)=5q(1q)v(q) = \text{Var}(X \mid Q = q) = mq(1-q) = 5q(1-q)

v=E[v(Q)]=5E[Q(1Q)]=5(E[Q]E[Q2])v = E[v(Q)] = 5\,E[Q(1-Q)] = 5\left(E[Q] - E[Q^2]\right)

Untuk Beta(α,β)\text{Beta}(\alpha, \beta): E[Q]=210=0,2E[Q] = \frac{2}{10} = 0{,}2 dan E[Q2]=α(α+1)(α+β)(α+β+1)=231011=6110E[Q^2] = \frac{\alpha(\alpha+1)}{(\alpha+\beta)(\alpha+\beta+1)} = \frac{2 \cdot 3}{10 \cdot 11} = \frac{6}{110}

v=5(0,26110)=5(226110)=5×16110=80110=811v = 5\left(0{,}2 - \frac{6}{110}\right) = 5\left(\frac{22 - 6}{110}\right) = 5 \times \frac{16}{110} = \frac{80}{110} = \frac{8}{11} a=Var(μ(Q))=Var(5Q)=25Var(Q)=25×αβ(α+β)2(α+β+1)=25×16100×11=4001100=411a = \text{Var}(\mu(Q)) = \text{Var}(5Q) = 25\,\text{Var}(Q) = 25 \times \frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)} = 25 \times \frac{16}{100 \times 11} = \frac{400}{1100} = \frac{4}{11} k=va=8/114/11=2k = \frac{v}{a} = \frac{8/11}{4/11} = 2 \checkmark

5. Verification k=2k = 2 konsisten dengan yang diperoleh dari formula kredibilitas di langkah (c). Mean posterior 5/3=1,6675/3 = 1{,}667 berada di antara mean prior 11 dan Xˉ=2\bar{X} = 2.

Hasil: Posterior Beta(10,20)\text{Beta}(10, 20); E[μ(Q)X]=5/3E[\mu(Q)\|\mathbf{X}] = 5/3; Z=2/3Z = 2/3; k=v/a=2k = v/a = 2 terbukti.

Exam Tips — Soal C

Target waktu: 6 menit. Common trap 1: Lupa bahwa β\beta^* menyerap nmxinm - \sum x_i (total eksposur dikurangi total klaim), bukan hanya nxin - \sum x_i. Common trap 2: Bekerja pada skala QQ tapi soal minta μ(Q)=mQ\mu(Q) = mQ — pastikan konversi ke skala yang benar. Shortcut: Hitung kk dari formula kredibilitas terlebih dahulu (Z=nm/(nm+α+β)Z = nm/(nm+\alpha+\beta)), lalu verifikasi dengan v/av/a.


Section 5 — Verifikasi & Sanity Check

Sanity Check 1 — Batas Kredibilitas

Faktor ZZ harus selalu berada di antara 0 dan 1. Ketika nn \to \infty, Z1Z \to 1 (data mendominasi). Ketika n=0n = 0, Z=0Z = 0 (estimator = mean prior). Jika ZZ Anda keluar dari [0,1][0, 1], ada kesalahan di perhitungan kk atau nn.

Sanity Check 2 — Mean Posterior di Antara Prior dan Data

Posterior mean selalu berada di antara mean prior μ\mu dan sample mean Xˉ\bar{X} (untuk kasus single-parameter conjugate dengan squared-error loss). Jika hasil Anda di luar rentang [min(μ,Xˉ),max(μ,Xˉ)][\min(\mu, \bar{X}),\, \max(\mu, \bar{X})], ada kesalahan numerik.

Sanity Check 3 — Konsistensi Bühlmann

Setelah menghitung ZZ dari posterior, hitung ulang dengan Z=n/(n+k)Z = n/(n+k) menggunakan k=v/ak = v/a dari distribusi prior. Kedua nilai harus identik untuk pasangan conjugate.

Metode Alternatif

Untuk pasangan Poisson–Gamma, terdapat cara cepat menghitung ZZ tanpa mengeksplisitkan posterior:

Z=nn+k,k=1βZ = \frac{n}{n + k}, \quad k = \frac{1}{\beta}

Cukup baca β\beta dari parameter prior Gamma dan langsung hitung k=1/βk = 1/\beta. Tidak perlu menghitung vv dan aa secara terpisah untuk exam.


Section 6 — Visualisasi Mental

Bayesian Update sebagai Pergerakan Distribusi:

Bayangkan distribusi prior π(θ)\pi(\theta) sebagai “kepercayaan awal” yang berbentuk kurva lonceng (atau kurva miring untuk Gamma). Setiap observasi data menggeser dan mempersempit kurva tersebut menuju nilai yang didukung data:

Prior π(θ):        lebar, terpusat di μ
                   ┌─────────────────────┐
                   │      ╭───╮          │
                   │    ╭─╯   ╰─╮        │
                   │ ╭──╯       ╰──╮     │
                   └─┴─────────────┴─────┘
                         μ (mean prior)

Posterior π(θ|X): lebih sempit, bergeser ke data
                   ┌─────────────────────┐
                   │           ╭─╮       │
                   │         ╭─╯ ╰─╮     │
                   │       ╭─╯     ╰─╮   │
                   └───────┴─────────┴───┘
                    μ    posterior  X̄
                         mean

Semakin banyak data (nn besar), semakin sempit posterior dan semakin dekat ke Xˉ\bar{X}.

Interpretasi ZZ sebagai bobot:

Z = 0 (n=0): ────────[μ]──────────────────[X̄]
             estimator = μ

Z = 0.5:     ────────[μ]────[E[μ(Θ)|X]]──[X̄]
             estimator di tengah

Z = 1 (n→∞): ────────[μ]──────────────────[X̄]
             estimator = X̄

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Rumus
Lebar kurva priora=Var(μ(Θ))a = \text{Var}(\mu(\Theta)) — VHM
Lebar kurva likelihoodv=E[v(Θ)]v = E[v(\Theta)] — EVPV
Posisi kurva posteriorZXˉ+(1Z)μZ\bar{X} + (1-Z)\mu
Pergeseran ke kanan (jika Xˉ>μ\bar{X} > \mu)Z(Xˉμ)Z(\bar{X} - \mu)
Penyempitan posteriorVar(μ(Θ)X)\text{Var}(\mu(\Theta)\|\mathbf{X}) mengecil seiring nn

Section 7 — Jebakan Umum

Kesalahan Parametrisasi — Gamma dan Beta

Distribusi Gamma memiliki dua konvensi parametrisasi: (shape,scale)(shape, scale) di mana mean =αβ= \alpha\beta, versus (shape,rate)(shape, rate) di mana mean =α/λ= \alpha/\lambda. Selalu periksa konvensi soal. Klugman menggunakan parametrisasi (α,θ)(\alpha, \theta) dengan mean =αθ= \alpha\theta (di mana θ\theta adalah scale, bukan rate).

Kesalahan Konseptual
  1. Mencampuradukkan nn (jumlah periode) dan xi\sum x_i (total klaim): Di update Poisson–Gamma, α=α+xi\alpha^* = \alpha + \sum x_i (bukan α+n\alpha + n), dan β\beta^* menggunakan nn bukan xi\sum x_i.
  2. Mengira Bühlmann selalu sama dengan Bayes: Bühlmann adalah aproksimasi linear terbaik dari estimator Bayesian. Hanya untuk pasangan conjugate keduanya persis sama.
  3. Menganggap posterior mean adalah parameter posterior: Posterior mean E[ΘX]E[\Theta \mid \mathbf{X}] berbeda dari parameter distribusi posterior α,β\alpha^*, \beta^*.
  4. Lupa bahwa μ=E[μ(Θ)]\mu = E[\mu(\Theta)], bukan mean sample: μ\mu adalah mean dari distribusi campuran (grand mean), yang dihitung dari distribusi prior — bukan dari data.
Kesalahan Interpretasi Soal

Jika soal minta “premi Bayesian untuk tahun depan”, yang dimaksud adalah E[Xn+1X]=E[μ(Θ)X]E[X_{n+1} \mid \mathbf{X}] = E[\mu(\Theta) \mid \mathbf{X}]. Ini berbeda dari E[ΘX]E[\Theta \mid \mathbf{X}] kecuali jika μ(θ)=θ\mu(\theta) = \theta (kasus Poisson di mana μ(θ)=θ\mu(\theta) = \theta).

Red Flags
  • Soal menyebut “distribusi prior” dan “data observasi” → langsung pikirkan Bayes update
  • Soal menyebut “Poisson” + “Gamma” → pasangan conjugate, gunakan formula update langsung
  • Soal minta verifikasi k=v/ak = v/a → hitung vv dan aa dari distribusi prior, bukan posterior
  • Soal berikan Xˉ\bar{X} dan mean prior → cari ZZ dari E[μ(Θ)X]=ZXˉ+(1Z)μE[\mu(\Theta)\|\mathbf{X}] = Z\bar{X} + (1-Z)\mu

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. Estimator Bayesian adalah mean posterior: E[μ(Θ)X]=μ(θ)π(θx)dθE[\mu(\Theta) \mid \mathbf{X}] = \int \mu(\theta)\,\pi(\theta \mid \mathbf{x})\,d\theta

  2. Update Poisson–Gamma: α=α+xi\alpha^* = \alpha + \sum x_i; β=β/(nβ+1)\beta^* = \beta/(n\beta+1); k=1/βk = 1/\beta

  3. Update Binomial–Beta: α=α+xi\alpha^* = \alpha + \sum x_i; β=β+nmxi\beta^* = \beta + nm - \sum x_i

  4. Identitas kunci Bühlmann–Bayes (hanya untuk conjugate): E[μ(Θ)X]=ZXˉ+(1Z)μ,Z=nn+k,k=vaE[\mu(\Theta) \mid \mathbf{X}] = Z\bar{X} + (1-Z)\mu, \quad Z = \frac{n}{n+k}, \quad k = \frac{v}{a}

  5. Verifikasi: posterior mean selalu di antara μ\mu dan Xˉ\bar{X}; Z[0,1]Z \in [0,1]

Kapan Digunakan

  • Soal menyebutkan distribusi prior atas parameter risiko
  • Soal minta premi prediktif berdasarkan data historis individu
  • Soal minta pembuktian bahwa estimator Bayesian sama dengan kredibilitas Bühlmann
  • Pasangan distribusi yang disebutkan adalah Poisson–Gamma, Binomial–Beta, atau Normal–Normal

Kapan TIDAK Boleh Digunakan

Quick Decision Tree

graph TD
    A["Soal minta estimasi kerugian masa depan?"] -->|"Ya"| B["Ada distribusi prior yang diberikan?"]
    A -->|"Tidak"| Z["Bukan topik kredibilitas"]
    B -->|"Ya"| C["Pasangan conjugate?<br>(Poisson-Gamma, Binomial-Beta, Normal-Normal)"]
    B -->|"Tidak"| D["Gunakan Bühlmann / Empirical Bayes<br>[[7.2]] atau [[7.4]]"]
    C -->|"Ya"| E["Update parameter posterior<br>Hitung mean posterior"]
    C -->|"Tidak"| F["Hitung numerik atau<br>aproksimasi Bühlmann"]
    E --> G["Nyatakan sebagai Z*Xbar + (1-Z)*mu<br>Verifikasi k = v/a"]

Follow-up Options
  1. “Berikan contoh soal variasi Normal–Normal conjugate untuk 7.3”
  2. “Jelaskan hubungan 7.3 Bayesian Credibility dengan 7.4 Empirical Bayesian Methods
  3. “Buat flashcard 1-halaman untuk rumus-rumus update conjugate di 7.3”

📖 Ref: Klugman et al. (2019) Bab 16–18; Tse (2009) Bab 6–9 | 🗓️ 2026-04-17 | #TA2 #BayesianCredibility #TeoriRisiko