Microsoft Word - 49_ro'fah - Revisi edit PENDUGA TIPE KERNEL.docx 458 ComTech Vol.2 No. 1 Juni 2011: 458-466 PENDUGAAN FUNGSI INTENSITAS PROSES POISSON PERIODIK DENGAN TREN FUNGSI PANGKAT MENGGUNAKAN METODE TIPE KERNEL Ro’fah Nur Rachmawati Jurusan Matematika, Fakultas Sains dan Teknologi, Binus University Jl. KH. Syahdan No. 9, Palmerah, Jakarta Barat 11480. rofah.nr@binus.ac.id ABSTRACT Stochastic process has an important role in many areas in everyday life, including the customer service process. The number of customers who come to a service center will be different for each particular time. A special form of stochastic process with continuous time and discrete state space is periodic Poisson process, which is a Poisson process with an intensity function of a periodic function. However, on the stochastic modeling of a phenomenon by a periodic Poisson process, the intensity function of the process is generally unknown. Therefore, a method is needed to infer the function. In this article, a Kernel estimator is formulated from a periodic Poisson process with a trend component in a rank function, which is divided into two cases; the identified rank function coefficient and the unidentified rank function coefficient. Keywords: stochastic process, periodic Poisson process, Kernel function, rank function trend ABSTRAK Proses stokastik memiliki peranan yang penting dalam berbagai bidang pada kehidupan sehari-hari, diantaranya pada proses pelayanan pelanggan. Banyaknya pelanggan yang datang pada suatu pusat servis akan berbeda untuk setiap waktu tertentu. Salah satu bentuk khusus dari proses stokastik dengan waktu kontinu dan ruang state diskret adalah proses Poisson periodik yaitu proses Poisson dengan fungsi intensitas berupa fungsi periodik. Namun pada pemodelan stokastik dari suatu fenomena yang dimodelkan dengan suatu proses Poisson periodik, fungsi intensitas dari proses tersebut umumnya tidak diketahui. Sehingga diperlukan suatu metode untuk menduga fungsi tersebut. Pada tulisan ini dirumuskan suatu penduga tipe Kernel dari suatu proses Poisson periodik dengan suatu komponen tren berbentuk fungsi pangkat, yang dibagi menjadi dua kasus, yaitu pada kasus pertama diasumsikan koefisien fungsi pangkat diketahui dan kasus kedua diasumsikan koefisien fungsi pangkat tidak diketahui. Kata kunci: Proses Stokastik, Proses Poisson Periodik, Fungsi Kernel, Tren Fungsi Pangkat Pendugaan Fungsi Intensitas… (Ro’fah Nur Rachmawati) 459 PENDAHULUAN Pemodelan stokastik melibatkan unsur peluang untuk menduga perilaku dari suatu sistem yang tidak diketahui dengan pasti pada masa yang akan datang. Sebagai contoh, banyaknya kendaraan yang melewati suatu ruas jalan raya akan berbeda untuk setiap waktu tertentu, dalam fenomena pelayanan pelanggan (costumer service) yaitu banyaknya pelanggan yang datang ke suatu pusat servis akan berbeda untuk setiap waktu tertentu. Salah satu proses stokastik yang menarik untuk dipelajari adalah proses Poisson Periodik, yaitu suatu proses Poisson dengan fungsi intensitas berupa fungsi periodik. Namun, pemodelan suatu fenomena yang dimodelkan dengan suatu proses Poisson periodik, fungsi intensitas dari proses tersebut umumnya tidak diketahui, sehingga diperlukan suatu metode untuk menduga fungsi tersebut. Tulisan ini mengkaji pendugaan tipe Kernel proses Poisson periodik yang memiliki peranan penting dalam berbagai bidang pada kehidupan sehari-hari. Pendugaan fungsi intensitas proses Poisson periodik tanpa melibatkan suatu komponen tren telah dilakukan kajiannya pada Mangku (2006), dan pendugaan dengan melibatkan suatu komponen tren linear telah dilakukan kajiannya pada Rachmawati (2008). Namun, jika fungsi intensitas pada proses Poisson periodik meningkat berdasarkan suatu fungsi pangkat terhadap waktu, model yang lebih tepat untuk digunakan adalah proses Poisson periodik dengan suatu komponen tren berbentuk fungsi pangkat. Dalam tulisan ini, pendugaan fungsi intensitas proses Poisson periodik dibagi menjadi dua kasus, yaitu untuk koefisien fungsi pangkat diketahui dan untuk koefisien fungsi pangkat yang tidak diketahui. Sehingga karya ilmiah ini bertujuan untuk menduga fungsi intensitas (lokal) di suatu titik 0, , dengan menggunakan sebuah realisasi tunggal dari proses Poisson periodik pada 0, (Mangku, 2006: 1) bagi kedua kasus. METODE Metode yang digunakan adalah studi pustaka, dengan dasar teori sebagai berikut: Misalkan , , … , adalah peubah acak dalam ruang peluang , , . Dikatakan bahwa barisan peubah acak konvergen dalam peluang ke , dinotasikan , jika untuk setiap 0, P | | 0 untuk ∞. (Grimmet and Stirzaker, 1992, p.274). Suatu statistik , , … , yang konvergen dalam peluang ke parameter , yaitu disebut penduga konsisten bagi (Hogg et al., 2005, p.204) Jika adalah peubah acak, maka untuk suatu 0, P | | E | | . (Ghahramani, 2005, p.476). Jika adalah peubah acak dengan nilai harapan dan ragam , untuk setiap 0 , P | | . (Ghahramani, 2005, p.477). Suatu titik s disebut titik Lebesgue dari suatu fungsi yang terintegralkan lokal, jika lim | | 0. (Wheeden RL, Zygmund A. 1997, p.239). Pendugaan fungsi intensitas proses Poisson periodik tanpa menggunakan suatu tren telah dikaji oleh Mangku (2006), yang menghasilkan: 460 ComTech Vol.2 No. 1 Juni 2011: 458-466 , 1 adalah penduga bagi di titik 0, dengan periode ( 0) diketahui (Mangku, 2006, p. 2). Pendugaan fungsi intensitas proses Poisson periodik dengan melibatkan suatu tren linear dengan menggunakan Kernel seragam telah dibahas oleh Rachmawati (2008), yang menghasilkan: , , 1 ⁄ 1 , 2 1 ⁄ adalah penduga bagi di titik 0, dengan periode ( 0) diketahui dan diketahui. Untuk kasus tidak diketahui diperoleh , , 1 ⁄ 1 , 2 , 1 ⁄ adalah penduga bagi di titik 0, dengan periode ( 0 ) diketahui, dengan , adalah penduga bagi yaitu , 0, . (Rachmawati, 2008, p.31). Misalkan adalah proses Poisson homogen pada inteval 0, ∞ dengan fungsi intensitas yang terintegralkan lokal dan terdiri atas dua komponen, yaitu komponen periodik dengan periode ( 0) diketahui dan suatu komponen tren berbentuk fungsi pangkat. Dengan kata lain untuk 0, ∞ fungsi intensitas dapat ditulis sebagai berikut dengan adalah komponen tren fungsi pangkat, 0 menyatakan kemiringan dari tren, dan tulisan ini dibatasi hanya untuk 0 1. Tulisan ini tidak mengasumsikan suatu bentuk parametrik dari , kecuali bahwa adalah periodik yaitu untuk setiap 0, ∞ dan berlaku (Mangku, 2006,p.1). Diasumsikan bahwa adalah titik Lebesgue bagi , yaitu lim | | 0. (Mangku, 2006, p.2). Misalkan untuk suatu Ω, hanya terdapat sebuah realisasi dari proses Poisson N yang terdefinisi pada suatu ruang peluang Ω, , P dengan fungsi intensitas seperti pada (1) yang diamati pada interval terbatas 0, 0, ∞ . Karena adalah fungsi periodik dengan periode , maka masalah menduga pada titik s dengan dapat direduksi menjadi masalah menduga pada titik s dengan 0, (Mangku, 2006,p.1). (1) (2) (3) Pendugaan Fungsi Intensitas… (Ro’fah Nur Rachmawati) 461 Misalkan : adalah suatu fungsi bernilai real, disebut Kernel, jika memenuhi sifat-sifat berikut: (K.1) K merupakan fungsi kepekatan peluang. (K.2) K terbatas. (K.3) K memiliki daerah definisi pada 1,1 . Dan adalah barisan dari bilangan real positif yang konvergen menuju nol, yaitu 0 untuk ∞ (Mangku, 2006, p. 2). Perlu diperhatikan bahwa, nilai asimtotik ∞ dengan nilai bukan menyatakan ukuran sample, namun menyatakan panjang interval dari peubah acak 0, yang menyatakan banyaknya kejadian (Mangku, 2006, p.2). Pada bahasan ini dikaji dua kasus, yaitu pertama diasumsikan bahwa kemiringan a diketahui sedangkan fungsi pada 0, tidak diketahui. Kedua diasumsikan fungsi intensitas global bagi diketahui, sedangkan kemiringan a dan fungsi pada 0, ∞ tidak diketahui. HASIL DAN PEMBAHASAN Hasil yang diperoleh pada pendugaan tipe Kernel dari pada 0, untuk kasus diketahui, dapat dirumuskan sebagai berikut: , , 1 , 1 , dengan , 1 , . Pembahasan mengenai penyusunan penduga tipe Kernel seperti pada persamaan (5) dapat diuraikan sebagai berikut: Pertama perlu diperhatikan bahwa, pendugaan tipe Kernel untuk fungsi intensitas proses Poisson periodik hanya didasarkan pada sebuah realisasi tunggal dari proses Poisson , dengan demikian diperlukan informasi yang cukup mengenai nilai pada interval 0, (Mangku, 2006, p.3). Dengan alasan ini, asumsi (2) memiliki peranan yang sangat penting dibalik penyusunan ide penduga , , untuk . Sehingga untuk setiap titik s dan dan dari (1) dan (2) maka . Dengan pemisalan seperti pada (5) maka persamaan di atas menjadi 1 , 1 1 , 1 , . (4) (5) (6) 462 ComTech Vol.2 No. 1 Juni 2011: 458-466 Pendekatan yang digunakan pada pendugaan fungsi intensitas lokal dari suatu proses Poisson di titik s adalah dengan menaksir nilai rata-rata dari banyaknya kejadian di sekitar titik s. Secara matematis, misalkan 0, ∞ dan 0, menyatakan banyaknya kejadian yang terjadi pada 0, , maka fungsi intensitas lokal di titik s dapat didekati dengan 1 2 , . Dengan demikian, nilai fungsi di titik s dapat didekati dengan nilai rataan dari banyaknya kejadian di sekitar s, yaitu pada interval , serta dengan menggunakan (3) maka (7) dapat ditulis , , 1 , 1 E , 2 , . Dengan mengganti E , dengan padanan stokastiknya yaitu , maka (8) dapat ditulis , , 1 , 1 , 2 , 1 , 1 1 2 I , , , 1 , 1 , dimana I , . Agar penduga yang diperoleh lebih umum maka digunakan fungsi Kernel umum K sehingga diperoleh penduga bagi fungsi intensitas proses Poisson periodik untuk kasus diketahui seperti pada persamaan (5). Selanjutnya dikaji untuk kasus tidak diketahui, pada kasus ini diperlukan penduga bagi dengan mengasumsikan bahwa fungsi intensitas global bagi diketahui. Pendekatan yang digunakan pada pendugaan fungsi intensitas global dari suatu proses Poisson adalah dengan menaksir nilai rata- rata dari banyaknya kejadian dalam selang 0, . Secara matematis penduga bagi fungsi intensitas global dapat dinyatakan dengan 0, . Dengan demikian fungsi intensitas global bagi yang merupakan nilai rata-rata dari pada 0, yaitu . Penduga dari a adalah: , 1 0, 1 . Sehingga dapat dihasilkan penduga tipe Kernel dari pada 0, untuk kasus tidak diketahui adalah (8) (9) (10) (7) Pendugaan Fungsi Intensitas… (Ro’fah Nur Rachmawati) 463 , , 1 , 1 , , dengan , ∑ dan . Dapat diperhatikan bahwa persamaan (10) sama seperti pada persamaan (5) dengan mengganti nilai menjadi penduganya , . Untuk mendapatkan penduga , cukup diperhatikan bahwa E 0, . Dengan mengganti E 0, dengan padanan stokastiknya yaitu 0, maka persamaan di atas menjadi 0, 1 atau dapat ditulis 1 0, 1 . Sehingga diperoleh penduga dari a yaitu , seperti pada (9). Lema: Misalkan fungsi intensitas seperti (1) dan terintegralkan lokal, dengan 1 menyatakan fungsi intensitas global bagi . Penduga , merupakan penduga yang konsisten bagi a, dengan Mean Square Error-nya adalah , 1 1 untuk ∞. Bukti: Berdasarkan definisi dari , maka , Bias , Var , dan Bias , E , . E , dapat dihitung sebagai berikut: E , E 1 0, 1 1 E 0, 1 . (11) (12) 464 ComTech Vol.2 No. 1 Juni 2011: 458-466 Karena fungsi intensitas seperti (1) maka persamaan diatas menjadi E , 1 1 1 1 1 1 . untuk ∞. Ragam dari , diperoleh dengan cara yang serupa yaitu: Var , Var 1 0, 1 1 Var 0, . Karena N adalah proses Poisson, maka persamaan di atas dapat ditulis Var , 1 E 0, 1 1 untuk ∞. Dengan mensubstitusikan persamaan (13) dan (14) ke (15) diperoleh , 1 1 untuk ∞ seperti pada persamaan (11). Dengan demikian Lema terbukti. Teorema (Kekonsistenan , ) Penduga , merupakan penduga konsisten bagi a, yaitu , untuk ∞. Bukti: Berdasarkan definisi akan ditunjukkan bahwa 0, P , 0 untuk ∞. Berdasarkan ketaksamaan segitiga, diperoleh (13) (14) (15) Pendugaan Fungsi Intensitas… (Ro’fah Nur Rachmawati) 465 , , E , E , . Berdasarkan persamaan (13) diperoleh, E , , berarti 0, ada N sehingga E , 2 , N. Dari persamaan di atas diperoleh P , P , E , 2 . Jadi untuk membuktikan (15) cukup ditunjukkan P , E , 2 0. Dengan ketaksamaan Chebychev, diperoleh P , E , 2 4Var , . Dari persamaan (14) diperoleh Var , 0 untuk ∞, sehingga diperoleh (16). Dengan demikian Teorema terbukti. Selanjutnya, pembahasan mengenai ide pembentukan penduga , , bagi untuk kasus tidak diketahui sama dengan kasus diketahui, yaitu dengan mengganti a pada kasus pertama dengan , . Sehingga hasil dari pendugaan , , pada kasus ini adalah: , , 1 , 1 , , dengan memisalkan , ∑ dan , , adalah penduga bagi a, K adalah suatu Kernel, dan adalah barisan bilangan real positif yang konvergen menuju nol, yaitu 0 untuk ∞. PENUTUP Dari kajian yang telah dilakukan, diperoleh: Penduga tipe Kernel bagi pada 0, untuk kasus diketahui adalah: , , 1 , 1 , dengan n adalah panjang interval pengamatan, , , ∑ , K adalah suatu Kernel, dan adalah barisan bilangan real positif yang konvergen menuju nol, yaitu 0 untuk ∞. Penduga tipe Kernel bagi pada 0, untuk kasus tidak diketahui adalah: (16) 466 ComTech Vol.2 No. 1 Juni 2011: 458-466 , , 1 , 1 , , dengan n adalah panjang interval pengamatan, , , ∑ , K adalah suatu Kernel, dan adalah barisan bilangan real positif yang konvergen menuju nol, yaitu 0 untuk ∞, serta , 1 0, 1 adalah penduga konsisten bagi a. Penduga , merupakan penduga yang konsisten bagi a, dengan Mean Square Error-nya adalah , 1 1 untuk ∞, dan fungsi intensitas seperti (1) dan terintegralkan lokal, dengan 1 menyatakan fungsi intensitas global bagi . DAFTAR PUSTAKA Ghahramani, S. (2005). Fundamental of Probability, (3rd ed.). New York: Prentice Hall. Hogg, R. V. & Craig, A. T. (2005). Introduction to Mathematical Statistics, (5th ed.). New Jersey: Prentice Hall, Engelwood Cliffs. Mangku, I. W. (2001). Estimating the Intensity of a Cyclic Poisson Process. Amsterdam: University of Amsterdam. Mangku, I. W. (2006). Asymptotic Normality of a Kernel-type Estimator for the Intensity of a Periodic Poisson Process. Journal of Mathematics and Its Applications, 5 (2), 13-22. Mangku, I. W. (2006). Weak and Strong Convergence of a Kernel-type Estimator for the Intensity of a Periodic Poisson Process. Journal of Mathematics and Its Applications, 5 (1), 1-12 Rachmawati, N. R. (2008). Sifat-sifat Statistika Penduga Fungsi Intensitas Proses Poisson Periodik dengan Tren Fungsi Pangkat. Skripsi tidak diterbitkan. Bogor: Fakultas Matematika dan Ilmu Pengetahuan Alam Institut Pertanian Bogor. Wheeden, R. L., Zygmund, A. (1997). Measure and Integral: An Introduction to Real Analysis. New York: Marcell Dekker.