AktuNotes
← Kembali
TA2 · Materi 1.1

Moment and Probability Generating Functions

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

📊 1.1 — Moment and Probability Generating Functions

Ringkasan Cepat

Topik: Moment and Probability Generating Functions | Bobot: ~5–10% | Difficulty: Medium Ref: Klugman et al. (2019) Loss Models 5th ed., Bab 3–5 | Prereq: None

Section 0 — Pemetaan Topik

Topik TA2Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Model Besar Klaim1.1Menghitung momen (mean, variance, skewness) dan probabilitas dari MGF dan PGF; mengidentifikasi distribusi dari bentuk MGF/PGF5–10%MediumNone1.2 Distribution Classes and Extreme Value, 1.4 Tail Characteristics, 2.1 Frequency MGF and PGFKPW (2019) Bab 3–5

Section 1 — Intuisi

Bayangkan kamu bekerja di bagian underwriting dan harus memutuskan apakah premi yang dikenakan cukup untuk menanggung klaim di masa depan. Untuk melakukan itu, kamu perlu memahami perilaku distribusi klaim — bukan hanya nilai rata-ratanya, tapi juga seberapa bervariasi klaim tersebut, seberapa sering klaim besar (extreme) bisa terjadi, dan bagaimana bentuk distribusinya secara keseluruhan. Semua informasi itu tersimpan dalam apa yang disebut momen.

Momen adalah ukuran ringkas yang mendeskripsikan suatu distribusi probabilitas. Momen pertama adalah nilai rata-rata (ekspektasi), momen kedua berkaitan dengan variansi, momen ketiga dengan kemiringan (skewness), dan seterusnya. Alih-alih menghitung setiap momen secara terpisah — yang bisa menjadi sangat tedious — para aktuaris menggunakan alat bantu matematika bernama fungsi pembangkit momen (MGF) dan fungsi pembangkit probabilitas (PGF). Fungsi-fungsi ini adalah semacam “kapsul waktu” yang mengemas semua informasi momen dan probabilitas dari suatu distribusi ke dalam satu fungsi yang elegan.

Dalam konteks pemodelan risiko, MGF sangat berguna ketika kita menjumlahkan banyak variabel acak independen — misalnya total klaim dari banyak pemegang polis. Karena MGF dari jumlah variabel acak independen adalah hasil kali MGF masing-masing variabel, pekerjaan analitis yang kompleks menjadi jauh lebih sederhana. PGF serupa dengan MGF, namun dirancang khusus untuk variabel acak diskrit seperti jumlah klaim (frequency), sehingga memungkinkan kita mengekstrak probabilitas klaim sebanyak kk kali langsung dari turunan fungsi tersebut.

Section 2 — Definisi Formal

Definisi Matematis

Moment Generating Function (MGF):

MX(t)=E[etX]=etxfX(x)dx(kontinu)M_X(t) = E[e^{tX}] = \int_{-\infty}^{\infty} e^{tx} f_X(x)\, dx \quad \text{(kontinu)} MX(t)=E[etX]=xetxpX(x)(diskrit)M_X(t) = E[e^{tX}] = \sum_{x} e^{tx} p_X(x) \quad \text{(diskrit)}

Probability Generating Function (PGF):

PX(z)=E[zX]=k=0zkpk,pk=P(X=k)P_X(z) = E[z^X] = \sum_{k=0}^{\infty} z^k \, p_k, \quad p_k = P(X = k)
SimbolMaknaCatatan
MX(t)M_X(t)MGF dari variabel acak XXTerdefinisi untuk tt di sekitar 0
ttParameter transformasiBiasanya t(h,h)t \in (-h, h) untuk suatu h>0h > 0
PX(z)P_X(z)PGF dari variabel acak diskrit XXHanya untuk X0X \geq 0 bilangan bulat
zzParameter transformasi PGFBiasanya z[0,1]z \in [0, 1] atau z1\|z\| \leq 1
pkp_kProbabilitas P(X=k)P(X = k)Koefisien ekspansi PGF
μk\mu_k'Momen ke-kk (raw moment)E[Xk]E[X^k]
μk\mu_kMomen sentral ke-kkE[(Xμ)k]E[(X - \mu)^k]
κk\kappa_kKumulan ke-kkDiturunkan dari log MGF

Rumus Utama

Menghitung momen ke-kk dari MGF:

E[Xk]=MX(k)(0)=dkdtkMX(t)t=0E[X^k] = M_X^{(k)}(0) = \left.\frac{d^k}{dt^k} M_X(t)\right|_{t=0}

