PENDUGAAN PARAMETER DISTRIBUSI GENERALIZED WEIBULL DENGAN

Download Metode pendugaan tersebut merupakan salah satu metode yang sangat populer ... Pendugaan Parameter Generalized Weibull dengan Menggunakan Me...

0 downloads 502 Views 98KB Size
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