ESTIMASI PARAMETER REGRESI ROBUST

DENGAN METODE ESTIMASI-S

PADA PRODUKSI JAGUNG DI INDONESIA TAHUN 2010Oleh Retno Jati Sahari

Fakultas Matematika dan Ilmu Pengetahuan Alam

Universitas Sebelas MaretABSTRAKProduksi jagung di Indonesia menunjukan peningkatan dari tahun ke tahun, tetapi Indonesia masih harus melakukan impor jagung. Ada 2 variabel penting yang mempengaruhi produksi jagung nasional di Indonesia yaitu luas lahan (Ha) dan produktivitas. Untuk menganalisis hubungan antara variabel prediktor dan variabel respon (produksi) diperlukan suatu metode, salah satunya yaitu analisis regresi. Melalui analisis regresi dapat ditentukan model prediksi nilai produksi jagung serta seberapa besar pengaruh variabel-variabel prediktor. Namun, dalam melakukan regresi persyaratan yang harus dipenuhi yaitu asumsi-asumsi klasik yang terkandung pada data. Ada 4 asumsi klasik yang harus dipenuhi yaitu normalitas, homokesdasitas, non autokorelasi dan non multikolinearitas. Apabila salah satu asumsi tidak dipenuhi maka model regresi yang dihasilkan tidak valid. Oleh karena itu diperlukan suatu metode untuk mengatasi data yang tidak memenuhi asumsi-asumsi. Pada penelitian ini digunakan data produksi jagung 33 provinsi di Indonesia. Regresi Robust estimasi-S salah satunya untuk mengatasi asumsi normalitas yang tidak dipenuhi data produksi jagung tahun 2010. Kelebihan dari estimasi-S dibanding estimasi M, LTS, GS, dan GM yaitu estimasi-S dapat mencapai breakdown point hingga 50% maka estimasi-S dapat mengatasi setengah dari pencilan dan memberikan pengaruh yang baik bagi pengamatan lainnya. Penggunaan regresi robust untuk mengestimasi nilai-nilai parameter pada model regresi produksi Jagung Indonesia tahun 2010 diharapkan dapat memberikan informasi yang valid bagi pemerintah Indonesia sehingga dapat dijadikan acuan untuk langkah peningkatan produksi di tahun-tahun berikutnya.Kata kunci : produksi, jagung, estimasi, regresi, RobustI. PENDAHULUAN

Jagung merupakan bahan makanan pokok yang sangat dibutukan karena setiap hari dikonsumsi oleh sebagian masyarakat di seluruh dunia. Kebutuhan bahan makanan pokok akan meningkat sejalan dengan pertambahan jumlah penduduk yang semakin menigkat.

Menurut data BPS, produksi jagung tahun 2010 (ATAP) sebesar 18,33 juta ton, meningkat sebanyak 697,89 ribu ton (3,96 persen) dibandingkan tahun 2009. Peningkatan produksi tersebut terjadi di Jawa sebesar 489,94 ribu ton, dan di luar Jawa sebesar 207,95 ribu ton. Angka ramalan I (Aram I) produksi jagung tahun 2011 sebesar 17,93 juta ton. Jumlah ini turun sekitar 438.960 ton atau 2,39 persen ketimbang produksi tahun lalu. Sebenarnya, kebutuhan jagung nasional hanya 16,3 juta ton. Dengan produksi jagung sebesar 18,33 juta ton di tahun 2010, seharusnya kebutuhan dalam negeri tercukupi.

Analisis regresi merupakan teknik statistika yang digunakan untuk menyelidiki dan nemodelkan hubungan antara variabel dependen dan variabel independen. Jika Y variabel dependen dan X1, X2, ... , XK variabel independen, maka model regresi linear secara umum dapat dinyatakan sebagai Dengan adalah parameter-parameter regresi dan adalah sisaan yang berdistribusi normal dengan mean nol dan variansi konstan (Sembiring, 2003). Permasalahan yang muncul dalam analisis regresi adalah menentukan estimator terbaik untuk menentukan . Dalam menetukan estimator terbaik sangat dipengaruhi oleh penggunaan metode. Metode yang biasa digunakan adalah Metode Kuadrat Terkecil (MKT).