Label: Turunan ke-kk MGF dievaluasi di t=0t=0 menghasilkan momen ke-kk.

Kumulan generating function (CGF) dan kumulan:

KX(t)=lnMX(t),κk=KX(k)(0)K_X(t) = \ln M_X(t), \quad \kappa_k = K_X^{(k)}(0)

Label: Kumulan pertama adalah mean, kumulan kedua adalah variansi.

Mean dan variansi dari kumulan:

κ1=E[X]=μ,κ2=Var(X)=σ2\kappa_1 = E[X] = \mu, \qquad \kappa_2 = \text{Var}(X) = \sigma^2

Label: Cukup hitung CGF sekali untuk mendapat mean dan variansi secara efisien.

Menghitung probabilitas dari PGF:

pk=P(X=k)=PX(k)(0)k!p_k = P(X = k) = \frac{P_X^{(k)}(0)}{k!}

Label: Koefisien Taylor dari PGF di z=0z=0 adalah probabilitas titik.

Hubungan MGF dan PGF:

PX(z)=MX(lnz),MX(t)=PX(et)P_X(z) = M_X(\ln z), \qquad M_X(t) = P_X(e^t)

Label: Dua sisi dari koin yang sama — pilih mana yang lebih mudah sesuai konteks.

MGF untuk jumlah variabel acak independen:

MX1+X2++Xn(t)=i=1nMXi(t)M_{X_1 + X_2 + \cdots + X_n}(t) = \prod_{i=1}^{n} M_{X_i}(t)

Label: Properti kunci untuk model agregat — total klaim = produk MGF individual.

PGF untuk jumlah variabel acak independen:

PX1++Xn(z)=i=1nPXi(z)P_{X_1 + \cdots + X_n}(z) = \prod_{i=1}^{n} P_{X_i}(z)

Label: Analog PGF dari properti additif MGF, berlaku untuk distribusi frekuensi.

Momen dari PGF — factorial moments:

E[X(X1)(X2)(Xk+1)]=PX(k)(1)E[X(X-1)(X-2)\cdots(X-k+1)] = P_X^{(k)}(1)

Label: Factorial moment ke-kk didapat dari turunan ke-kk PGF dievaluasi di z=1z=1.

Hubungan factorial moment dan momen biasa:

E[X]=PX(1),E[X2]=PX(1)+PX(1),Var(X)=PX(1)+PX(1)[PX(1)]2E[X] = P_X'(1), \qquad E[X^2] = P_X''(1) + P_X'(1), \qquad \text{Var}(X) = P_X''(1) + P_X'(1) - [P_X'(1)]^2

Label: Cara efisien menghitung momen distribusi frekuensi tanpa MGF.

Asumsi Eksplisit

  1. MGF MX(t)M_X(t) terdefinisi (finite) pada suatu interval terbuka yang mengandung t=0t = 0.
  2. PGF hanya berlaku untuk variabel acak diskrit non-negatif dengan nilai integer (X{0,1,2,}X \in \{0, 1, 2, \ldots\}).
  3. Jika MGF terdefinisi, ia menentukan distribusi secara unik (theorema identifikasi distribusi).
  4. Turunan dan integral dapat dipertukarkan dengan ekspektasi (dominated convergence theorem terpenuhi).
  5. Untuk properti perkalian, variabel-variabel acak harus saling independen.

Section 3 — Jembatan Logika

Dari Definisi ke Rumus

MGF muncul secara alami dari ekspansi Taylor etXe^{tX} di sekitar t=0t = 0:

etX=1+tX+t2X22!+t3X33!+e^{tX} = 1 + tX + \frac{t^2 X^2}{2!} + \frac{t^3 X^3}{3!} + \cdots

Ambil ekspektasi kedua sisi:

MX(t)=E[etX]=1+tE[X]+t22!E[X2]+t33!E[X3]+M_X(t) = E[e^{tX}] = 1 + t\,E[X] + \frac{t^2}{2!}\,E[X^2] + \frac{t^3}{3!}\,E[X^3] + \cdots

Artinya, momen ke-kk adalah koefisien tkk!\frac{t^k}{k!} dalam ekspansi MGF. Jika kita turunkan kk kali dan set t=0t=0, semua suku berpangkat 1\geq 1 gugur, menyisakan E[Xk]E[X^k]. Inilah mengapa MX(k)(0)=E[Xk]M_X^{(k)}(0) = E[X^k].

Support dan Domain

