ANALISIS PENGARUH LIMBAH SUSU SAPI CAIR MT.GKSI BOYOLALI TERHADAP PERTUMBUHAN DAN HASIL TANAM TANAMAN BAWANG MERAH ( Allium cepa var.ascallonicum )

## Edy Sulistiyawan
Program Studi Statistika Universitas PGRI Adi Buana Surabaya Alamat e-mail : edy.sulistiyawan@gmail.com

## ABSTRAK
Untuk mengurangi pencemaran lingkungan akibat pengolahan limbah yang kurang tepat dapat dilakukan pemanfaatan limbah tersebut sebagai pupuk, terutama limbah organik yang masih banyak mengandung unsur Nitrogen, Fosfat dan Kalium fosfat . MT.GKSI Boyolali merupakan satu-satunya pabrik pengolahan susu sapi cair ( cilling ) yang telah melakukan uji kualitas air susu sebelum masuk ke pabrik pengolahan susu. Limbah susu yang dihasilkan berwarna putih dan apabila keberadaannya terlalu banyak tentu dapat merusak lingkungan karena dapat menimbulkan bau yang menyengat. Dalam penelitian ini dapat diketahui pengaruh limbah susu sapi cair dan kadar limbah susu sapi cair yang optimal bagi pertumbuhan dan hasil tanam tanaman bawang merah. Dari hasil pengolahan data dengan menggunakan analisis varians multivariate dapat disimpulkan bahwa konsentrasi limbah susu sapi cair yang berbeda-beda ternyata memberikan pengaruh yang berarti terhadap hasil tanam tanaman bawang merah. Dan secara univariate efek konsentrasi limbah susu sapi cair memberikan pengaruh yang berarti terhadap ketiga respon hasil tanam yaitu berat basah, berat kering dan jumlah umbi perumpun. Pengujian perbedaan rata-rata masing-masing perlakukan, dimana metode yang digunakan adalah Newman-Keuls didapatkan bahwa limbah susu sapi cair akan memberikan pengaruh yang optimal pada konsentrasi sebesar 60%. Sedangkan untuk data yang pengamatannya berdasarkan tiap minggunya, menunjukkan bahwa ada pengaruh pengamatan pada waktu yang berbeda terhadap variabel respon pertumbuhan apabila diperlakukan pada konsentrasi limbah yang berbeda-beda. Demikian pula secara univariate juga ada pengaruh pengamatan pada waktu yang berbeda pada respon tinggi tanaman dan jumlah daun apabila diperlukan pada konsentrasi yang berbeda-beda.
Kata Kunci : Analisis Varians Multivarite, Limbah Susu Sapi Cair, Newman-Keuls , Tanaman Bawang Merah.

