PENGENDALIAN KUALITAS “X” MENGGUNAKAN GRAFIK PENGENDALI HOTELLING T2 UNIVARIAT DAN MULTIVARIAT Lellie Sulistyawati Darmawan, Adi Setiawan, Lilik Linawati Program Studi Matematika, Fakultas Sains dan Matematika Universitas Kristen Satya Wacana Salatiga Jl. Diponegoro 52-60 Salatiga 50711
Abstrak Pengendalian kualitas memiliki peranan yang penting khususnya dalam meningkatkan volume penjualan produk. Salah satu metode statistik dalam mengendalikan produk adalah dengan grafik pengendali. Kualitas suatu produk biasanya ditentukan oleh lebih dari satu karakteristik. Grafik pengendali untuk beberapa karakteristik kualitas disebut grafik pengendali multivariat. Salah satu grafik pengendali multivariat yaitu grafik pengendali Hotelling T2. Pada penelitian ini akan diaplikasikan grafik pengendali Hotelling T2 univariat dan multivariat untuk ketiga karakteristik data proses produksi ”X” selama bulan Januari-Agustus 2010. Pada grafik Hotelling T2 multivariat dihitung batas pengendali atas sebesar 14,6064. Dari batas tersebut, terdeteksi dua sampel yang berada di luar control yaitu sampel 50 dan 101. Selanjutnya, dari ketiga karakteristik kualitas tersebut akan dianalisis karakteristik yang paling dominan mempengaruhi sampel out of control dengan menggunakan boxplot untuk menganalisis sampel outliers dan diperoleh karakteristik yang paling dominan mempengaruhi adalah brix (kadar gula). Dari studi simulasi dapat dibangkitkan sampel dengan ukuran berbeda- beda dan diperoleh hasil proporsi sampel out-of-control mendekati nilai batas kesalahan (level of significance) α = 0,0027 . Setelah sampel out-of-control terdeteksi, diasumsikan telah dilakukan penelusuran penyebab proses tidak terkontrol untuk selanjutnya dihitung nilai indeks kemampuan proses baik secara univariat maupun multivariat. Nilai indeks kemampuan proses menunjukkan nilai lebih besar dari 1 sehingga dapat disimpulkan bahwa proses produksi telah berjalan dengan baik. Kata kunci : grafik pengendali, Hotelling T2, indeks kemampuan proses
1. PENDAHULUAN Latar Belakang Di era globalisasi yang semakin kompetitif ini, para pelaku bisnis tentu menginginkan agar produknya diterima oleh konsumen dan mampu bersaing di pasaran. Salah satu faktor yang mempengaruhi keputusan konsumen dalam memilih suatu produk adalah kualitas produk tersebut. Konsumen akan merasa puas apabila kualitas produk yang mereka pilih sesuai dengan harapan mereka. Tingkat kepuasan konsumen dapat tercermin pada keputusan untuk membeli produk dan melakukan pembelian ulang terhadap produk tersebut. Selanjutnya, kepuasan konsumen terhadap produk diharapkan dapat meningkatkan volume penjualan produk dan akhirnya berpengaruh pada keuntungan perusahaan. Oleh sebab itu, masalah kualitas menjadi hal yang penting dan perlu mendapat perhatian perusahaan. Pengendalian kualitas produk dalam proses produksi perlu dilakukan secara terus menerus untuk dapat mendeteksi ketidaknormalan secara cepat sehingga dapat dilakukan tindakan perbaikan yang diperlukan sebelum timbul banyak produk yang tidak sesuai dengan spesifikasi perusahaan.
PT. Y merupakan salah satu perusahaan yang memproduksi minuman dalam kemasan dengan orientasi kualitas ekspor. Oleh sebab itu, kualitas produk menjadi perhatian utama perusahaan untuk menjaga loyalitas konsumen terhadap perusahaan dan dengan demikian meningkatkan dominasi pasar. Untuk menjaga kualitas produk, perusahaan secara terus menerus melakukan pengendalian kualitas terhadap produk yang dihasilkan. Sehubungan dengan efisiensi, maka pengamatan kualitas produk dilakukan dengan ukuran sampel satu. Pengamatan seperti ini disebut pengamatan individual. Salah satu produk yang dihasilkan oleh PT. Y adalah minuman teh dalam kemasan dengan merk ”X”. Kualitas produk ”X” ditentukan oleh pengukuran empat karakteristik antara lain temperatur (suhu), brix (kadar gula), pH (kadar keasaman), dan vacuum. Perusahaan telah menentukan batas spesifikasi untuk masing masing karakteristik kualitas tersebut. Produk dianggap cacat atau tidak memenuhi syarat jika terdapat setidaknya satu pengukuran karakteristik yang tidak berada dalam interval spesifikasi yang telah ditentukan perusahaan. Pada tugas akhir ini, akan diaplikasikan grafik pengendali Hotteling T2 untuk memonitor pergeseran mean proses. Selanjutnya akan dihitung indeks kemampuan proses (IKP) untuk mengetahui apakah proses produksi yang berlangsung telah memenuhi spesifikasi yang ditentukan perusahaan. Perumusan Masalah Berdasarkan latar belakang tersebut, maka dapat dirumuskan permasalahan : 1. Bagaimana mengaplikasikan grafik pengendali Hotelling T2 (Hotelling T2 Control Chart) untuk mengendalikan kualitas produk ”X” yang terdiri dari beberapa karakteristik kualitas? 2. Berapakah nilai indeks kemampuan proses (IKP) pada proses produksi ”X” ? Tujuan Penelitian Tujuan dari penelitian ini adalah : 1. Mengaplikasikan metode Hotelling T2 pada produk ”X” yang terdiri dari beberapa karakteristik. 2. Memperoleh nilai indeks kemampuan proses pada produksi “X”. Batasan Masalah Data yang digunakan dalam penelitian ini adalah data yang diperoleh dari PT. Y. Jenis produk yang akan dianalisis adalah produk ”X” selama kurun waktu produksi bulan Januari 2010 sampai Agustus 2010. Uji distribusi normal untuk keempat karakteristik produk menunjukkan data tidak berdistribusi normal untuk karakteristik pertama (suhu) sehingga pada penelitian ini akan dianalisis tiga karakteristik antara lain : brix (kadar gula), kadar keasaman (pH), dan vacuum. Komputasi dilakukan dengan menggunakan bantuan software Matlab 6.5 dan Maple 9.5. 2. DASAR TEORI Pengendalian Kualitas Statistik (Statistical Quality Control) Tiap produk mempunyai sejumlah unsur yang menggambarkan kecocokan penggunaannya. Unsur – unsur ini biasanya dinamakan ciri kualitas. Terdapat beberapa jenis ciri kualitas antara lain : fisik (panjang, berat, kekentalan), indera (rasa, warna), dan orientasi waktu
(ketahanan, dapatnya dirawat). Kualitas adalah keadaan fisik, fungsi dan sifat suatu produk yang dapat memenuhi selera dan kebutuhan konsumen dengan memuaskan sesuai nilai uang yang dikeluarkan. Pengendalian kualitas adalah aktivitas keteknikan dan manajemen yang dengan aktivitas tersebut diukur ciri – ciri kualitas produk, membandingkannya dengan spesifikasi dan mengambil tindakan penyehatan yang sesuai apabila ada perbedaan antara penampilan yang sebenarnya dengan yang standar (Montgomery, 1990:3). Oleh karena itu, dalam pengendalian kualitas diupayakan sedini mungkin dapat mendeteksi adanya ketidakstabilan sehingga segera dapat dilakukan tindakan pembetulan sebelum timbul banyak unit yang tidak sesuai dengan spesifikasi yang perusahaan. Statistik adalah pengambilan keputusan tentang suatu proses atau populasi berdasarkan pada suatu analisis informasi yang terkandung dalam suatu sampel dari populasi tersebut. Metode statistik juga memainkan peranan penting dalam jaminan kualitas. Metode statistik memberikan cara – cara pokok dalam pengambilan sampel produk, pengujian serta evaluasinya, dan informasi di dalam data itu digunakan untuk mengendalikan dan meningkatkan proses pembuatan (Montgomery, 1990 : 27). Salah satu alat dalam pengendalian kualitas statistik yang digunakan untuk memonitor stabilitas proses adalah grafik pengendali (control chart). Grafik Pengendali Grafik pengendali adalah teknik pengendali proses pada jalur yang digunakan secara luas yang biasanya digunakan untuk menaksir parameter suatu proses produksi menentukan kemampuan dan memberikan informasi yang berguna dalam meningkatkan proses itu (Montgomery, 1990 : 120). Grafik pengendali merupakan suatu grafik karakteristik kualitas yang dibentuk oleh data karakteristik kualitas sampel dan waktu sampel. Grafik pengendali terdiri dari nilai karakteristik kualitas sampel, garis batas pengendali atas (BPA) dan garis batas pengendali bawah (BPB). Berdasarkan banyaknya karakteristik kualitas yang diukur, grafik pengendali dibedakan menjadi 2 jenis yaitu grafik pengendali univariat dan grafik pengendali multivariat. Grafik pengendali univariat digunakan jika hanya ada satu karakteristik kualitas yang diukur, sedangkan grafik pengendali multivariat digunakan jika diperlukan pengendalian dua atau lebih karakteristik kualitas yang berhubungan secara bersama sama. Secara umum, grafik pengendali dibedakan atas dua macam, yaitu grafik pengendali untuk variabel dan grafik pengendali untuk atribut. a. Grafik Pengendali Variabel Grafik pengendali dengan pemeriksaan secara variabel digunakan sebagai alat untuk mengontrol proses dengan mengumpulkan dan menganalisis data dari pemeriksaan kualitas suatu produk yang ditentukan satu atau lebih parameter karakteristik kualitas yang mempunyai harga ukur tertentu. Pengendalian rata – rata proses biasanya dilakukan dengan menggunakan grafik pengendali x sedangkan pengendalian variabilitas (pemencaran) dilakukan dengan menggunakan grafik pengendali untuk deviasi standar (grafik S) atau grafik pengendali untuk rentang yaitu grafik R (Montgomery 1990 : 204). b. Grafik Pengendali Atribut Seringkali karakteristik kualitas suatu produk tidak dapat dengan mudah dinyatakan secara numerik. Dalam kasus seperti itu, biasanya produk yang diperiksa diklasifikasikan sesuai dengan spesifikasi pada karakteristik kualitas tersebut. Istilah ”cacat” dan ”tidak cacat” sering digunakan untuk mengidentifikasi kedua klasifikasi produk tersebut. Karakteristik kualitas
seperti ini dinamakan sifat (atribut). Terdapat beberapa grafik pengendali sifat yang sering digunakan yaitu yang berhubungan dengan bagian yang tidak sesuai dan grafik pengendali untuk ketidaksesuaian. Grafik untuk bagian yang tidak sesuai dinamakan grafik p sedangkan grafik untuk ketidaksesuaian dinamakan grafik c atau u. Grafik Pengendali Hotelling T2 Univariat Salah satu pendekatan yang digunakan dalam memantau kualitas produk pada kasus multivariat adalah dengan menggunakan metode grafik pengendali Hotelling T2. Grafik pengendali Hotelling T2 diperkenalkan oleh Harold Hotelling pada tahun 1947 dengan menggunakan data pembidik bom selama Perang Dunia II. Hotelling T2 merupakan generalisasi dari distribusi-t. Jika x adalah mean dari sampel acak berukuran n dari suatu distribusi normal dengan mean μ dan variansi s maka
t=
x−μ
s/ n akan mempunyai distribusi-t dengan derajat bebas n-1. Selanjutnya nilai
(x − μ ) =
(
2
t
2
)( ) (x − μ )
= n x − μ s2
−1
'
(2.1) s2 / n akan mengikuti distribusi χ 2 dengan derajat bebas 1 (Young, 1999). Sejalan dengan persamaan (2.1) nilai T 2 pada sampel ke-j ( T j2 ) dapat dinyatakan sebagai berikut :
Tj
2
(x =
j
−x
)
2
s2 / n
dengan batas pengendali atas
(
)( ) (x
= n x j − x s2
−1
j
−x
)
'
(2.2)
n −1 2 χn . n
Grafik Pengendali Hotelling T2 untuk Pengamatan Subgrup Diketahui sampel bedistribusi normal yang terdiri dari p-karakteristik kualitas, dengan m menggambarkan banyaknya sampel, dan masing – masing sampel berukuran n. Didefinisikan vektor rataan sampel X j = X 1 j X 2 j K X pj j = 1,2,..., m
[
]
dengan X ij mewakili rataan sampel dari karakteristik mutu ke-i untuk sampel ke-j dan n
X ij =
∑X k =1
ijk
, i = 1, 2, . . ., p dan j = 1, 2, . . ., m. n Selanjutnya, vektor rataan nominal untuk tiap – tiap karakteristik mutu dinotasikan dengan X = [ X 1 X 2 K X p ] dengan X i merupakan vektor rataan nominal untuk karakteristik mutu ke i dan 1 m X i = ∑ X ik i = 1, 2, . . ., p. m k =1 Variansi sampel untuk karakteristik mutu ke-i dalam sampel ke-j dihitung dari
n 1 ( X ijk − X ij ) 2 i = 1, 2, . . ., p dan j = 1, 2, . . ., m ∑ (n − 1) k =1 sedangkan kovariansi antara karakteristik mutu i dan h dalam sampel ke-j adalah n 1 S ihj = ∑ ( X ijk − X ij )( X hjk − X hj ) j = 1, 2, . . ., m dan i ≠ h. (n − 1) k =1
S ij2 =
Selanjutnya dihitung rata – rata S ij2 dan S ihj meliputi seluruh m sampel untuk memperoleh m
S i2 =
1 m
∑ S ij2
S ih =
1 m
∑ S ihj
i = 1, 2, . . ., p
j =1
dan m
i ≠ h.
j =1
Matriks kovarian S dapat disusun sebagai berikut ⎡ S12 S12 L S1 p ⎤ ⎢ ⎥ S 22 L S 2 p ⎥ ⎢ S= . ⎢ M M ⎥ ⎢ ⎥ S p2 ⎦⎥ ⎣⎢ Nilai T2 untuk masing – masing sampel adalah T j = n( X j − X ) S −1 ( X j − X ) ' dengan j = 1, 2, . . ., m. (2.3) -1 dengan n adalah ukuran masing masing sampel dan S merupakan invers dari matriks kovariansi S sedangkan batas grafik pengendali dapat ditentukan dari persamaan berikut : ⎛ mnp − mp − np + p ⎞ ⎟⎟ Fα , p ,( mn − m − p +1) BPA = ⎜⎜ (2.4) ⎝ mn − m − p + 1 ⎠ dengan m menggambarkan banyak sampel, masing – masing berukuran n. Jika nilai T2 untuk sampel ke-j, yaitu T j2 > BPA, hal ini menunjukkan sampel ke -j di luar kendali (Young, 1999). 2
Grafik Pengendali Hotelling T2 Untuk Pengamatan Individual Diketahui sampel berukuran m yang terdiri dari p-karakteristik dengan matriks data digambarkan sebagai berikut : ⎡ x11 x12 L x1 p ⎤ ⎢x x 22 L x 2 p ⎥⎥ 21 ⎢ X = . ⎢ M M O M ⎥ ⎢ ⎥ ⎣⎢ x m1 x m 2 L x mp ⎦⎥ Vektor rataan sampel didefinisikan sebagai X = X1 X 2 K X p
[
dengan
]
m
∑ X ji Xi =
j =1
i = 1, 2, . . ., p. m Kovariansi dihitung dengan menggunakan persamaan
(
)(
1 m S ik = ∑ xij − xi xkj − xk m − 1 j =1
)
(2.5)
dengan i = 1,2,..p dan k = 1,2,..p. Statistik Hotelling T2 dihitung dengan menggunakan persamaan T j2 = ( X j − X ) S −1 ( X j − X ) ' . (2.6) Batas grafik pengendali dapat ditentukan dari persamaan berikut ⎛ p(m + 1)(m − 1) ⎞ ⎟⎟ Fα , p,( m− p ) BPA = ⎜⎜ (2.7) ⎝ m(m − p ) ⎠ dengan m menggambarkan banyaknya sampel. Jika nilai T2 untuk sampel ke-j, yaitu T j2 > BPA, hal ini menunjukkan sampel ke-j di luar kendali (Young, 1999). Indeks Kemampuan Proses Univariat Kemampuan / kapabilitas proses adalah alat ukur suatu proses, apakah proses produksi tersebut sudah memenuhi spesifikasi yang telah ditentukan oleh perusahaan. Analisis kemampuan proses dilakukan berdasarkan asumsi bahwa data karakteristik produk berdistribusi normal dan proses berada di bawah kendali. Analisis kemampuan proses dilakukan dengan membandingkan data proses yang berada di bawah kendali dengan batas spesifikasi yang telah ditentukan oleh perusahaan. Proses dikatakan capable apabila penyebaran proses aktual lebih kecil daripada sebaran yang diijinkan perusahaan. Perbandingan antara penyebaran spesifikasi produk dengan penyebaran proses sebenarnya disebut indeks kemampuan proses. Metode perhitungan indeks kemampuan proses yang banyak digunakan dalam dunia industri yaitu Cp, Cpk dan Cpm. Kane (1986) mengembangkan Indeks Kemampuan Proses Cp dan Cpk yang umum digunakan dalam dunia industri untuk mengevaluasi proses yang melibatkan satu karakteristik kualitas. Indeks kemampuan proses ini mengasumsikan data berdistribusi normal dan dapat ditulis dalam bentuk USL − LSL Cp = 6σ ⎛ USL − μ μ − LSL ⎞ C pk = min⎜ , ⎟ 3σ ⎠ ⎝ 3σ dengan USL merupakan batas spesifikasi atas dan LSL batas spesifikasi bawah. Mean dan standar deviasi proses berturut – turut ditunjukkan oleh μ dan σ . Selanjutnya, nilai estimasi sampel untuk masing masing indeks kemampuan proses tersebut adalah USL − LSL , (2.8) Cˆ p = 6s ⎛ USL − x x − LSL ⎞ (2.9) , Cˆ pk = min⎜ ⎟. 3s ⎠ ⎝ 3s
Nilai indeks kemampuan proses Cp dan Cpk tidak memperhitungkan selisih antara mean proses dengan nilai mean target. Chan (1988) mengembangkan nilai indeks kemampuan proses Cpm yang dapat ditentukan dengan persamaan USL − LSL C pm = 6 σ 2 + (μ − T ) 2 dengan T menunjukkan nilai mean target yang diharapkan perusahaan. Selanjutnya nilai estimasi sampel ditentukan dengan persamaan USL − LSL . (2.10) Cˆ pm = 2 2 6 s + (x − T ) Indeks Kemampuan Proses Multivariat ( Multivariate Capability Process ) Indeks Kemampuan Proses Multivariat (Multivariate Capability Process) adalah suatu indeks proses yang menunjukkan nilai rasio antara penyebaran (variabilitas) spesifikasi produk yang diijinkan dan penyebaran proses aktual yang melibatkan lebih dari satu variabel. Ada beberapa macam metode perhitungan indeks kemampuan proses, salah satunya adalah metode indeks kemampuan proses MCpm (Zahid, 2008). Perhitungan nilai indeks kemampuan proses MCpm ini didefinisikan sebagai rasio dari dua volume yaitu vol ( R1 ) MC pm = vol ( R2 ) dengan R1 merupakan daerah toleransi modifikasi, sedangkan R2 merupakan daerah proses 99.73 %. Jika data berdistribusi normal multivariat maka R2 berbentuk ellipsoid sedangkan R1 merupakan ellipsoid terbesar yang berada dalam daerah spesifikasi dan berpusat pada target dengan volume R1 adalah p
vol ( R1 ) =
2∏ Ti i =1
π p/2
(2.11) p Γ( p / 2 ) dengan Ti merupakan nilai tengah spesifikasi ke-i (i=1,2,...p). Selanjutnya, indeks kemampuan proses multivariat dapat dituliskan sebagai Vol ( R1 ) MC pm = Vol ( x − μ ) ∑ T−1 ( x − μ ) T ≤ K ( p)
(
)
dengan K(p) merupakan kuantil 99.73 % dari distribusi χ 2 dengan derajat bebas p dan
∑T−1 = E[( x − T )( x − T ) ' ] . Bagian penyebut dapat dituliskan dalam bentuk vol ( R2 ) = ∑ MCˆ pm
1/ 2
[
(πK ( p )) p / 2 [Γ( p / 2 + 1)] × 1 + ( x − μ ) ∑ −1 ( x − μ ) ' −1
]
1/ 2
.
Nilai estimasi indeks MCpm ditentukan dengan rumus vol ( R1 ) 1 = 1/ 2 × 1/ 2 −1 p/2 T S (πK ( p )) [Γ( p / 2 + 1)] n ⎤ ⎡ −1 1 x T S x T + − − ⎥ ⎢ n −1 ⎦ ⎣
(
)
(
)
atau MCˆ pm =
Cˆ p Dˆ
(2.12)
dengan Cˆ p =
vol ( R1 )
S
1/ 2
(πK ( p)) p / 2 [Γ( p / 2 + 1)]
−1
(2.13)
dan
(
) (
)
1/ 2
n ⎡ ⎤ (2.14) x − T ' S −1 x − T ⎥ . Dˆ = ⎢1 + ⎣ n −1 ⎦ Jika nilai indeks lebih dari 1 maka proses mempunyai variasi lebih kecil dibandingkan dengan batas spesifikasi sehingga dapat dikatakan proses produksi telah berjalan dengan baik. Sebaliknya, jika indeks bernilai kurang dari 1 hal tersebut menunjukkan variasi proses lebih besar daripada batas spesifikasi perusahaan. Artinya proses tersebut banyak menghasilkan produk yang tidak sesuai dengan spesifikasi. Untuk memperjelas konsep IKP multivariat, akan dihitung IKP untuk data produk keramik Can-cup S-9 (Zahid, 2008) pada contoh berikut. Diambil 500 data produk keramik yang terdiri dari dua karakteristik yaitu panjang dan berat. Spesifikasi untuk masing – masing karakteristik didefinisikan sebagai interval panjang 4 – 6 cm dan berat pada kisaran 190 – 242 gr. Batas – batas spesifikasi tersebut membentuk daerah toleransi (persegi panjang). Pusat spesifikasi tersebut, yaitu 5 dan 216, diasumsikan sebagai target spesifikasi. Nilai mean sampel untuk tiap karakteristik berturut turut yaitu 5,55 dan 219,73. Matriks kovariansi diperoleh ⎡0,45 0,46 ⎤ S=⎢ ⎥. ⎣0,46 62,70⎦ Selanjutnya, akan ditunjukkan daerah toleransi spesifikasi, ellips terbesar yang memenuhi daerah tersebut dan daerah proses dalam Gambar 2.1.
Gambar 2.1. Gambaran geometris untuk perbandingan antara daerah toleransi spesifikasi dan daerah proses yang sebenarnya
Dari Gambar 2.1. dapat dilihat bahwa sebagian daerah ellips proses berada di luar spesifikasi yang telah ditentukan. Selain itu, ellips proses yang sebenarnya lebih besar dari ellips toleransi spesifikasi sehingga volume R2 lebih besar daripada volume R1. Akibatnya nilai IKP akan lebih kecil dari 1 sehingga kita dapat memprediksi bahwa nilai IKP data produksi keramik tersebut lebih kecil dari 1. Secara matematis, nilai Indeks Kemampuan Proses diperoleh sebagai berikut Cˆ p 0,72 MCˆ pm = = = 0,39 . 1,827 Dˆ Nilai IKP menunjukkan nilai kurang dari 1, maka dapat dikatakan proses produksi banyak menghasilkan produk yang berada di luar spesifikasi. Hasil tersebut analog dengan gambaran yang diperoleh dari Gambar 2.1 yang menunjukkan banyak titik sampel yang diproduksi berada di luar persegi panjang merah sebagai area spesifikasi perusahaan. Dari nilai tersebut, dapat disimpulkan bahwa proses produksi belum berjalan dengan baik. 3. METODE PENELITIAN Sumber Data Penelitian Data yang digunakan dalam penelitian ini adalah data sekunder pada proses produksi ”X” kemasan botol. Produk ”X” kemasan botol merupakan salah satu produk minuman yang diproduksi oleh PT. Y. Data yang akan diteliti merupakan data produksi minuman ”X” kemasan botol selama periode bulan Januari 2010 sampai Agustus 2010. Namun, proses produksi tidak dilakukan setiap hari sepanjang bulan tersebut. Proses produksi dilakukan sesuai jadwal dan persediaan produk. Adapun karakteristik kualitas produk ”X” kemasan botol yang digunakan dalam penelitian ini antara lain : a. Brix atau kadar gula dalam larutan minuman, dengan batas spesifikasi yang ditetapkan oleh perusahaan 8,25 – 8,55. b. Ph atau kadar keasaman larutan minuman ”X”, dengan batas spesifikasi 6,1 – 7,1. c. Vacuum, atau tekanan pada saat penyegelan kemasan, dengan batas spesifikasi 25 – 45 (batas yang dianjurkan). Metode Analisis Data Langkah langkah dalam analisis data dijabarkan sebagai berikut : 1. Pengujian data berdistribusi normal multivariat. Pengujian distribusi normal multivariat dilakukan dengan menggunakan metode Chi-Square plot. 2. Membuat program grafik pengendali Hotelling T2 pengamatan individual dengan menggunakan software Matlab 6.5 untuk kasus univariat dan multivariat. 3. Penerapan grafik pengendali Hotelling T2 individual untuk ketiga karakteristik produk ”X”. 4. Jika grafik pengendali sudah menunjukkan bahwa semua sampel berada di bawah kendali maka dapat dilanjutkan dengan menghitung indeks kemampuan proses. Jika tidak, maka harus dilakukan penelusuran faktor penyebab ketidakstabilan proses. Setelah faktor penyebab diketahui, dilanjutkan dengan melakukan pengontrolan sampai semua data berada di bawah kendali. Pada penelitian ini, untuk tiap sampel yang terdeteksi out of control diasumsikan telah dilakukan penelusuran penyebab
ketidakstabilan proses sehingga dapat dihilangkan untuk menghitung indeks kemampuan proses. 5. Menghitung indeks kemampuan proses dari data yang sudah terkontrol. Mulai Pengumpulan data karakteristik produk ”X” di PT. Y
Uji distribusi normal data karakteristik produk ”X”
Penerapan Grafik Pengendali Hotelling T2 pada data karakteristik kualitas produk ”X”
Out of control?
ya
Penelusuran penyebab ketidakstabilan
tidak Hitung IKP
Proses telah berjalan dengan baik
tidak
IKP < 1 ?
ya
Proses belum berjalan dengan baik
Selesai
Gambar 3.1. Diagram alir penelitian
4. ANALISIS DATA DAN PEMBAHASAN
Pada bab ini akan dianalisis data karakteristik kualitas produk ”X” kemasan botol selama proses produksi bulan Januari sampai Agustus 2010 (Lampiran 1). Sebelum data tersebut dianalisis dengan menggunakan Grafik Pengendali Hotelling T2 individual, terlebih dahulu dilakukan pengujian distribusi normal multivariat. Pengujian distribusi normal multivariat ini dengan menggunakan metode chi-square dan QQ-plot. Hasil pengujian distribusi normal multivariat terhadap ketiga karakteristik produk ”X” menunjukkan bahwa data tersebut mengikuti distribusi normal multivariat. Setelah asumsi distribusi normal multivariat dipenuhi maka analisis data dengan menggunakan grafik pengendali Hotelling T2 dapat dilakukan. Spesifikasi Perusahaan Perusahaan telah menetapkan standar spesifikasi atau batasan nilai untuk masing masing karakteristik kualitas produk. Produk dianggap ”cacat” jika tidak memenuhi batas spesifikasi yang telah ditentukan oleh perusahaan sehingga dapat dilihat apakah ada data yang tidak memenuhi batas spesifikasi. Dari data diperoleh satu sampel yang tidak memenuhi spesifikasi yang telah ditentukan perusahaan khususnya pada pengukuran variabel pertama yaitu kadar gula (Brix). Sampel tersebut yaitu sampel ke 101 yang menunjukkan nilai 8,21. Nilai tersebut tidak memenuhi batas bawah variabel pertama yaitu 8,25. Selanjutnya, data produksi “X” disajikan pada Gambar 4.1 dengan batas spesifikasi untuk masing masing karakteristik yang telah ditentukan perusahaan. Penerapan Grafik Pengendali Hotelling T2 Univariat Data produksi ”X” akan diaplikasikan pada grafik pengendali Hotteling univariat untuk masing masing karakteristik kualitas produk. Brix (Kadar Gula) Rataan untuk data brix (kadar gula) adalah 8,3767 dengan variansi s 2 = 0,0013 . Untuk masing masing sampel dihitung nilai T2 dengan persamaan (2.2) sedangkan nilai BPA diperoleh 9,0290. Selanjutnya, nilai T2 untuk tiap sampel disajikan dalam grafik pengendali Hotelling univariat pada Gambar 4.2. Dari Gambar 4.2. diperoleh dua sampel yang berada di luar batas pengendali yaitu sampel ke 50 dan 101. Nilai T2 untuk masing masing sampel tersebut adalah 12,43 dan 21,51. Jika dibandingkan dengan batas spesifikasi perusahaan, sampel ke 101 juga berada di luar batasan yang telah ditetapkan perusahaan. Pada grafik spesifikasi perusahaan, sampel ke 50 tidak teridentifikasi sebagai sampel out of control, melainkan berada tepat pada batas bawah spesifikasi perusahaan sebesar 8,25. pH (Kadar Keasaman) Data kadar keasaman mempunyai rataan x = 6,5632 dengan variansi s 2 = 0,0223 . Nilai T2 untuk masing masing sampel dihitung dengan persamaan (2.2) lalu dibandingkan dengan batas pengendali atas. Dari ke 309 sampel diperoleh tidak ada satupun sampel yang berada di luar BPA. Gambar 4.3. menunjukkan nilai T2 untuk tiap sampel dalam grafik pengendali Hotelling univariat.
Gambar 4.1. Grafik tiap karakteristik kualitas produk ”X” dibandingkan dengan spesifikasi perusahaan.
Gambar 4.2. Grafik pengendali Hotelling T2 univariat untuk karakteristik brix (kadar gula).
Gambar 4.3. Grafik pengendali Hotelling T2 univariat untuk karakteristik pH (kadar keasaman).
Vacuum Nilai rataan data vacuum yaitu x = 34,2654 dengan variansi s 2 = 4,3449 . Untuk masing masing sampel dihitung nilai T2 dengan persamaan 2.2. Dari 309 sampel diperoleh 1 sampel yang berada di luar BPA yaitu sampel ke 57. Nilai T2 untuk sampel tersebut menunjukkan nilai 9,0346 yang nilainya lebih besar dari BPA = 9,0290. Jika dibandingkan dengan batas spesifikasi perusahaan sampel ke 57 berada di dalam batas spesifikasi. Namun, jika dibandingkan dengan sampel yang lain, sampel ke-57 menunjukkan nilai paling minimum diantara semua sampel lainnya. Untuk lebih jelasnya, ditunjukkan grafik pengendali Hotelling univariat pada Gambar 4.4.
Gambar 4.4. Grafik pengendali Hotelling T2 univariat untuk karakteristik vacuum.
Penerapan Grafik Pengendali Hotelling T2 Multivariat Setelah syarat data berdistribusi normal multivariat dipenuhi, data dianalisis dengan grafik pengendali Hotelling T2 berdasarkan teori yang telah diuraikan pada bab II. Rataan untuk masing masing karakteristik kualitas produk adalah X = [8,3767 6,5632 34,2654] . Nilai matriks kovariansi dihitung dengan menggunakan persamaan 2.5.sehingga diperoleh nilai ⎡ 0,0013 − 0,0006 0,0105⎤ S = ⎢⎢− 0,0006 0,0223 0,0048⎥⎥ . ⎢⎣ 0,0105 0,0048 4,3449⎥⎦
Selanjutnya dihitung nilai T2 untuk 309 sampel dengan menggunakan persamaan 2.6. Sebagai contoh, akan dihitung nilai T2 untuk sampel pertama. Data karakteristik produk untuk sampel pertama dinyatakan dalam X 1 = [8,36 6,25 30] sehingga nilai X 1 − X = [8,36 6,25 30] − [8,3767 6,5632 34,2654] = [− 0,0167 − 0,3132 − 4,2654] . Dengan menggunakan persamaan 2.6 dihitung nilai 2 T1 = ( X 1 − X ) S −1 ( X − X ) ' −1 ⎡ 0,0013 − 0,0006 0,0105⎤ ⎡ − 0,0167 ⎤ = [− 0,0167 − 0,3132 − 4,2654] ⎢⎢− 0,0006 0,0223 0,0048⎥⎥ ⎢⎢ − 0,3132 ⎥⎥ ⎢⎣ 0,0105 0,0048 4,3449⎥⎦ ⎢⎣− 4,2654⎥⎦ = 8,6276 sehingga diperoleh nilai T12 = 8,6276 . Masing masing nilai T2 kemudian dibandingkan dengan nilai BPA sebesar 14,6064. Grafik pengendali Hotelling T2 disajikan dalam Gambar 4.5. Gambar 4.5. menunjukkan bahwa dari data karakteristrik produk ”X” terdeteksi sebanyak dua sampel yang out of control yaitu sampel ke 50 dan 101. Nilai T2 kedua sampel tersebut berturut turut adalah 17,4784 dan 22,7456. Hasil ini bersesuaian dengan hasil yang diperoleh pada grafik Hotelling univariat untuk karakteristik brix dan pH.
Gambar 4.5. Grafik pengendali Hotelling T2 multivariat untuk 3 karakteristik.
Analisis Penyebab Proses di Luar Kendali Pada bagian ini akan diselidiki dari ketiga karakteristik produk tersebut, apakah ada variabel dominan yang mempengaruhi sehingga proses berada di luar kendali. Analisis dilakukan dengan menyusun boxplot data yang telah distandardisasi dengan menggunakan persamaan x ji − xi dengan i = 1,2,...,p dan j = 1,2,...,m (4.1) z ji = si Selanjutnya dari boxplot untuk masing masing karakteristik dapat dianalisis apakah sampel outliers bersesuaian dengan sampel yang teridentifikasi sebagai sampel out of control pada Gambar 4.5. Jika sampel outliers bersesuaian dengan sampel out of control maka karakteristik tersebut merupakan karakteristik dominan yang menyebabkan sampel out of control pada grafik pengendali. Sebaliknya, jika sampel out of control bukan sampel outliers maka dapat disimpulkan bahwa faktor yang menyebabkan sampel out of control merupakan akumulasi dari nilai ketiga karakteristik kualitas tersebut. Gambar 4.6 menunjukkan boxplot dari data proses produksi yang telah distandarisasi dengan menggunakan persamaan 4.1. Selanjutnya diselidiki apakah sampel ke 50 dan 101 yang teridentifikasi sebagai sampel out of control pada grafik pengendali Hotelling T2 multivariat merupakan salah satu sampel outliers pada boxplot tersebut. Diperoleh bahwa kedua sampel tersebut merupakan outliers untuk karakteristik pertama, sedangkan untuk karakteristik lain, kedua sampel tersebut bukan merupakan sampel outliers. Dengan demikian dapat disimpulkan bahwa karakteristik pertama yaitu brix (kadar gula) merupakan faktor penyebab dominan yang mempengaruhi out of control pada produk. Hal ini juga didukung hasil pengamatan grafik pengendali univariat untuk karakteristik pertama (Gambar 4.2). Dari grafik univariat tersebut teridentifikasi dua sampel out of control yaitu sampel ke 50 dan 101 yang juga merupakan out of control pada pengamatan multivariat.
Gambar 4.6 Boxplot data terstandardisasi untuk masing masing karakteristik kualitas produk.
Studi Simulasi Dengan menggunakan nilai rata rata X = [8,3767 6,5632 34,2654] dan kovariansi sampel ⎡ 0,0013 − 0,0006 0,0105⎤ S = ⎢⎢− 0,0006 0,0223 0,0048⎥⎥ , ⎢⎣ 0,0105 0,0048 4,3449⎥⎦ akan dibangkitkan data acak berdistribusi normal multivariat. Dari data acak tersebut akan diaplikasikan grafik pengendali Hotelling T2 multivariat. Simulasi dilakukan beberapa kali dengan menggunakan jumlah sampel yang berbeda beda. Diambil contoh untuk banyak sampel 1000 dan dihitung nilai batas pengendali atas 14,2929. Gambar 4.6 menunjukkan grafik Hotelling untuk banyak sampel 1000, dideteksi tiga sampel yang berada di luar batas yaitu sampel ke-491, ke-584 dan ke-667. Hasil tersebut mendekati nilai teoritis untuk nilai α = 0,0027 untuk 1000 sampel maka akan diperoleh sampel out-of-control sebanyak 0,0027 × 1000 = 2,7 ≈ 3 .
Gambar 4.7. Grafik pengendali Hotelling T2 multivariat untuk data simulasi berukuran 1000 Ringkasan banyaknya sampel (n), nilai batas pengendali atas, serta sampel yang out-ofcontrol disajikan dalam tabel 4.1. Dari tabel tersebut kita dapat membandingkan proporsi out of control dengan jumlah teoritis untuk α = 0,0027 .
n 1000
BPA 14,2929
Jumlah observasi di luar BPA 3
Proporsi out of control 3 1000
= 0,003
= 0,003 = 0,0023
2000
14,2243
6
6 2000
3000
14,2016
7
5000
14,1834
15
7 3000 15 5000
= 0,003
Tabel 4.1. Simulasi untuk n sampel, BPA, dan jumlah observasi di luar BPA
Indeks Kemampuan Proses Univariat Indeks kemampuan proses merupakan indeks yang dapat menggambarkan apakah proses memenuhi spesifikasi atau tidak. Oleh karena itu, data yang digunakan harus terkontrol. Dalam hal ini terkontrolnya data adalah karena telah dilakukan penelusuran penyebab ketidakstabilan proses dan penyebab tersebut telah diperbaiki sehingga dapat diperoleh data yang telah terkontrol baik. Selanjutnya akan dilakukan perhitungan indeks kemampuan proses univariat untuk masing masing karakteristik produk. Sebagai contoh untuk variabel pertama yaitu kadar gula, ditemukan dua sampel yang berada di luar BPA yaitu sampel ke 50 dan 101. Sampel ini selanjutnya dihilangkan dengan asumsi sudah ditelusuri penyebab ketidakstabilan proses. Proses penghilangan sampel out-of-control tersebut diulang hingga semua sampel berada di bawah kendali seperti pada Gambar 4.8. Diperoleh 300 sampel yang berada di bawah batas pengendali atas yang baru sebesar 9,0299. Nilai mean dan standar deviasi sampel berturut turut x = 8,3766 dan s = 0,0316 . Selanjutnya akan dihitung indeks kemampuan proses masing masing untuk Cp, Cpk dan Cpm. USL − LSL Cˆ p = 6s =
8.55 − 8.2 = 1.6009 6 × 0.0312
⎛ USL − x x − LSL ⎞ , Cˆ pk = min⎜ ⎟ 3s ⎠ ⎝ 3s ⎛ 8,55 − 8,3770 8,3770 − 8,25 ⎞ = min⎜ , ⎟ 3 × 0,0312 ⎠ ⎝ 3 × 0,0312 = min(1,8460 ;1,3588) = 1,3588 Cˆ pm = =
USL − LSL 6 s2 + (x − T )2 8 ,55 − 8 , 25 6 ( 0 , 0316 ) 2 + (8,3770 − 8, 4 ) 2
= 1, 2897
Dari ketiga metode perhitungan indeks kemampuan proses univariat tersebut diperoleh hasil indeks kemampuan proses lebih dari 1 menunjukkan bahwa secara univariat kadar gula masih berada di dalam batas spesifikasi yang telah ditentukan perusahaan. Selanjutnya, dilakukan prosedur serupa untuk karakteristik kualitas yang lain. Ringkasan nilai indeks kemampuan proses untuk ketiga karakteristik dapat dilihat pada Tabel 4.2. Karakteristik kualitas Kadar gula (brix) Kadar keasaman (Ph) Vacuum
Cp 1.8485 1.1209 1.6357
Cpk 1.8314 1.0359 1.5217
Cpm 1.8461 1.0861 1.5476
Tabel 4.2. Indeks Kemampuan proses univariat untuk ketiga karakteristik produk
Gambar 4.8. Grafik pengendali Hotelling T2 untuk karakteristik kadar gula setelah sampel out of control dihilangkan Indeks Kemampuan Proses Bivariat Pada bagian ini akan ditunjukkan gambaran geometris indeks kemampuan proses untuk tiap pasangan karakteristik dan dihitung nilai indeks kemampuan proses bivariat. Sebelumnya dilakukan uji korelasi untuk mengetahui apakah masing masing karakteristik saling independen satu sama lain. Uji korelasi menunjukkan hasil yang signifikan untuk karakteristik pertama dan ketiga sehingga dapat disimpulkan kedua karakteristik tersebut saling berkorelasi. Sebaliknya, untuk pasangan karakteristik yang lain tidak berkorelasi sehingga diasumsikan nilai kovariansi s12 = 0. Selanjutnya, akan dianalisis untuk karakteristik pertama (Brix) dan kedua (pH). Batas spesifikasi yang telah ditentukan perusahaan untuk kedua karakteristik tersebut berturut turut
adalah 8,25 – 8,55 dan 6,1 – 7,1. Ellips terbesar yang memenuhi batas spesifikasi tersebut adalah ellips dengan persamaan 2 2 ( ( x − 8,4) y − 6,6) R1 ≡ + = 1. 0,0225 0,25 Daerah R2 merupakan ellips dengan persamaan 2 2 1 (4.1) R2 ≡ 2 2 s 22 x1 − x 1 + s12 x 2 − x 2 = T 2 . s1 s 2 Persamaan 4.1. dapat diubah ke dalam bentuk (x − x1 )2 (x2 − x2 )2 2 + =T (4.2) s12 s 22
[
(
)
(
)]
dengan x1 , x2 merupakan rata rata dan s12 , s 22 variansi berturut turut untuk variabel pertama dan kedua. Persamaan ellips R2 untuk karakteristik pertama dan kedua adalah (x − 8,38)2 + (x2 − 6,5605)2 = 12,0639 . 0,0011 0,0219 2 Nilai T ditentukan dengan menggunakan persamaan ⎛ p(m + 1)(m − 1) ⎞ ⎟⎟ Fα , p,( m− p ) T 2 = ⎜⎜ ⎝ m(m − p ) ⎠ dengan p banyaknya karakteristik dan m banyaknya sampel yang diobservasi sedangkan α dipilih 0,0027. Gambar 4.9 menunjukkan indeks kemampuan proses tersebut secara geometris. Pada gambar 4.9 persegi panjang merah menunjukkan batas spesifikasi yang telah ditentukan perusahaan sedangkan daerah R1 merupakan ellips berwarna biru yang merupakan ellips terbesar yang dibatasi batas spesifikasi untuk kedua karakteristik.Titik pusat ellips tersebut merupakan nilai target (titik tengah batas spesifikasi). R2 merupakan daerah yang dibatasi ellips berwarna hijau menunjukkan daerah proses yang sebenarnya. Dari Gambar 4.9 kita dapat memperkirakan bahwa nilai indeks kemampuan proses lebih besar dari 1. Hal ini diindikasikan dari R2 yang hampir seluruh daerahnya berada dalam R1. Selain itu, semua sampel jatuh di dalam daerah R1 yang menunjukkan bahwa semua sampel memenuhi spesifikasi perusahaan. Secara matematis, nilai indeks kemampuan proses bivariat dihitung dengan menggunakan persamaan (2.12) dengan Cˆ p = 1,2785 dan Dˆ = 1,2312 Cˆ p 1,2785 MCˆ pm = = = 1,0385 . 1,2312 Dˆ Nilai IKP untuk kedua karakteristik tersebut menunjukkan nilai lebih besar dari 1 sehingga dapat disimpulkan bahwa dilihat dari kedua karakteristik tersebut proses sudah berjalan dengan baik.
Gambar 4.9. Gambaran geometris untuk perbandingan antara daerah toleransi spesifikasi dan daerah proses yang sebenarnya untuk karakteristik pertama dan kedua Selanjutnya akan dihitung nilai IKP untuk karakteristik pH dan vacuum dengan spesifikasi perusahaan masing masing 6,1 – 7,1 dan 25 – 45. Ellips terbesar yang memenuhi batas spesifikasi tersebut adalah ellips dengan persamaan (x − 6,6 )2 + ( y − 35)2 = 1 . 0,25 100 Dengan menggunakan persamaan (4.2) kita dapat menentukan daerah R2 yaitu ellips dengan persamaan (x − 6,5605)2 + ( y − 34,3191)2 = 12,0639 . 0,0219 4,1916 Gambar 4.10. menunjukkan area yang dibatasi oleh spesifikasi perusahaan dan daerah proses yang sebenarnya. Dari gambar 4.10 kita dapat memperkirakan bahwa nilai indeks kemampuan proses lebih besar dari 1. Hal ini diindikasikan dari R2 yang hampir seluruh daerahnya berada dalam R1. Selain itu, semua sampel jatuh di dalam daerah R1 yang menunjukkan bahwa semua sampel memenuhi spesifikasi perusahaan. Secara matematis, nilai indeks kemampuan proses bivariat dihitung dengan menggunakan persamaan 2.12 dengan Cˆ p = 1,3950 dan Dˆ = 1,0840 .
Gambar 4.10. Gambaran geometris untuk perbandingan antara daerah toleransi spesifikasi dan daerah proses yang sebenarnya untuk karakteristik kedua dan ketiga. Cˆ p 1,3950 = = 1,2870 1,0840 Dˆ Selanjutnya akan dihitung nilai IKP untuk karakteristik brix dan vacuum dengan spesifikasi perusahaan masing masing 8,25 – 8,55 dan 25 – 45. Ellips terbesar yang memenuhi batas spesifikasi tersebut adalah ellips dengan persamaan (x − 8,4 )2 + ( y − 35)2 = 1 . 0,0225 100 Daerah R2 merupakan ellips dengan persamaan 1 ⎡s 2 x − x 2 + s 2 y − x 2 − 2 s 2 x − x y − x ⎤ = T 2 (4.3) 1 2 1 2 1 12 2 ⎢ ⎥⎦ ⎣ 2 s12 s 22 − s12 sehingga R2 untuk karakteristik brix dan vacuum dapat ditentukan sebagai berikut 2 2 4,1916 ( x1 − 8,38) + 0,0011( x 2 − 34,32) − 0,0116(x1 − 8,38)(x 2 − 34,32) = −0,0132 . Dengan cara serupa diperoleh nilai Cˆ p 2,1637 MCˆ pm = = = 2,0390 . 1,0612 Dˆ Dari nilai IKP tersebut kita dapat menyimpulkan bahwa untuk karakteristik brix dan vacuum proses berada dalam spesifikasi perusahaan. Ketiga hasil IKP bivariat tersebut menunjukkan nilai lebih besar dari 1 sehingga dapat disimpulkan bahwa proses telah berjalan dengan baik. MCˆ pm =
(
)
(
)
(
)(
)
Indeks Kemampuan Proses Multivariat Pada penelitian ini, data yang digunakan untuk menghitung indeks kemampuan proses diperoleh dengan menghilangkan sampel yang terdeteksi out of control. Hal ini dilakukan dengan asumsi sudah dilakukan penelusuran penyebab ketidakstabilan proses dan prosesnya telah diperbaiki. Dari Gambar 4.5 terdeteksi adanya sinyal out of control sebanyak 2 sampel yaitu sampel ke 50 dan 101. Kemudian 2 sampel out of control tersebut dihilangkan dan selanjutnya dilakukan pengontrolan yang baru menggunakan grafik pengendali Hotelling T2. Proses ini diulang hingga diperoleh semua nilai sampel berada di dalam batas pengendali grafik. Proses pengendalian, data sampel out of control dan batas pengendali yang baru dapat dilihat pada Tabel 4.3.
Pengontrolan ke 1 2 3 4
Banyaknya out of control 2 2 1 0
Nomor sampel out of control 50,101 100, 143 142 -
BPA 14,6064 14,6094 14,6125 14,6140
Tabel 4.3. Proses penghilangan sampel out of control Setelah dilakukan pengontrolan keempat kalinya, diperoleh 304 sampel yang berada di bawah grafik pengendali. Dari ke 304 data sampel tersebut akan dilakukan perhitungan indeks kemampuan proses multivariat. Ellipsoid terbesar yang memenuhi daerah spesifikasi adalah ellipsoid dengan persamaan 2 2 2 ( ( ( x − 8,4) y − 6,6) z − 35) R1 ≡ + + =1 0,0225 0,25 100 sedangkan ellipsoid proses yang sebenarnya ditentukan oleh persamaan (x − 8,38)2 + ( y − 6,5605)2 + (z − 34,32)2 = 14,4703 . 0,0011 0,0219 4,1916 Gambar 4.10 menunjukkan dua ellipsoid yang mewakili daerah batas spesifikasi dan daerah 99,73 % proses. Daerah R1 (batas spesifikasi) ditunjukkan oleh ellipsoid wiregrid berwarna sedangkan daerah proses R2 ditunjukkan oleh ellipsoid patch berwarna abu abu. Dari Gambar 4.10 terlihat daerah R2 hampir seluruhnya berada di dalam ellipsoid R1. Melihat fakta tersebut, kita dapat memperkirakan nilai indeks kemampuan proses lebih besar dari 1. Secara matematis, dengan menggunakan persamaan (2.11) diperoleh nilai p
vol ( R1 ) = =
2∏ Ti i =1
p
π p/2
Γ( p / 2 )
2 × 0,75 π 3 / 2 3 0,8862
= 3,1416.
Selanjutnya nilai Cˆ p dihitung sebagai berikut
Cˆ p = =
vol ( R1 )
S
1/ 2
(πK ( p)) p / 2 [Γ( p / 2 + 1)]
−1
3,1416
0,0101(π × 14,1563)
3/ 2
0,7523
3,1416 = 1,3935 2,2544 Dengan menggunakan persamaan (2.12) dihitung nilai indeks kemampuan proses Cˆ p 1,3935 = = 1,1066 . MCˆ pm = 1,2593 Dˆ Terlihat bahwa nilai indeks kemampuan proses diperoleh lebih dari 1. Hal ini menunjukkan variasi proses lebih kecil daripada batas spesifikasi perusahaan sehingga dapat disimpulkan bahwa proses produksi ”X” sudah dalam keadaan baik. =
Gambar 4.10. Gambaran geometris untuk perbandingan antara daerah toleransi spesifikasi dan daerah proses yang sebenarnya untuk ketiga karakteristik kualitas.
5. SIMPULAN DAN SARAN SIMPULAN 1. Aplikasi grafik pengendali univariat untuk masing masing karakteristik, menunjukkan proses terkendali untuk karakteristik pH, sedangkan untuk karakteristik brix terdapat dua sampel yang di luar kendali yaitu sampel ke 50 dan 101. Grafik univariat untuk karakteristik vacuum menunjukkan satu sampel yang di luar kendali yaitu sampel ke 57.
2. Penerapan grafik pengendali Hotelling T2 Multivariat untuk proses produksi ”X” selama bulan Januari 2010 – Agustus 2010 menunjukkan terdapat 2 sampel yang berada di luar kontrol dengan batas pengendali atas 14,6064. 3. Dengan menggunakan rata-rata dan matriks kovariansi data produksi, dapat dibangkitkan data acak berdistribusi normal multivariat untuk ukuran sampel yang berbeda dan diperoleh proporsi banyaknya sampel out of control mendekati nilai α = 0,0027 . 4. Analisis sinyal out of control dilakukan dengan membuat boxplot untuk tiap karakteristik, selanjutnya dari analisis sampel outliers diperoleh karakteristik yang secara dominan mempengaruhi ketidaknormalan proses adalah karakteristik brix (kadar gula). 5. Nilai indeks kemampuan proses univariat menunjukkan nilai IKP > 1 sehingga secara univariat proses produksi telah berjalan dengan baik. Demikian pula, untuk pengamatan multivariat, nilai indeks kemampuan proses multivariat MCpm untuk proses produksi ”X” adalah sebesar 1,106. Hal ini menunjukkan proses telah memenuhi spesifikasi perusahaan dengan asumsi telah dilakukan penyelidikan terhadap penyebab ketidakstabilan proses dan dilakukan perbaikan. SARAN Dalam menghitung nilai indeks kemampuan proses sebaiknya benar – benar dilakukan penelusuran terhadap penyebab ketidakstabilan proses selanjutnya diperbaiki sehingga nilai indeks kemampuan proses yang diperoleh lebih mencerminkan keadaan proses yang sebenarnya, bukan sekedar menghilangkan sampel-sampel yang out of control. 6. DAFTAR PUSTAKA
Johnson, Richard. Dean Wichern. 2007. Applied Multivariate Statistical Analysis, 6th ed. New Jersey : Prentice Hall. Linn, Richard. 2006. Supplier Selection Based on Process Capability and Price Analysis. Quality Engineering Vol 18 : 123 – 129. M. Abdolshah. 2009. Overall Processes Capability Index for Assembly Production Lines. Journal of Applied Science 9 (20) : 3764 – 3769. Montgomery, D.C. 1990. Pengantar Pengendalian Kualitas Statitik. Alih bahasa : Zanzawi Soejoeti. Yogyakarta : Universitas Gadjah Mada. Young, Timothy. 1999. Multivariate Control Charts of MDF and OSB Vertical Density Profile Attributes. Forest Products Journal.Vol 49 : 79 – 86. Zahid, Abu. Arifa Sultana. 2008. Assesment and Comparison of Multivariate Process Capability Indices in Ceramic Industry. Journal of Mechanical Engineering Vol. ME39 : 18 – 25. Web 1 : Hotelling Control Charts. http://www.itl.nist.gov/div898/handbook/pmc/section3 diakses tanggal 15 Oktober 2010. Web 2 : What is Process Capability?. http://www.itl.nist.gov/div898/handbook/pmc/section1/pmc16.htm diakses tanggal 19 Desember 2010.