AktuNotes
← Kembali
TA1 · Materi 1.3

Curtate Future Lifetime

Medium Bobot: 15–25% (bersama Topik 1) London (1997), Survival Models and Their Estimation, Bab 1–4; Dickson et al. (2009), Bab 2
TA1AnalisisSurvivalCurateFutureLifetimeVariabelAcakDiskrit

📊 1.3 — Curtate Future Lifetime

Ringkasan Cepat

Topik: Curtate Future Lifetime | Bobot: ~15–25% (Topik 1) | Difficulty: Medium Ref: London (1997) Bab 1–4; Dickson et al. (2009) Bab 2 | Prereq: 1.1 Survival and Lifetime Variables, 1.2 Survival and Hazard Functions


Section 0 — Pemetaan Topik

Topik TA1Sub-topik IDSkill DiujiBobotDifficultyPrerequisiteConnected TopicsReferensi
Analisis Survival1.3Mendefinisikan KxK_x; menghitung PMF, CDF, dan momen KxK_x; menghubungkan KxK_x dengan TxT_x15–25%Medium1.1 Survival and Lifetime Variables, 1.2 Survival and Hazard Functions1.4 Parametric Survival Models, 1.5 Censoring and Non-Parametric EstimationLondon (1997), Bab 1–4

Section 1 — Intuisi

Bayangkan sebuah perusahaan asuransi jiwa yang menjual polis kepada seseorang berusia 40 tahun. Aktuaria di perusahaan itu perlu menghitung berbagai nilai polis — premi, cadangan, dan manfaat — yang semuanya bergantung pada kapan tertanggung meninggal dunia. Dalam banyak produk asuransi jiwa tradisional, manfaat kematian dibayarkan di akhir tahun meninggal, bukan tepat pada saat kematian terjadi. Perbedaan ini bukan sekadar soal administrasi — ia mengubah cara kita memodelkan risiko secara matematis.

Di sinilah curtate future lifetime muncul sebagai konsep kunci. Jika sisa hidup seseorang dinyatakan dalam variabel TxT_x (waktu kontinu hingga kematian), maka KxK_x adalah versi “dibulatkan ke bawah” dari TxT_x — yaitu bilangan bulat yang menunjukkan berapa tahun penuh seseorang itu masih hidup setelah usia xx. Misalnya jika Tx=7,3T_x = 7{,}3 tahun, maka Kx=7K_x = 7: orang tersebut meninggal di tahun kedelapan, setelah menyelesaikan 7 tahun penuh.

Mengapa versi diskrit ini begitu berguna? Dalam praktiknya, banyak tabel mortalitas — seperti tabel mortalitas yang diterbitkan PAI atau TMI (Tabel Mortalitas Indonesia) — disajikan dalam interval satu tahun. Semua probabilitas disajikan per-tahun, bukan per-detik. Memodelkan KxK_x memungkinkan kita bekerja langsung dengan tabel mortalitas diskrit ini secara alami, tanpa perlu mengintegrasikan fungsi kontinu yang kompleks. Ini adalah jembatan antara teori survival kontinu dan aplikasi aktuaria yang berorientasi kalender tahunan.


Section 2 — Definisi Formal

Definisi Matematis

Diberikan variabel acak sisa hidup kontinu Tx0T_x \geq 0 untuk seseorang berusia xx, maka curtate future lifetime didefinisikan sebagai:

Kx=TxK_x = \lfloor T_x \rfloor

yaitu bagian bilangan bulat (floor) dari TxT_x. KxK_x adalah variabel acak diskrit dengan support k=0,1,2,3,k = 0, 1, 2, 3, \ldots

SimbolMaknaCatatan
TxT_xSisa hidup kontinu untuk seseorang berusia xxTx0T_x \geq 0, variabel acak kontinu
KxK_xCurtate future lifetime — sisa hidup diskritKx=Tx{0,1,2,}K_x = \lfloor T_x \rfloor \in \{0, 1, 2, \ldots\}
kkNilai realisasi dari KxK_xBilangan bulat non-negatif
kpx{}_{k}p_xProbabilitas seseorang usia xx masih hidup setelah kk tahunkpx=Sx(k)=P(Tx>k){}_{k}p_x = S_x(k) = P(T_x > k)
qx+kq_{x+k}Probabilitas seseorang usia x+kx+k meninggal dalam satu tahun ke depanqx+k=1px+kq_{x+k} = 1 - p_{x+k}
kqx{}_{k\|}q_xDeferred mortality probability — prob. meninggal tepat di tahun ke-(k+1)(k+1)kqx=kpxqx+k{}_{k\|}q_x = {}_{k}p_x \cdot q_{x+k}
exe_xCurtate expectation of life — ekspektasi KxK_xex=E[Kx]e_x = E[K_x]
2ex{}^2 e_xSecond moment pendukung perhitungan Var(Kx)\text{Var}(K_x)Dihitung dari E[Kx2]E[K_x^2]