MGF tidak selalu ada untuk semua distribusi. Distribusi heavy-tail seperti Pareto atau Lognormal tidak memiliki MGF yang terdefinisi di seluruh bilangan real (hanya terdefinisi untuk t0t \leq 0 atau bahkan hanya untuk t=0t = 0). Dalam praktik aktuaria, ini penting: distribusi klaim seringkali heavy-tailed, sehingga PGF lebih aman untuk frekuensi (distribusi diskrit) dan analisis momen langsung lebih umum digunakan untuk distribusi severity.

Derivasi: Mean dan Variansi dari CGF

Langkah 1 — Definisikan CGF:

KX(t)=lnMX(t)K_X(t) = \ln M_X(t)

Langkah 2 — Turunkan terhadap tt:

KX(t)=MX(t)MX(t)K_X'(t) = \frac{M_X'(t)}{M_X(t)}

Langkah 3 — Evaluasi di t=0t = 0: karena MX(0)=E[e0]=1M_X(0) = E[e^0] = 1 dan MX(0)=E[X]M_X'(0) = E[X]:

KX(0)=E[X]1=E[X]=μK_X'(0) = \frac{E[X]}{1} = E[X] = \mu

Langkah 4 — Turunan kedua CGF:

KX(t)=MX(t)MX(t)[MX(t)]2[MX(t)]2K_X''(t) = \frac{M_X''(t) M_X(t) - [M_X'(t)]^2}{[M_X(t)]^2}

Langkah 5 — Evaluasi di t=0t = 0:

KX(0)=E[X2]1[E[X]]21=E[X2][E[X]]2=Var(X)K_X''(0) = \frac{E[X^2] \cdot 1 - [E[X]]^2}{1} = E[X^2] - [E[X]]^2 = \text{Var}(X)

Kesimpulan: Kumulan pertama = mean, kumulan kedua = variansi. CGF lebih efisien karena tidak perlu menghitung E[X2]E[X^2] secara terpisah.

Derivasi: Menghitung Probabilitas dari PGF

Langkah 1 — Tulis definisi PGF:

PX(z)=k=0pkzk=p0+p1z+p2z2+P_X(z) = \sum_{k=0}^{\infty} p_k z^k = p_0 + p_1 z + p_2 z^2 + \cdots

Langkah 2 — Turunkan kk kali terhadap zz:

PX(k)(z)=k!pk+(k+1)k2pk+1z+P_X^{(k)}(z) = k!\, p_k + (k+1)k\cdots 2\, p_{k+1} z + \cdots

Langkah 3 — Evaluasi di z=0z = 0: semua suku mengandung zz gugur:

PX(k)(0)=k!pkP_X^{(k)}(0) = k!\, p_k

Langkah 4 — Isolasi pkp_k:

pk=PX(k)(0)k!p_k = \frac{P_X^{(k)}(0)}{k!}

Ini adalah formula ekuivalen dengan koefisien seri Taylor — probabilitas adalah koefisien ekspansi.

Dilarang
  1. Jangan evaluasi MGF di t0t \neq 0 untuk mendapatkan momen — momen selalu dari turunan dievaluasi di t=0t = 0, bukan di titik lain.
  2. Jangan gunakan PGF untuk variabel kontinu — PGF hanya valid untuk distribusi diskrit integer non-negatif. Untuk distribusi kontinu, gunakan MGF atau Laplace transform.
  3. Jangan asumsikan MGF selalu ada — distribusi heavy-tail khas klaim (Pareto, Lognormal) sering tidak memiliki MGF yang finite. Selalu verifikasi domain tt sebelum menggunakan MGF.

Section 4 — Contoh Soal

Soal A — Fundamental

Variabel acak XX mengikuti distribusi Eksponensial dengan rata-rata θ=500\theta = 500. Hitunglah E[X]E[X], E[X2]E[X^2], dan Var(X)\text{Var}(X) menggunakan MGF.

Solusi Soal A

Pendekatan: Turunkan MGF Eksponensial dua kali, evaluasi di t=0t=0.

1. Identifikasi Variabel

  • Distribusi: XEksponensial(θ)X \sim \text{Eksponensial}(\theta) dengan θ=500\theta = 500
  • fX(x)=1θex/θf_X(x) = \frac{1}{\theta} e^{-x/\theta} untuk x>0x > 0

2. Identifikasi Distribusi / Model Distribusi Eksponensial digunakan sebagai model severity sederhana (waktu antar klaim / besar klaim kecil). MGF-nya memiliki bentuk tertutup yang mudah diturunkan.

3. Setup Persamaan

MX(t)=E[etX]=0etx1θex/θdx=1θ0ex(1/θt)dxM_X(t) = E[e^{tX}] = \int_0^\infty e^{tx} \frac{1}{\theta} e^{-x/\theta}\, dx = \frac{1}{\theta} \int_0^\infty e^{-x(1/\theta - t)}\, dx

4. Eksekusi Aljabar

Integral konvergen jika t<1/θt < 1/\theta, sehingga:

MX(t)=1θ11/θt=11θtM_X(t) = \frac{1}{\theta} \cdot \frac{1}{1/\theta - t} = \frac{1}{1 - \theta t}

Turunan pertama:

MX(t)=θ(1θt)2M_X'(t) = \frac{\theta}{(1 - \theta t)^2} E[X]=MX(0)=θ1=500E[X] = M_X'(0) = \frac{\theta}{1} = 500

Turunan kedua:

MX(t)=2θ2(1θt)3M_X''(t) = \frac{2\theta^2}{(1 - \theta t)^3} E[X2]=MX(0)=2θ2=2×5002=500,000E[X^2] = M_X''(0) = 2\theta^2 = 2 \times 500^2 = 500{,}000 Var(X)=E[X2][E[X]]2=500,000250,000=250,000\text{Var}(X) = E[X^2] - [E[X]]^2 = 500{,}000 - 250{,}000 = 250{,}000

5. Verification Distribusi Eksponensial memiliki E[X]=θE[X] = \theta dan Var(X)=θ2\text{Var}(X) = \theta^2 — cocok dengan hasil di atas. Coefficient of variation =σ/μ=500/500=1= \sigma/\mu = 500/500 = 1, konsisten dengan sifat memoryless Eksponensial.

Hasil: E[X]=500E[X] = 500, E[X2]=500,000E[X^2] = 500{,}000, Var(X)=250,000\text{Var}(X) = 250{,}000.

Exam Tips — Soal A

Target waktu: 2–3 menit. Common trap: Lupa batas domain t<1/θt < 1/\theta — jika soal meminta MX(t)M_X(t) untuk tt tertentu, verifikasi dulu apakah tt valid. Shortcut: Hafalkan MGF distribusi standar: Eksponensial =(1θt)1= (1-\theta t)^{-1}, Gamma =(1θt)α= (1-\theta t)^{-\alpha}, Normal =eμt+σ2t2/2= e^{\mu t + \sigma^2 t^2/2}.

Soal B — Exam-Typical

Variabel acak diskrit NN merepresentasikan jumlah klaim suatu polis dalam setahun, dengan PGF:

PN(z)=0.410.6zP_N(z) = \frac{0.4}{1 - 0.6z}

Hitunglah P(N=0)P(N = 0), P(N=1)P(N = 1), P(N=2)P(N = 2), E[N]E[N], dan Var(N)\text{Var}(N).

Solusi Soal B

Pendekatan: Identifikasi distribusi dari bentuk PGF, lalu gunakan factorial moments untuk momen.

1. Identifikasi Variabel

  • PGF: PN(z)=0.410.6zP_N(z) = \frac{0.4}{1 - 0.6z}
  • Perlu: p0,p1,p2p_0, p_1, p_2, E[N]E[N], Var(N)\text{Var}(N)

2. Identifikasi Distribusi / Model PGF Geometrik dengan parameter pp adalah PN(z)=p1(1p)zP_N(z) = \frac{p}{1 - (1-p)z}. Bandingkan: p=0.4p = 0.4, (1p)=0.6(1-p) = 0.6. Jadi NGeometrik(p=0.4)N \sim \text{Geometrik}(p = 0.4) dengan P(N=k)=0.4×0.6kP(N = k) = 0.4 \times 0.6^k untuk k=0,1,2,k = 0, 1, 2, \ldots

3. Setup Persamaan

Probabilitas dari ekspansi series:

PN(z)=0.4k=0(0.6z)k=k=00.4×0.6kzkP_N(z) = 0.4 \sum_{k=0}^{\infty} (0.6z)^k = \sum_{k=0}^{\infty} 0.4 \times 0.6^k \cdot z^k

Momen dari PGF:

E[N]=PN(1),Var(N)=PN(1)+PN(1)[PN(1)]2E[N] = P_N'(1), \quad \text{Var}(N) = P_N''(1) + P_N'(1) - [P_N'(1)]^2

4. Eksekusi Aljabar

Probabilitas langsung:

P(N=0)=0.4×0.60=0.4P(N=0) = 0.4 \times 0.6^0 = 0.4 P(N=1)=0.4×0.61=0.24P(N=1) = 0.4 \times 0.6^1 = 0.24 P(N=2)=0.4×0.62=0.144P(N=2) = 0.4 \times 0.6^2 = 0.144

Turunan PGF:

PN(z)=0.4×0.6(10.6z)2=0.24(10.6z)2P_N'(z) = \frac{0.4 \times 0.6}{(1-0.6z)^2} = \frac{0.24}{(1-0.6z)^2} E[N]=PN(1)=0.24(0.4)2=0.240.16=1.5E[N] = P_N'(1) = \frac{0.24}{(0.4)^2} = \frac{0.24}{0.16} = 1.5 PN(z)=2×0.24×0.6(10.6z)3=0.288(10.6z)3P_N''(z) = \frac{2 \times 0.24 \times 0.6}{(1-0.6z)^3} = \frac{0.288}{(1-0.6z)^3} PN(1)=0.288(0.4)3=0.2880.064=4.5P_N''(1) = \frac{0.288}{(0.4)^3} = \frac{0.288}{0.064} = 4.5 Var(N)=4.5+1.5(1.5)2=62.25=3.75\text{Var}(N) = 4.5 + 1.5 - (1.5)^2 = 6 - 2.25 = 3.75

5. Verification Geometrik: E[N]=(1p)/p=0.6/0.4=1.5E[N] = (1-p)/p = 0.6/0.4 = 1.5 ✓ dan Var(N)=(1p)/p2=0.6/0.16=3.75\text{Var}(N) = (1-p)/p^2 = 0.6/0.16 = 3.75 ✓.

Hasil: P(N=0)=0.4P(N=0) = 0.4, P(N=1)=0.24P(N=1) = 0.24, P(N=2)=0.144P(N=2) = 0.144, E[N]=1.5E[N] = 1.5, Var(N)=3.75\text{Var}(N) = 3.75.

Exam Tips — Soal B

Target waktu: 3–4 menit. Common trap: Gunakan PN(1)P_N'(1) untuk E[N]E[N], bukan PN(0)P_N'(0) — kesalahan ini sangat umum! Shortcut: Kenali bentuk PGF standar (Geometrik, Poisson, Binomial) — identifikasi distribusi langsung dari PGF menghemat waktu daripada menghitung turunan manual.

Soal C — Challenging

Variabel acak XX dan YY independen dengan XGamma(α=2,θ=1000)X \sim \text{Gamma}(\alpha = 2, \theta = 1000) dan YEksponensial(θ=500)Y \sim \text{Eksponensial}(\theta = 500). Misalkan S=X+YS = X + Y merepresentasikan total klaim dari dua sumber berbeda. Gunakan MGF untuk menentukan E[S]E[S], Var(S)\text{Var}(S), dan E[S3]E[S^3].

Solusi Soal C

Pendekatan: Gunakan additivity MGF untuk jumlah independen, lalu turunkan MGF gabungan.

1. Identifikasi Variabel

  • XGamma(α=2,θ=1000)X \sim \text{Gamma}(\alpha=2, \theta=1000): MX(t)=(11000t)2M_X(t) = (1-1000t)^{-2}, berlaku untuk t<1/1000t < 1/1000
  • YEksponensial(θ=500)Y \sim \text{Eksponensial}(\theta=500): MY(t)=(1500t)1M_Y(t) = (1-500t)^{-1}, berlaku untuk t<1/500t < 1/500
  • S=X+YS = X + Y, independen

2. Identifikasi Distribusi / Model Karena XX dan YY independen, MGF SS adalah produk MGF individual. Domain tt: irisan keduanya, yaitu t<1/1000=0.001t < 1/1000 = 0.001.

3. Setup Persamaan

MS(t)=MX(t)MY(t)=(11000t)2(1500t)1M_S(t) = M_X(t) \cdot M_Y(t) = (1-1000t)^{-2} \cdot (1-500t)^{-1} KS(t)=lnMS(t)=2ln(11000t)ln(1500t)K_S(t) = \ln M_S(t) = -2\ln(1-1000t) - \ln(1-500t)

4. Eksekusi Aljabar

Dari Gamma: E[X]=αθ=2000E[X] = \alpha\theta = 2000, Var(X)=αθ2=2,000,000\text{Var}(X) = \alpha\theta^2 = 2{,}000{,}000

Dari Eksponensial: E[Y]=500E[Y] = 500, Var(Y)=250,000\text{Var}(Y) = 250{,}000

Karena independen:

E[S]=E[X]+E[Y]=2000+500=2500E[S] = E[X] + E[Y] = 2000 + 500 = 2500 Var(S)=Var(X)+Var(Y)=2,000,000+250,000=2,250,000\text{Var}(S) = \text{Var}(X) + \text{Var}(Y) = 2{,}000{,}000 + 250{,}000 = 2{,}250{,}000

Untuk E[S3]E[S^3], gunakan relasi momen-kumulan. Kumulan ketiga κ3=E[(Sμ)3]\kappa_3 = E[(S-\mu)^3] dan:

E[S3]=μ3+3μσ2+κ3E[S^3] = \mu^3 + 3\mu\sigma^2 + \kappa_3

Kumulan ketiga additif: κ3(S)=κ3(X)+κ3(Y)\kappa_3(S) = \kappa_3(X) + \kappa_3(Y).

Untuk Gamma: κ3(X)=2αθ3=2×2×109=4×109\kappa_3(X) = 2\alpha\theta^3 = 2 \times 2 \times 10^9 = 4 \times 10^9

Untuk Eksponensial: κ3(Y)=2θ3=2×5003=250,000,000=2.5×108\kappa_3(Y) = 2\theta^3 = 2 \times 500^3 = 250{,}000{,}000 = 2.5 \times 10^8

κ3(S)=4×109+2.5×108=4.25×109\kappa_3(S) = 4 \times 10^9 + 2.5 \times 10^8 = 4.25 \times 10^9 E[S3]=25003+3×2500×2,250,000+4.25×109E[S^3] = 2500^3 + 3 \times 2500 \times 2{,}250{,}000 + 4.25 \times 10^9 =15,625,000,000+16,875,000,000+4,250,000,000= 15{,}625{,}000{,}000 + 16{,}875{,}000{,}000 + 4{,}250{,}000{,}000 =36,750,000,000= 36{,}750{,}000{,}000

5. Verification E[S]=2500E[S] = 2500 masuk akal sebagai rata-rata total klaim dari dua sumber. Gunakan sanity check: E[S3]>(E[S])3=15.625×109E[S^3] > (E[S])^3 = 15.625 \times 10^9 ✓ karena distribusi right-skewed membuat momen ke-3 selalu melebihi kubik mean.

Hasil: E[S]=2500E[S] = 2500, Var(S)=2,250,000\text{Var}(S) = 2{,}250{,}000, E[S3]=36,750,000,000E[S^3] = 36{,}750{,}000{,}000.

Exam Tips — Soal C

Target waktu: 5–6 menit. Common trap: Lupa bahwa kumulan (bukan momen mentah) yang bersifat additif — κk(X+Y)=κk(X)+κk(Y)\kappa_k(X+Y) = \kappa_k(X) + \kappa_k(Y) untuk independen, tetapi E[(X+Y)k]E[Xk]+E[Yk]E[(X+Y)^k] \neq E[X^k] + E[Y^k]. Shortcut: Untuk Gamma(α,θ)(\alpha, \theta): κk=(k1)!αθk\kappa_k = (k-1)!\,\alpha\theta^k. Hafalkan rumus kumulan distribusi standar.

Section 5 — Verifikasi & Sanity Check

Cross-check MGF dengan Momen Langsung

Untuk distribusi Eksponensial(θ)(\theta), verifikasi selalu bisa dilakukan:

E[Xk]=MX(k)(0)=k!θkE[X^k] = M_X^{(k)}(0) = k!\,\theta^k

Ini konsisten dengan hasil integrasi langsung 0xk1θex/θdx=k!θk\int_0^\infty x^k \frac{1}{\theta} e^{-x/\theta}\, dx = k!\,\theta^k. Jika kedua cara tidak cocok, ada kesalahan dalam penurunan MGF.

Cross-check PGF: Normalisasi Probabilitas

Evaluasi PGF di z=1z = 1 harus selalu menghasilkan 1 (jumlah seluruh probabilitas = 1):

PN(1)=k=0pk=1P_N(1) = \sum_{k=0}^{\infty} p_k = 1

Gunakan ini sebagai sanity check pertama sebelum menghitung momen. Jika PN(1)1P_N(1) \neq 1, PGF yang diberikan tidak valid atau ada kesalahan rekognisi distribusi.

Metode Alternatif

Untuk distribusi yang memiliki formula rekursif (kelas (a,b,0)(a,b,0) — lihat 2.2 (a,b,0) and (a,b,1) Distribution Classes), momen dapat dihitung secara rekursif tanpa menurunkan PGF secara eksplisit. Hal ini lebih efisien ketika distribusi frekuensi memenuhi relasi Panjer. Untuk distribusi dengan MGF bentuk tertutup, turunan langsung seringkali lebih cepat dari pendekatan CGF, kecuali untuk momen tingkat tinggi.

Section 6 — Visualisasi Mental

Diagram Konseptual MGF:

                  t kecil (dekat 0)
                       |
f_X(x) ──── E[e^{tX}] ──── M_X(t)
              "kapsul"          |
                           d^k/dt^k |_{t=0}
                                |
                           E[X^k]   ← momen ke-k

Diagram Konseptual PGF:

p_k = P(N=k)          z ∈ [0,1]
    |                      |
    └──── P_N(z) = Σ p_k z^k ────── PGF
                   |
         P_N^(k)(0) / k!      P_N^(k)(1)
              |                    |
           p_k                  factorial moment ke-k

Keterangan visual:

  • MGF “mengemas” distribusi kontinu: sumbu horizontal adalah tt, MGF naik cepat untuk distribusi heavy-tail
  • PGF adalah polinomial / power series di z[0,1]z \in [0,1]: nilai di z=1z=1 selalu 1 (normalisasi)
  • Kemiringan (slope) PGF di z=1z=1 adalah E[N]E[N] — secara grafis terlihat sebagai gradien kurva PGF tepat di ujung kanan

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Rumus
”Kapsul” MGF mengemas semua E[Xk]E[X^k]MX(t)=k=0E[Xk]k!tkM_X(t) = \sum_{k=0}^{\infty} \frac{E[X^k]}{k!} t^k
Turunan di t=0t=0 “membuka kapsul”E[Xk]=MX(k)(0)E[X^k] = M_X^{(k)}(0)
Slope PGF di z=1z=1E[N]=PN(1)E[N] = P_N'(1)
Kurvatura PGF di z=1z=1Terkait Var(N)\text{Var}(N) via PN(1)+PN(1)[PN(1)]2P_N''(1) + P_N'(1) - [P_N'(1)]^2
PN(1)=1P_N(1) = 1Normalisasi probabilitas

Section 7 — Jebakan Umum

Kesalahan Parametrisasi

Konvensi parameter yang sering tertukar:

  • Distribusi Gamma: ada yang menulis Gamma(α,β)\text{Gamma}(\alpha, \beta) dengan β=1/θ\beta = 1/\theta (rate), ada yang menulis Gamma(α,θ)\text{Gamma}(\alpha, \theta) dengan θ\theta sebagai scale. MGF berbeda: (1θt)α(1-\theta t)^{-\alpha} vs (1+t/β)α(1 + t/\beta)^{-\alpha}. Selalu verifikasi konvensi sebelum menulis MGF.
  • Distribusi Normal: MX(t)=eμt+12σ2t2M_X(t) = e^{\mu t + \frac{1}{2}\sigma^2 t^2}σ2\sigma^2 adalah variansi (bukan standar deviasi). Kesalahan umum: menulis eμt+σt2/2e^{\mu t + \sigma t^2/2} dengan σ\sigma sebagai SD — ini sudah benar jika σ=SD\sigma = \text{SD}, tapi pastikan definisi konsisten dengan soal.
Kesalahan Konseptual
  1. Evaluasi di titik salah: Momen diperoleh dari turunan MGF di t=0t = 0, bukan di t=1t = 1 atau nilai lain. Factorial moments diperoleh dari turunan PGF di z=1z = 1, bukan z=0z = 0.
  2. MGF vs PGF untuk frekuensi: Kedua fungsi bisa digunakan untuk distribusi diskrit, tapi PGF lebih natural. Jangan bingung antara MN(t)M_N(t) dan PN(z)P_N(z) — keduanya berbeda fungsi, meskipun terhubung via PN(z)=MN(lnz)P_N(z) = M_N(\ln z).
  3. Additivity: Yang additif adalah MGF (produk, bukan jumlah) untuk variabel independen: MX+Y(t)=MX(t)MY(t)M_{X+Y}(t) = M_X(t) \cdot M_Y(t). Bukan MX+Y(t)=MX(t)+MY(t)M_{X+Y}(t) = M_X(t) + M_Y(t).
  4. Eksistensi MGF: Pareto dan Lognormal tidak memiliki MGF yang finite untuk t>0t > 0. Jangan gunakan MGF untuk distribusi ini dalam analisis ekor.
Kesalahan Interpretasi Soal
  • Kata “momen ke-k” bisa berarti raw moment E[Xk]E[X^k] atau central moment E[(Xμ)k]E[(X-\mu)^k] — bedakan! Raw moment dari MGF, central moment butuh konversi.
  • Kata “factorial moment ke-k” spesifik berarti E[X(X1)(Xk+1)]E[X(X-1)\cdots(X-k+1)] — ini dari PGF, bukan raw moment biasa.
  • “Identifikasi distribusi dari MGF/PGF”: ini meminta pengenalan bentuk fungsi, bukan kalkulasi momen — hafalkan bentuk MGF/PGF distribusi standar.
Red Flags
  • Soal menyebut “compound distribution” atau “aggregate claims” → properti perkalian MGF kemungkinan digunakan (lihat 4.2 Compound Distributions)
  • Soal menyebut “characteristic function” → mirip MGF tapi menggunakan bilangan kompleks, bukan real — berbeda prosedur
  • Soal memberikan MGF dan meminta distribusinya → hafalkan bentuk kanonik: Eksponensial (1θt)1(1-\theta t)^{-1}, Gamma (1θt)α(1-\theta t)^{-\alpha}, Normal eμt+σ2t2/2e^{\mu t + \sigma^2 t^2/2}, Poisson eλ(et1)e^{\lambda(e^t - 1)}, Binomial (q+pet)n(q + pe^t)^n
  • Soal menyebut “cumulant generating function”KX(t)=lnMX(t)K_X(t) = \ln M_X(t), bukan MGF sendiri

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. MGF mendefinisikan momen:
E[Xk]=MX(k)(0)E[X^k] = M_X^{(k)}(0)
  1. CGF mendefinisikan kumulan (lebih efisien):
KX(t)=lnMX(t),κ1=E[X],κ2=Var(X)K_X(t) = \ln M_X(t), \quad \kappa_1 = E[X], \quad \kappa_2 = \text{Var}(X)
  1. MGF independen bersifat multiplikatif:
MX+Y(t)=MX(t)MY(t)(X,Y independen)M_{X+Y}(t) = M_X(t) \cdot M_Y(t) \quad (X, Y \text{ independen})
  1. PGF untuk probabilitas distribusi diskrit:
pk=PX(k)(0)k!,E[N]=PN(1),Var(N)=PN(1)+PN(1)[PN(1)]2p_k = \frac{P_X^{(k)}(0)}{k!}, \quad E[N] = P_N'(1), \quad \text{Var}(N) = P_N''(1) + P_N'(1) - [P_N'(1)]^2
  1. Hubungan MGF–PGF:
PX(z)=MX(lnz)P_X(z) = M_X(\ln z)

Kapan Digunakan

  • Menghitung momen distribusi klaim (mean, variance, skewness) secara efisien
  • Mengidentifikasi distribusi dari bentuk MGF/PGF yang diberikan
  • Menghitung probabilitas titik distribusi frekuensi dari PGF
  • Menentukan distribusi jumlah variabel acak independen (model agregat dasar)
  • Membuktikan bahwa sumlah Gamma/Eksponensial independen tetap dalam keluarga Gamma

Kapan TIDAK Boleh Digunakan

  • MGF tidak digunakan untuk distribusi heavy-tail khas klaim (Pareto, Lognormal) — gunakan momen langsung atau Laplace transform
  • PGF tidak digunakan untuk variabel acak kontinu
  • Jika distribusi tidak diketahui bentuknya, MGF tidak bisa langsung diidentifikasi tanpa tabel referensi

Quick Decision Tree

graph TD
    A["Perlu menghitung momen atau probabilitas?"] -->|"Momen"| B["Distribusi kontinu atau diskrit?"]
    A -->|"Probabilitas titik"| C["Distribusi diskrit?"]
    B -->|"Kontinu"| D["Apakah MGF ada? (bukan heavy-tail)"]
    B -->|"Diskrit"| E["Gunakan MGF atau PGF (keduanya valid)"]
    D -->|"Ya"| F["Turunkan MGF k kali, evaluasi t=0<br>E[X^k] = M(k)(0)"]
    D -->|"Tidak (Pareto, Lognormal)"| G["Hitung E[X^k] dari integral langsung"]
    E --> H["Gunakan PGF: turunkan k kali, evaluasi z=1<br>untuk factorial moments"]
    C -->|"Ya"| I["Gunakan PGF: p_k = P^(k)(0) / k!"]
    C -->|"Tidak"| J["Distribusi kontinu: tidak ada PGF"]

Follow-up Options
  1. “Berikan contoh soal variasi identifikasi distribusi dari bentuk MGF”
  2. “Jelaskan hubungan 1.1 Moment and Probability Generating Functions dengan 2.1 Frequency MGF and PGF
  3. “Buat flashcard 1-halaman untuk MGF distribusi standar yang sering keluar di TA2”

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