Dalam kasus model regresi linear, dimungkinkan terdapat data outlier (pencilan) yaitu pengamatan dengan nilai mutlak sisaan jauh lebih besar daripada sisaan-sisaan lain sehingga akan mempengaruhi model regresi yang terbentuk. Data pencilan tersebut tidak boleh dibuang begitu saja karena akan mempengaruhi model prediksi serta menghasilkan estimasi parameter yang kurang tepat. Untuk menyelesaikan masalah tersebut diperlukan adanya metode yang bersifat robust dimana nilai estimasinya tidak boleh dipengaruhi perubahan kecil dalam data.

Regresi Robust merupakan metode regresi yang digunakan ketika distribui dari sisaan tidak normalatau adanya beberapa pncilan yang berpengaruh pada model. Dalam regresi robust terdapat beberapa metode estimsi seperti estimasi-M, estimasi Least Median Square (LSM), estimasi Least Trimmed Squarre (LTS), estimasi-S, estimasi-MM (Chen, 2002)

Dalam penelitian ini penulis membahas dengan metode estimasi-S karena metode ini mempunyai kelebihan yaitu bisa digunakan untuk pencilan dengan proporsi hingga 50% serta digunakan ketik variabel dependen dan variabel independn terdapat pencilan.

Metode estimasi-S prtama kali dikembangkan oleh Rousseeuw dan Yohai (1984) dimana metode ini merupakan keluarga high breakdown point yaitu ukuran umum proporsi dari data pencilan yang dapat ditangani sebelum pengamatan tersebut mempengaruhi model prediksi. Disebut estimasi-S karena mengestimasi berdasarkan skala. Skala yang digunakan adalah standart deviasi sisaan. II. TUJUAN DAN MANFAAT PENELITIAN

a. Tujuan penelitian

Tujuan yang ingin dicapai dalam penelitian ini adalah menentukan estimasi produksi jagung di Indonesia tahun 2010 menggunakan metode regresi robust estimasi-S

b. Manfaat penelitian

Manfaat yang diharapkan dalam penelitian ini adalah dapat mengembangkan ilmu pengetahuan dalam bidang statistika dan industri. Pada bidang statistika, metode estimasi-S dapat diaplikasikan terhadap data yang mengandung pencilan pada variabel dependen dan independennya. Sedangkan pada bidang industri dapat memberikan masukan dalam meningkatkan produksi jagung di Indonesia.

III. METODOLOGIMetode penelitian yang digunakan dalam penulisan makalah ini adalah studi kasus, yaitu melakukan estimasi regresi robust pada model produksi jagung di 33 provinsi di Indonesia tahun 2010. Data yang digunakan adalah data sekunder yang diambil dari Badan Pusat Statistik (BPS) Indonesia.1. Pengujian Asumsi Analisis Regresi

Pada model regresi, perlu dilakukan uji asumsi analisis regresi untuk mengetahui apakah model memenuhi asumsi atau tidak. Uji asumsi yang dilakukan pada model regresi adalah

1. Normalitas

Analisis regresi linier mengasumsikan bahwa sisaan berdistribusi normal. Menurut Gujarati (1978) pada regresi linier klasik diasumsikan bahwa tiap didistribusikan secara random dengan Salah satu cara untuk menguji asumsi kenormalan adalah dengan uji Kolmogorov-Smirnov. Uji ini didasarkan pada nilai D dengan

Dengan adalah fungsi distribusi frekuensi kumulatifrelatif dari distribusi teoritis dibawah . adalah distribusi frekuensi kumulatif pengamatan sebanyak sampel. adalah sisaan berdistribusi normal. Selanjutnya nilai D ini dibandingkan dengan nilai D kritis dengan signifikansi (tabel Kolmogorov-Smirnov). Apabila nilai atau , maka asumsi kenormalan tidak dipenuhi.