Rumus Utama

PMF (Probability Mass Function) dari KxK_x:

P(Kx=k)=kpxqx+k=kqx,k=0,1,2,P(K_x = k) = {}_{k}p_x \cdot q_{x+k} = {}_{k\|}q_x, \quad k = 0, 1, 2, \ldots

Label: Probabilitas meninggal tepat di tahun ke-(k+1)(k+1), setelah bertahan kk tahun penuh.

CDF dari KxK_x:

P(Kxk)=1k+1px=k+1qx,k=0,1,2,P(K_x \leq k) = 1 - {}_{k+1}p_x = {}_{k+1}q_x, \quad k = 0, 1, 2, \ldots

Label: Probabilitas meninggal dalam k+1k+1 tahun pertama.

Survival function diskrit (tail probability):

P(Kx>k)=P(Kxk+1)=k+1px,k=0,1,2,P(K_x > k) = P(K_x \geq k+1) = {}_{k+1}p_x, \quad k = 0, 1, 2, \ldots

Label: Probabilitas masih hidup pada ulang tahun ke-(k+1)(k+1).

Ekspektasi (Curtate Expectation of Life):

ex=E[Kx]=k=0kkqx=k=1kpxe_x = E[K_x] = \sum_{k=0}^{\infty} k \cdot {}_{k\|}q_x = \sum_{k=1}^{\infty} {}_{k}p_x

Label: Rata-rata jumlah tahun penuh yang dijalani seseorang usia xx sebelum meninggal.

Momen kedua dan Variansi:

E[Kx2]=k=1(2k1)kpxE[K_x^2] = \sum_{k=1}^{\infty} (2k - 1) \cdot {}_{k}p_x Var(Kx)=E[Kx2](ex)2=k=1(2k1)kpx(k=1kpx)2\text{Var}(K_x) = E[K_x^2] - (e_x)^2 = \sum_{k=1}^{\infty} (2k-1) \cdot {}_{k}p_x - \left(\sum_{k=1}^{\infty} {}_{k}p_x\right)^2

Label: Alternatif: Var(Kx)=2k=1kkpxexex2\text{Var}(K_x) = 2\sum_{k=1}^{\infty} k \cdot {}_{k}p_x - e_x - e_x^2.

Hubungan KxK_x dengan TxT_x:

KxTx<Kx+1K_x \leq T_x < K_x + 1

Label: KxK_x adalah tahun penuh yang diselesaikan; fraksi tahun yang tersisa ada di interval [Kx,Kx+1)[K_x, K_x+1).

Asumsi Eksplisit

  1. Tabel mortalitas lengkap tersedia: Nilai kpx{}_{k}p_x atau qx+kq_{x+k} diketahui untuk semua k0k \geq 0.
  2. Support tak terbatas (ω = ∞) atau terbatas: Jika tabel mortalitas memiliki batas akhir ω\omega, maka support Kx{0,1,,ωx1}K_x \in \{0, 1, \ldots, \omega - x - 1\}.
  3. Tidak ada tambahan informasi (time-homogeneous): Probabilitas transisi hanya bergantung pada usia saat ini, bukan pada riwayat masa lalu.
  4. Pembulatan ke bawah (floor function): Kx=TxK_x = \lfloor T_x \rfloor, bukan pembulatan biasa.
  5. Independensi dari waktu kalender: Model mengasumsikan tabel mortalitas stasioner, tidak bergantung pada tahun kalender.

Section 3 — Jembatan Logika

Dari Definisi ke Rumus

Mengapa PMF KxK_x berbentuk kpxqx+k{}_{k}p_x \cdot q_{x+k}? Logikanya sederhana: agar Kx=kK_x = k terjadi, diperlukan dua hal sekaligus — (1) orang tersebut harus bertahan hidup selama kk tahun penuh (probabilitas kpx{}_{k}p_x), dan (2) orang tersebut harus meninggal dalam tahun ke-(k+1)(k+1), yaitu antara ulang tahun ke-kk dan ke-(k+1)(k+1) (probabilitas qx+kq_{x+k}). Karena kedua kejadian ini bersyarat dan sekuensial, probabilitasnya adalah perkalian dari keduanya.

Support dan Domain

