J. Sains MIPA, Edisi Khusus Tahun 2008, Vol. 14, No. 1, Hal.: 17 - 22 ISSN 1978-1873
PENDUGAAN PARAMETER DISTRIBUSI GENERALIZED WEIBULL DENGAN MENGGUNAKAN METODE KEMUNGKINAN MAKSIMUM Rani Sari Hermita, Warsono dan Dian Kurniasari Jurusan Matematika, FMIPA Universitas Lampung Jl. S. Brojonegoro No.1 Bandar Lampung 35145, Indonesia Diterima 28 Agustus 2007, perbaikan 10 Desember 2007, disetujui untuk diterbitkan 27 Desember 2007
ABSTRACT In this paper we discuss procedures of a maximum likelihood method in estimating parameters of generalized Weibull distribution. To assess unbiased properties we develop Monte Carlo Simulation. We show that how Newton-Raphson iterations can be utilized in the maximum likelihood method. We also demonstrate that the bias of maximum likelihood estimates of the parameters is smaller as the sample size increases, and the confident intervals is shorter as the sample size increases. Keywords:generalized Weibull distribution, maximum likelihood method, Monte Carlo Simulation, Newton-Raphson iterations
1. PENDAHULUAN Memilih model peluang terbaik dalam data kelangsungan hidup bukanlah sesuatu hal yang mudah untuk dilakukan. Satu pendekatan untuk mengatasi masalah ini adalah dengan menggunakan model-model umum (general models). Salah satu model umum yang dapat digunakan adalah model distribusi generalized Weibull karena memiliki potensi yang bagus untuk mencocokkan data kelangsungan hidup. Misalkan X adalah peubah acak dari distribusi generalized Weibull dengan tiga parameter, maka fungsi kepekatan peluang dari peubah acak tersebut adalah 1)
x
f ( x)
1 x
e
;
0, > 0, >0
(1)
dengan X = peubah acak yang didefinisikan sebagai waktu mati/rusak/gagal (failure time). = parameter lokasi (treshold) yang menujukkan lokasi waktu, dimana pada saat lokasi waktu tersebut belum ada obyek pengamatan yang mati/rusak/gagal maupun hilang. = parameter skala yang menunjukkan besarnya keragaman data distribusi generalized Weibull. = parameter bentuk yang menunjukkan laju kematian/kerusakan data distribusi generalized Weibull. Dimana fungsi kumulatif dari distribusi generalized Weibull sebagai berikut: x F ( x) 1 e
(2)
Untuk menduga parameter distribusi generalized Weibull pada penelitian ini digunakan metode kemungkinan maksimum (maximum likelihood method). Metode pendugaan tersebut merupakan salah satu metode yang sangat populer dalam pendugaan parameter suatu distribusi 2, 3). Metode kemungkinan maksimum didasarkan pada teori data yang berukuran besar. Tujuan penelitian ini adalah untuk mengkaji prosedur metode kemungkinan maksimum dalam menduga parameter distribusi generalized Weibull dan mengkaji sifat ketakbiasan dari dugaan yang dihasilkan dari metode kemungkinan maksimum.
2. METODE PENELITIAN Metode yang digunakan dalam penelitian ini mengikuti langkah-langkah sebagai berikut: 1. Menduga parameter distribusi generalized Weibull dengan metode kemungkinan maksimum sebagai berikut:
2008 FMIPA Universitas Lampung
17
Rani Sari Hermita dkk Pendugaan Parameter Distribusi Generalized Weibull
a. b. c.
2.
Membentuk fungsi kemungkinan yang berasal dari fungsi distribusi generaliized Weibull. Menurunkan fungsi kemungkinan distribusi generaliized Weibull dengan fungsi ln. Mencari turunan pertama dari ln fungsi kemungkinan terhadap parameter , , dan yang hendak di duga dan menyamaknnya dengan nol. d. Apabila solusi dari persamaan yang dihasilkan dari langkah c, prosedur pendugaan dilanjutkan dengan iterasi Newton-Raphson. Memverifikasi kinerja metode kemungkinan maksimum dalam menduga parameter distribusi generalized Weibull, melalui studi simulasi Monte Carlo4). Skenario simulasi Monte Carlo yang dilakukan adalah dengan mengambil kombinasi nilai parameter =3, ß=2 dan =1 dan ukuran sampel sebanyak n = 10, n = 30, dan n = 100 yang masingmasing diulang sebanyak 100 kali, dengan nilai toleransi masing-masing sebesar 0.1, 0.01, 0.001, dan 0.0001.
3. HASIL DAN PEMBAHASAN 3.1. Pendugaan Parameter Generalized Weibull dengan Menggunakan Metode Kemungkinan Maksimum Fungsi kemungkinan Persamaan (1) dari sebuah sampel acak dari n observasi pada distribusi generalized Weibull dapat dilihat sebagai berikut: n
Lx , ,
n
xi
1
xi
n n
e
n
n
1
1
xi
e
xi
i 1
i 1
i 1
kemudian dimaksimumkan dengan log atau dengan ln L, sehingga menjadi n
n
ln L x , ,
ln
n
n
n
1
1
xi
e
xi
i 1
i 1
sehingga diperoleh logaritma fungsi kemungkinan maksimum distribusi generalized Weibull n
ln L x , ,
n ln
n ln
n(
1) ln
1
1
ln xi i 1
n
(3)
xi i 1
Setelah dimaksimumkan dengan menggunakan ln L , maka diperoleh pendugaan parameter dari metode kemungkinan maksimum yaitu dengan mencari turunan pertama dari logaritma parameter-parameter yang akan diduga dan menyamakannya dengan nol, sehingga diperoleh: Untuk parameter n
1 i 1
1
n
1
xi
xi
(4)
0
i 1
Untuk parameter
1 n
1
n
(5)
xi i 1
Untuk parameter n
n
n
ln xi
n ln
i 1
xi
ln
ln xi
i 1
n
xi
0
(6)
i 1
Karena solusi Persamaan (6) di atas tidak dapat diselesaikan secara analitik, maka diperlukan pendekatan numerik. Dalam tulisan ini untuk menyelesaikan persamaan tersebut menggunakan metode iterasi Newton-Raphson, yaitu dengan mensubtitusikan Persamaan (4) ke dalam Persamaan (6) sehingga diperoleh:
g ,
0 n
g
,
n
ln xi i 1
n n
xi
n
xi
ln xi
0
(7)
i 1
i 1
18
2008 FMIPA Universitas Lampung
J. Sains MIPA, Edisi Khusus Tahun 2008, Vol. 14, No. 1
Seperti yang telah dijabarkan di atas, nilai penduga parameter Oleh karena itu sebelum dicari nilai penduga parameter
dan
dan
ditentukan melalui penduga parameter
.
, terlebih dahulu dihitung nilai penduga parameter .
Langkah-langkah metode Newton-Raphson bagi distribusi generalized Weibull adalah sebagai berikut 5): 1.
menentukan nilai awal
2.
menentukan persamaan fungsi g
0
yang ditentukan mendekati
sebagai fungsi skor Fisher yang diduga pada nilai
.
yang diperloleh dari persamaan (6) n
n
n
g
n
xi
ln xi
i 1
ln xi
(8) n
i 1
xi i 1
dan turunan pertamanya g '
g
g'
n
n
n
n
xi
ln xi
n
i1
ln xi i1
xi
2
n
i1
xi
n
n
xi
ln xi
i1
xi
ln xi
2
i1 n
xi i1
i1
diperoleh Persamaan (9) 3.
masukkan persamaan fungsi g dan turunan pertamanya ke dalam rumus metode Newton-Raphson sampai dengan error yang mana berdasarkan prosedur, berlaku:
g
g
g
(10)
0
0
0 0
sehingga dari Persamaan (9) dan (10) diperoleh dugaan parameter bagi
0
g'
1 0
g
g
0
0
g
0
,
sebagai berikut: (11)
0
0 0
selanjutnya untuk menyederhanakan perhitungan nilai penduga parameter kemungkinan maksimum bagi dituliskan sebagai berikut
g r 1
r
r
g
,
dapat
(12) r
r r
2008 FMIPA Universitas Lampung
19
Rani Sari Hermita dkk Pendugaan Parameter Distribusi Generalized Weibull
dengan nilai r selalu bertambah satu dari satu iterasi ke iterasi berikutnya. Dengan memberikan nilai awal r pada persamaan (12) maka akan diperoleh sebuah pendugaan yang lebih baik hinga perbedaan antara pendugaan dengan parameter yang sebenarnya mendekati nol. 3.2. Sifat Ketakbiasan Penduga Kemungkinan Maksimum pada Distribusi Generalized Weibull Hasil yang diperoleh dari simulasi adalah sebagai berikut. Seperti terlihat pada Tabel 1 untuk nilai toleransi sebesar 0.1, semakin besar ukuran sampel, nilai dugaan semakin mendekati nilai parameter sebenarnya. Dengan kata lain, bias dugaan semakin kecil dengan semakin besarnya ukuran sampel. Sementara itu, selang kepercayaan bagi dugaan parameter , , yang diperoleh dengan menggunakan metode kemungkinan maksimum cenderung semakin pendek dengan semakin meningkatnya ukuran sampel. Hasil yang serupa diperoleh untuk nilai toleransi masing-masing sebesar 0.01, 0.001, dan 0.0001 (Tabel 2, 3, dan 4. Bias dari dugaan parameter semakin kecil dengan semakin meningkatnya ukuran sampel dan selang kepercayaan yang dihasilkan semakin sempit dengan semakin besarnya ukuran sampel. Tabel 1. Nilai penduga parameter =3, =2 dan =1 distribusi generalized Weibull untuk ukuran sampel n = (10, 30, dan 100) dan nilai toleransi = 0.1 Penduga Parameter
Ukuran Sampel
n = 10
n = 30
n = 100
Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan(95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam
Lokasi ( )
Skala ( )
Bentuk ( )
3.2002 [2.8105 3.5898] 0.2002 0.0395 3.0565 [2.9324 3.1805] 0.0565 0.0040 3.0174 [2.9856 3.0491] 0.0174 2.6238e-004
1.5384 [0.3624 2.7144] -0.4616 0.3600 1.7984 [1.1508 2.4459] -0.2016 0.1091 1.9473 [1.5281 2.3666] -0.0527 0.0458
0.7329 [0.4793 0.9865] -0.2671 0.0167 0.8996 [0.6942 1.1051] -0.1004 0.0110 0.9663 [0.8173 1.1153] -0.0337 0.0058
Tabel 2. Nilai penduga parameter =3, =2 dan =1 distribusi generalized Weibull untuk ukuran sampel n = (10, 30, dan 100) dan nilai toleransi = 0.01
Ukuran Sampel
N = 10
N = 30
n = 100
20
Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan(95%) Bias Ragam
Lokasi ( ) 3.2002 [2.8105 3.5898] 0.2002 0.0395 3.0491 [2.9487 3.1495] 0.0491 0.0026 3.0175 [2.9847 3.0503] 0.0175 2.7998e-004
Penduga Parameter Skala ( ) 1.5373 [0.3617 2.7129] -0.4627 0.3598 1.8439 [1.1813 2.5065] -0.1561 0.1143 1.9347 [1.5198 2.3495] -0.0653 0.0448
Bentuk ( ) 0.7317 [0.4783 0.9850] -0.2683 0.0167 0.9052 [0.6740 1.1363] -0.0948 0.0139 0.9634 [0.8183 1.1086] -0.0366 0.0055
2008 FMIPA Universitas Lampung
J. Sains MIPA, Edisi Khusus Tahun 2008, Vol. 14, No. 1
Tabel 3. Nilai penduga parameter =3, =2 dan =1 distribusi generalized Weibull untuk ukuran sampel n = (10, 30, dan 100) dan nilai toleransi = 0.001. Ukuran Sampel
n = 10
n = 30
n = 100
Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam
Penduga Parameter Skala ( ) 1.5460 [0.5529 2.5392] -0.4540 0.2568 1.7872 [1.1244 2.4500] -0.2128 0.1144 1.9478 [1.5345 2.3610] -0.0522 0.0445
Lokasi ( ) 3.2113 [2.7996 3.6230] 0.2113 0.0441 3.0603 [2.9359 3.1847] 0.0603 0.0040 3.0209 [2.9776 3.0641] 0.0209 4.8731e-004
Bentuk ( ) 0.7474 [0.4822 1.0126] -0.2526 0.0183 0.8835 [0.6959 1.0711] -0.1165 0.0092 0.9633 [ 0.8274 1.0993] -0.0367 0.0048
Tabel 4. Nilai penduga parameter =3, =2 dan =1 distribusi generalized Weibull untuk ukuran sampel n = (10, 30, dan 100) dan nilai toleransi = 0.0001.
Ukuran Sampel
n = 10
n = 30
n = 100
Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam Nilai Tengah Selang Kepercayaan (95%) Bias Ragam
Penduga Parameter Skala ( ) 1.5048 [0.5586 2.4510] -0.4952 0.2330 1.8336 [1.0894 2.5777] -0.1664 0.1441 1.9240 [1.5211 2.3269] -0.0760 0.0422
Lokasi ( ) 3.2076 [2.8087 3.6065] 0.2076 0.0414 3.0672 [2.9308 3.2036] 0.0672 0.0048 3.0187 [2.9851 3.0522] 0.0187 2.9349e-004
Bentuk ( ) 0.7401 [0.4912 0.9890] -0.2599 0.0161 0.8987 [0.6886 1.1088] -0.1013 0.0115 0.9641 [0.8183 1.1100] -0.0359 0.0055
4. KESIMPULAN Berdasarkan hasil yang diperoleh, maka dapat diambil beberapa kesimpulan sebagai berikut: n
ln L x , ,
n ln
n ln
n(
1) ln
1
ln xi i 1
1
n
xi i 1
Karena logaritma fungsi kemungkinan distribusi generalized Weibull tidak dapat diselesaikan secara analitik, maka diatasi dengan menggunakan metode iterasi Newton-Raphson. Penduga parameter distribusi generalized Weibull dengan menggunakan metode kemungkinan maksimum mempunyai sifat ketakbiasan yang diperoleh untuk parameter ( , , ) mempunyai bias yang kecil untuk ukuran sampel besar, dengan kata lain semakin besar jumlah sampel yang digunakan, maka hasil pendugaannya semakin tak bias.
DAFTAR PUSTAKA 1.
Jhonshon, N.L. and Kotz, S. 1970. Continous Univariate Distribution. John Wiley, New York.
2.
Hogg, R.V. and Craig, A.T. 1995. Introduction to Mathematical Statistics. Fifth Edition. Prentice-hall Inc., New jersey.
2008 FMIPA Universitas Lampung
21
Rani Sari Hermita dkk Pendugaan Parameter Distribusi Generalized Weibull
3.
Hogg, R.V. and Tanis, E.A. 1997. Probability and Statistical Inference. Sixth Edition. Prentice-hall Inc., New jersey.
4.
Ross, S. M. 1997. Simulation. Second Edition. Academic Press, California.
5.
Munir, R.. 2003. Metode Numerik. Informatika, Bandung.
22
2008 FMIPA Universitas Lampung