MAKALAH ELASTISITAS PERMINTAAN DAN PENAWARAN BERAS DI LUAR JAWA

ABSTRAK
Tujuan utama penelitian ini adalah menganalis penawaran dan permintaan beras luar Jawa dan prospek kawasan ini dalam mendukung swasembada beras dengan menggunakan model persamaan simultan. Model ini terdiri atas sub model; produksi, konsumsi dan perdagangan. Hasil Penelitian Menunjukkan:
Pertama, Perilaku areal panen padi di luar Jawa ternnyata hanya dipengaruhi oleh harga padi. Walaupun demikian elastisitas areal panen terhadap harga padi adalah inelastis. Fenomena ini menunjukkan harga padi akan mendorong petani meningkatkan produksi padi melalui peningkatan areal (ekstensifikasi), bukan melalui peningkatan produktivitas (intensifikasi) karena harga padi tidak siknifikan pengaruhnya terhadap produktivitas padi.
Kedua, Produksi beras luar Jawa tidak siknifikan pengaruhnya terhadap impor beras menunjukkan produksi beras di luar Jawa belum mampu menjadi kontributor yang siknifikan dalam mengurangi impor beras nasional.
Hal ini menunjukkan luar Jawa belum berperan besar sebagai pensuplai beras nasional.
Ketiga, Permintaan beras di luar Jawa tidak dipengaruhi oleh harga beras
tetapi sangat ditentukan oleh jumlah penduduknya. Hal ini menunjukkan permintaan beras luar Jawa di masa mendatang akan semakin meningkat sejalan
dengan semakin bertambahnya jumlah penduduk luar Jawa.
Keempat, Harga padi di luar Jawa sangat ditentukan oleh harga dasar namun respon (elastisitas) harga padi terhadap harga dasar adalah inelastis (kurang dari satu), dan kelima, harga beras eceran luar Jawa dipengaruhi oleh
harga dasar dan harga padi dengan nilai elastisitas harga beras eceran terhadap
harga dasar dan harga padi itu adalah inelastis.
Hasil penelitian ini menyarankan kebijakan yang dapat menolong harga padi di tingkat petani seperti harga dasar dan subsidi input adalah penting untuk tetap diterapkan untuk memacu produksi beras di luar Jawa. Hal ini karena harga dasar ini akan mempengaruhi harga padi dan selanjutnya harga padi akan memacu petani untuk meningkatkan produksi melalui ekstensifikasi

Pendahuluan
Pasca penghargaan swasembada, ada kesan di kalangan pengambil kebijakan masalah pangan khususnya beras dianggap telah tuntas. Pemerintah terlena dengan penghargaan FAO atas keberhasilan mewujudkan swasembada pangan di tahun 1984, setidaknya hal ini diindikasikan dengan semakin menyusutnya lahan-lahan sawah subur di Pulau Jawa sejak tahun 1984 tersebut untuk berbagai kepentingan industri dan perumahan.
Walaupun selanjutnya ada Keppres No 32 tahun 1992 tentang larangan
pengalihan fungsi lahan irigasi teknis di Pulau Jawa, tetapi gagal mencegah proses
konversi lahan-lahan irigasi di Jawa. Akibatnya produksi beras nasional turun drastis, terbukti di tahun 1989 kita telah mengimpor beras sebesar 464.449 ton bahkan 10 tahun kemudian kita dikejutkan dengan jumlah impor yang sangat spektakuler yakni sebesar 5,8 juta ton di tahun 1998. Penyusutan lahan persawahan di Jawa disebabkan oleh desakan pertambahan penduduk, perkembangan sektor industri, konversi lahan produktif menjadi real estate, daerah wisata dan peruntukan lainnya yang saling tumpang tindih (Kasryno,
1996). Hal ini dapat dilihat misalnya dari laju konversi lahan pertanian (sawah) yang cepat. Menurut Hermanto, dalam dekade terakhir rata-rata konversi lahan sawah di Jawa berkisar 13.400 sampai 87.600 hektar per tahun (Irawan, 1997b). Pada masa mendatang trend konversi sawah di Jawa diperkirakan masih akan terjadi sehingga beban wilayah ini sebagai penghasil beras nasional akan semakin berat. Tabel 1. Rata-rata Luas Lahan Sawah yang Dikonversikan di Pulau Jawa Referensi Jenis lahan Ha/tahun Keterangan
JICA, 1988
Deli Hydarulic, 1991
World Bank, 1988
Sumaryanto et.al, 1995
BCEOM, 1988
Nasoetion dan Winoto
Sawah irigasi 20.000 –
Sawah irigasi 22.500 1990-2000
Sawah 20.000 –
Sawah 22.637 1981-1993
Sawah irigasi 13.400 1981-1985
Sawah irigasi 27.633 1981-1986

Sumber: Irawan (1997)
Selama ini kecenderungan konversi lahan yang tinggi pada lahan-lahan pertanian disekitar sentral pertumbuhan ekonomi dan industri yang umumnya adalah kota-kota besar di Jawa seperti wilayah Jabotabek, Bandung, Semarang, Yogyakarta, Surabaya dan Malang (Dahuri dan Saefuddin, 1996). Akibat problem konversi lahan di pulau Jawa, maka potensi produksi gabah hilang sekitar 7,5 ton per tahun. Jika konversi lahan dengan laju yang begitu cepat sampai tahun 2020, maka potensi kehilangan gabah di Jawa sekitar 82 juta ton per tahun, setara dengan pemenuhan kebutuhan beras bagi seratus juta penduduk pulau Jawa tahun 2020, sementara hingga saat ini 63 persen suplai beras nasional masih bersumber dari pulau Jawa (Irawan, 1998).
Penyusutan lahan di Jawa seperti dikemukakan di atas tak pelak lagi menjadi kontribusi utama turunnya produksi beras nasional. Di tahun 1998, pemerintah terpaksa mengimpor beras dalam jumlah yang sangat mengejutkan yakni sebesar 5,8 juta ton.
Sementara itu sentra produksi beras di luar Pulau Jawa sampai sekarang belum juga mampu menyamai prestasi lahan dan petani di Pulau Jawa, dimana produktivitas padi sawah dan ladang tahun 1996-2000 di Jawa adalah rata-rata 50,14 kuintal per hektar lebih tinggi 43% dibanding produktivitas luar Jawa yang rata-rata hanya sebesar 35,05 kuintal per hektar (Irawan, 2000).
Lahan sawah di Pulau Jawa sendiri tampaknya akan terus bergeser menjadi lahan untuk industri dan jasa sehingga dalam era selanjutnya kawasan luar Jawa diharapkan akan menjadi tumpuan harapan untuk berperan besar. Untuk mengetahui apakah memang kawasan luar Jawa mampu menjadi subtitusi penting sebagai pensuplai beras nasional dikaitkan dengan kondisi permintaan (konsumsi) beras dan sejumlah kebijakan perberasan di Indonesia, maka sangat diperlukan informasi tentang perilaku penawaran dan permintaan beras luar Jawa.
Informasi dari keragaan pasar padi/beras luar Jawa ini diharapkan dapat
memberikan suatu saran kebijakan sehubungan dengan kawasan luar Jawa di masa mendatang sebagai pendukung ketersedian pangan (beras) nasional.

