MediumBobot: 5–10%Klugman, Panjer & Willmot (2019), Loss Models 5th ed., Bab 6
TA2FrekuensiKlaimMGFPGFTeoriRisiko
📊 2.1 — Frequency MGF and PGF
≡Ringkasan Cepat›
Topik: Fungsi Pembangkit Momen & Probabilitas untuk Model Frekuensi | Bobot: ~5–10% | Difficulty: Medium
Ref: Klugman et al. (2019), Bab 6 | Prereq:1.1 Moment and Probability Generating Functions
Section 0 — Pemetaan Topik
Topik TA2
Sub-topik ID
Skill Diuji
Bobot
Difficulty
Prerequisite
Connected Topics
Referensi
Model Frekuensi Klaim
2.1
Menentukan MGF dan PGF untuk distribusi frekuensi; menghitung momen dan probabilitas dari kedua fungsi
Bayangkan Anda adalah seorang aktuaris di perusahaan asuransi kendaraan bermotor. Setiap bulan, ratusan nasabah mengajukan klaim. Anda perlu memodelkan berapa banyak klaim yang masuk — bukan seberapa besar nilai klaim tersebut, melainkan jumlahnya. Inilah yang disebut model frekuensi klaim: suatu distribusi probabilitas atas bilangan cacah N=0,1,2,3,… yang merepresentasikan banyaknya klaim dalam satu periode.
Masalahnya, distribusi frekuensi memiliki banyak bentuk — Poisson, Binomial, Negatif Binomial, dan lainnya. Untuk masing-masing, kita perlu bisa menghitung rata-rata klaim, variansi, bahkan peluang tepat k klaim terjadi. Menghitung ini langsung dari distribusi probabilitas bisa sangat melelahkan. Di sinilah Fungsi Pembangkit menjadi alat yang ampuh.
Terdapat dua fungsi pembangkit utama: MGF (Moment Generating Function) dan PGF (Probability Generating Function). Keduanya adalah cara untuk “mengemas” seluruh distribusi ke dalam satu fungsi yang kompak. Ketika kita turunkan (differentiate) fungsi tersebut, informasi tentang momen atau probabilitas “keluar” dengan sendirinya. Khusus untuk distribusi frekuensi (diskrit, bernilai non-negatif), PGF sangat elegan: turunan ke-k di titik nol langsung memberi probabilitas P(N=k), dikali faktorial.
Section 2 — Definisi Formal
ℹDefinisi Matematis›
Misalkan N adalah variabel acak diskrit non-negatif (banyaknya klaim) dengan fungsi massa probabilitas pk=P(N=k) untuk k=0,1,2,…
MGF:MN(t)=E[etN]=∑k=0∞etkpk
PGF:PN(z)=E[zN]=∑k=0∞zkpk
Simbol
Makna
Catatan
N
Variabel acak frekuensi klaim
Diskrit, bernilai 0,1,2,…
pk
P(N=k), fungsi massa probabilitas
∑k=0∞pk=1
MN(t)
MGF dari N
Terdefinisi di sekitar t=0
PN(z)
PGF dari N
PN(1)=1; konvergen untuk ∥z∥≤1
E[N]
Mean frekuensi klaim
Rata-rata jumlah klaim per periode
Var(N)
Variansi frekuensi klaim
Ukuran sebaran jumlah klaim
Rumus Utama
Hubungan MGF dan PGF:
PN(z)=MN(lnz)atau ekuivalenMN(t)=PN(et)
Label: Substitusi z=et menghubungkan kedua fungsi secara langsung.
Momen dari MGF:
E[Nr]=MN(r)(0)=dtrdrMN(t)t=0
Label: Turunan ke-r dari MGF di t=0 menghasilkan momen ke-r.
N adalah variabel acak diskrit non-negatif (N∈{0,1,2,…}).
∑k=0∞pk=1 — distribusi terdefinisi dengan baik.
MGF diasumsikan terdefinisi di sekitar t=0 (ada interval terbuka yang memuat t=0).
PGF konvergen setidaknya untuk ∣z∣≤1, sehingga PN(1)=1 selalu valid.
Jika distribusi memiliki support terbatas (mis. Binomial), MGF dan PGF terdefinisi di seluruh R dan C.
Section 3 — Jembatan Logika
◈Dari Definisi ke Rumus›
PGF lahir dari ekspansi deret pangkat (power series). Karena N hanya bernilai 0,1,2,…, kita dapat menulis E[zN]=∑k=0∞zkpk — ini persis deret pangkat dengan koefisien pk. Akibatnya, pk dapat “diekstrak kembali” dengan menurunkan k kali lalu evaluasi di z=0: pk=PN(k)(0)/k!. Inilah kekuatan utama PGF untuk frekuensi.
◈Support dan Domain›
PGF selalu konvergen untuk ∣z∣≤1 karena ∑∣z∣kpk≤∑pk=1.
MGF mungkin tidak ada untuk semua distribusi, tetapi PGF selalu ada untuk distribusi diskrit non-negatif.
Di titik z=1: PN(1)=∑pk=1 selalu terpenuhi — gunakan ini untuk sanity check.
Di titik z=0: PN(0)=p0=P(N=0) — PGF di nol adalah peluang tidak ada klaim!
Verifikasi via formula langsung: p2=2!e−3⋅32=29e−3 ✓
5. Verification
Untuk Poisson, E[N]=Var(N)=λ=3. Kedua nilai sama. ✓ Magnitude p2≈0.224 masuk akal sebagai mode distribusi Poisson(λ=3).
Hasil:E[N]=3, Var(N)=3, P(N=2)=29e−3≈0.2240
▲Exam Tips — Soal A›
Target waktu: 2–3 menit. Common trap: Lupa faktor k!=2 saat menghitung p2 dari turunan kedua PGF. Shortcut: Untuk Poisson, mean = variansi = λ — konfirmasi ini sebagai sanity check instan.
Soal B — Exam-Typical
Jumlah klaim bulanan N berdistribusi Negatif Binomial dengan parameter r=2 dan β=1.5.
(a) Tentukan PGF dari N.
(b) Hitung E[N], E[N2], dan Var(N) menggunakan PGF.
(c) Hitung P(N=0) dan P(N=1) dari PGF.
✓Solusi Soal B›
Pendekatan: Gunakan PGF Negatif Binomial, turunkan secara analitik untuk momen, evaluasi di z=0 untuk probabilitas.
1. Identifikasi Variabel
N∼NegBin(r=2,β=1.5)
Mean teoritis: rβ=2×1.5=3
Variansi teoritis: rβ(1+β)=2×1.5×2.5=7.5
2. Identifikasi Distribusi / Model
Negatif Binomial — dipilih karena variansi > mean (overdispersion), cocok untuk data klaim yang mengelompok (contagious).
Target waktu: 3–4 menit. Common trap: Kesalahan saat menerapkan chain rule pada turunan (2.5−1.5z)−2 — jangan lupa mengalikan dengan −1.5 dari turunan bagian dalam. Shortcut: Evaluasi di z=1 selalu memberi (2.5−1.5)−k=1−k=1 — kalkulasi jauh lebih sederhana.
Soal C — Challenging
Diketahui PGF dari suatu distribusi frekuensi N adalah:
PN(z)=1−0.6z0.4
(a) Identifikasi distribusi N beserta parameternya.
(b) Hitung E[N], Var(N).
(c) Misalkan S=X1+X2+⋯+XN adalah model agregat dengan klaim individual Xi i.i.d. berdistribusi Eksponensial dengan mean 1000, independen dari N. Tentukan E[S] dan Var(S).
✓Solusi Soal C›
Pendekatan: Kenali bentuk PGF → identifikasi distribusi → hitung momen → gunakan formula model agregat kolektif.
1. Identifikasi Variabel
PN(z)=1−0.6z0.4 — bentuk geometri
Xi∼Exp(θ=1000), jadi E[X]=1000, E[X2]=2θ2=2,000,000, Var(X)=θ2=1,000,000
2. Identifikasi Distribusi / Model
PGF Geometric dengan parameter β: PN(z)=1+β1⋅1−1+ββz1.
Bandingkan: 1−0.6z0.4. Maka 1+β1=0.4⇒β=1.5. Cek: 1+ββ=2.51.5=0.6 ✓.
Jadi N∼Geometric(β=1.5) (kasus khusus NegBin dengan r=1).
5. VerificationVar(S)>E[N]⋅Var(X) karena variansi frekuensi juga berkontribusi. Rasio Var(S)/E[S]2=5,250,000/2,250,000=2.33>1 — menunjukkan dispersi tinggi yang konsisten dengan Geometric distribution.
Target waktu: 4–6 menit. Common trap: Soal C menghubungkan PGF dengan model agregat — jangan lupa formula Var(S)=E[N]Var(X)+Var(N)(E[X])2, bukan hanya E[N]Var(X). Shortcut: Kenali bentuk PGF langsung: 1−dzc selalu Geometric; eλ(z−1) selalu Poisson.
Section 5 — Verifikasi & Sanity Check
✓Sanity Check 1 — PGF di z=1
Untuk semua distribusi diskrit, PN(1)=1 tanpa kecuali.
Gunakan ini untuk verifikasi PGF Anda: substitusi z=1, pastikan hasilnya 1.
Contoh: Poisson: eλ(1−1)=e0=1 ✓; NegBin: [1/(1+β⋅0)]r=1 ✓.›
✓Sanity Check 2 — PGF di z=0 adalah p0PN(0)=∑k=0∞0kpk=p0⋅1+0+0+⋯=p0.
Jika soal memberi PGF, evaluasi di z=0 untuk mendapat P(N=0) secara instan.
Contoh: Poisson: PN(0)=eλ(0−1)=e−λ=p0 ✓.›
✓Sanity Check 3 — Relasi MGF–PGF›
MN(lnz)=PN(z). Substitusi t=lnz (atau z=et) harus menghasilkan ekspresi yang sama.
Contoh Poisson: MN(t)=eλ(et−1), substitusi t=lnz: eλ(z−1)=PN(z) ✓.
Metode Alternatif
Untuk menghitung momen, dapat juga menggunakan MGF langsung:
PGF lebih natural untuk distribusi frekuensi karena:
Koefisiennya langsung berupa probabilitas.
PGF selalu konvergen untuk ∣z∣≤1.
PGF distribusi komposit (compound) memiliki bentuk yang elegan: PS(z)=PN(PX(z)).
Section 6 — Visualisasi Mental
Kurva PGF Poisson(λ=3): Fungsi PN(z)=e3(z−1) adalah kurva eksponensial yang melewati titik (0,e−3)≈(0,0.050) di sumbu-y dan (1,1) di kanan. Fungsi ini cekung ke atas (convex), meningkat dari kiri ke kanan. Nilai di z=0 adalah p0; kemiringan di z=1 adalah mean λ.
NegBin konvensi berbeda-beda: Klugman menggunakan (r,β) di mana mean =rβ, variansi =rβ(1+β). Beberapa referensi lain menggunakan (r,p) di mana p=1/(1+β). Cek selalu: “mean berapa?” sebelum menulis PGF.
MGF vs PGF — mana yang digunakan? Turunan MGF di t=0 → momen biasa. Turunan PGF di z=1 → momen faktorial. Turunan PGF di z=0 dibagi k! → probabilitas pk.
Lupa faktor k!:pk=PN(k)(0)/k!, bukanPN(k)(0). Ini bukan kesalahan kecil — akan memberi jawaban salah besar.
Momen faktorial vs momen biasa:PN′′(1)=E[N(N−1)]=E[N2]. Untuk mendapat E[N2]: tambahkan E[N], yaitu E[N2]=PN′′(1)+PN′(1).
Asumsi konvergensi: Jangan gunakan PGF di ∣z∣>1 tanpa memverifikasi konvergensinya terlebih dahulu.
⬡Kesalahan Interpretasi Soal›
Soal bertanya “factorial moment kedua”: ini adalah PN′′(1), bukanE[N2].
Soal bertanya “momen kedua”: ini E[N2]=PN′′(1)+PN′(1).
“PGF dievaluasi di z=0” → jawaban adalah p0, bukan 0.
Jika soal memberikan MGF dan meminta PGF, ingat PN(z)=MN(lnz).
▲Red Flags›
Kata “tentukan probabilitas P(N=k) dari PGF” → turunkan k kali, evaluasi di z=0, bagi k!.
Kata “tentukan mean menggunakan PGF” → turunkan sekali, evaluasi di z=1.
Bentuk PGF diberikan dan diminta identifikasi distribusi → kenali: eλ(z−1)=Poisson, (1−q+qz)m=Binomial, [1/(1+β(1−z))]r=NegBin.
Soal melibatkan model kolektif S=∑i=1NXi → butuh E[N] dan Var(N) dari PGF, lalu gunakan conditional variance formula.
Untuk distribusi kontinu (besar klaim) — gunakan MGF biasa, bukan PGF.
Untuk momen bersyarat tanpa kondisi yang jelas — PGF hanya untuk distribusi marginal.
Jika distribusi tidak memiliki support di {0,1,2,…} — PGF tidak terdefinisi.
Quick Decision Tree
graph TD A["Soal melibatkan distribusi frekuensi N"] --> B{"Apa yang diminta?"} B --> |"Mean atau Var(N)"| C["Turunkan PGF<br>Evaluasi di z=1<br>E[N] = P'(1)"] B --> |"P(N=k)"| D["Turunkan k kali<br>Evaluasi di z=0<br>p_k = P^(k)(0)/k!"] B --> |"Identifikasi distribusi"| E{"Bentuk PGF?"} B --> |"Momen faktorial ke-r"| F["Turunkan r kali<br>Evaluasi di z=1<br>P^(r)(1)"] E --> |"exp(lambda*(z-1))"| G["Poisson(lambda)"] E --> |"(1-q+qz)^m"| H["Binomial(m,q)"] E --> |"[1/(1+beta*(1-z))]^r"| I["NegBin(r,beta)"] C --> J["Var(N) = P''(1)+P'(1)-[P'(1)]^2"]
❝Follow-up Options›
“Berikan contoh soal variasi tentang PGF distribusi campuran (mixed Poisson)”