Dalam penelitian ini asumsi yang digunakan adalah asumsi dari sisaan tidak berdistribusi normal, sehingga MKT tidak layak untuk digunakan (Draper dan Smith, 1998).

2. Homoskedastisitas

Salah satu asumsi penting dalam analisis regresi adalah variasi sisaan pada setiap variabel independen adalah homoskedastisitas. Asumsi ini dapat ditulis sebagai berikut

Salah satu cara menguji kesamaan variansi yaitu dengan melihat pola tebaran sisaan terhadap nilai estimasi Y. Jika tebaran sisaan bersifat acak (tidak membentuk pola tertentu), maka dikatakan bahwa variansi sisaan homogen (Draper dan Smith, 1998).

Untuk lebih tepatnya, menurut Gujarati (1978) salah satu cara untuk mendeteksi heteroskedastisitas adalah dengan pengujian korelasi rank Spearman yang didefinisikan sebagai berikut

Dengan perbedaan dalam rank yang ditempatkan pada dua karakteristik yang berbeda dari individual atau fenomena ke-i dan n adalah banyaknya individual yang dirank. Koefisien rank korelasi tersebut dapat digunakan untuk mendeteksi heterokedastisitas dengan mengasumsikan . Adapun tahapnnya dalah sebagai berikut

1. Mencocokkan regresi terhadap data mengenai Y dan X dan mendapatkan sisaan .

2. Dengan mengabaikan tanda dari , yaitu dengan mengambil nilai mutlaknya |, meranking baik harga mutlak || dan sesuai dengan urutan yang meningkat atau menurun dan menghitung koefisien rank korelasi Spearman yang telah diberikan sebelumnya.

3. Dengan mengasumsikan bahwa koefisien rank korelasi populasi adalah nol dan, signifikan dari yang disampel dapat diuji dengan pengujian t sebagai berikut :

Jika nilai t yang dihitung melebihi nilai t kritis maka H0 ditolak, artinya asumsi homoskedastitas tidak dipenuhi. Jika model regresi meliputi lebih dari satu variabel X, dapat dihitung antara dan tiap-tiap variabel X secara terpisah dan dapat di uji untuk tingkat penting secara statistik dengan pengujian t yang diberikan di atas.

3. Non autokorelasi

Salah satu asumsi penting dari regresi linear adalah bawa tidak ada autokrelasi antara serangkaian pegamatan yang diurutkan menurut waktu. Adanya kebebasan antar sisaan dapat dideteksi secara grafis dan empiris. Pendeteksian autokorelasi secara grafis yaitu denan melihat pola tebaran sisaan terhadap urutan waktu. Jika tebaran sisaan terhadap urutan waktu tidak membentuk suatu pola tertentu atau bersifat acak maka dapat disimpulkan tidak ada autokorelasi antar sisaan (Draper dan Smith, 1998)

Pengujian secara empiris dilakukan dengan menggunakan statistik uji Durbin-Watson. Hipotesis yang diuji adalah:

H0: Tidak terdapat autokorelasi antar sisaan

H1: Terdapat autokorelasi antar sisaan

Adapun rumusan matematis uji Durbin-Watson adalah:

Kaidah keputusan dalam uji Durbin-Watson adalah:

1. Jika atau , maka H0 ditolak berarti bahwa terdapat autokorelasi antar sisaan.

2. Jika , maka H0 tidak ditolak yang berarti bahwa asumsi non autokorelasi terpenuhi.

3. Jika atau maka tidak dapat diputudkan apakah H0 diterima atau ditolak, sehingga tidak dapat disimpulan ada atau tidak adanya autokorelasi.

4. Untuk statistik dari Durbin-Watson dapat dilihat pada tabel4. Non Multikolinearitas