KERANGKA PEMIKIRAN
Fungsi Produksi dan Penawaran Beras
Fungsi produksi padi dapat dirumuskan sebagai berikut:
Qs = q (A, F, L V) …………………(2.1)
Dimana: Qs = Jumlah produksi padi (unit)
A = Luas Areal padi (unit)
F = Jumlah pemakaian pupuk (unit)
L = Jumlah tenaga kerja (unit)
V = Faktor produksi lainnya (unit)
Untuk memaksimumkan produksi padi dibutuhkan biaya tertentu. Perumusan biaya dalam bentuk anggaran total adalah sebagai berikut:
B = Bo + Pa A + Pf F + Pl*L + Pv*V ………… (2.2)
Dimana:
B = Biaya total (Rp)
Bo = Biaya Peubah (Rp)
Pa = Harga lahan (Rp/unit)
Pf = Harga pupuk (Rp/unit)
Pl = Harga tenaga kerja (Rp/unit)
Pv = Harga faktor produksi lainnya (Rp/unit)
Sehingga fungsi keuntungan produksi padi dapat dirumuskan sebagai berikut:
= Ps * Qs – B …………………………………………(2.3)
= Ps * q(A, F, L, V) – (Bo + Pa * A + Pf * F + Pl *l + Pv* V) … .(2.4)
dimana:
= Keuntungan (Rp)
Ps = Harga padi (Rp/unit)
Fungsi keuntungan diperoleh jika turunan pertama sama dengan nol dan turunan
kedua mempunyai nilai Hessian Determinan lebih besar nol.
Turunan pertamanya adalah:
/A = Ps * A’ – Pa = 0 atau Ps * A’= Pa ………(2.5)
/F = Ps * F’ –Pf = 0 atau Ps * F’= Pf ………(2.6)
/L = Ps * L’ – Pl = 0 atau Ps * L’= Pl ………..(2.7)
/V = Ps * V’ – Pv = 0 atau Pp * V’= pv ……….(2.8)
Dimana A’ F’ L’ dan V’ adalah produk marginal masing-masing produksi. Oleh
sebab itu keuntungan maksimum diperoleh jika produk marginal sama dengan rasio harga faktor terhadap harga produk. Dari persamaan di atas diketahui bahwa seluruh peubah harga merupakan peubah eksogen sedangkan lainnya adalah endogen.
Fungsi permintaan faktor produksi oleh petani dapat dirumuskan sebagai berikut:
A = a(Ps, Pa, Pf, Pl, Pv) ……………………………………..(2.9)
F = f (Ps, Pa, Pf, Pl, Pv) ……………………………………..(2.10)
L = l(Ps, Pa, Pf, Pl, Pv) ……………………………………..(2.11)
V = v (Ps, Pa, Pf, Pl, Pv) ……………………………………..(2.12)
Peningkatan atau penurunan harga padi akan meningkatkan atau menurunkan jumlah produksi pada dan meningkatkan atau menurunkan jumlah permintaan faktor terhadap faktor produksi. Dengan mensubtitusikan persamaan (2.9), (2.10), (2.11),
(2.12), ke persamaan (2.1) maka fungsi penawaran dapat dirumuskan sebagai berikut:
Qs = q(Ps, Pa, Pf, Pl, Pv) …………………….(2.13)
Dengan mengembangkan fungsi penawaran di atas, Dolan menunjukkan selain harga komoditi dan harga input produksinya, peubah penting lainnya yang
mempengaruhi penawaran suatu komoditi adalah harga komoditi lainnya, biaya faktor produksi, tujuan perusahaan, tingkat teknologi, pupuk, subsidi, harapan harga dan keadaan alam (Hutauruk, 1996).

