Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
OPTIMASI MULTIRESPON DENGAN MENGGUNAKAN METODE GABUNGAN DATA ENVELOPMENT ANALYSIS (DEA) AGGRESSIVE DAN RESPONSE SURFACE (STUDI KASUS: PT.PHILLIPS INDONESIA) Rendra Erdkhadifa1, Sony Sunaryo2dan Muhammad Sjahid Akbar3 1)Program Studi Magister Statistika, Institut Teknologi Sepuluh Nopember e-mail:
[email protected] 2,3) Jurusan Statistika, Institut Teknologi Sepuluh Nopember Jl Raya ITS - Surabaya, 60111, Indonesia ABSTRAK Data Envelopment Analysis (DEA) merupakan teknik pemrograman matematika nonparametrik yang digunakan untuk mengevaluasi efisiensi relatif dari kehomogenan Decision Making Units (DMUs) dengan multiple input dan multiple output. Secara kontras, formula aggresive meningkatkan perbedaan efisiensi dengan menghendaki bahwa efisiensi dapat bernilai lebih dari satu dan memperhatikan “peer-evaluation” yang ditunjukkan dalam nilai cross-efficiencies. Berdasarkan nilai efisiensi tersebut kemudian dilakukan optimasi dengan menggunakan metode response surface hingga didapatkan level optimum. Kedua metode tersebut akan diterapkan pada kasus di PT.Phillips Indonesia dimana pool distance dan constant mounting length diduga berpengaruh terhadap efisiensi. Efisiensi diperoleh dari hasil perhitungan yang melibatkan variabel lumen, wattage, dan life time. Berdasarkan hasil analisis didapatkan bahwa terdapat 3 DMUs yang dapat dikatakan efisien dan memiliki lebih dari 1. Jika dibandingkan dengan penelitian acuan, untuk kriteria respon larger the better DEA aggressive menghasilkan nilai yang lebih baik. Namun untuk kriteria respon smaller the better penelitian acuan menghasilkan nilai yang lebih baik. Kata kunci : Cross-efficiencies, Data Envelopment Analysis Aggressive, Response Surface PENDAHULUAN PT. Phillips Indonesia merupakan perusahaan terbesar yang memproduksi lampu pijar dan neon di Indonesia. Salah satu tujuan dari perusahaan ini adalah mengembangkan jalur distribusi yang merata di seluruh Indonesia dengan mengimplementasikan pemasaran untuk mendukung kampanye 'kampung terang hemat energi' (Ramitha, 2008). Pada tahun 2008 perusahaan tersebut juga memproduksi 'mini tornado' yang merupakan lampu hemat energi berkualitas atau compact fluorescent (CFL). Salah satu bentuk kemajuan teknologi yang digunakan dalam perusahaan ini adalah satu mesin yang digunakan dalam kecepatan rendah mampu mempro-duksi lampu pijar sebanyak 2000 lampu/jam. Permasalahan yang terjadi di PT Philips Indonesia yang dipaparkan dalam penelitian (Fatima,2012) adalah pengukuran kedua variabel tersebut yaitu constant mounting length (CML) dan pool distance (PD) yang memiliki nilai Cpk rendah Hal ini dikarenakan untuk pengukuran kedua variabel memang tidak mudah karena berukuran mm. Ketika proses pembentukan lampu untuk mengukur lumen, wattage, dan life time. Berdasarkan pengukuran tersebut diperoleh bahwa diperoleh bahwa life time membe-rikan hasil yang bervariasi. Hal ini disebabkan oleh pengukuran variabel PD dan CML yang tidak tepat. Variabel PD dan CML yang diduga berpengaruh terhadap efisiensi yang dibentuk dari variabel lumen, wattage, dan life time. ISBN : 978-602-97491-7-5 A-3-1
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
Beberapa penelitian yang telah menerapkan dan mengembangkan metode Data Envelopment Analysis pada beberapa kasus untuk mendapatkan nilai efisiensi [Ali, Lerme, dan Seifore (1995), Angiz, Mustafa, dan Kamali (2013), Azadeh dkk (2011), Charnes, Cooper, dan Rhodes (1978), Ertay dan Ruan (2005), Hong dkk (1999), Mousavi-Avval (2011), Ramadany (2012), Ramanathan (2003), Rochma (2012)]. Beberapa penelitian yang juga menerapkan dan pengembangan metode response surface untuk mendapatkan nilai optimasi [Brandvik dan Daling (1998), Fatima (2012), Puri, Beg, dan Gupta (2001), Wong, Hobbs, dan Onof (2005), Ying (2012)]. Adanya suatu pengembangan DEA CCR dalam kesalahan pengaturan pembobotan yang berpengaruh terhadap nilai efisiensi. Oleh karena itu adanya suatu pengembangan metode baru yang menghendaki nilai efisiensi bisa bernilai lebih dari satu disebut dengan DEA dengan formulasi aggressive (Al-Refaie dan Li, 2008). Penelitian mengenai gabungan antara response surface dan data envelopment analysis CCR guna menjadikan multiple responses ke dalam bentuk single respon efisiensi untuk dioptimalkan (Tsai, Tong, dan Wang, 2010). Berdasarkan penelitian-penelitian tersebut, maka pada makalah ini adalah untuk penga-jian metode DEA aggressive. Pada makalah ini juga menerapkan kasus di PT. Phillips Indo-nesia, dimana mendapatkan nilai efisiensi hasil pengubahan multirespon menjadi single respon untuk setiap decision making units (DMUs) dan nilai optimasi dari hasil pengaruh coil terhadap efisiensi. Kemudian juga dilakukan perbandingan nilai lumen, wattage, dan life time dari hasil penelitian acuan (Fatima, 2012) dengan penelitian dengan metode yang disarankan. METODE Metode Data Envelopment Analysis CCR Data Envelopment Analysis (DEA) merupakan tek-nik pemrograman matematika nonparametrik yang di-gunakan untuk mengevaluasi efisiensi relatif dari kehomogenan Decision Making Units (DMUs) dengan multiple input dan multiple output (Al-Refaie dan Li, 2008). Model CCR menggunakan ‘virtual multiplier’ menggabungkan multiple inputs dan multiple outputs ke dalam indeks tunggal. ‘Virtual multiplier’ yang digunakan adalah membangkitkan dari jumlahan pembobot output yang dibagi dengan jumlahan pembobot input. Efisiensi relatif dari setiap DMUs ke-t misal dapat disimbolkan dengan t , dan untuk mendapatkan nilai efisiensi dapat menyelesaikan persamaan linier sebagai bentuk pengaturan input (Al-Refaie dan Li, 2008). n mak t v j y jt j 1
(1) p
Dengan fungsi kendala u k x kt 1
k 1 p n v j y ji uk xki 0 ; j 1 k 1
u1 , u 2 , , u p 0 v1, v2 , , vn 0
dimana i 1,2, , m
j 1,2, , n k 1,2, , p
Metode Data Envelopment Analysis Aggressive Konsep utama dari DEA aggressive yang merupakan penyempurnaan dari DEA CCR dimana akan mengoptimalkan model pada CCR yang meminimumkan cross-eficiencies DMUs lainnya. Secara umum, berdasarkan hasil analisis, nantinya didapatkan nilai DMUs ISBN : 978-602-97491-7-5 A-3-2
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
yang memiliki efisien. Sehingga untuk persoalan tersebut langkah yang perlu dilakukan adalah dengan “cross-evaluation” yang diaplikasikan ranking “best performer”. Pada penelitian ini, kriteria-kriteria tersebut digunakan untuk menentukan respon mana yang menjadi output dan input dengan ketentuan sebagai berikut (Al-Refaie dan Li, 2008). 1. Jika semua variabel respon memiliki kriteria STB, maka diantara respon tersebut dipilih salah satu untuk dijadikan sebagai input. 2. Jika semua variabel respon memiliki kriteria NTB, maka menghitung quality loss berdasarkan rumus berikut. s2 Li c i 2 dimana c merupakan nilai koefisien quality loss dan yi , si merupakan rata y i
rata dan standard deviasi dari perulangan respons setiap DMUs. Kemudian Li dijadikan sebagai input dan satu respons untuk output untuk semua DMUs. 3. Jika respons memiliki ketiga kriteria, maka respon yang memiliki kriteria STB dan Lj untuk respon kri-teria NTB dijadikan sebagai input. Respon dengan kriteria LTB dijadikan sebagai output. 4. Jika respons memiliki kriteria STB, dan LTB, untuk respon STB dijadikan sebagai input dan respon dengan kriteria LTB dijadikan sebagai output Berikut ini merupakan model dari konsep DEA aggressive pada DMUs ke-t (Al-Refaie dan Li, 2008).
n p min v jt y ji ukt xki k 1 t i j 1 t i
p
(2)
Dengan fungsi kendala ukt xki 1 k 1 t i p n v jt y ji u kt x ki ; j 0 j 1 k 1 p n v jt y jt Ett ukt xkt 0 j 1 k 1 v jt , ukt 0
Cross-efficiencies DMUs dapat dihitung dengan pembobot yang telah optimal. Perhitungan nilai cross-eficiencies disimbolkan dengan Eti yang dinyatakan sebagai berikut. n v jt y ji j 1 Eti p ukt xki k 1
(3)
Kemudian dilakukan estimasi terhadap nilai mean dari cross-eficiencies yang merupakan ratarata cross-eficiencies tiap DMUs untuk membandingkan performance dari banyaknya DMUs. E ti t ei i n 1
(4)
Metode Response Surface Response surface methodology (RSM) merupakan kumpulan teknik statistika yang digunakan untuk memodelkan dan menganalisis persoalan yang mana satu atau lebih dari satu ISBN : 978-602-97491-7-5 A-3-3
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
respon yang dipengaruhi oleh beberapa variabel dan secara objektif mendapatkan hubungan antara respon dengan beberapa variabel serta untuk mendapatkan nilai optimasi (Park, 1996). Model dengan persamaan polynomial dari derajat yang lebih tinggi atau disebut dengan model “orde kedua” yang meliputi kuadratik dan interaksi faktor. Berikut ini merupakan pendugaan model orde kedua. p p p p y 0 k xk kk xk 2 kk *xk xk * k 1 k 1 k 1 k* 1 k k*
(5)
Metode yang digunakan dalam penaksiran parameter metode response surface salah satunya adalah dengan MLE (Maximum Likelihood Estimation).
1 X TY
βˆ X T X
T Y Xβˆ Y Xβˆ 2 ˆ m
(6)
(7)
Berdasarkan model orde kedua dapat ditentukan titik stasioner. Titik stasioner digunakan untuk menentukan jenis dari karakteristik respon. 1 xs B 1βˆ 2
(8)
1 ys ˆ0 xs ' βˆ 2
(9)
Di samping titik stasioner, analisis yang dapat digunakan untuk menduga jenis karakteristik respon dengan anali-sis kanonik Pengujian residual perlu dilakukan, dimana residual memiliki beberapa asumsi yang harus dipenuhi yaitu identik, independen, dan berdistribusi normal. Metode yang digunakan untuk mendeteksi asumsi identik adalah uji glejser. Untuk melihat asumsi independen dapat dilihat melalui plot ACF (Autocorrelation Function) (Wei, 2006). Asumsi residual yang harus dipenuhi adalah berdistribusi normal dan pengujiannya adalah uji kolmogorov-smirnov HASIL DAN PEMBAHASAN Langkah awal yang dilakukan sebelum dilakukan analisis adalah menentukan variabel yang dijadikan sebagai variabel output dan input. Pada analisis ini, ditentukan variabel wattage dijadikan sebagai variabel input karena memiliki kriteria respon “smaller the better”, sedangkan lumen dan life time sebagai output dengan kriteria respon “larger the better”. Setelah ditentukan, kemudian dilakukan analisis DEA CCR dan DEA aggressive yang ditunjukkan pada Tabel 1.
ISBN : 978-602-97491-7-5 A-3-4
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013 Tabel 1 Nilai Effisiensi DEA
DMUs 1 2 3 4 5 6 7 8 9 10 11 12 13
Efisiensi DEA CCR 0,9563182 1 0,9268453 0,9028295 0,8814115 0,9419899 0,9399003 0,9386219 0,9358197 0,9365087 0,9562703 0,9677844 1
Efisiensi DEA Aggressive 0,702 1,218 0,853 0,841 1,062 0,843 0,847 0,858 0,612 0,696 0,702 0,545 1,225
Best Performer DEA Aggressive 5 12 9 6 11 7 8 10 2 3 4 1 13
Tabel 1 menginformasikan untuk DEA CCR, DMUs yang ke-2 dan ke-13 dapat dikatakan efisien karena memiliki nilai efisien sebesar 1. Sedangkan untuk nilai efisiensi terrendah yaitu DMUs ke-5 dengan nilai 0,8814. Secara umum nilai efisiensi DMUs tersebut bernilai 0,9. Sehingga tidak bisa dilihat DMUs mana yang lebih baik. Untuk mengatasi kondisi tersebut maka perlu dilakukan analisis DEA aggressive untuk melihat DMUs yang paling baik. Berdasarkan hasil analisis DEA formulasi aggressive diperoleh bahwa yang memiliki nilai efisiensi terbesar adalah DMUs ke-12 yaitu sebesar 1,225 sedangkan untuk nilai terkecil adalah DMUs ke-12 sebesar 0,545. Berdasarkan hasil analisis juga menginformasikan ranking efisiensi yang digunakan untuk melihat ‘best performer’. Pemeringkatan tersebut didapatkan dari nilai rata-rata yang terkecil sampai yang terbesar. Terdapat 3 DMUs yang memiliki nilai rata-rata efisiensi yang lebih dari 1 yaitu DMUs 2,5, dan 12. Nilai efisiensi DEA CCR dan DEA aggressive, kemudian dilakukan optimasi dengan menggunakan response surface. Setelah dilakukan analisis optimasi, didapatkan persamaan sebagai berikut. yˆ 0,958896 0,00123x1 0,0294 x2 0,0136 x12 0,009077 x22 0,01554 x1x2
(10)
yˆ 0,86254 0,03818x1 0,0226 x2 0,10544 x12 0,07919 x22 0,13552 x1x2
(11)
Persamaan 10 merupakan model efisiensi orde kedua untuk DEA CCR dengan nilai koefisien determinasi sebesar 72,46%. Sedangkan untuk persamaan 14 adalah model efisiensi DEA aggressive orde kedua. Nilai koefisien determinasi yang dimiliki sebesar 42,26%. Kedua model tersebut melibatkan interaksi antar faktor dan kuadratik faktor. Hasil optimasi efisiensi DEA CCR didapatkan bahwa optimasi efisiensi DEA CCR sebesar 1,0079 dengan nilai x1 sebesar 1,9032 dan nilai x2 sebesar -3,2503. Jika dikonversikan ke dalam nilai sebenarnya, maka didapatkan bahwa pool distance sebesar 26,9032 dan constant mounting length sebesar 24,49952. Sedangkan hasil analisis DEA aggressive didapatkan nilai bahwa untuk mendapatkan efisiensi secara optimal sebesar 1,113384 dengan nilai faktor pool distance sebesar -0,497573 dan nilai dari constant mounting length sebesar 1,323787. Pengkonversian pada nilai sebenarnya maka pool distance sebesar 24,502427 dan ISBN : 978-602-97491-7-5 A-3-5
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
nilai constant mounting length sebesar 27,5819408. Jika dibandingkan dengan hasil penelitian acuan maka dapat dilihat pada Tabel 2. Tabel 2 Hasil Perbandingan Analisis
Hasil Penelitian Fatima (2012) DEA CCR DEA Aggressive
Response Faktor Efisiensi Y1 Y2 Y3 X1 X2 1381,204 993,962 1428,18 24,1 28,5 1587,32 105,667 339,043 269,032 244,995 1,0079 1413,597 100,052 1755,99 245,024 275,819 1,1134
Tabel 2 menginformasikan bahwa untuk variabel respon lumen dengan menggunakan metode DEA CCR menghasilkan nilai yang lebih baik sesuai dengan kriteria respon LTB dimana meningkatkan nilai sebesar 14,923%. Kriteria respon LTB lainnya yaitu life time dengan menggunakan DEA aggressive menghasilkan nilai yang lebih baik dengan meningkatkan nilai sebesar 22,95%. Namun untuk nilai respon dengan kriteria STB menghasilkan nilai yang lebih baik dari hasil penelitian Fatima (2012). Secara keseluruhan, hasil analisis dengan menggunakan DEA aggressive menghasilkan nilai yang lebih baik karena terjadi peningkatan dibandingkan penelitian acuan untuk kriteria LTB. KESIMPULAN DAN SARAN Hasil analisa menunjukkan bahwa untuk hasil DEA CCR terdapat 2 DMUs yang efisien. Sedangkan untuk DEA aggressive terdapat 3 DMUs. ‘Best performer’ pada DEA aggressive adalah DMUs yang ke-13 dengan nilai efisiensi sebesar 1,225. Pada hasil optimasi dengan menggunakan DEA CCR efisiensi optimal pada 1,0079. Nilai tersebut dicapai ketika nilai pool distance sebesar 26,903 dan constant mounting length sebesar 24,499. Sedangkan DEA aggressive sebesar 1,1134 ketika nilai pool distance sebesar 24,502 dan nilai constant mount-ing length sebesar 27,582. Hasil perbadingan antara penelitian acuan, hasil analisis DEA CCR dan DEA aggressive menunjukkan bahwa secara umum aggressive memberikan hasil yang lebih baik untuk respon dengan kriteria ‘larger the better’. Penelitian acuan memberikan hasil yang lebih baik untuk respon dengan kriteria ‘smaller the better’. Untuk penelitian selanjutnya, sebaiknya perlu dimasukkan variabel lain yang diduga berpengaruh secara signifikan terhadap efisiensi baik secara parsial, interaksi dengan faktor lain, maupun kuadratik. DAFTAR PUSTAKA Al-Refaie, A., Li, M. H. (2008). Solving the Multiresponse Problem in Taguchi Method by Aggressive Formulation in DEA. London : Proceedings of the World Congress on Engineering Vol II Ali, A. A, Lerme, C. S., Seiford, L. M. (1995). Components of Efficiency Evaluation In Data Envelopment Analysis. European Journal of Operational Research, 80, 462-473. Angiz, M. Z., Mustafa, A., & Kamali, M. J. (2013). Cross-ranking of Decision Making Units in Data Envelopment Analysis. Applied Mathematical Modelling, 37, 398–405.
ISBN : 978-602-97491-7-5 A-3-6
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
Azadeh, A., Ghaderi, S. F., Mirjalili, M., & Moghaddam, M. (2011). Integration of Analytic Hierarchy Process and Data Envelopment Analysis for Assessment and Optimization of Personnel Productivity In A Large Industrial Bank. Expert Systems with Applications, 38, 5212–5225 Brandvik, P. J., Daling, P. S. (1998). Optimisation of Oil Spill Dispersant Composition By Mixture Design and Response Surface Methods. Chemometrics and Intelligent Laboratory Systems 42, 63–72. Charnes, A., Cooper, W. W., Rhodes, E. (1978). Measuring The Efficiency Of Decision Making Units. North-Holland : European Journal of Operational Research, 2, 429444 Daniel, W. (1989). Statistika Nonparametrik Terapan. Jakarta: PT. Gramedia. Draper, N., R., Smith, H. (1998). Applied Regression Analysis : Third Edition. Canada : John Wiley & Sons, Inc. Ertay, T., Ruan, D. (2005). Data Envelopment Analysis Based Decision Model For Optimal Operator Allocation In CMS. European Journal of Operational Research, 164, 800– 810. Fatima, A. (2012). Optimasi Lumen, Wattage Dan Life Time Pada Lampu GLS Dengan Menggunakan Pendekatan Fungsi Desirability pada Metode Multiresponse Surface (Studi Kasus : PT. Philips Indonesia). Surabaya: Program Sarjana, Institut Teknologi Sepuluh Nopember. Hong, H. K., H, S. H., Shin, C. K., Park, S. C., & Kim, S. H. (1999). “Evaluating The Efficiency of System Integration Projects Using Data Envelopment Analysis (DEA) and Machine Learning”. Expert Systems with Applications, 16, 283–296. Hutasoit, G. (2012). Lampu Pijar. http://gressellahutasoit.blogspot.-com/2012/12/ lampu-pijar. html. [diakses tanggal, 20 Februari 2013, pukul 13:25] Khuri, A. I., & Cornell, J. A. (1996). Response Surfaces Designs and Analysis. New York : Marcel Dekker. Montgomery, D. C. (1997). Design and Analysis of Expeiments 5th Edition. New York: John Wiley, Sons. Mousavi-Avval, S. H., Rafiee, S., Jafari, A., & Mohammadi, A. (2011). Improving Energy Use Efficiency of Canola Production Using Data Envelopment Analysis (DEA) Approach. Karaj : Department of Agricultural Machinery Engineering, Faculty of Agricultural Engineering and Technology, University of Tehran Park, S. H. (1996). Robust Design and Analysis For Quality Engineering. London : Chapman & Hall, 2-6 Boundary Row, London SE1 8HN, UK. Puri, S, Beg, Q. K., & Gupta, R. (2001). Optimization of Alkaline Protease Production from Bacillus sp. by Response Surface Methodology. New Delhi : Department of Microbiology, University of Delhi South Campus.
ISBN : 978-602-97491-7-5 A-3-7
Prosiding Seminar Nasional Manajemen Teknologi XVIII Program Studi MMT-ITS, Surabaya 27 Juli 2013
Ramadany, R. (2012). Analisis Tingkat Efisiensi Pelayanan Kesehatan di Tiap Kabupaten/ Kota Se-Jawa Timur Dengan Metode Data Envelopment Analysis (DEA). Surabaya: Program Sarjana, Institut Teknologi Sepuluh Nopember. Ramanathan, R. (2003). An Introduction to Data Envelopment Analy-sis : A Tool for Performance Measurement. New Delhi : Sage Publi-cations India Pvt Ltd. Ramitha, V. (2008). Pelopor Lampu Hemat Energi. http://m.inilah.com/read/detail/23797/pelopor-lampu-hemat-energi.[diakses tanggal, 20 Februari 2013, pukul 13:20] Rochma, A. (2012). Optimasi Multirespon Pada Metode Taguchi Dengan Pendekatan Data Envelopment Analysis (Dea) Formula Agre-sif (Studi Kasus Mesin Wire-EDM). Surabaya: Program Sarjana, Institut Teknologi Sepuluh Nopember. Tsai, C., Tong, L, Wang, C. (2010). Optimization od Multiple Responses Using Data Envelopment Analysis and Response Surface Methodology. Tamkang Journal of Science and Engineering, Vol. 12, No. 2, pp. 197-203. Wei,W.W.(2006).Time Series Analysis.USA: PearsonEducation Wikipedia (2013). Lampu Pijar. http://id.wikipedia.org/wiki/Lampu_pijar, [diakses tanggal, 20 Februari 2013, pukul 13:30] Wong, S. M., Hobbs, R. E., Onof, C. (2005). An Adaptive Response Surface Method For Reliability Analysis of Structures With Multiple Loading Sequences. Structural Safety, 27, 287–308. Ying, L. (2012). Application of Stochastic Response Surface Method in the Structural Reliability. Procedia Engineering, 28, 661 – 664.
ISBN : 978-602-97491-7-5 A-3-8