Menurut Montgomery dan Peck (1992), kolinearitas terjadi karena terdapat korelasi yang cukup tinggi di antara variabel independen. VIF (Variance Inflation Factor) merupakan salah satu cara untuk mengukur besar kolineritas dan didefinisikan sebagai berikut

Dengan m= 1,2,...,p dan p adalah banyaknya variabel independen adalah koefisien determinasi yang dihasilkan dari regresi variabel independen dengan variabel independen lain . nilai VIF menjadi semakin besar jika terdapat korelasi yang semakin besar diantara variabel independen. Jika VIF lebih dari 10, multikolinearitas memberikan pengaruh yang serius pada pendugaan metode kuadrat kecil.2. Pencilan

Pada beberapa kasus dimungkinkan adanya data yang jauh dari pola kumpulan dan keseluruhan, yang lazim didefinisikan sebagai data pencilan. Keberadaan dari pencilan akan menyebabkan kesulitan dalam proses analisis data dan perlu untuk dihindari. Permasalahan yang uncul akibat adanya pencilan antara lain:

1. Sisaan yang besar dari model yan terbentuk 2. Variansi dari data akan menjadi lebih besar

3. Estimasi interval akan memiliki rentang yag lebih besar

Menurut Drape dan smith (1998) metode yang digunakan dalam mengidentifikasi pencilan terhadap variabel Y adalah Studientized Deleted Residual (TRES) yang didfinisikan sebagai:

Dimana: = 1,2,.....,n

= = = simpangan baku beda = =p+1

= banyaknya pengamatan

Hipotesis untuk menguji adanya pencilan adalah:

H0: pengamatan ke-i bukan pencilan

H1: pengamatan ke-i merupakan pencilan

adalah stastistik uji untuk mengetahui pencilan terhadap Kriteria pengujian yang melandasi keputusan adalah:

Metode yang diunakan dalam mengidentifikasi pencilan terhadap variabel adalah nlai pengaruh (Leverage Point). Nilai pengaruh dari penamatan menunjukan besarnya peranan terhadap dan didefinisikan sebagai:

=Dimana i: 1,2,...,n

Xi= adalah vektor baris yang berisi nilai-nilai dari peubah variabel independen dalam pengamatan ke-i. Nilai berada diantara 0 dan 1 dengan k=p+1. Jika lebih besar dari dengan

Maka pengamatan ke-i dikatakan pencilan terhadap X.

3. Estimasi-S

Estimasi-S pertama kali diperkenalkan oleh Rousseeuw dan Yohai (1984) merupakan estimasi robust yang dapat mencapai breakdown point hingga 50%. Breakdown point adalah ukuran umum proporsi dari pencilan yang dapat ditangani sebelum pengamatan tersebut mempengaruhi model. Karena estimasi-S dapat mencapai breakdown point hingga 50% maka estimasi-S dapat mengatasi setengah dari pencilan dan memberikan pengaruh yang baik bagi pengamatan lainnya.

Estimasi-S didefinisikan

dengan menentukan nilai estimator skala robust yang minimum dan memenuhi

dengan

merupakan fungsi pembobot Tukeys biweight

Penyelesaian persamaan (2.8) adalah dengan cara menurunkannya terhadap sehingga diperoleh

disebut fungsi pengaruh yang merupakan turunan dari . Sehingga bias dituliskan yaitu

Dengan merupakan fungsi pembobot IRLS dimana dan c = 1,547.

Sisaan awal yang digunakan pada estimasi-S adalah sisaan yang diperoleh dari metode kuadrat terkecil. Persamaan (2.9) dapat diselesaikan dengan MKT terboboti secara iterasi yang disebut Iteratively Reweighted Least Square (IRLS) hingga mencapai konvergen.

Tabel 2.1. Fungsi objektif dan fungsi pembobot untuk MKT dan Tukeys biweight

MetodeFungsi objektifFungsi pembobotInterval