## PENDAHULUAN
Indonesia sebagai Negara yang masih dalam taraf berkembang, dituntut untuk lebih aktif mengadakan pembangunan di berbagai sektor, untuk mengejar
ketinggalalan dari bangsa yang lebih maju. Hal ini dimaksudkan untuk mewujudkan Tujuan Pembangunan Nasional yang tertuang dalam GBHN Tap MPR No. II/MPR/1988 butir 1 yang berbunyi: “ Mewujudkan suatu
masyarakat yang makmur yang merata materiil spirituil berdasarkan pancasila”. Pembangunan jangka panjang di Indonesia dititik beratkan pada bidang ekonomi, hal ini dikarenakan, dengan peningkatan hasil- hasil dalam bidang ekonomi akan tersedia sumber-sumber pembangunan akan lebih luas bagi bidang yang lain. Sasaran dalam bidang ekonomi adalah struktur ekonomi yang seimbang dengan titik berat kekuatan dalam bidang industri
didukung oleh bidang pertanian yang kuat.
Dari hasil pembangunan dibidang ekonomi khususnya sektor industri telah banyak kita rasakan manfaatnya akan tetapi pembangunan di sektor industri tidak hanya menghasilkan manfaat, melainkan juga membawa resiko. Resiko yang ditimbulkan oleh industrialisasi adalah berupa akses-akses yang kurang baik bagi lingkungan, yaitu turunnya kualitas lingkungan yang dapat menjadi bumerang bagi kehidupan manusia itu sendiri. Dengan turunnya kualitas lingkungan, maka kebutuhan dasar untuk hidup yang wajar tidak terpenuhi. Dalam era industrialisasi saat ini perlu adanya perencanaan dan usaha untuk melestarikan alam atau sekurang-
kurangnya menurunkan akses negatif terhadap lingkungan, antara lain berupa pengendalian limbah sebelum memasuki suatu komponen lingkungan. Pengendalian tersebut dapat dilakukan secara tradisional, teknik madya atau teknik mesin. Untuk mengurangi pencemaran lingkungan akibat pengolahan limbah yang kurang tepat dapat dilakukan pemanfaatan limbah tersebut semaksimal mungkin sebagai sumber energi yang dapat diperbaharui. Pemanfaatan limbah tersebut antara lain digunakan sebagai pupuk terutama limbah organik yang masih banyak mengandung unsur Nitrogen, Fosfat dan Kalium Fosfat . MT. GKSI Boyolali merupakan satu- satunya pabrik pengolahan susu sapi cair ( cilling ) yang telah melakukan uji kualitas air susu sebelum masuk pabrik pengolahan susu, yang berada di jawa tengah. Limbah susu cairnya bewarna putih, berasal dari sisa atau pembuangan air susu yang tidak lolos seleksi seleksi dalam uji kualitas. Limbah ini apabila keberadaannya terlalu banyak tentu akan dapat merusak lingkungan karena dapat menimbulkan bau yang menyengat. Hal ini disebabkan limbah air susu tersebut
masih mengandung unsur kamia organik dalam berbagai kosentrasi.
Masalah yang harus diselesaikan dalam penelitian ini adalah diketahuinya pengaruh limbah susu sapi cair MT. GKSI Boyolali tersebut sebagai pupuk dalam budidaya tanaman bawang merah ( Alium cepa ) dan pada konsentrasi berapa limbah tersebut dapat menghasilkan atau memberikan pengaruh yang optimum pada tanaman bawang merah. Penelitian ini bertujuan untuk mengetahui pengaruh limbah susu sapi cair MT. GKSI Boyolali terhadap pertumbuhan dan hasil tanam tanaman bawang merah ( Allium cepa ) dan untuk mengetahui kadar limbah susu sapi cair yang potensi bagi pertumbuhan dan hasil tanam tanaman bawang merah.
Manfaat yang diharapkan agar dapat diketahui pengaruh limbah susu sapi cair terhadap pertumbuhan tanam tanaman bawang merah, sehingga nantinya dapat memberikan tambahan informasi bagi para petani bawang merah. Dalam penelitian ini, peneliti membatasi permasalahan pada penelitian lapangan yang dilakukan di rumah plastik, Ds. Poncol, kecamatan Poncol Magetan.