KxK_x adalah variabel acak diskrit non-negatif. Nilai Kx=0K_x = 0 berarti orang tersebut meninggal sebelum ulang tahun berikutnya (meninggal di tahun pertama). Nilai Kx=kK_x = k berarti orang tersebut merayakan ulang tahun ke-kk setelah usia xx, tetapi meninggal sebelum ulang tahun ke-(k+1)(k+1).

Derivasi Rumus Ekspektasi Alternatif:

Langkah 1 — Mulai dari definisi:

ex=E[Kx]=k=0kP(Kx=k)=k=1kkqxe_x = E[K_x] = \sum_{k=0}^{\infty} k \cdot P(K_x = k) = \sum_{k=1}^{\infty} k \cdot {}_{k\|}q_x

Langkah 2 — Substitusi kqx=kpxk+1px{}_{k\|}q_x = {}_{k}p_x - {}_{k+1}p_x (karena kqx=kpxqx+k=kpxk+1px{}_{k\|}q_x = {}_{k}p_x \cdot q_{x+k} = {}_{k}p_x - {}_{k+1}p_x):

ex=k=1k(kpxk+1px)e_x = \sum_{k=1}^{\infty} k \cdot ({}_{k}p_x - {}_{k+1}p_x)

Langkah 3 — Lakukan Abel’s summation (summation by parts):

ex=k=1kkpxk=1kk+1pxe_x = \sum_{k=1}^{\infty} k \cdot {}_{k}p_x - \sum_{k=1}^{\infty} k \cdot {}_{k+1}p_x

Langkah 4 — Re-indeks jumlah kedua: misalkan j=k+1j = k+1, sehingga k=j1k = j-1:

k=1kk+1px=j=2(j1)jpx\sum_{k=1}^{\infty} k \cdot {}_{k+1}p_x = \sum_{j=2}^{\infty} (j-1) \cdot {}_{j}p_x

Langkah 5 — Gabungkan:

ex=k=1kkpxj=2(j1)jpx=1px+k=2[k(k1)]kpx=k=1kpxe_x = \sum_{k=1}^{\infty} k \cdot {}_{k}p_x - \sum_{j=2}^{\infty} (j-1) \cdot {}_{j}p_x = {}_{1}p_x + \sum_{k=2}^{\infty} [k - (k-1)] \cdot {}_{k}p_x = \sum_{k=1}^{\infty} {}_{k}p_x

Hasil: ex=k=1kpxe_x = \sum_{k=1}^{\infty} {}_{k}p_x — formula ini jauh lebih mudah dihitung dari tabel mortalitas.

Dilarang
  1. Jangan tulis P(Kx=k)=kqxP(K_x = k) = {}_{k}q_x — ini salah! kqx=P(Txk){}_{k}q_x = P(T_x \leq k) adalah CDF kontinu, bukan PMF diskrit.
  2. Jangan gunakan ex=E[Tx]e_x = E[T_x] — itu adalah complete expectation of life ex\overset{\circ}{e}_x, bukan curtate. Curtate exexe_x \leq \overset{\circ}{e}_x.
  3. Jangan abaikan perbedaan P(Kxk)=k+1qxP(K_x \leq k) = {}_{k+1}q_x vs P(Kx<k)=kqxP(K_x < k) = {}_{k}q_x — beda satu indeks bisa membatalkan seluruh perhitungan.

Section 4 — Contoh Soal

Soal A — Fundamental

Soal: Diketahui bahwa untuk seseorang berusia xx: px=0,97p_x = 0{,}97, px+1=0,95p_{x+1} = 0{,}95, px+2=0,92p_{x+2} = 0{,}92. Hitung P(Kx=0)P(K_x = 0), P(Kx=1)P(K_x = 1), P(Kx=2)P(K_x = 2), dan P(Kx3)P(K_x \geq 3).

Solusi Soal A

Pendekatan: Gunakan PMF P(Kx=k)=kpxqx+kP(K_x = k) = {}_{k}p_x \cdot q_{x+k} langsung.

1. Identifikasi Variabel

  • px=0,97p_x = 0{,}97, sehingga qx=0,03q_x = 0{,}03
  • px+1=0,95p_{x+1} = 0{,}95, sehingga qx+1=0,05q_{x+1} = 0{,}05
  • px+2=0,92p_{x+2} = 0{,}92, sehingga qx+2=0,08q_{x+2} = 0{,}08
  • 2px=pxpx+1=0,97×0,95=0,9215{}_{2}p_x = p_x \cdot p_{x+1} = 0{,}97 \times 0{,}95 = 0{,}9215
  • 3px=2pxpx+2=0,9215×0,92=0,847780{}_{3}p_x = {}_{2}p_x \cdot p_{x+2} = 0{,}9215 \times 0{,}92 = 0{,}847780