MKT

Tukeys biweight

4. Langkah-langkah yang dilakukan dalam mengestimasiu parameter pada regresi robust estimasi-S adalah

a. Menduga koefisien regresi dengan MKT (Metode Kuadrat Terkecil)

b. Menguji asumsi klasik analisis regresi linear

c. Mendeteksi adanya pencilan pada data dengan metode TRES dan hiid. Menduga koefisien regresi dengan estimasi-S

Langkah-langkah metode estimasi-S :

a. Menghitung sisaan awal yang diperoleh dari MKT

b. Menghitung standar deviasi sisaan untuk mendapat nilai c. Menghitung nilai pembobot d. Menghitung MKT terbobot untuk mendapatkan penduga kuadrat terkecil terbobot

e. Menjadikan sisaan langkah (c) sebagai sisaan awal langkah (b) sehingga diperoleh nilai dan pembobot yang baru

f. Melakukan pengulangan iterasi sampai didapatkan kekonvergenan sehingga diperoleh yang merupakan estimasi-S

IV. HASIL DAN PEMBAHASANPada bab ini akan disajikan hasil analisis data sekunder produksi jagung di Indonesia tahun 2010 yang diperoleh dari BPS. Data tersebut meliputi produksi jagung sebagai variabel dependen (Y) sedangkan luas lahan dan produktivitas sebagai variabel independen (X). Data tersaji pada Tabel 4.1 berikut

NoProvinsiLuas Panen(Ha)Produktivitas(Ku/Ha)Produksi(Ton)

1Nanggroe Aceh Darussalam16709038.0743885

2Sumatera utara137771850.13274822

3Sumatera barat35426259.2459801

4Riau4186223.218044

5Kepulauan Riau96121.17454

6Jambi3069137.078280

7Sumatera Selatan12579637.2533769

8Kepulauan Bangka Belitung105530.94341

9Bengkulu7433136.2320516

10Lampung212657147.52447509

11DKI Jakarta3120.6715

12Jawa Barat92396260.08153778

13Banten2855732.848697

14Jawa Tengah305871048.41631816

15daerah Istimewa Yogyakarta34557639.886837

16Jawa Timur558731844.421257721

17Bali6635524.8526706

18Nusa Tenggara Barat24900540.4361593

19Nusa Tenggara Timur65362026.72244583

20Kalimantan Barat16827337.3845014

21Kalimantan Tengah934528.783247

22Kalimantan Selatan11644951.5622584

23Kalimantan Timur1199325.564693

24Sulawesi Utara44614436.59121930

25Gorontalo67916747.22143833

26Sulawesi Tengah16230637.9742747

27Sulawesi Selatan134304444.27303375

28Sulawesi Barat5802043.613308

29Sulawesi Tenggara7484025.2829607

30Maluku1527324.276293

31Maluku Utara205461910813

32Papua683417.513903

33Papua Barat193116.621162

4.1 Metode Kuadrat Terkecil

Model regresi ganda dengan metode kuadrat terkecil adalah

dengan

: Produksi jagung provinsi di Indonesia tahun 2010 (Ton)

: Produktivitas (Ku/Ha)

: Luas Panen (Ha)

Hasil di atas merupakan output dari Software Minitab 16

Selanjutnya dilakukan uji asuimsi klasik untuk melihat apakah model regresi yang diperoleh memenuhi asumsi klasik atau tidak. Berikut merupakan hasil uji asumsi klasik tersebut

4.1.1 Uji Normalitas

Pengujian kenormalan digunakan untuk mengetahui apakah sistem berdistribusi normal atau tidak. Plot kenormalan untuk sisaan dari model Produksi jagung Indonesia tahun 2010 sebagai berikut

Gambar 4.1 Plot probabilitas dari sistemGambar 4.1 terlihat bahwa pola penyebaran sisaan mengikuti garis lurus, ini

