Model Multi-State dalam Matematika
Model Multi-State dalam Matematika
Faihatuz Zuhairoh
DEFINISI
T ⊆ {(i, j)|i ̸= j; i, j ∈ S }
MODEL HIDUP-MATI
0:Hidup 1:Meninggal
−µ01 µ01 00 01
x+t x+t t px t px
Qx+t = t Px =
0 0 0 1
Untuk model hidup-mati, peluang (x) untuk tetap berada pada state 0
(hidup) selama periode t adalah
Z t
t px = exp − µx+s ds
0
1:Meninggal-Kecelakaan
0:Hidup
0: Sehat 1: Cacat
2: Meninggal
Hal penting dari model ini adalah bahwa kecacatan terjadi secara per-
manen maka tidak ada panah dari state 1 kembali ke state 0.
Matriks intensitas transisi dan probabilitas transisi diberikan sebagai ber-
ikut:
−µ01 02 µ01 µ02 00 01 02
x+t − µx+t x+t x+t t px t px t px
Qx+t 0 −µ12
x+t
12
µx+t t Px = 0 t x p 11 p 12
t x
0 0 0 0 0 1
Faihatuz Zuhairoh Model Multi-state 9 / 53
Jenis-jenis Model Multi-state
µ01
x+t
0: Sehat 1: Sakit
µ10
x+t
µ02
x+t
µ12
x+t
2: Meninggal
meninggalkan state 0.
Sehingga force of mortality-nya dapat dituliskan sebagai berikut.
Probabilita transisi
Z t Z t
(1) (τ ) (1) (2) (τ ) (2)
t qx = s px µx+s ds, t qx = s px µx+s ds
0 0
PROBABILITAS TRANSISI
Untuk rantai markov waktu diskrit, diasumsikan bahwa proses pada state
i pada waktu t , akan berpindah ke state j pada waktu t + 1 ditulis dengan
ptij .
SIFAT MARKOV
Untuk status terserap (absorbing state ), semua elemen pada baris adalah
0, kecuali diagonal utama. Diagonal utama sama dengan 1.
Contoh
Hitunglah P61
Jawab
Dari soal, jelas bahwa status 1 dan 2 menyerap (absorbing ). Untuk usia
(61) probabilitas dia tetap hidup setelah 1 tahun adalah
560(τ )
= 0.7
p61 =
800
Probabilitas bahwa dia akan meninggal karena sebab 1 dan sebab 2 dalam
tahun itu adalah
P61 = 0 1.00 0
0 0 1.00
PERSAMAAN CHAPMAN-KOLMOGOROV
Persamaan Chapman-Kolmogorov
t+s Px = t Px s Px+t
Contoh
Untuk model multi status dengan 3 status, Sehat (0), cacat (1), dan
Meninggal (2)
(i) Untuk k = 0, 1;
(00) (01) (10) (12)
px+k = 0.7, px+k = 0.2, px+k = 0.1, px+k = 0.25
(ii)Terdapat 100 orang di awal, semuanya sehat. Status masa depan
mereka independen.
Hitunglah variansi dari 100 orang tersebut yang meninggal dalam dua
tahun pertama.
Jawab
0.7 0.2 ?
Jawab
2 D0 ∼ B(100, 0.22)
Jawab
Ada dua jenis arus kas yang terkait dengan rantai Markov waktu diskrit,
yaitu:
1 Arus kas yang timbul karena berada dalam state tertentu.
Contoh
0.4 0.6
Terdapat rantai Markov waktu diskrit homogen dengan P = .
0.8 0.2
Ruang state-nya adalah E = {0, 1}. Subjek sekarang dalam state 0, dan
i = 5%.
a. Misalkan terdapat arus kas 1 untuk state 0 sekarang dan awal dua
tahun ke depan. Tentukan APV dari aliran arus kas.
b. Misalkan terdapat arus kas 1 jika transisi dari state 1 ke state 0
terjadi pada akhir tahun pertama, kedua dan ketiga. Tentukan
APV dari aliran arus kas.
Jawab
00
Pr(Y2 = 0|Y0 = 0) = 2 p0
= Pr(0 → 0 → 0) + Pr(0 → 1 → 0)
= 0.4 × 0.4 + 0.6 × 0.8
= 0.64
Jawab
b. Pr(Y0 = 1, Y1 = 0|Y0 = 0) = 0
Sama seperti dalam multi decrement, dalam input Rantai Markov Waktu
Kontinu adalah intensitas transisi. Misalkan ada n+ 1 status yang mung-
kin (n > 0) sedemikian sehingga ruang state E adalah 0, 1, 2, . . . , n.
Untuk i, j ∈ E dan t ≥ 0.
ij
h px
µijx+t = limh→0+ laju transisi sesaat-t dari state i ke j untuk i ̸= j .
h
Contoh
00 ,
Hitunglah 10 p60 01 02
10 p60 dan 10 p60 untuk model cacat permanen jika di-
ketahui µx = 0.028, µx = 0.023 dan µ12
01 02
x = 0.032
Jawab
Pada kasus ini, semua intensitas transisi konstan. Matriks Q nya adalah:
Z 10
00
10 p60 = exp − 0.051dt
0
= e −0.51
= 0.60050
Jawab
Z t
01 00 01 11
t p60 = s p60 µ60+s t−s p60+s ds
Z0 t
= e 0.051s 0.028e 0.032(t−s) ds
0
1 − e −0.019t
= 0.028e −0.032t
0.019
01 −0.32 1 − e −0.19
t p60 = 0.028e
0.019
= 0.18517
Karena 00 01 02 02
10 p60 +10 p60 +10 p60 =1, maka 10 p60 =0.21433
p12 p23
1: Susceptible 2: Infected 3: Recovered
Transisi setiap pasien diamati serta total waktu pasien berada di setiap
state juga diamati. Berdasarkan Gambar model multi-state SIR yang
perlu diketahui adalah
Sehingga dapat dihitung nilai µ12 dan µ23 dengan persamaan berikut.
N11 + N12 N22 + N23
µ12 = dan µ23 =
t1 t2
sehingga diperoleh matriks intensitas transisi (Q) berikut.
−0.01321 0.01321 0
MODEL SIRD
µ31
3: Sembuh
µ23
µ12
1: Rentan 2: Terinfeksi
µ24
µ14
4: Meninggal
p11 p12 p13 p14
p21 p22 p23 p24
Pij =
p31 p32 p33 p34
0 0 0 1
0
−(µ12 + µ14 ) µ12 µ14
0 −(µ23 + µ24 ) µ23 µ24
Qij =
0 −µ31 0
µ31
0 0 0 0
Dalam model multi-state, arus kas dapat muncul dari berbagai situasi.
Misalnya, dalam program asuransi cacat permanen, beberapa atau semua
hal berikut mungkin disediakan:
anuitas sementara tertanggung dinonaktifkan (state 1)
jumlah sekaligus pada saat menjadi cacat (transisi dari state 0 ke
state 1)
jumlah sekaligus pada saat kematian (transisi ke state 2)
Secara umum, terdapat arus kas berikut dalam model multi status:
Benets cash ows:
1. Arus kas yang berasal dari transisi tertentu (mis. Transisi dari
status 1 ke status 2)
2. Arus kas yang timbul karena meninggalkan status tertentu (mis.
Setiap kali ketika Y meninggalkan status 0)
3. Arus kas yang muncul saat memasuki status tertentu (mis. Setiap
kali ketika Y memasuki status 1)
Premium cash ows:
4. Aliran uang yang muncul ketika Y berada dalam status tertentu
(mis. Ketika Y berada pada status 0 pada tahun pertama atau
pada awal setiap tahun)
Berikut diasumsikan orang yang berusia (x) dalam status i pada waktu
0.
1. Misalkan bt dibayarkan pada waktu t jika transisi dari status k ke
status l terjadi pada waktu t dengan t ≤ n. Probabilitas untuk
mendapatkan bt dalam (t, t + dt) adalah probabilitas bahwa ada
transisi dari status k ke status l dalam (t, t + dt), yaitu, t pxik µkl
x+t .
Jadi APV-nya adalah
Z n
bt v t t pxik µkl
x+t dt
0
Demikian pula, jika anuitas dibayarkan pada awal setiap tahun, di-
peroleh:
n−1
ij
X
äx:n = v t t pxij
t=0
Contoh Soal
Terdapat model disability income dimana laju transisi konstan pada se-
mua usia dengan µ01 = 5%, µ02 = 3%, µ10 = 1% dan µ12 = 7, 5%.
Asumsikan bahwa suku bunga adalah 4%.
Hitung APV dari manfaat cacat sebesar 200 yang dibayarkan per tahun
secara kontinu untuk orang yang tepat berusia 40 tahun dan cacat, se-
lama periode cacat ini, dibayarkan hingga usia maksimum 60 tahun.
Jawab
Intensitas transisi yang keluar dari state cacat adalah 0.085. Probabilitas
untuk pembayaran adalah t px11 ¯
= e −0.085t Jadi, diperoleh
Z 20
APV = 200 e −0.04t e −0.085t dt
0
200
= (1 − e −2.5 )
0.125
= 1468.664
SOAL
SOAL
µ02
µ03 µ12
x+t
x+t x+t
3: Withdrawn 2: Meninggal
SOAL
P= 0 0.6 0.4
0 0 1
JAWABAN
1. 0.2590
2. 0.112994 dan 0.122699
px01 = 0 s px00 µ01 11
Rt
3. a. t x+s t−s px+s ds
03 t 00 03
R
b. t px = 0 s px µx+s ds
01 00 01 03
c. t px = 1 − t px − t px − t px
4. a. 280.5888
b. 762.3716
5. 215.0459