2. Identifikasi Distribusi / Model KxK_x berdistribusi diskrit dengan PMF kqx=kpxqx+k{}_{k\|}q_x = {}_{k}p_x \cdot q_{x+k}.

3. Setup Persamaan

P(Kx=k)=kpxqx+kP(K_x = k) = {}_{k}p_x \cdot q_{x+k}

4. Eksekusi Aljabar

P(Kx=0)=0pxqx=1×0,03=0,0300P(K_x = 0) = {}_{0}p_x \cdot q_x = 1 \times 0{,}03 = 0{,}0300 P(Kx=1)=1pxqx+1=0,97×0,05=0,0485P(K_x = 1) = {}_{1}p_x \cdot q_{x+1} = 0{,}97 \times 0{,}05 = 0{,}0485 P(Kx=2)=2pxqx+2=0,9215×0,08=0,073720P(K_x = 2) = {}_{2}p_x \cdot q_{x+2} = 0{,}9215 \times 0{,}08 = 0{,}073720 P(Kx3)=3px=0,847780P(K_x \geq 3) = {}_{3}p_x = 0{,}847780

5. Verification Jumlah: 0,0300+0,0485+0,073720+0,847780=1,0000{,}0300 + 0{,}0485 + 0{,}073720 + 0{,}847780 = 1{,}000 ✓ PMF valid.

Hasil: P(Kx=0)=0,030P(K_x=0) = 0{,}030; P(Kx=1)=0,0485P(K_x=1) = 0{,}0485; P(Kx=2)=0,07372P(K_x=2) = 0{,}07372; P(Kx3)=0,84778P(K_x \geq 3) = 0{,}84778.

Exam Tips — Soal A

Target waktu: 2 menit. Common trap: Lupa bahwa 0px=1{}_{0}p_x = 1 (selalu), sehingga P(Kx=0)=qxP(K_x=0) = q_x langsung. Shortcut: Hitung kpx{}_{k}p_x secara kumulatif sebagai produk berurutan dari px,px+1,p_{x}, p_{x+1}, \ldots


Soal B — Exam-Typical

Soal: Untuk seseorang berusia 60 tahun, diketahui q60=0,020q_{60} = 0{,}020, q61=0,022q_{61} = 0{,}022, q62=0,025q_{62} = 0{,}025, q63=0,028q_{63} = 0{,}028. Hitung curtate expectation of life e60e_{60} menggunakan data mortalitas ini, dengan asumsi kp60=0{}_{k}p_{60} = 0 untuk k5k \geq 5 (populasi kecil, disederhanakan). Juga hitung Var(K60)\text{Var}(K_{60}).

Solusi Soal B

Pendekatan: Hitung terlebih dahulu semua kp60{}_{k}p_{60}, lalu gunakan formula e60=k=14kp60e_{60} = \sum_{k=1}^{4} {}_{k}p_{60} dan Var(K60)=2k=14kkp60e60e602\text{Var}(K_{60}) = 2\sum_{k=1}^{4} k \cdot {}_{k}p_{60} - e_{60} - e_{60}^2.

1. Identifikasi Variabel

  • q60=0,020q_{60}=0{,}020, q61=0,022q_{61}=0{,}022, q62=0,025q_{62}=0{,}025, q63=0,028q_{63}=0{,}028
  • p60=0,980p_{60}=0{,}980, p61=0,978p_{61}=0{,}978, p62=0,975p_{62}=0{,}975, p63=0,972p_{63}=0{,}972

2. Identifikasi Distribusi / Model K60K_{60} diskrit, support {0,1,2,3,4}\{0,1,2,3,4\} (karena 5p60=0{}_{5}p_{60} = 0). Gunakan formula e60=k=14kp60e_{60} = \sum_{k=1}^{4} {}_{k}p_{60}.

3. Setup Persamaan

kp60=j=0k1p60+j{}_{k}p_{60} = \prod_{j=0}^{k-1} p_{60+j} e60=k=14kp60,Var(K60)=2k=14kkp60e60e602e_{60} = \sum_{k=1}^{4} {}_{k}p_{60}, \quad \text{Var}(K_{60}) = 2\sum_{k=1}^{4} k \cdot {}_{k}p_{60} - e_{60} - e_{60}^2

4. Eksekusi Aljabar