Berarti asumsi kenormalan tidak terpenuhi karena gambar plot terdapat pencilan. Untuk menguji kenormalan dapat juga digunakan uji Kolmogorof-Smirnof sebagai berikut

i. : sisaan berdistribusi normal: sisaan tidak berdistribusi normal

ii. Pilih iii. Daerah kritis: ditolak jika p-value < iv. Statistik ujiBerdasarkan output software Minitab 16, diperoleh hasil output pada gambar 4.1 dengan v. KesimpulanBerdasarkan hasil regresi dapat dilihat p-value < 0.05 maka ditolak artinya sisaan tidak berdistribusi normal.

Dengan demikian asumsi kenormalan pada data produksi jagung ke-33 provinsi di Indonesia tahun 2010 tidak dapat dipenuhi.

4.1.2 Uji Homoskedastisitas

Untuk mendeteksi homoskedaktisitas dapat dilakukan dengan metode plot. Plot kesamaan variansi untuk data sisaan pada model produksi jagung di indonesia tahun 2010 adalah sebagai berikut

Gambar 4.2. Plot sisaan dengan Pada gambar 4.2 tampak bahwa variansi sisaan dari satu pengamatan ke pengamatan lain berpola acak yang mengindikasikan nahwa variansi sisaan konstan sehingga dapat diindikasikan asumsi homoskedastisitas dipenuhi. Dari hasil tyersebut dapat diambil kesimpulan bahwa asumsi homoskesdastisitas dapat dipenuhi.

4.1.3 Uji Non Autokorelasi

Autokorelasi diartikan sebagai korelasi antara anggota serangkaian observasi yang diurutkan menurut waktu. Uji non autokorelasi dapat dideteksi dengan rumus Durbin-Watson.

Uji Durbin-Watson (Uji DW)

i. , artinya tidak ada autokorelasi, artinya ada autokorelasi

ii. Pilih iii. Daerah kritis

Pada dan serta diperoleh nilai dan sehingga dan ditolak jika iv. Statistik uji

Dari perhitungan dengan bantuan software Minitab 16 diperoleh v. KesimpulanBerdasarkan hasil regresi dapat diperoleh bahwa maka tidak ditolak artinya asumsi non autokorelasi pada model; produksi jagung Indonesia tahun 2010 dapat dipenuhi

4.1.4 Uji Non Multikolinearitas

Pengujian multikolinearitas bertujuan untuk mengetahui ada tidaknya hubungan linear antara variabel independen. Untuk mendeteksi adanya multikolinearitas dapat dilakukan dengan berbagai uji. Salah satu deteksi ada tidaknya multikolinearitas adalah dengan melihat pada nilai VIF. Nilai VIF diperoleh dengan melakukan regresi secara parsial dan kemudian menghitung nilai VIF. Dengan bantuan software Minitab 16, diperoleh hasil output sebagai berikut

Tabel 4.2. Hasil output uji multikolinearitasVariabel IndependenVIFKeterangan

(Produktivitas)1.173 < 10Tidak terdapat multikolinearitas

(Luas Panen)1.173 < 10Tidak terdapat multikolinearitas

Berdasarkan hasil output pada Tabel 4.2. dapat dilihat bahwa nilai VIF untuk semua variabel independen, baik variabel produktivitas () maupun Luas panen () adalah lebih kecil dari 10, sehingga dapat disimpulkan bahwa asumsi non multikolinearitas dipenuhi.

Berdasarkan pengujian asumsi klasik pada model produksi jagung di Indonesia tahun 2010 menggunakan analisis regresi diperoleh bahwa semua asumsi klasik terpenuhi.

4.2 Deteksi Pencilan

Berdasarkan statistik uji untuk mengetahui pencilan terhadap Y yaitu TRES dengan menarik kesimpulan menolak apabila nilai maka diperoleh kesimpulan bahwa terdapat pencilan data ke-14