## METODE PENELITIAN
Penelitian limbah susu sapi cair
MT.GKSI Boyolali dimana lingkup penelitin ini meliputi pengamatan
terhadap variabel respon. Adapun variabel-variabel respon tersebut yaitu : tinggi tanaman, jumlah daun, berat basah, berat kering, jumlah umbi perumpun. Pengamatan terhadap variabel respon pertumbuhan tanaman bawang merah dilakukan selama 9 minggu sedangkan pengukuran variabel respon hasil tanam tanaman bawang merah dilakukan setelah panen. Dalam penelitian ini pengolahan data dilakukan dengan menggunakan paket progam SPSS dan Minitab.
Untuk mengetahui pengaruh limbah susu sapi cair terhadap hasil tanam tanaman bawang merah yaitu pada respon berat basah, berat kering, dan jumlah umbi perumpun digunakan metode manova satu arah. Dimana terdapat faktor yang masing-masing
terdiri dari enam level kosentrasi dengan enam kali pengulangan, yaitu : Kosentrasi limbah 0%, 20%, 40%, 60%, 80%, 100%.
Tabel 1. Pengorganisasian Data Rancangan Manova Perlakuan (Kosentrasi) 0% 20% 40% 60% 80% 100% peng amat an ke-1 Y 111 Y 211 Y 311 Y 411 Y 511 Y 611 Y 112 Y 212 Y 312 Y 412 Y 512 Y 612
Y 113 Y 212 Y 313 Y 413
Y 513 Y 613 peng amat an ke-2 Y 121 Y 221 Y 321 Y 421 Y 521 Y 621 Y 122 Y 222 Y 322 Y 422 Y 522 Y 622 Y 123 Y 223 Y 323 Y 423 Y 523 Y 623 …
… … …
… … …
peng amat an ke-5
Y 151 Y 251 Y 351 Y 451 Y 551 Y 651 Y 152 Y 252 Y 352 Y 452 Y 552 Y 652 Y 153 Y 253 Y 353 Y 453 Y 553 Y 653
dimana :
Y ij1 = berat basah tanaman bawang merah pada kosentrasi ke-i pada pengamatan ke-j. Y ij2 = berat kering tanaman bawang merah pada kosentrasi ke-i pada pengamatan ke-j. Y ij3 = jumlah umbi perumpun tanaman bawang merah pada kosentrasi ke-i pada pengamatan ke-j
Sedangkan untuk mengetahui pada kosentrasi berapa limbah susu sapi cair akan memberikan pengaruh optimum pada pertumbuhan dan hasil tanam tanaman bawang merah dilakukan dengan menggunakan metode Newman –
Keuls .
Kemudian untuk mengetahui pengaruh limbah susu sapi cair terhadap pertumbuhan tanaman bawang merah dilakukan penelitian dengan menggunakan metode anacova satu arah dengan satu variabel covariate yaitu waktu pengamatan tiap minggunya dengan enam faktor masing-masing terdiri dari 2 variabel respon, yaitu tinggi tanaman dan jumlah daun.
Tabel 2. Managemen Data Rancangan Manacova Va
ria
bel Co ver iat e Perlakuan (FAKTOR) 1 2 . . . 6 1 2 . . .
9
Y ଵଵଵ ଵ Y ଵଵଵ ଶ Y ଶଵଵ ଵ Y ଶଵଵ ଶ . . . Y ଺ଵଵ ଵ Y ଺ଵଵ ଶ Y ଵଵଶ ଵ Y ଵଵଶ ଶ Y ଶଵଶ ଵ Y ଶଵଵ ଶ . . . Y ଺ଵଶ ଵ Y ଺ଵଶ ଶ . . . . . . . . . . . . . . . . . . . . . Y ଵଵହ ଵ Y ଵଵହ ଶ Y ଶଵହ ଵ Y ଶଵହ ଶ . . . Y ଺ଵହ ଵ Y ଺ଵହ ଶ Y ଵଶଵ ଵ Y ଵଶଵ ଶ Y ଶଶଵ ଵ Y ଶଶଵ ଶ . . . Y ଺ଶଵ ଵ Y ଺ଶଵ ଶ
Y ଵଶଶ ଵ Y ଵଶଶ ଶ Y ଶଶଶ ଵ Y ଶଶଶ ଶ . . . Y ଺ଶଶ ଵ Y ଺ଶଶ ଶ . . . . . . . . . . . . . . . . . . . . . Y ଵଶହ ଵ Y ଵଶହ ଶ Y ଶଶହ ଵ Y ଶଶହ ଶ . . . Y ଺ଶହ ଵ Y ଺ଶହ ଶ . . . . . . . . . . . . . . . . . . . . . Y ଵଽଵ ଵ Y ଵଽଵ ଶ Y ଶଽଵ ଵ Y ଶଽଵ ଶ . . . Y ଺ଽଵ ଵ Y ଺ଶହ ଶ Y ଵଽଶ ଵ Y ଵଽଶ ଶ Y ଶଽଶ ଵ Y ଶଽଵ ଶ . . . Y ଺ଽଶ ଵ Y ଺ଽଶ ଶ . . . . . . . . . . . . . . . . . . . . . Y ଵଽହ ଵ Y ଵଽହ ଶ Y ଶଽହ ଵ Y ଶଽହ ଶ . . . Y ଺ଽହ ଵ Y ଺ଽହ ଶ Dimana :
Y ୧୨୩ ଵ = Tinggi Tanaman kosentrasi ke-I pada variabel covariate ke-j, pengamatan je-k. Y ୧୨୩ ଶ = Jumlah Daun kosentrasi ke-I pada variabel Covariate ke-j, pengamatan ke-k.

## HASIL PENELITIAN