1p60=0,980{}_{1}p_{60} = 0{,}980 2p60=0,980×0,978=0,958440{}_{2}p_{60} = 0{,}980 \times 0{,}978 = 0{,}958440 3p60=0,958440×0,975=0,934479{}_{3}p_{60} = 0{,}958440 \times 0{,}975 = 0{,}934479 4p60=0,934479×0,972=0,908314{}_{4}p_{60} = 0{,}934479 \times 0{,}972 = 0{,}908314 e60=0,980000+0,958440+0,934479+0,908314=3,781233e_{60} = 0{,}980000 + 0{,}958440 + 0{,}934479 + 0{,}908314 = 3{,}781233

Untuk variansi:

2k=14kkp60=2[1(0,980)+2(0,95844)+3(0,934479)+4(0,908314)]2\sum_{k=1}^{4} k \cdot {}_{k}p_{60} = 2[1(0{,}980) + 2(0{,}95844) + 3(0{,}934479) + 4(0{,}908314)] =2[0,980+1,91688+2,803437+3,633256]=2×9,333573=18,667146= 2[0{,}980 + 1{,}91688 + 2{,}803437 + 3{,}633256] = 2 \times 9{,}333573 = 18{,}667146 Var(K60)=18,6671463,781233(3,781233)2=18,6671463,78123314,297723=0,588190\text{Var}(K_{60}) = 18{,}667146 - 3{,}781233 - (3{,}781233)^2 = 18{,}667146 - 3{,}781233 - 14{,}297723 = 0{,}588190

5. Verification e603,78<4e_{60} \approx 3{,}78 < 4: masuk akal untuk support maksimum 4. Var(K60)>0\text{Var}(K_{60}) > 0 ✓. Juga cek: e60max(K60)=4e_{60} \leq \max(K_{60}) = 4 ✓.

Hasil: e603,781e_{60} \approx 3{,}781 tahun; Var(K60)0,588\text{Var}(K_{60}) \approx 0{,}588.

Exam Tips — Soal B

Target waktu: 4 menit. Common trap: Formula variansi Var(Kx)=2kkpxexex2\text{Var}(K_x) = 2\sum k \cdot {}_{k}p_x - e_x - e_x^2 — jangan lupa kurangi exe_x (bukan hanya ex2e_x^2). Shortcut: Hitung kpx{}_{k}p_x secara kumulatif satu kali, simpan hasilnya, lalu gunakan untuk exe_x dan komponen variansi sekaligus.


Soal C — Challenging

Soal: Asumsikan De Moivre dengan ω=100\omega = 100 dan seseorang berusia x=40x = 40. Di bawah model ini, kp40=60k60{}_{k}p_{40} = \frac{60 - k}{60} untuk k=0,1,,60k = 0, 1, \ldots, 60. Hitung: (a) e40e_{40}, (b) Var(K40)\text{Var}(K_{40}), dan (c) P(K40>e40)P(K_{40} > e_{40}).

Solusi Soal C

Pendekatan: Gunakan formula deret aritmetik untuk k=16060k60\sum_{k=1}^{60} \frac{60-k}{60}, lalu hitung momen kedua, dan evaluasi tail probability.

1. Identifikasi Variabel

  • Model De Moivre: T40Uniform[0,60]T_{40} \sim \text{Uniform}[0, 60]
  • kp40=60k60{}_{k}p_{40} = \frac{60-k}{60}, untuk k=1,2,,59k = 1, 2, \ldots, 59 (untuk k=60k=60: 60p40=0{}_{60}p_{40} = 0)

2. Identifikasi Distribusi / Model K40DiskritK_{40} \sim \text{Diskrit}, support {0,1,,59}\{0, 1, \ldots, 59\}. Karena T40U[0,60]T_{40} \sim U[0,60], PMF K40K_{40} seragam: P(K40=k)=160P(K_{40} = k) = \frac{1}{60} untuk k=0,1,,59k = 0, 1, \ldots, 59.

3. Setup Persamaan

e40=k=159kp40=k=15960k60e_{40} = \sum_{k=1}^{59} {}_{k}p_{40} = \sum_{k=1}^{59} \frac{60-k}{60}

4. Eksekusi Aljabar

(a) Curtate Expectation:

e40=160k=159(60k)=160j=159j=16059×602=592=29,5e_{40} = \frac{1}{60}\sum_{k=1}^{59}(60-k) = \frac{1}{60}\sum_{j=1}^{59} j = \frac{1}{60} \cdot \frac{59 \times 60}{2} = \frac{59}{2} = 29{,}5

(b) Variansi — hitung E[K402]E[K_{40}^2] terlebih dahulu:

Karena K40Diskrit Seragam{0,1,,59}K_{40} \sim \text{Diskrit Seragam}\{0, 1, \ldots, 59\}:

E[K402]=160k=059k2=16059×60×1196=59×1196=702161170,167E[K_{40}^2] = \frac{1}{60}\sum_{k=0}^{59} k^2 = \frac{1}{60} \cdot \frac{59 \times 60 \times 119}{6} = \frac{59 \times 119}{6} = \frac{7021}{6} \approx 1170{,}167 Var(K40)=E[K402](e40)2=1170,167(29,5)2=1170,167870,25=299,917\text{Var}(K_{40}) = E[K_{40}^2] - (e_{40})^2 = 1170{,}167 - (29{,}5)^2 = 1170{,}167 - 870{,}25 = 299{,}917

Catatan: Untuk distribusi diskrit seragam {0,,n1}\{0,\ldots,n-1\}: Var=n2112\text{Var} = \frac{n^2-1}{12}. Dengan n=60n = 60: 3600112=359912299,917\frac{3600-1}{12} = \frac{3599}{12} \approx 299{,}917 ✓.

(c) P(K40>29,5)=P(K4030)P(K_{40} > 29{,}5) = P(K_{40} \geq 30):

P(K4030)=30p40=603060=3060=0,5P(K_{40} \geq 30) = {}_{30}p_{40} = \frac{60 - 30}{60} = \frac{30}{60} = 0{,}5

5. Verification e40=29,5=ωx12=592e_{40} = 29{,}5 = \frac{\omega - x - 1}{2} = \frac{59}{2} ✓ (formula curtate expectation De Moivre). P(K40>e40)=0,5P(K_{40} > e_{40}) = 0{,}5: masuk akal untuk distribusi simetris seragam — setengah dari populasi melewati ekspektasi. Bandingkan dengan complete expectation: e40=ωx2=30>e40=29,5\overset{\circ}{e}_{40} = \frac{\omega - x}{2} = 30 > e_{40} = 29{,}5 ✓ (curtate selalu lebih kecil).

Hasil: e40=29,5e_{40} = 29{,}5; Var(K40)299,92\text{Var}(K_{40}) \approx 299{,}92; P(K40>e40)=0,5P(K_{40} > e_{40}) = 0{,}5.

Exam Tips — Soal C

Target waktu: 5 menit. Common trap: Pada De Moivre, K40K_{40} adalah diskrit seragam {0,,59}\{0,\ldots,59\} (bukan {1,,60}\{1,\ldots,60\}) — perhatikan batas bawah selalu 0. Shortcut: Kenali langsung e40=ωx12e_{40} = \frac{\omega - x - 1}{2} untuk De Moivre — ini formula jadi yang harus dihafalkan.


Section 5 — Verifikasi & Sanity Check

Check 1 — Jumlah PMF = 1

Untuk support terbatas {0,1,,m}\{0, 1, \ldots, m\}:

k=0mP(Kx=k)=k=0mkqx=1\sum_{k=0}^{m} P(K_x = k) = \sum_{k=0}^{m} {}_{k\|}q_x = 1

Ini setara dengan 0pxm+1px=10=1{}_{0}p_x - {}_{m+1}p_x = 1 - 0 = 1 (jika m+1px=0{}_{m+1}p_x = 0). Selalu verifikasi ini setelah menghitung seluruh PMF.

Check 2 — Hubungan Curtate vs Complete Expectation

Selalu berlaku:

exexex+1e_x \leq \overset{\circ}{e}_x \leq e_x + 1

Di bawah UDD (Uniform Distribution of Deaths dalam satu tahun): ex=ex+12\overset{\circ}{e}_x = e_x + \frac{1}{2}. Jika hasil exe_x lebih besar dari ex\overset{\circ}{e}_x, ada kesalahan.

Metode Alternatif

Selain formula ex=k=1kpxe_x = \sum_{k=1}^{\infty} {}_{k}p_x, ekspektasi dapat dihitung langsung via PMF:

ex=k=0kkqx=0qx+1pxqx+1+22pxqx+2+e_x = \sum_{k=0}^{\infty} k \cdot {}_{k\|}q_x = 0 \cdot q_x + 1 \cdot p_x q_{x+1} + 2 \cdot {}_{2}p_x q_{x+2} + \cdots

Keduanya equivalen. Formula kpx\sum {}_{k}p_x lebih efisien jika tabel kpx{}_{k}p_x sudah tersedia, sedangkan formula PMF berguna jika distribusi mortalitas tahunan yang diketahui.


Section 6 — Visualisasi Mental

Garis Waktu Diskrit KxK_x:

Usia x        x+1       x+2       x+3       ...
  |-----+-----|---------|---------|-------->
  0     T?    1         2         3
  
Jika T_x = 2.7:   K_x = floor(2.7) = 2
                   Meninggal di tahun ke-3, setelah 2 tahun penuh

K_x = 0:  Meninggal di tahun pertama (antara usia x dan x+1)
K_x = 1:  Meninggal di tahun kedua  (antara usia x+1 dan x+2)
K_x = k:  Meninggal di tahun ke-(k+1) (antara usia x+k dan x+k+1)

PMF K_x (untuk data hipotetis):

P(K_x=k)
  |
  |  ***
  | *   *
  |*     *
  |       *
  |        *
  |         *
  +-+--+--+--+--+---> k
    0  1  2  3  4  5

Bentuk: biasanya unimodal atau monoton turun (bergantung pada
hazard function). Untuk populasi muda, monoton menurun ringan.

Distribusi kontinu TxT_x vs diskrit KxK_x:

f(t) — densitas T_x (kontinu):
  |     ____
  |    /    \
  |   /      \
  |  /        \___
  +-+--+--+--+--+-> t
    0  1  2  3  4

P(K_x=k) — PMF K_x (diskrit): area di bawah f(t) pada [k, k+1)
  |
  |  []  []
  |  []  []  []
  |  []  []  []  []
  +-+--+--+--+--+--> k
    0   1   2   3

Hubungan Visual ↔ Rumus

Elemen VisualKomponen Rumus
Panjang interval [k,k+1)[k, k+1)Satu tahun polis / tabel mortalitas
Area di bawah f(t)f(t) pada [k,k+1)[k, k+1)P(Kx=k)=kpxqx+kP(K_x = k) = {}_{k}p_x \cdot q_{x+k}
Luas kumulatif [0,k+1)[0, k+1)P(Kxk)=k+1qxP(K_x \leq k) = {}_{k+1}q_x
Luas sisa [k+1,)[k+1, \infty)P(Kx>k)=k+1pxP(K_x > k) = {}_{k+1}p_x
Pusat massa distribusi PMFex=E[Kx]e_x = E[K_x]

Section 7 — Jebakan Umum

Kesalahan Parametrisasi — Indeks pada CDF

Salah: P(Kxk)=kqxP(K_x \leq k) = {}_{k}q_x Benar: P(Kxk)=k+1qxP(K_x \leq k) = {}_{k+1}q_x

Penjelasan: KxkK_x \leq k berarti meninggal paling lambat di tahun ke-(k+1)(k+1), jadi melibatkan k+1k+1 periode, bukan kk. Ini sering menyebabkan kesalahan off-by-one yang fatal di soal.

Kesalahan Konseptual
  1. KxK_x vs TxT_x: ex=E[Kx]e_x = E[K_x] adalah curtate expectation, berbeda dengan complete expectation ex=E[Tx]\overset{\circ}{e}_x = E[T_x]. Di bawah UDD: ex=ex+0,5\overset{\circ}{e}_x = e_x + 0{,}5.
  2. Floor vs Round: Kx=TxK_x = \lfloor T_x \rfloor menggunakan floor (bulatkan ke bawah), bukan pembulatan biasa. Jika Tx=2,9T_x = 2{,}9, maka Kx=2K_x = 2, bukan 3.
  3. Support dimulai dari 0: Kx=0K_x = 0 adalah nilai yang valid dan berarti (meninggal di tahun pertama). Jangan mulai penjumlahan dari k=1k=1.
  4. PMF vs Tail: kpx=P(Kxk){}_{k}p_x = P(K_x \geq k), bukan P(Kx=k)P(K_x = k). PMF adalah kpxk+1px=kqx{}_{k}p_x - {}_{k+1}p_x = {}_{k\|}q_x.
Kesalahan Interpretasi Soal
  • “Probabilitas meninggal sebelum ulang tahun ke-(k+1)(k+1)P(Kxk)=k+1qxP(K_x \leq k) = {}_{k+1}q_x (inklusif: meninggal tepat di tahun ke-(k+1)(k+1) termasuk).
  • “Probabilitas meninggal tepat di tahun ke-(k+1)(k+1)P(Kx=k)=kqxP(K_x = k) = {}_{k\|}q_x.
  • “Probabilitas masih hidup setelah kk tahun”kpx=P(Kxk){}_{k}p_x = P(K_x \geq k)strictly greater than perlu P(Kx>k)=k+1pxP(K_x > k) = {}_{k+1}p_x.