Fungsi Permintaan
Fungsi permintaan beras diturunkan dari fungsi utilitas konsumen. Fungsi utilitas dapat dirumuskan sebagai berikut:
U = u(Qd, R) ……………………………………(2.14)
Dimana:
U = Total utilitas dari beras (unit)
Qd = Jumlah beras yang dikonsumsi (unit)
R = Jumlah komoditi lain yang dikonsumsi (unit)
Konsumen yang rasional akan memaksimumkan kepuasannya dari konsumsi suatu komoditi pada tingkat harga yang berlaku dan tingkat pendapatan tertentu.
Dengan demikian sebagai kendala untuk memaksimumkan fungsi utilitas adalah sebagai berikut:
Y = Pb * Qd + Pr * R ……………………………. (2.15)
Dimana:
Y = Tingkat Pendapatan (Rp)
Pb = Harga beras (Rp/unit)
Pr = Harga komoditi lain (Rp/unit)
Dari persamaan (2.14 dan 2.15) dapat dirumuskan fungsi kepuasan yang akan
dimaksimumkan yakni sebagai berikut:
Z = u(Qd, R) + (Y – Pb * Q – Pr * R) …………………(2.16)
Dimana adalah Lagrange Multiplier, jika syarat pertama dan kedua dipenuhi maka fungsi utilitas dapat dimaksimukan sebagai berikut:
Qd’ – Pb = 0 atau Qd’ = Pb …………………….(2.17)
R’ – Pr = 0 atau R’ =Pr ………………………(2.18)
Y – Pb * Qd – Pr * R = 0 ……………………………………(2.19)
Dimana Qd’ dan R’ adalah utilitas marginal dari komoditi Q dan R sehingga
= Q’/Pq = R’/Pr ……………………………….(2.20)
Persamaan (2.20) menunjukkan bahwa kepuasan maksimum konsumen tercapai jika utilitas marginal dibagi dengan harga harus sama bagi kedua komoditi tersebut dan juga harus sama dengan utilitas marginal dari pendapatan.
Dari persamaan (2.19) dan (2.20) diketahui bahwa Pb, Pr dan Y merupakan peubah eksogen yang mempengaruhi permintaan beras.
Dengan demikian fungsi permintaan beras dapat dirumuskan sebagai berikut:
Qd = d (Pb, Pr, dan Y) …………………………………(2.21)
Persamaan 2.21 menunjukkan bahwa jumlah permintaan beras merupakan fungsi dari harga beras (Pb), harga komoditi lain (Pr) dan Pendapatan (Y). Selanjutnya Koutsoyiannis (1982) memberi perluasan teori permintaan ini yang menyatakan bahwa permintaan suatu barang dipengaruhi oleh harga barang tersebut, harga barang lain, selera, pendapatan, distribusi pendapatan, jumlah penduduk dan harapan harga.
Respon Beda Kala Komoditas Pertanian
Salah satu karakteristik utama produk pertanian adalah tenggang waktu antara menanam dengan memanen (gestation period). Dengan demikian hasil yang diperoleh petani didasarkan pada prakiraan-prakiraan di masa mendatang serta pengalaman di masa lalu.
Pada kebanyakan komoditas pertanian, harga output tidak dapat dipastikan saat produk itu ditanam. Petani harus mengambil keputusan-keputusan produksi berdasarkan pengalaman masa lampau. Perkiraan petani itu dilakukan berdasarkan pengalaman di masa lampau.
Hal ini mengacu pada adanya beda kala (lag) antara dua periode, yaitu saat
menanam dan panen. Respon petani terjadi setelah beda kala sebagai dampak perubahan harga input, output dan kebijakan pemerintah. Jika terjadi peningkatan harga beras (gabah) maka tidak segera diikuti oleh peningkatan produktivitas dan areal karena keputusan alokasi sumberdaya telah ditetapkan pada saat sebelumnya, hal ini menurut Tomek dan Robinson (1990) disebut asset fixity (kekakuan aset). Jika petani memperkirakan bahwa peningkatan harga ini akan berlangsung terus pada periode berikutnya barulah petani mengubah komposisi inputnya pada masa tanam mendatang, sehingga pengaruh kenaikan harga tersebut baru terlihat pada periode tanam berikutnya. Jika praduga adanya ekspetasi demikian dapat diterima maka hubungan-hubungan yang spesifik antara harga harapan dengan harga masa lalu dapat dibuat (Colman dan Young, 1990), sehingga model dapat dikembangkan menjadi dinamik yang dirintis antara lain oleh Nerlove melalui model penyesuaian parsial.
Nerlove (1958) berpendapat para petani di setiap periode merevisi dugaan mereka terhadap apa yang mereka anggap sebagai proporsi yang normal terhadap perbedaan yang terjadi dengan yang sebelumnya dianggap normal. Dengan kata lain, petani menyesuaikan perkiraan harga di masa mendatang dalam bentuk proporsi dari selisih antara perkiraan dengan kenyataan. Petani tidak akan dapat segera menyesuaikan kegiatan produksi mereka sebagai respon setelah adanya stimulus pasar.
Hal ini terjadi karena alasan-alasan sebagai berikut:
1. Secara psikologis, adanya hambatan untuk segera melakukan perubahan karena telah terbiasa (habit) dengan perilaku lama. Di sini muncul faktor kelembaman (inersia) dalam menyesuaikan diri terhadap perubahan. Misalnya saja orang tidak akan segera meningkatkan konsumsinya begitu terjadi kenaikan pendapatan mereka karena penyesuaian untuk meningkatkan konsumsi ini memerlukan waktu. Dalam bidang pertanian kelembaman ini juga terjadi, misalnya jika terdapat perubahan yang melibatkan adopsi teknologi baru yang secara tradisional tidak diusahakan.
2. Perlunya penyesuaian parsial karena adanya kendala kelembagaan seperti kuota produksi dan ketersediaan sarana pendukung berupa kredit usahatani. Jika terjadi perubahan harga faktor produksi maka petani memerlukan waktu untuk melakukan subtitusi input dan hal ini membutuhkan tenggang waktu (gestation period).
3. Adanya kendala kelembagaan seperti adanya kontrak/perjanjian, maka selama masa kontrak pihak yang terlibat harus mentaatinya, dalam hal pertanian alokasi sumberdaya baru dapat dilakukan setelah kontrak selesai.