Berdasar statistik uji untuk mengetahui pencilan terhadap X yaitu hii yang dengan menarik kesimpulan bahwa pengamatan menolak apabila nilai maka diperoleh kesimpulan bahwa terdapat pencilan data ke- 16.

Tabel 4.3. Hasil perhitungan TRES dan hii untuk mendeteksi pencilan

No

1-0.399270.037748

20.938610.077901

3-0.263270.193402

40.234210.064494

50.739280.076610

6-0.121760.040633

7-0.319880.037425

80.240860.039656

9-0.191640.037549

100.872070.092328

110.777800.079843

121.445330.180429

130.046400.037273

142.277750.163457

15-0.617390.037328

16-0.908130.716534

170.003930.056058

18-0.485160.041738

19-6.351260.069429

20-0.404610.036318

210.300590.043303

22-0.640950.118668

230.424200.053630

24-1.044480.030598

25-0.153560.063677

26-0.391380.037680

27-0.484150.054003

28-0.399490.063002

29-0.063630.054053

300.450100.059276

310.569990.091749

320.839330.103368

330.972530.110840

4.3 Model Regresi Robust dengan Estimasi- S

Proses perhitungan estimasi-S yang iteratif dimulai dengan menentukan estimasi awal koefisien regresi, yang diperoleh dari MKT yaitu kemudian berdasarkan algoritma estimasi-S, dihitung nilai dan sisa . Proses iterasi menggunakan MKT terboboti dilanjutkan dengan menghitung sisaan dan pembobot yang baru dan dilakukan pendugaan parameter secara berulang-ulang sampai konvergen. Kekonvergenan tercapai jika koefisien regresi sudah sama dengan koefisien regresi sebelumnya. (Salibian dan Yohai,2006).Tabel 4.4. Nilai dan dan pada estimasi-S

No

1533475326582.00.011680.9999

2943910535-1095.6-0.156060.9797

318344799-2965.2-0.422370.8565

458306206-375.7-0.053510.9976

56-19291935.40.275680.9375

653205966-645.9-0.092000.9929

71166412879-1215.4-0.173120.9751

852-20172068.90.294690.9287

927191924795.00.113240.9893

1073258087-761.9-0.108530.9902

110-99609959.61.418650.0253

1255823529802843.00.404960.8676

131166212786-1124.0-0.160100.9787

1418799216096427028.03.8498926.9696

153824445591-7347.1-1.046530.2942

16339491342405-2914.2-0.415100.8612

1755545947-393.0-0.055980.9974

1893122118457-25335.5-3.6088019.7300

1917805971183.00.168500.9764

2034774527-1050.3-0.149610.9814

2127642604160.50.022860.9996

22380927091100.20.156720.9796

2322042894-690.4-0.098340.9919

2476278636-1008.7-0.143690.9828

2534033495-91.8-0.013070.9999

2635554128-573.0-0.081620.9944

273571134846865.50.123280.9873

2831955226-2030.8-0.289260.9313

2932033378-174.6-0.024870.9995

3011831009174.50.024850.9995

31944747197.30.028110.9993

324152409161.10.008700.9999

33600-7391339.20.190750.9698

Berikut ini merupakan nilai dan pada iterasi demi iterasi yang didapat dari software Minitab 16

Tabel 4.5. Tabel nilai dan Iterasi

1.93766.4(- 90002 ; 2406 ; 4.45)

2.97868.1(-76027 ; 2021 ; 4.44)

3.99671.0(-58936 ; 1433 ; 4.44)

4.102097(- 54142 ; 1260 ; 4.44)

5.102919(- 52627 ; 1204 ; 4.44)

6.103189(- 52136 ; 1185 ; 4.44)

7.103277(- 51976 ; 1179 ; 4.44)

8.103305(- 51925 ; 1178 ; 4.44)

9.103314(- 51908 ; 1177 ; 4.44)

10.103317(- 51903 ; 1177 ; 4.44)

11.103318(- 51901 ; 1177 ; 4.44)