Red Flags
  • Soal menyebut “manfaat dibayar di akhir tahun kematian” → gunakan KxK_x, bukan TxT_x.
  • Soal memberikan tabel mortalitas tahunan (qxq_x, qx+1q_{x+1}, …) → ini adalah konteks diskrit, pakai KxK_x.
  • Soal menyebut De Moivre atau TxU[0,ωx]T_x \sim U[0, \omega-x] → gunakan formula ex=ωx12e_x = \frac{\omega-x-1}{2} langsung.
  • Soal meminta variansi KxK_x → ingat formula Var=2kkpxexex2\text{Var} = 2\sum k \cdot {}_{k}p_x - e_x - e_x^2, jangan hanya hitung E[Kx2](E[Kx])2E[K_x^2] - (E[K_x])^2 tanpa mengetahui formula efisiennya.

Section 8 — Ringkasan Eksekutif

Must-Remember
  1. Definisi: Kx=TxK_x = \lfloor T_x \rfloor — bagian bilangan bulat dari sisa hidup kontinu.

  2. PMF:

P(Kx=k)=kpxqx+k=kqxP(K_x = k) = {}_{k}p_x \cdot q_{x+k} = {}_{k\|}q_x
  1. CDF: P(Kxk)=k+1qxP(K_x \leq k) = {}_{k+1}q_x; Tail: P(Kx>k)=k+1pxP(K_x > k) = {}_{k+1}p_x

  2. Curtate expectation (dua formula ekuivalen):

ex=k=1kpx=k=0kkqxe_x = \sum_{k=1}^{\infty} {}_{k}p_x = \sum_{k=0}^{\infty} k \cdot {}_{k\|}q_x
  1. Variansi:
Var(Kx)=2k=1kkpxex(ex)2\text{Var}(K_x) = 2\sum_{k=1}^{\infty} k \cdot {}_{k}p_x - e_x - (e_x)^2
  1. De Moivre shortcut: ex=ωx12e_x = \dfrac{\omega - x - 1}{2}; KxDiskrit Seragam{0,,ωx1}K_x \sim \text{Diskrit Seragam}\{0,\ldots,\omega-x-1\}

Kapan Digunakan

  • Soal menyebut manfaat/premi dibayar di akhir tahun (whole life, term life — annual payment)
  • Data mortalitas diberikan dalam tabel tahunan (qxq_x, pxp_x)
  • Pertanyaan tentang ekspektasi hidup curtate atau variansi sisa hidup diskrit
  • Model survival De Moivre dengan batas usia ω\omega
  • Sebagai langkah awal sebelum menghitung APV (Actuarial Present Value) asuransi jiwa

Kapan TIDAK Boleh Digunakan

  • Ketika manfaat dibayar saat kematian (tepat waktu) → gunakan TxT_x kontinu dan ex\overset{\circ}{e}_x
  • Ketika dibutuhkan model hazard rate kontinu μx\mu_x → lihat 1.2 Survival and Hazard Functions
  • Ketika data berbentuk kontinu atau diketahui distribusi parametrik lengkap → pertimbangkan TxT_x langsung

Quick Decision Tree

graph TD
    A["Soal melibatkan sisa hidup seseorang usia x"] --> B{"Manfaat/pembayaran\nberbasis waktu?"}
    B -->|"Akhir tahun kematian"| C["Gunakan K_x = floor(T_x)"]
    B -->|"Tepat saat kematian"| D["Gunakan T_x kontinu"]
    C --> E{"Apa yang dicari?"}
    E -->|"P(K_x = k)"| F["PMF: kpx * q(x+k)"]
    E -->|"P(K_x <= k)"| G["CDF: (k+1)qx"]
    E -->|"E[K_x]"| H["ex = sum kpx dari k=1"]
    E -->|"Var(K_x)"| I["Var = 2*sum(k*kpx) - ex - ex^2"]
    H --> J{"Model De Moivre?"}
    J -->|"Ya, omega diketahui"| K["ex = (omega - x - 1) / 2"]
    J -->|"Tidak"| L["Hitung kumulatif kpx dari tabel"]

Follow-up Options
  1. “Berikan contoh soal variasi dengan tabel mortalitas TMI (Tabel Mortalitas Indonesia)”
  2. “Jelaskan hubungan 1.3 Curtate Future Lifetime dengan 1.4 Parametric Survival Models
  3. “Buat flashcard 1-halaman untuk topik ini”

📖 Ref: London (1997), Survival Models and Their Estimation, Bab 1–4; Dickson et al. (2009), Bab 2 | 🗓️ 2026-04-19 | #TA1 #AnalisisSurvival #CurateFutureLifetime