Model Distribusi Beda Kala
Nerlove (1958) berpendapat tidak mudah menghitung elastisitas penawaran jangka pendek karena sebenarnya merupakan elastisitas titik sehingga nilainya berubah pada titik yang berbeda. Sedang elastisitas jangka panjang sukar dihitung secara langsung. Ia menawarkan cara baru dengan model distribusi beda kala penyesuaian parsial yang kemudian terbukti dapat menjelaskan fenomena penawaran dengan lebih baik.
Elastisitas jangka pendek selalu lebih kecil dari elastisitas jangka panjang karena dalam jangka panjang memungkinkan bergesernya fungsi penawaran dan penyesuaian sumberdaya. Ia berpendapat masalah formulasi hubungan-hubungan ekonomi yang memasukkan distribusi beda kala sesungguhnya terletak pada estimasi elastisitas jangka panjang.
Estimasi distribusi beda kala dapat dilakukan dalam tiga kelompok pendekatan (Nerlove (1958):
1. Metode penaksiran khusus. Merupakan suatu pendekatan dengan tidak mengemukakan suatu asumsi yang spesifik tentang bentuk distribusinya dengan cara menghentikan proses model regresi jika koefisien pada peubah beda kala yang dimasukkan mulai tidak nyata secara statistik atau koefisien regresinya berubah tanda saat peubah bebasnya ditambah.
2. Dengan mengasumsikan suatu bentuk umum distribusi beda kala dan mengestimasi parameter yang menegaskan distribusinya dengan tepat. Dalam kelompok ini ada cara dari Koyck, Friedman, Fisher dan Cagan.
3. Metode penyesuaian parsial dengan mengembangkan suatu model dinamik yang eksplisit mengenai perilaku produsen, yang antara lain dikembangkan oleh (Nerlove, 1958). Pendekatan pertama dikritik karena memiliki beberapa kelemahan seperti tidak adanya petunjuk apriori mengenai panjang maksimum beda kala, masalah pada derajat bebas; jika peubahnya banyak menyebabkan derajat bebasnya kecil serta ada kecenderungan multikolinearity yang tinggi.

Model Penyesuaian Parsial
Model distribusi beda kala penyesuaian parsial yang dikembangkan Nerlove merupakan model yang populer digunakan dalam studi-studi respon penawaran.
Dalam bentuk yang paling sederhana yaitu model penyesuaian parsial sederhana berderajat satu, misalnya dalam konteks respon areal padi terhadap harga. Areal panen padi yang diinginkan (A*t) dipengaruhi oleh tingkat harga komoditas, maka persamaannya menjadi:
A*t = bo + b1 X t + ut ………………………………………… (2.22)
dimana A*t = areal panen yang diinginkan
X t = harga beras
Luas areal yang diharapkan tidak dapat diamati secara langsung sehingga untuk mengatasinya didalilkan suatu hipotesis yang merupakan hipotesis perilaku penyesuaian parsial:
A t – At-1 = d(A*t – A t-1 ) ………………………………………………..(2.23)
Perubahan areal yang sebenarnya terjadi (A t – At-1 ) merupakan proporsi tertentu dari perubahan yang diinginkan d(A*t – A t-1 ).
Proporsi tertentu ini disebut koefisien penyesuaian parsial. Nilai d ini terletak di antara dua nilai ekstrim 0 dan 1.
jika d = 0, maka tidak ada perubahan apapun dalam areal d = 1, maka areal yang diharapkan sama dengan yang dicapai sehingga penyesuaiannya seketika.
Persamaan (2.23) dapat diatur kembali sehingga dituliskan:
A t = d A*t + (1-d) A t-1 ……………………………………..(2.24)
Areal panen padi yang diamati pada periode tertentu dipengaruhi oleh luas areal yang diinginkan dan luas areal yang ada pada permulaan periode sebelumnya. Jika (2.22) disubtitusikan ke (2.24) dan menyusunnya kembali maka diperoleh:
Y t = d (bo + b1 X t + u t ) + (1-d) Y t-1
= d bo + db1 X t + (1-d) Y t-1+ du t………………………………..(2.25)

Model Respon Penawaran Padi
Model empiris penawaran padi yang digunakan dalam penelitian ini pada
dasarnya menggunakan model penyesuaian parsial Nerlove. Untuk memperoleh dugaan respon penawaran maka dilakukan pendugaan tak langsung. Ghatak dan Ingersen dalam Irawan (1999) mengatakan bahwa hubungan antara luas areal panen, produktivitas dan out put, dalam bentuk yang sederhana adalah output (Q) dispesifikasikan sebagai perkalian antara luas areal panen (A), produktivitas (Y) dan peubah teknis dan ekonomi lainnya (Z), sehingga dapat dituliskan sebagai berikut:
Q = A.Y ……………………………………………… (2.26)
Diasumsikan areal panen (A) dan produktivitas (Y) merupakan fungsi dari harga, sehingga respon terhadap perubahan harga (P) adalah sebagai berikut:
A = A(P,Z…)
Y = Y(P, A, Z, …)
Maka jika (2.26) didiferensiasikan total terhadap harga adalah sebagai berikut:
dQ/dP = Y (dA/dP) + A (dY/dP) ………………………………(2.27)
Karena produktivitas respon terhadap perubahan areal, maka pengaruh perubahan areal terhadap produktivitas akibat adanya perubahan harga dapat ditulis sebagai berikut:
dY/dP = Y/P (Y/ A . A / P) ……………………………..(2.28)
Artinya perubahan produktivitas karena perubahan harga terdiri atas perubahan produktivitas secara parsial terhadap harga padi dan perubahan produktivitas karena terjadinya perubahan areal akibat perubahan harga padi. Jika persamaan (2.28) disubtitusikan ke (2.26) maka:
dQ/dP = Y (dA/dP) + A [(Y/P + (Y/ A . A / P)] ………………….(2.29)
kedua ruas dikalikan P/Q maka diperoleh:
(dQ/dP) P/Q = (P/Q )Y (dA/dP) + A(P/Q) [(Y/P + (Y/ A . A / P)] ……..(2.30)
(dQ/dP) P/Q = A/A (P/Q )Y (dA/dP) + A(P/Q) [(Y/P + (Y/ A . A / P)] ……(2.31)
karena A.Y = Q maka:
(dQ/dP) P/Q = P/A (dA/dP) + P/Y (Y/P) + (Y/ A (A/Y).((P/A) (A / P)) ….(2.32)
Jika dinyatakan dalam elastisitas maka:
E(P) = E(Y,P) +[ E (A,P) (1 + E(Y,A))] ……………………..(2.33)
dimana : E(P) = Elastisitas (respon) penawaran padi
E(Y,P) = Elastisitas (respon) produktivitas terhadap harga
E (A,P) = Elastisitas (respon) Areal terhadap harga
E(Y,A) = Elastisitas (respon) produktivitas terhadap areal
Dengan demikian maka elastisitas (respon) penawaran agregat E(P) dapat diduga secara tidak langsung dengan melakukan pendugaan lebih dahulu terhadap E(Y,P), E (A,P), dan E(Y,A).

KONTRUKSI MODEL DAN PROSEDUR ANALISIS
Konstruksi Model
Secara teoritis petani yang bertujuan memaksimumkan keuntungan akan
memutuskan berapa output yang akan diproduksi untuk memaksimumkan pendapatan mereka.
Produksi beras luar Jawa yang merupakan komponen dari penawaran beras, tergantung pada luas panen dan tingkat produktivitas padi, jadi produksi merupakan hasil perkalian dari kedua peubah tersebut. Penawaran beras luar Jawa pada periode tertentu berasal dari tiga sumber yaitu jumlah produksi domestik, stok yang dilepas ke pasar, dan jumlah impor, sedangkan jumlah beras yang tersedia ini akan digunakan untuk konsumsi pangan. Selanjutnya, hubungan keterkaitan dari masing-masing peubah dalam ekonomi beras luar Jawa dapat dilihat pada model ekonometrika berikut.
Model dalam studi ini dirumuskan melalui persamaan-persamaan yang terdiri dari peubah-peubah penjelas (explanatory variables) yang digunakan dalam model ekonomi beras luar Jawa ini. Terdapat 10 peubah endogen yang menjelaskan perilaku luas areal panen, produktivitas padi, produksi padi dan beras, impor beras nasional, permintaan beras luar Jawa, perubahan stok luar Jawa, harga padi dan harga beras.
Berdasarkan kriteria ekonomi, statistik dan ekonometrika, maka model struktural yang disajikan dalam studi ini adalah model yang dinyatakan yang paling memuaskan.
Model sistem persamaan simultan tersebut adalah sebagai berikut:

a. Persamaan Luas Areal Panen
Luas areal panen padi merupakan fungsi dari harga padi, harga dasar rill, harga tanaman yang berkompetisi dengan padi dalam menggunakan lahan (dalam hal ini adalah jagung), tingkat bunga, curah hujan dan luas areal panen tahun sebelumnya.
Dengan demikian luas areal padi di Luar Jawa adalah sebagai berikut:
ALJW = a0 + a1 HPDT + a2 HDPTRIL + a3 HKJT + a4 SUKT + a5 CHLT + a6 ALJW(_1) + u1 …………………………..(3.1)
Dimana:
ALJW = Luas areal panen padi (sawah+ladang) luar Jawa (000 ha)
HPDT = Harga gabah (Rp/kg)
HDPTRIL = Harga dasar gabah riil (Rp/kg)
HKJT = Harga jagung (Rp/kg)
SUKT = Suku bunga pinjaman (%)
CHLT = Curah hujan (mm/thn)
ALJW(_1) = areal panen satu tahun sebelumnya
u1 = peubah pengganggu
Tanda parameter yang diharapkan:
a1, a2, a5 > 0, a2, a4 >0 dan 0 < a6 < 1 b. Persamaan Produktivitas Produktivitas padi di luar Jawa adalah merupakan fungsi dari Harga Padi, persentase Lahan Irigasi Jawa, Areal Panen luar Jawa dan Produktivitas padi satu tahun sebelumnya. Persamaan produktivitas padi di luar Jawa dapat dirumuskan sebagai berikut: YLPT = b0 + b1 HPDT + b2 ILT + b3 ALJW + b4 YLPT(_1)+ u2 ……..(3.2) Dimana: YLPT = Produktivitas (padi sawah + lading) di luar Jawa HPDT = Harga gabah (Rp/kg) ILT = luas lahan irigasi di Jawa (000 ha) ALJW = Luas areal panen padi (sawah+ladang) luar Jawa (000 ha) YLPT(_1) = produktivitas padi satu tahun sebelumnya u2 = peubah pengganggu Tanda parameter dugaan yang diharapkan: b1, b2, b5 > 0, dan 0 < b4 < 1
Jumlah produksi padi di luar Jawa selama tahun t merupakan perkalian dari Luas Areal dan Produktivitas, ditunjukkan pada persamaan di bawah ini:
QPLT = ALJW * YLPT ……………………(3.3)
QPLT = Produksi beras Luar Jawa
Selanjutnya total produksi beras di luar Jawa didalam penelitian ini diperoleh dengan mengalikan suatu faktor konversi dengan total produksi padi di luar Jawa. Angka konversi ini adalah 65 persen gabah akan menjadi beras melalui proses penggilingan, dengan demikian angka konversi yang digunakan adalah 0.65. Jadi total produksi beras luar Jawa dapat dirumuskan sebagai berikut:
QBLT = 0.65 * QPLT …………………..(3.4)
Dimana : QBLT = Produksi beras luar Jawa

c. Persamaan Impor Beras Nasional
Setiap tahun pemerintah menentukan sejumlah beras yang akan diimpor yang mana jumlah Impor Beras Nasional dipengaruhi oleh Harga Impor Beras, Harga Beras Eceran, Produksi Beras luar Jawa, Nilai Tukar dan Jumlah Impor Beras satu tahun sebelumnya, sehingga persamaan jumlah impor beras dapat dirumuskan sebagai berikut:
IMNAS = c0 + c1 HIMT + c2 HBRT + c3 QBLT + c4 NTRT +
+ c5 IMNAS(_1) + u3 …………………………..(3.5)
IMNAS = Impor beras nasional (000 ton)
HIMT = Harga impor beras (US$/ton)
HBRT = Harga eceran beras (Rp/kg)
QBLT = Produksi beras luar Jawa (000 ton)
NTRT = Nilai Tukar (Rp/US$)
IMNAS(_1) = Impor beras nasional satu tahun sebelumnya (000 ton)
Tanda parameter dugaan yang diharapkan:
c1, c3, c4< 0, c2 > 0 dan 0< c5 < 1

d. Perubahan Stok Beras
Perubahan stok beras luar Jawa ini merupakan pengurangan dari jumlah stok akhir tahun sebelumnya terhadap jumlah stok akhir tahun. Data perubahan stok beras luar Jawa ini didapat dari mengalikan stok beras nasional dengan persentase penduduk Jawa dibanding populasi penduduk nasional (32 persen). Peubah perubahan stok beras ini dipengaruhi oleh peubah-peubah; Harga Eceran Beras, Impor Beras Nasional, Produksi Beras luar Jawa dan Stok Beras tahun sebelumnya.
Dengan demikian persamaan struktural yang diperoleh adalah sebagai berikut:
STLJT = d0 + d1 HBRT + d2 IMNAS + d3 QBLT + d4 STLJT(_1) + u4…(3.5)
Dimana:
STLJT = Stok beras luar Jawa (000 ton)
HBRT = Harga beras eceran (Rp/kg)
IMNAS = Impor Beras Nasional (000 ton)
QBLT = Produksi beras luar Jawa (000 ton)
STLJT(_1) = Stok beras luar Jawa satu tahun sebelumnya (000 ton)
Parameter dugaan yang diharapkan: d1 < 0, d2, d3 > 0 dan 0< d4 < 1

e. Permintaan Beras Luar Jawa
Jumlah permintaan terhadap beras untuk konsumsi domestik meliputi konsumsi penduduk, makanan ternak, dan penggunaan lainnya. Tetapi karena kesulitan memperoleh data untuk masing-masing komponen tersebut maka peubah permintaan beras ini hanya mencakup data konsumsi penduduk, karena dianggap paling dominan dalam konsumsi beras domestik sehingga peubah lain diabaikan.
Dalam studi ini persamaan permintaan beras luar Jawa didekati dengan mengalikan jumlah permintaan beras nasional dengan persentase penduduk luar Jawa (32 persen). Dengan demikian permintaan beras luar Jawa dapat dirumuskan sebagai berikut:
QDBT = e0 + e1 HBRT + e2 POPT + e3 QDBT(_1) + u5…(3.6)
QDBT = Jumlah permintaan beras luar Jawa (000 ton)
HBRT = Harga eceran beras (Rp/kg)
POPT = Populasi penduduk luar Jawa (000 orang)
QDBT(_1) = Jumlah permintaan beras luar Jawa tahun sebelumnya (000 ton)
Parameter dugaan yang diharapkan: e1 < 0, e2 > 0 dan 0< e3 < 1

f. Integrasi Pasar
Keterkaitan masing-masing segmen pasar dapat dijelaskan melalui hubunganhubungan harga padi dan harga beras. Harga beras yang akan selalu lebih tinggi dari harga padi karena adanya biaya harga baik berupa margin keuntungan maupun margin yang disebabkan oleh pertambahan biaya transportasi, penyimpanan dan lain-lain.
Dalam studi ini persamaan harga dapat dirumuskan sebagai berikut:
HPDT = f0 + f1 HBRT + f2 HDPT + u6…(3.7)
HBRT = g0 + g1 HDPT + g2 HPDT + u7…(3.8)
Dimana:
HPDT = Harga gabah (Rp/kg)
HBRT = Harga eceran beras (Rp/kg)
HDPT = Harga dasar gabah (Rp/kg)
u6 dan u7 = peubah pengganggu

Identifikasi Model
Dalam identifikasi model struktural ada beberapa unsur yang harus diperhatikan (Koutsoyiannis, 1977) antara lain:
1. Jumlah current endogenous variables dalam model =G
2. Jumlah current endogenous variables yang terdapat pada setiap persamaan =g
3. Jumlah predetermined variables dalam model = K
4. Jumlah predetermined variables yang terdapat pada setiap persamaan = k
Menurut order condition, suatu persamaan dapat diidentifikasikan jika jumlah peubah yang tercakup dalam persamaan lebih besar atau sama dengan jumlah seluruh peubah endogen dikurangi satu. Dengan demikian rumus identifikasi model struktural menurut order condition adalah sebagai berikut:
(G-g) + (K-k) (G-1) atau (K-k) (g-1)
Jika (K-k) sama dengan (g-1) maka persamaan dalam model dikatakan exactly
identified, jika (K-k) lebih kecil dari (g-1), maka persamaan dalam model dikatakan unidentified, dan jika (K-k) lebih besar dari (g-1) maka persamaan dalam model dikatakan over identified.
Dari model struktural yang dirumuskan di atas diketahui bahwa jumlah current endogenous variables sebanyak 10 buah, current exogenous variables sebanyak 10 dan lagged endogenous variables sebanyak 5 buah.
Dengan demikian jumlah seluruh peubah yang tercakup dalam model sebanyak 25 buah terdiri dari 15 peubah predetermined. Berdasarkan kondisi order, maka setiap persamaan struktural adalah over identified.

Motode Pendugaan Model
Jika persamaan dalam model struktural semuanya over identified maka persamaan ini dapat diduga dengan metode LIML (Limited Information Maximum Likelihood), FIML (Full Information Maximum Likelihood), 2SLS (Two Stage Least Squares) atau 3SLS (Three Stage Least Squares). Pemilihan metode di atas disesuaikan dengan tujuan penelitian yaitu untuk mendapatkan koefisien persamaan struktural secara simultan. Pendugaan parameter secara simultan akan membantu simulasi kebijakan secara tepat dan efisien.
Dalam studi ini metode yang akan digunakan dalam menduga parameter
struktural adalah Two Stage Least Squares (2SLS). Metode estimasi 2SLS dibentuk dengan asumsi:
1. Syarat-syarat gangguan harus memenuhi asumsi stochastic sama dengan nol, varians konstan dan kovarians sama dengan nol.
2. Spesifikasi model struktural adalah tepat sekali sejauh yang menyangkut peubah predetermine. Hal ini menyangkut untuk semua peubah predetermine dalam model sudah persis diketahui.
3. Jumlah pengamatan sampel adalah lebih besar dari jumlah peubah predetermine dalam model.
4. Peubah penjelas tidak mengalami kolinearitas sempurna.
Dengan memperhatikan asumsi di atas ada masalah praktis yang timbul yaitu nilai t hitung dan Durbin-Watson statistik tidak valid untuk menduga persamaan struktural dari model persamaan simultan, terutama dengan adanya peubah endogen beda kala pada right hand side dari persamaan struktural.

Sumber Data
Data yang digunakan adalah data series tahun 1969-1998 berasal dari berbagai sumber yakni: Biro Pusat Statistik, Badan Urusan Logistik, IMF, IRRI dan Departemen Pertanian.

HASIL DAN PEMBAHASAN
Adapun hasil pendugaan parameter peubah-peubah dalam Model Pasar Beras di Luar Jawa adalah seperti yang tercantum pada Tabel 4 dengan rincian sebagai berikut:
Areal Panen dan Produktivitas Padi di Luar Jawa
Koefisien determinasi (R2) dari persamaan areal dan produktivitas cukup tinggi yakni masing-masing 0,8597 dan 0,8719. Nilai R2 dari persamaan areal sebesar 0,8597 berarti 85,97 persen dari variasi dari areal panen dapat dijelaskan oleh variasi peubahpeubah penjelasnya. Sedangkan nilai R2 dari persamaan produktivitas sebesar 0,8719 menunjukkan bahwa 87,19 persen variasi peubah produktivitas dapat dijelaskan oleh peubah peubah penjelasnya.
Hanya satu peubah penjelas (harga padi) pada persamaan areal yang siknifikan pada taraf = 0,10. Sedangkan peubah-peubah penjelas lainnya (harga dasar ril, harga komoditi jagung, tingkat bunga pinjaman, curah hujan, luas areal panen satu tahun sebelumnya) tidak siknifikan pengaruhnya terhadap areal panen. Begitu juga pada persamaan produktivitas hanya satu peubah penjelas yang siknifikan pada taraf = 0,01 (produktivitas padi satu tahun sebelumnya).
Informasi di atas menunjukkan bahwa yang mempengaruhi perilaku produksi padi di Luar Jawa adalah harga padi, dimana petani merespon harga yang tinggi dengan meningkatkan produksi padi dengan cara meningkatkan luas areal (ekstensifikasi).
Adapun nilai elastisitas luas areal terhadap harga adalah inelastis yakni 0,4725 yang berarti jika harga padi (gabah) naik sebesar 10 persen maka akan meningkatkan areal padi sebesar 4,725 persen.
Stok Beras Luar Jawa
Koefisien determinasi (R2) persamaan stok beras luar Jawa adalah 0,5358 yang menunjukkan 53,58 persen variasi peubah stok beras luar Jawa dapat dijelaskan oleh variasi peubah-peubah penjelasnya.
Harga beras eceran (tingkat rumah tangga) berpengaruh siknifikan terhadap stok beras luar Jawa pada taraf = 0,10. Nilai elastisitas stok beras luar Jawa terhadap harga beras eceran adalah inelastis dalam jangka pendek (-0,771) dan elastis dalam jangka panjang (-2,3914). Artinya jika harga beras eceran naik sebesar 10 persen maka akan menurunkan stok beras luar Jawa sebesar 7,71 persen dalam jangka pendek dan 23,914 persen dalam jangka panjang. Penurunan stok beras yang diakibatkan oleh naiknya harga beras eceran, karena stok beras yang ada digunakan untuk operasi pasar.
Tabel 4. Hasil Pendugaan Parameter Peubah-peubah dalam Model Pasar Beras
di Luar Jawa.
No Peubah Parameter Dugaan t-hitung Elastisitas
Jangka Pendek Jangka Panjang
01. Intersep 2295, 3286 – – –
HPDT 11,702 2,072 (C) 0,4725 –
HDPTRIL 298,709 0,620 (ns) – –
HKJT -10,831 -1,319 (ns) – –
SUKT 17.57 0,809 (ns) – –
CHLT 0,0935 0,439 (ns) – –
ALJW_ 1 0,1182 0,638 (ns) – –
R square = 0.8597 DW = 2,406 Koefisien adjustment = –
02. Intersep 3,445 – – –
HPDT 0,001 1,298 (ns) – –
ILT -0,001 -0,690 (ns) – –
ALJW 0,00002 0,143 (ns) – –
YLPT_ 1 0,665 4,454 (A) – –
R square = 0.8719 DW = 1,912 Koefisien adjustment = –

03 Stok Beras luar Jawa
Intersep -742,796
HBRT -0,989 -1,678 (C) -0,771 -2,3914
IMNAS 0,1122 1,709 (C) 0,2248 0,6972
QBLT 0,1212 2,168 (B) 2,66 8,25
STLJT_ 1 0,6776 3,439 (A) – –
R square = 0.5358 DW = 1,953 Koefisien adjustment = 0,3224
04 Impor Beras Nasional
Intersep 2582,364
Intersep 2582,364 – – –
HIMT -3,1022 -1,585 (D) -1,089 -2,1195
HBRT 7,8268 2,870 (A) 3,045 5,926
QBLT -0,1449 -0,928 (ns) – –
NTRT -2,1654 -2,327 (B) -2,869 -5,58
IMNAS_ 1 0,4862 2,186 (B) – –
R square = 0.6102 DW = 1,644 Koefisien adjustment = 0,5138

Keterangan: Huruf dibelakang nilai t-hitung menunjukkan taraf nyata () yaitu:
A berarti berbeda nyata dengan nol pada taraf = 0,01
B berarti berbeda nyata dengan nol pada taraf = 0,05
C berarti berbeda nyata dengan nol pada taraf = 0,10
D berarti berbeda nyata dengan nol pada taraf = 0,15
E berarti berbeda nyata dengan nol pada taraf = 0,20
Impor nasional juga berpengaruh siknifikan terhadap stok beras luar Jawa pada taraf = 0,10. Elastisitas stok beras luar Jawa terhadap impor nasional ini adalah inelastis baik dalam jangka pendek (0,2248) dan jangka panjang (0,6972).
Hal ini menunjukkan jika impor beras nasional naik sebesar 10 persen maka akan meningkatkan stok beras luar Jawa sebesar 2,248 persen dalam jangka pendek dan 6,972 dalam jangka panjang.
Produksi beras luar Jawa juga siknifikan pengaruhnya terhadap stok beras luar Jawa pada taraf = 0,05. Adapun nilai elastisitas stok beras luar Jawa terhadap produksi beras luar Jawa adalah elastis baik dalam jangka pendek (2,66) dan jangka panjang (8,25) yang menunjukkan jika produksi beras luar Jawa naik sebesar 10 persen maka akan meningkatkan stok beras luar Jawa sebesar 26,6 persen dalam jangka pendek dan 82,5 persen dalam jangka panjang.

Impor Beras Nasional
Koefisien determinasi (R2) persamaan impor beras nasional adalah 0,6102 yang berarti 61,02 persen variasi peubah impor beras nasional dapat dijelaskan oleh variasi peubah-peubah penjelasnya.
Harga impor berpengaruh siknifikan terhadap impor beras nasional pada taraf = 0,15. Adapun nilai elastisitas impor beras nasional terhadap harga impor adalah elastis baik dalam jangka pendek (-1,089) dan jangka panjang (-2,1195). Artinya jika harga impor naik sebesar 1 persen maka akan menurunkan impor sebesar 1,089 persen dalam jangka pendek dan 2,1195 persen dalam jangka panjang. Sedangkan harga beras eceran berpengaruh positip dan siknifikan terhadap impor beras nasional pada taraf =0,01 yang menunjukkan semakin tinggi harga beras eceran akan cenderung semakin meningkatkan impor beras nasional.
Nilai elastisitas impor nasional terhadap harga beras eceran adalah elastis baik dalam jangka pendek (3,045) dan jangka panjang (5,926) yang menunjukkan jika harga beras eceran naik sebesar 10 persen maka akan meningkatkan impor beras nasional sebesar 30,45 persen dalam jangka pendek dan 59,26 dalam jangka panjang.
Adapun produksi beras luar Jawa tidak siknifikan berpengaruh nyata terhadap impor beras nasional. Walaupun demikian nilai koefisien parameter dugaan peubah produksi beras luar Jawa yang negatif menunjukkan bahwa peningkatan produksi beras luar Jawa akan cenderung menurunkan impor beras nasional.
Nilai tukar ternyata berpengaruh siknifikan terhadap impor beras nasional pada taraf =0,05. Sedangkan nilai elastisitas impor beras nasional terhadap nilai tukar (Rp per US dolar) adalah elastis baik dalam jangka pendek (-2,869) dan jangka panjang (-5,58). Hal ini menunjukkan jika rupiah mengalami depresiasi sebesar 1 persen maka akan menurunkan impor sebesar 2,869 persen dalam jangka pendek dan 5,58 dalam jangka panjang.

Permintaan Beras Luar Jawa
Koefisien determinasi (R2) persamaan permintaan beras luar Jawa adalah 0,9892 yang berarti 98,92 persen variasi peubah permintaan beras luar Jawa dapat dijelaskan oleh variasi peubah-peubah penjelasnya.
Harga beras eceran berpengaruh tidak siknifikan terhadap permintaan beras luar Jawa, namun tanda koefisien harga beras eceran yang negatif telah sesuai dengan yang diharapkan yakni jika harga eceran beras naik maka jumlah beras yang diminta cenderung turun.
Populasi berpengaruh positip dan siknifikan terhadap jumlah permintaan beras pada taraf =0,01. Sedangkan nilai elastisitas permintaan beras terhadap populasi penduduk luar Jawa ini adalah elastis baik dalam jangka pendek (1,388) dan jangka panjang (1,913). Ini berarti jika penduduk bertambah 1 persen maka akan meningkatkan jumlah beras yang diminta sebesar 1,388 persen dalam jangka pendek dan 1,913 persen dalam jangka panjang.

Harga Padi dan Beras
Koefisien determinasi (R2) persamaan harga padi adalah 0,9939 yang berarti 99,39 persen variasi peubah harga padi dapat dijelaskan oleh variasi peubah-peubah penjelasnya. Sedangkan koefisien determinasi (R2) persamaan harga beras adalah 0,9954 menunjukkan 99,54 persen variasi peubah harga beras dapat dijelaskan oleh variasi peubah-peubah penjelasnya.
Dari persamaan harga padi tampak bahwa harga padi lebih ditentukan oleh harga dasar dibanding harga eceran beras. Hal ini ditunjukkan oleh peubah harga dasar yang siknifikan pengaruhnya terhadap harga padi pada taraf =0,01, sedangkan harga beras eceran tidak siknifikan pengaruhnya terhadap harga padi.
Elastisitas harga padi terhadap harga dasar adalah inelastis (0,8436) yang berarti jika harga dasar naik 10 persen maka akan menaikkan harga padi sebesar 8,436 persen.
Dari persamaan harga beras eceran dapat diketahui bahwa harga dasar dan harga padi berpengaruh siknifikan terhadap harga beras eceran pada taraf =0,01. Adapun nilai elastisitas harga beras eceran terhadap harga dasar dan harga padi adalah inelastis yakni masing-masing sebesar 0,412 dan 0,55. Hal ini berarti jika harga dasar naik sebesar 10 persen maka harga beras eceran naik sebesar 4,12 persen, sedangkan kenaikan harga padi sebesar 10 persen maka akan menaikkan harga beras eceran sebesar 5,5 persen.

KESIMPULAN DAN IMPLIKASI KEBIJAKAN
Kesimpulan
Dari uraian di atas dapat disimpulkan hal-hal sebagai berikut:
1. Perilaku areal panen padi di luar Jawa ternnyata hanya dipengaruhi oleh harga padi. Walaupun demikian elastisitas areal panen terhadap harga padi adalah inelastis. Fenomena ini menunjukkan harga padi akan mendorong petani meningkatkan produksi padi melalui peningkatan areal (ekstensifikasi), bukan melalui peningkatan produktivitas (intensifikasi) karena harga padi tidak siknifikan pengaruhnya terhadap produktivitas padi.
2. Produksi beras luar Jawa tidak siknifikan pengaruhnya terhadap impor beras menunjukkan produksi beras di luar Jawa belum mampu menjadi kontributor yang siknifikan dalam mengurangi impor beras nasional. Hal ini menunjukkan luar Jawa belum berperan besar sebagai pensuplai beras nasional.
3. Permintaan beras di luar Jawa tidak dipengaruhi oleh harga beras tetapi sangat ditentukan oleh jumlah penduduknya. Hal ini menunjukkan permintaan beras luar Jawa di masa mendatang akan semakin meningkat sejalan dengan semakin bertambahnya jumlah penduduk luar Jawa.
4. Harga padi di luar Jawa sangat ditentukan oleh harga dasar namun respon (elastisitas) harga padi terhadap harga dasar adalah inelastis (kurang dari satu).
5. Harga beras eceran luar Jawa dipengaruhi oleh harga dasar dan harga padi dengan nilai elastisitas harga beras eceran terhadap harga dasar dan harga padi itu adalah inelastis.

Implikasi Kebijakan
Kebijakan yang dapat menolong harga padi di tingkat petani seperti harga dasar dan subsidi input adalah penting untuk tetap diterapkan untuk memacu produksi beras di luar Jawa. Hal ini karena harga dasar ini akan mempengaruhi harga padi dan selanjutnya harga padi akan memacu petani untuk meningkatkan produksi melalui ekstensifikasi.

DAFTAR PUSTAKA
Colman, D. and T. Young. 1990. Principles of Agricultural Economic Market and
Prices

in Less Developed Countries. Cambridge University Press. Cambridge.
Hutauruk, J. 1996. Analisis Dampak Kebijakan Harga Dasar Padi dan Subsidi
Pupuk

Terhadap Permintaan dan Penawaran beras di Indonesia. Tesis Magister Sains
Institut Pertanian Bogor.

Intrigator, M.D. 1978. Econometric Models, Techniques and Application.
Prentice-Hall International. New Delhi.

Irawan, A. 1997 (a). Kebijakan Harga dan Keberlanjutan Produksi Padi. Ekonomi
Dan Keuangan Indonesia (EKI): 15 (4): 579-586.
——–. 1997(b). Kebijakan Harga dan Keberlanjutan Produksi Padi di luar Jawa.
Kompas 13/2/ 1997.

——– . 1998. Analisis Respon Penawaran Padi Sawah dan Ladang di Jawa dan
Luar Jawa. Studi Respon Penawaran. Tesis Magister Sains Institut
Pertanian Bogor.

——–. 1999. Analisis Respon Penawaran Padi Sawah dan Ladang di Jawa dan
Luar Jawa. Ekonomi dan Keuangan Indonesia (EKI): 17 (1):19-31.

——–. 2000. Perilaku Suplai Padi Indonesia dan Implikasinya terhadap
Peningkatan Produksi dalam Prosiding Perspektif Pembangunan Pertanian dan Kehutanan tahun 2001 ke Depan. ISBN:979-8094-68-9.

Kasyrino, F. 1996. “Arah Pengembangan Agribisnis di Pulau Jawa pada Abad
21”.