12.103319(- 51901 ; 1177 ; 4.44)

Berdasarkan Tabel 4.5. terlihat bahwa koefisien regresi sudah konvergen di iterasi ke-12 dengan model

Interpretasi model yaitu sebesar 99.9% produksi jagung di Indonesia pada tahun 2010 dapat diterangkan oleh variabel luas panen dan produktivitas tiap provinsi di Indonesia, sedangkan 0.1% diterangkan oleh variabel lain. Setiap peningkatan satu hektar (ha) luas panen maka produksi jagung di Indonesia akan betambah sebesar 4.44 ton, setiap peningkatan produktivitas sebesar satu kuintal/hektar(1 ku/ha) maka produksi jagung di Indonesia akan bertambah sebesar 1177 ton.

Untuk mengetahui variabel independen yang berpengaruh dilakukan uji signifikansi model regresi robust estimasi-S

i. (Luas panen dan produktivitas tidak berpengaruh signifikan terhadap produksi jagung di Indonesia tahun 2010)

untuk suatu (paling tidak terdapat salah satu luas panen atau produktivitas yang berpengaruh secara signifikan terhadap produksi jagung di Indonesia tahun 2010.

ii. Pilih iii. Daerah kritis : ditolak jika

iv. Statistik uji

Berdasarkan output Minitab 16 diperoleh nilai v. Kesimpulan

Karena maka ditolak artinya paling tidak ada salah satu diantara luas panen atau produktivitas yang berpengaruh secara signifikan terhadap produksi jagung di Indonesia tahun 2010.

Selanjutnya dilakukan uji parsial untuk mengetahui signifikansi atau pengaruh masing-masing variabel terhadap model regresi yang dihasilkan.

Tabel 4.6. Hasil uji pada estimasi S

Variabel IndependenP-ValueKesimpulan

Luas Panen0.00 < 0.05Signifikan

Produktivitas0.014 < 0.05Signifikan

Berdasarkan Tabel 4.5. dapat disimpulkan bahwa luas panen dan produktivitas ke-33 provinsi di Indonesia berpengaruh secara signifikan terhadap produksi jagung di Indonesia tahun 2010.V. KESIMPULAN DAN SARAN1. Kesimpulan

Berdasarkan hasil analisis, maka dapat disimpulkan bahwa

1. Hasil estimasi produksi jagung di Indonesia pada tahun 2010 dengan metode regresi robust estimasi-S diperoleh sebagai berikut

Interpretasi model yaitu sebesar 99.9% produksi jagung di Indonesia pada tahun 2010 dapat diterangkan oleh variabel luas panen dan produktivitas tiap provinsi di Indonesia, sedangkan 0.1% diterangkan oleh variabel lain. Setiap peningkatan sati hektar (ha) luas panen maka produksi jagung di Indonesia akan betambah sebesar 4.44 ton, setiap peningkatan produktivitas sebesar satu kuintal/hektar(1 ku/ha) maka produksi jagung di Indonesia akan bertambah sebesar 1177 ton.

2. Variabel independen yang berpengaruh dalam estimasi produksi jagung di Indonesia pada tahun 2010 dengan metode robust estimasi-S adalah luas panen dan produktivitas.2. Saran

Bagi para pembaca untuk menganalisis data yang lebih valid maka dapat digunakan metode lain selain estimasi-S yaitu estimasi-MM, estimasi LTS dan estimasi LMS untuk menyelesaikan masalah yang ada.VI. DAFTAR PUSTAKAhttp:// www.bps.go.id. Tanggal akses 17 Oktober 2012Artiana, Griya. (2012). Skripsi Estimasi Parameter Regresi Robust Dengan Metode Estimasi-S pada Penjualan Energi Listrik di Jawa Tengah Tahun 2009.(2.8)

(2.9)

EMBED MtbGraph.Document.16

EMBED MtbGraph.Document.16

_1412959585

_1412959587