## Deskriptif data
Tabel 3. Deskriptif Data Pertumbuhan
variabel Mean Std Dev Minimum Maximum V3 29.55 6.53 16.30 47.50 V4 12.19 3.34 7.66 21.12 V5 9.53 1.91 6.00 14.00 Dari Tabel 3 dapat diketahui bahwa untuk variabel berat basah pola pemusatan datanya terletak pada nilai mean sebesar 29.55 dan pola penyebaran data yang tidak lebar dimana nilai standar deviasinya 6.53 jadi pola simetrinya. Varibel berat kering mempunyai pola pemusatan data terletak pada nilai mean sebesar 12.19 pola penyebaran data yang tidak lebar pada nilai standar deviasi sebesar 3.34 dan polanya simetris. Sedangkan untuk variabel jumlah umbi perumpun mempunyai pola pemusatan dengan nilai rata-rata yang terkecil dari lainnya yaitu sebesar 9.53 dengan pola penyebaran datanya lebih kecil dari lainnya pula yaitu dengan nilai standar
deviasi sebesar 1.91 dan pola data simetris.

## Pengujian Secara Serentak
Analisis Manova digunakan untuk
mengetahui efek perlakuan kosentrasi limbah yang berbeda terhadap pertumbuhan tanaman bawang merah. Adapun model analisis varians
multivariate satu arah dalah sebagai berikut :
Y ୧୨୩ = μ ୩ + ߬ ୧୩ + ߳ ୧୨୩
Untuk : i = 1, 2, 3, 4, 5, 6 j = 1, 2, 3, 4, 5 k = 1, 2, 3
dengan : i = indeks yang menunjukkan kosentrasi limbah
j = indeks yng menunjukkan banyaknya ulangan k = indeks yang menunjukkan
respon yang diamati
dimana :
Y ijk = respon ke-k yang di amati pada kosentrasi limbah yang ke-i dan ulangan ke-j μ k = rata-rata dari respon ke-k т ij = efek dari setiap kosentrasi limbah pada setiap pengamatan terhadap respon
Tabel 4. Pengujian Statistik secara serentak untuk respon hasil tanam Test Name Value Approx. F Hypoth. DF Error DF Sig. of F Piilais Hotellings Wilks Roys 1.25049 7.25595 .08396 .87134 3.43089 9.99709 5.92283 15.00 15.00 15.00 72.00 62.00 61.13 .000 .000 .000 Hipotesa : Ho = tidak ada pengaruh penigkatn
kosentrasi limbah susu sapi cair trhadap ketiga variabel respon.
H1 = ada pengaruh peningkatan kosentrasi limbah susu sapi cair terhadap ketiga variabel respon. Statistik uji : Wilks = 0.08396 f = 5.92283 F tabel = F (5,24,5%) = 2.62000 Karena F-hitung > F-tabel, maka tolak H0
Artinya : ada efek kosentrasi limbah susu sapi cair terhadap pertumbuhan tanamn bawang merah
a. Pengujian Secara Individu
Model : Y ijk = μ k + т ik + ϵ ijk
Tabel 5. Pengujian Statistik Secara Individu Untuk Setiap Respon Hail Tanam
Variab el Hypoth . SS Error SS Hypoth . MS Error MS F Sig. of F V3 V4 V5 946.78 246.15 39.867 290.22 78.143 65.600 189.35 49.230 7.9733 12.092 3.255 2.733 15.66 15.12 2.917 .000 .000 .034 Hipotesa :
H0 : т 1k = т 2k = т 3k = . . . = т 6k = 0 H1 : paling tidak ada satu nilai T ik tidak sama dengan nol
Statistik uji :
Untuk variabel respon V3 (berat bash)
F hitung = 15.65884
F tabel = 2.62
Untuk variabel respon V4 (berat kering)
F hitung = 15.11998
F tabel = 2.62 Untuk variabel respon V5 (jumlah umbi perumpun)
F hitung = 2.91707
F tabel = 2.62
Kesimpulan :
1. Untuk variabel respon berat basah F hitung > F tabel maka tolak H0 2. Untuk variabel respon berat kering F hitung > F tabel maka tolak H0 3. Untuk variabel respon jumlah umbi perumpun F hitung > F tabel maka tolak H0 Interprestasi : Dari pengujian yang telah dilakukan maka dapat diambil kesimpulan bahwa efek kosentrasi limbah susu sapi cair memberikan pengaruh berarti terhadap ketiga variabel respon hasil tanam tanaman bawang merah. Pengujian Perbandingn rata-rata antar perlakuan dengan Metode Newman – Keuls
Setelah diketahui adanya efek kosentrasi limbah terhadap semua respon, maka dengan menggunakan analisis
Newman – Keuls akan diketahui perbandingan rata-rata antar kosentrasi untuk masing-masin respon. Sehingga dapaat diketahhui pada kosentrasi berapakah yang memberikan hasil yang optimum. Adapun hasil dari output didapatkan : a. Untuk respon tinggi tanaman
didapatkan bahwa kosentrasi 100% yang paling optimal.
b. Untuk respon jumlah
daun
didapatkan bahwa kosentrasi 60% yang optimal.
c. Untuk respon berat basah
didapatkan bahwa kosentrasi 60% yang optimal.
d. Untuk respon berat kering bahwa hasil limbah pada kosentrsi 40% berbeda denga hasil limbah pada kosentrasi 0%, sedangkan pada hasil limbah pada kosentrasi 60% berbeda dengan hasil limbah pada kosentrasi 0%, 100%, 20%, 80%, dan 40%.
e. Untuk respon jumlah umbi permpun didapatkan bahea kosentrasi 0% yang optimum.

## Kesimpulan secara keseluruhan
Hasil dari pengujian Newman – Keuls dapat diambil satu kesimpulan bahwa pada kosentrasi 60% limbah susu sapi cair akan memberikan pengaruh yang optimum terhadap pertumbuhan dan hasil tanam tanaman bawang merah.
Pengujian Efek Kosentrasi Terhadap Respon Pertumbuhan Apabila Diamati Pada Waktu Yang Berbeda Menggunakan Analisis Covarian Multivariate Deskriptif data Tabel 6. deskriptif data pertumbuhan Variable Mean Std Dev Minimum Maximum V1 V2 30.61 27.62 14.19 14.53 2.26 1.00 47.28 50.00
Dari Tabel 6 diatas dapat diketahui bahwa untuk variabel tinggi tanman pola pemusatan datanya terletak pada nilai sebesar 30.61 dan pola penyebaran data tidak lebar damana pola standar deviasinya 14.19 jadi polanya simetri. Variabel jumlah daun pola pemusatan datanya terletak pada mean sebesar 27.62 pelebaran data yang tidak lebar pada nilai standar deviasi sebesar 14.53 dan polanya simetris. Data penelitian untuk hasil tanam pengukurannya tidak dipengaruhi oleh waktu yang berbeda, sedangkan untuk data pertumbuhan pengukurannya berdasarkan waktu. Jadi untuk menguji kosentrasi limbah terhadap respon tinggi tanaman dan jumlah daun apabila diamati pada lama waktu yang berbeda
digunakan metode analisis covariate multivariate.

## Pengujian Secara Serentak
Tabel 7. Pengujian Stastik secara serentak untuk respon pertumbuhan Test Name Value Approx. F Hypoth. DF Error DF Sig. of F Pillais Hotellings Wilks Roys .3065 .4019 .7051 .2623 9.52184 10.48890 10.00518 10.0 10.0 10.0 526 522 524 .000 .000 .000 Untuk mengetahui apakah pengaruh terhadap variabel covariate terhadap
10ariable yang diamati, dilakukan pengujian hipotesisi sebagai berikut: Hipotesis : H 0 = Γ = 0 , berarti semua koefisien regresi bernilai nol, sehingga variabel covariate tidak berpengaruh pada respon pertumbuhan.
H 1 = Γ = 0 , berarti terdapat sedikitnya satu keofisien regresi yang bernilai tidak nol, sehinggavariabel covariate berpengaruh pada respon pertumbuhan.
Statistik uji :
Dari hasil pengolahan diperoleh nilai statistik uji :
Wilks = 0.70505 F = 10.00518
F tabel = F ( ୒ି୲ , ୒ି୲ିଵ , α ) = F ( ହ , ଶ଺ଷ , ଴ . ଴ହ ) = 2.21
F ୦୧୲୳୬୥ > F ୲ୟୠୣ୪ maka tolak H0. Karena F rasio > F tabel maka dapat disimpulkan bahwa ada pengaruh pengamatan pada waktu yang berbeda terhadap variabel respon apabila diperlakukan pada kosentrasi limbah yang berbeda beda.
Variabel respon disini adalah variabel pertumbuhan yaitu tinggi tanaman dan jumlah daun. Jadi apabila dilakukan pengamatan
setelah periode waktu
tertentu maka pertumbuhan veriabel respon bertambah sampai batas tertentu.

## Pengujian Secara Individu
Tabel 8. Pengujian Statistik secara individu untuk respon pertumbuhan
Vari able Hypoth. SS Error SS Hypot. MS Error MS F Sig. of F V1 V2 628.372 222.666 9988.732 1831.506 125.674 44.533 37.979 6.963 3.31 6.39 .006 .000 Hipotesa :
H0 : variabel covariate tidak berpengaruh pada variabel respon.
H1 : variabel covariate berpengaruh pada variabel respon.
Statistik uji : Untuk variabel respon V1 (tinggi tanaman) F hitung = 3.30897 F tabel = 2.210
Untuk variabel respon V2 (jumlah daun) F hitung = 6.39488
F tabel = 2.210 Kesimpulan :
1. Untuk variabel respon tinggi
tanaman F hitung > F tabel maka tolak H0 2. Untuk variabel respon jumlah daun F hitung > F tabel maka tolak H0
Dari pengujian yang telah dilakukan dapat disimpulkan bahwa ada pengaruh pengamatan pada waktu yang berbeda terhadap respon tinggi tanaman dan
jumlah daun aabila dilakukan perlakuan pada kosentrasi yang berbeda.
Uji Kesamaan matriks varians kovarians Asumsi yang haris terpenuhi dalam analisis covarians multivariate adalah matriks varians kovarians harus homogen. Pengujian ini dadasarkan atas hipotesa awal bahwa matris varians kovarians sampel ke-1 sampai dengan matriks varians kovarians sampel ke-2 dan seterusnya. Dan disimpulkan bahwa matriks varians kovarians sampel homogen.

## KESIMPULAN
Berdasarkan hasil dan pembahasan yang telah diuraikan diatas, maka dapat diambil beberapa kesimpulan antara lain: 1. Kosentrasi limbah susu sapi cair yang berbeda-beda ternyata memberikan pengaruh yang berarti terhadap pertumbuhan dan hasil tanam tanaman bawang merah.
2. Efek kosentrasi limbah susu sapi cair memberikan pengaruh yang berarti secara univariate terhadap variabel respon hasil tanam tanaman bawang merah.
3. Dari hasil analisis Newman – Keuls didapatkan bahwa pemanfaatan limbah susu sapi cair sebagai pupuk akan menghasilkan pertumbuhan dan hasil tanam yang optimum pada kosentrasi 60% (60% disini merupakan campuran antara 600cc limbah susu sapi cair dengan 400cc air murni/aqua).
4. Ada pengaruh pengamatan pada waktu yang berbeda terhadap variabel respon pertumbuhan tanaman bawang merah apabila diperlakukan pada kosentrasi limbah yang berbeda-beda.
5. Secara univariate juga ada pengaruh pengamatan pada waktu yang berbeda pada respon tinggi tanaman
dan jumlah daun
apabila diperlakukan pada kosentrasi yang berbeda-beda.

## DAFTAR PUSTAKA
[1] Draper, Norman R and Smith Harry,
1992, Analisis Regresi Terapan , PT. gramedia pustaka utama, jakarta,
edisi ke dua.
[2] Budi Samadi, Ir. Bambang Cahyono, 1996, Intensifikasi Budidaya Bawang Merah , Kanisius, Yogyakarta.
[3] Johnson, Richard. A and Wichern, Dean. W., 1988, Applied Multivariate Statistical Analysis , Prentice Hall Inc. Engglewood Cliffs, New Jersey.
[4] J. Supranto, 1983, Ekonometrika , buku satu, Lembaga Pebit Fakultas Ekonomi Universitas Indonesia.
[5] Katz, M., 1971, Effeck Of Pollution Upon Aquatic Life , dalam buku
‘Weste And Water Pollution Hand Book’, vol. 1, L. L. Ciaccio ed,
Marcel Dekker. Inc, New York.
[6] Montgomery, Douglas C, 1970, Design Analysis Of Experiments ,
Georgia Institute Of Technology, John Willey and Sons. [7] Morrison, Donald F, 1976, Multivariate Statistical Method ,
Second Edition, Mc Graw Hill. Inc.
[8] Sudjana, 1995, Desain dan Analisis Esperimen , Tarsito, Bandung.
[9] Sudjana, 1975, Metode Statistika , Edisi ke 5, Tarsito, Bandung.
[10] Stock, H and Mark W.W., 2003, Introduction
to Econometrics ,
Addison Wesley, New York.