A multiplicative compound criterion is used to generate a unique design which can simultaneously fulfil two important

objectives: parameter estimation and response prediction. Thegene ralised variance criterion is chosen for the first objectivec while the variance function or the integrated variance criterion is fOr the second

Haron, Kassim

Universiti Putra Malaysia Institutional Repository

Pertanika 15(2),151-157 (1992)Kriterium Majmuk bagi Penganggaran Parameter dan Peramalan SambutanKASSIM HARONjabalan MatemalikUniversiti Pertanian Malaysia43400 UPM, Serdang, SelangorDamlEhsan, Malaysia.Kata kunci: objektif, kriterium, rekabentuk, varians, kecekapan.ABSTRAKK17lerium majmuk berdaya darab digunakan untuk menjana sualu rekabentuk unik yang dapat memenuhi dua ob-jektifpenting serenlok: penganggaran parameter dan peramalan sambutan. Kriterium varians teritlak dipilih untuhobjektifpertama manakala kriterium fungsi varians atau kriteriwn uarians terkamir untuk yang kedua.ABSTRACTA multiplicative compound criterion is used to generate a unique design which can simultaneouslyfulfil two importantol!jeclives: parameter estimation and responseprediction. Thegeneralised variance criterion is chosenfor thefirst objectivewhile the variance function or the integrated variance criterion is fOT the second.,:1;. I,. ,(1.2)(l.I)dengan mo,;e ';1,m;'2,~e ~1, L.., ,manakala kriterium berdaya darab pula ditulissebagaimengambil kira objektif-objektif seseorang penye-lidik secara serentak perlu difikirkan. Pengop-timuman yang sesuai, jika dilakukan ke atas kri-terium seperti itu dapat menghasilkan suatu re-kabentuk unik yang mampu digunakan bagi men-capai kesemua objektif. Setakat ini kriteriummajmuk yang pernah dibincangkan oleh parapenyelidik di dalam bidang rekabentuk ujikaji dapatdipecahkan kepada dua kategori iaitu kriteriumberdaya tambah (Hill et aL 1968; Lauter 1974;Rasch & Herrendorfer 1980; Cook & Nachtsheim1982) dan kriterium berdaya darab (Lauter 1974;Atkinson & Cox 1974; Smith 1987).Secara matematiknya, kriterium majrnuk ber-daya tambah boleh ditulis sebagaiKebanyakan penyelidikan yang dijalankan ber-hubung dengan rekabentuk permukaan sambutansetakat ini lebih tertumpu kepada masalah pen-janaan rekabentuk terbaik bagi memenuhi objektiftunggal. Dan itu, pertimbangan yang sew~arnyaperlu diben kepada masalah yang dihadapi olehpenyelidik yang mempunyai beberapa objektiftertentu untuk dicapai daripada sesuatu ujikaji.Dalam hal ini, suatu prosedur penjanaan rekaben~tuk yang lebih sesuai dan praktikal perlu dicari.Menghasilkan sesuatu set data dengan meng-gunakan rekabentuk yang berbeza bagi setiap ob-jektif sudah tentil merugikan kerana prosessedemikian memeriukan masa dan kos perbe-lanjaan yang lebih. Harus diingat bahawa rekaben-tuk yang berbeza diperlukan kerana rekabentukyang cekap bagi memenuhi suatu objektifmungkinkurang cekap bagi memenuhi objektif yang lain.Misalnya, rekabentuk yang cekap bagi objektifmendiskriminasikan beberapa model tak linearmungkin kurang cekap jika digunakan pula bagitujuan menganggar parameter model (Hill et al.1968) .Bagi mengatasi masalah objektif berganda,penggunaan suatu kriterium majmuk yang dapatl.PENDAHULUANKASSIM HARaNdan(2.3),,' = 1 (2.1)andaikan(2.4)r:::; 1. 2•...• 5n<Nx ,Z,]Z,Z,n p. X~-+lZr:::;N~::' J Jjzl (~:+~~Xj+~:XD4dandenganmin GV ( x,~) =min{ I (F F)I},~ -- 1! --~m!n{I(~' ~)Ir'Matriks varians-kovarian boleh diperoleh denganterlebih dahulu melakukan penglinearan fungsimodel di sekitar nilai anggaran awal 13*, cukuphampir dengan p. Seterusnya. deng;;'n meng-gunakan kembangan siri Taylor hingga ke pering-kat pertama, kita peroleh sambutan hampiran2. LATARBELAKANGTEORImaka matriks maklumatnya boleh ditulis sebagaiJika ditakrifkanF = Ifiji ; i = 0,1, 2,j = 1, 2, ... , Ndengan2.1 Penganggaran ParameterSeperti yang telah dinyatakan. kriterium variansteritlak, ditandakan dengan GV, diperlukan bagimemastikan bahawa dengan pemilihan rekaben-tuk yang sesuai, setiap parameter dapat dianggardengan kepersisan yang tinggi secara serentak.Kriterium ini menyarankankan kita supayameminimumkan varians teritlak iaitu penentubagimatriksvarians-kovarians (~' ~)-la2 terhadaparas-aras faktor ~ :::; (>11' X2' ...• Xn )' • Secara mate-rnatiknya, masalah ditulis sebagai(1.3)Xj"O;~o'~,>Oj~ 1,2, .... ,Nxy ~ c:----=,.-'-'---=---;;j ~ + ~ x + ~ x'o 1 J 2 JE(y) ~,Dalam kes iniC; mewakili kriterium ke-i8, mewakili parameter pemberat yang sepa-dan dengan kriterium ke-i. Nilai yang diberikepada 8 j akan mencerminkan tahap keuta-maan objektifke-i kepada penyelidik.Di dalam kajian penulis, kriterium majmuk ber-claya darab menjadi pilihan kerana ia bukan sahaJamudah dikendalikan tetapijuga rekabentuk multi-tujuan I?M yang akan dihasilkan adalah tak varianterhadap perubahan dalam nilai parameter bagimodel y ~ f (x; P) (Smith 1987). Di peringkat inijuga. hanya dua kriterium yang mewakili ciuaobjektif berlainan sahaja dipertimbangkan. Kri-terium-kriterium yang dipilih akan mewakili ciuaobjektif penting yang biasanya menjadi tumpuankepada para penyelidik. Objektif-objektif tersebutialah:(a) menganggar parameter model dengan ke-persisan yangtinggi. lui penting kerana parameter-parameter yang terdapat di dalam sesuatu modelbiasanya mewakili fenomena tertentu dalam ujikajiyang dijalankan.(b) meramal sambutan, juga dengan kepersi-san yang tinggi. pada aras tertentu sesuatu faktoratau pada aras-arasyang terkandung di dalam suaturantau ujikaji X.Untuk objektif (a), kriterium yang diketahuipaling sesuai digunakan ialah kriterium variansteritlak (Box & Lucas 1959) manakala kriterium-kriterium fungsi varians atau varians terkamir pulasesuai bagi mencapai objektif (b).Di sepanjang kajian ini. polinomial songsangkuadratik satu faktor dengan asalan sifar, ditan-dakan dengan IFQIP(Po, PI' ]3,) dipilih sebagaimodelnya memandangkan kesuaian model sepertiini dalam memperihalkan banyak ujikaji dalambidang pertanian, biologi dan industri (NeIder1966) .Jika ralat ujikaji diandaikan bertabur secaranormal dengan min sifar dan varians malar 0'2. makasambutan min kej, boleh ditulis sebagai152 PERTANlKA VOL. 15 NO.2, 1992KRITERIUM MAJMUK BAGI PENGANGGARAN PARAMETER DAN PERAMAlAN SAMBUTAN2.2 Peramalan Sambutanx; Var B, + 2{ x, kov(B o'~,) +x; kov(B),J + x; kov(B ,J,) } ]2.2.1 Kriterium fungsi variansSekiranya peramalan sambutan di suatu arastertentu Xo dengan varians minimum menjadimatlamat seseorang penyelidik. maka penggunaanfungsi varians V(xo) diperlukan. Rumus bagi V(xo)diterbitkan seperti berikut:Dari (2.2), (2.3) dan (2.4) kiGl perolehVar(y) '" vaJtB ,E,,]l,~ndengan-, fn = dx o'"i ~o, I, 2r~O,I;s~I,2r"s2.2.2 Kriterium varians terkarnirAda ketikanya pula seseorang penyelidik itu bermi-nat untuk membuat ramalan terhadap nilai-nilaisambutan pacta aras faktor yang terkandung didalam suatu ran tau ujikaji X. Bagi objektif sepertiini, kriterium varians terkamir IV(xo) dirasakan pa-ling sesuai digunakan kerana ia dapat memastikanbahawa, secara puratanya, varians bagi sambutansuaian yadalah minimum didalam X,. Varians terka-mir, N(x,) diGlkrifkan sebagaiN( x,) ~ NQo" fVar y(x,) dx,v'denganb.. ~Var(~ ).. ,b.. ~ kov (~ "B.)(2.5)[Var B + x' Var B +0, '(f~: x;]'L.ovar e2=denganvar B ~ (z z -z')" . 0 2;, If' ¥Ivar 81 = (Z]Z5-Z:)0'If' ¥IJadi, daTi takrif fungsi varians (Box & Hunter1957), makaV(xo)~No" Var[y(x,»)3. REKABENTUK DWITUJUAN DM DAN DM •SERTA PERBINCANGAN -danVar y(x,) seperti yang boleh didapati dalam (2.5)Dalam bahagian ini tumpuan diberikan dahulukepada masalah penjanaan rekabentuk tunggalyang mengambil kira objektif-objektif:(i) menganggar parameter ~"~, dan~,(ii)meramalkan sambutan di suatu aras/titiktertentu xo.Bersesuaian dengan model yang dipilih iaitu kua-dratiksongsang lFQIP(~"~,,~,) yang seringdikaitkandengan pemodelan hasil tanaman/aras baja dalambidang pertanian, maka peramalan hasil di arasbaja maksimum Xm.k~ (::J dengan kepersisanyangmemuaskan merupakan antara matlamat parapenyelidik dalam bidang tersebut.Dalam hal ini, dicadangkan penggunaan kri·terium majrnuk berdaya darabC,,:(GV)"'(V(x,»)l-O; o,;e,;1 (31)dengan q mewakili bilangan parameter dalammodel. Dalam kes ini q ~ 3Tidaklah begitu sukar untuk menunjukkan, secaraalgebra. bahawa peminirnuman eM akan meng-(2 1 Z, -Z:)Ir ¥IA A _ (Z,Z.-Z,Z,) 0'(p"p,)- iI::' ¥I(z,Z.-Z~)(B"B,) ~ iI::' ¥I 0';A A ~(Z,Z,-Z,Z,) 0'(p "p,) IF' ¥IkovkovkovPERTANIKA VOL. 15 0.2,1992 153KASSIM HARONJADUALIRekabenluk dwitujuan D apabila X o~ 1M_ .\1 V J-1 0 /1J 2hasilkan penyelesaian optimum yang mempunyaisifa(~sifatberikut:(a) tiga litik berbeza xM = (xl' x2• x,) = xG\'dengan XLV mewakili titik D-optimum (Box& Lucas 1959).(b) nishah cerapan pada litik-ritik Xl' 'S dan x,masing-masing bersamaan p, = p', p,~ 1-2p'dan p, = p'; 0 p' 1.D = X[X, x, x,l- M P PI p! r,Jdengan Pi mewakili nisbah bilangan cerapan yangdiambil pacta titik ke-x j • i = 1, 2, 3, dijana untukmodeIIFQIP(9, 1,1) padasuatujujukan nilai pa-rameter pemberat e dengan bantuan rutin peng-optimuman yang sesuai yang terdapaldalam perpus-takaan NAG. Jelas kepada kita rekabentukI? M yang terhasil mempunyai sifat-sifatseperti yangdinyatakan di alas.KaJi3n seterusnya mendapati bahawa untukseliap niIai enishah bilangan cerapan pacta ketiga-tiga [i{ik yang berpadanan tidak berubah untuknilai ~ * yang lain.Sementara itu, Jadual 2 pula dimuatkan de-ngan darjah kebagusan setiap rekabentuk yangModellFQIP e D_ .\l(~" ~" ~,)1.0 x [0.90 3.00 9.99]P 1/3 1/3 1/30.8 x [0.90 3.00 9.99]P 0.27 0.46 0.270.7 x [0.90 3.00 999]P 0.23 0.54 0.23(9, I, I) 0.5 x [0.90 3.00 9.99]P 0.17 0.66 0.170.2 x [0.90 3.00 9.99]P 0.07 0.86 0.070.1 x [0.90 3.00 999]P 0.Q3 0.94 0.030.0 x [0.90 3.00 9.99]P 0.0 1.0 0.0(3.2)c,46.7053.3066.7086.7093.00100.00Kecekapan E (%)c,100.089.5578.4649.9210.353.66e1.00.80.70.50.20.10.0Perkara yang menarik yang dapat dirumuskan dari-padaJadual2 ialab peratus kecekapan E(%) yangagak peka terhadap perubahan dalam nilai 8.Bagaimanapun, kepekaan yang diperlihatkan ti-daklah memeranjatkan kerana kita sebenarnya telahmelakukan suatu kompromi dengan mengga-bungkan dua kriterium yang masing-masing akanmenghasilkan rekabentuk terbaikyang berlawananbentuk. Rekabentuk D-optimum iaitu rekabentukterbaik bagi penganggaran parameter memerlukancerapan yang sama banyak diambil pada tiga acasltitik yang berbeza. Sebaliknya, kriterium fungsivarians V(xo ) pula menghasilkan rekabentuk ter-baikdengan kesemuacerapan perludi ambil hanyapada titik x, itu 'endiri (,ila lihatJadual 1). Inilahpunca yang nyara kenapa pencapaian E(%) bag! D"tidak begitu terserlah. Sungguhpun begitu, nila(S= 0.7 masih mampu memberikan peratus kecekapanyangagak baikiaitu 78%jikadigunakan bag! maksudrnenganggar parameter dan 53% bagi ramalansambutan di XIIlU' Perlu ditegaskan bahawa ke-lemahan dari segi pencapaian peratus kecekapanseperti ini boleh dielakkan sekiranya kompromiyang melibatkan kriterium-kriterium yang lebihsesuai digunakan. Ini akan dijelaskan dalam per-bincangan berikutnya. Sementara itu kecekapanJADUAL2Kecekapan E (%) hagi rekabenruk I?IoldenganC;, i = 1, 2 mewakili nilai kriterium ke-i1?~ mewakili rekabentuk optimum bagi kri-terium ke-id!.jana dengan ~. ~ (9, I, 1) dalam memenuhio bjektif yang sepadan. Kebagusan setiap rekaben-[uk tersebut diukur berdasarkan kepada peraruskecekapan E(%) yang ditakrifkan sebagai(c,Il~:)E (%) = (C,"..\) x 100D- ..Jadual I memaparkan rekabentuk dwitujuanyang ditakrifkan sebagai154 PERTAI'olIKA VOL. 15 NO.2, 1992KRITERIUM MAJMUK BAGI PENGANGGARAN PARAMETER DAN PERAMAlAN SAMBUTAND. _ X[0.90 3.00 9.99]-" - p 0.00 1.00 0.00JADUAL3Rekabentuk dwitujuan n' dalam rantau ujikaji x:-"(0,12)Selain itu. kitajuga dapati bahawa tidak terdapatperubahan yang banyak di dalam rekabentuk J? ~lapabila nilai-nilai 8 beralih dari °ke I.1.0 x [0.90 3.00 9.99]p 1/3 1/3 1/30.8 x [0.91 3.02 9.88]P 0.32 0.33 0.350.6 x [0.91 3.05 9.74]P 0.30 0.33 0.37(9, I, 1) 0.4 x [0.92 3.07 9.58]P 0.28 0.33 0.390.2 x [0.93 3.09 9.49]P 0.25 0.33 0.410.0 x [0.95 3.12 9.10]P 0.22 0.32 0.454. PENGGUNAAN DAN ILUSTRASIJusteru itu, rekabentukX[0.92 3.07 9.58]p 0.28 0.33 0.39yang dijana dengan e ~ 0.4 misalnya,jika diguna-kan bagi menganggar parameter dan meramalsambutan, masing-masing sebagai objektiftunggal,di dalam ramau X:(O, 12), dijangka dapat meng-hasilkan kepersisan yang tinggi. Sebenarnya, de-ngan menggunakan sukatan kecekapanseperti yangdiberikan dalam (3.2), masing-masingdidapati men-catatkan kepersisan yang melebihi 97%.Hasil ujikaji yang dijalankan oleh Morrison et aL(1980) dipertimbangkan. Dalam kajian tersebut,rekabentuk sarna ruang dengan enam aras yangberbeza telah digunakan untuk mengutip data bagimemperihalkan hubungan antara hasil 'ryegrass'dan aras baja nitrogen. Set-set data yang dikutipdaripada 21 kaw3san pertanian digunakan untukpenyuaian model. Polinomial songsang kuadratikdidapati model yang sesuai.Nilai parameter f3 bagi model yang mewakiliKawasan 6 (Kassim f989) dipilih sebagai nilaipenganggar awal ~ *. Sekiranya menganggar pa-rameter model ~ sebagai objektif utama dan me-ramal hasil makSimulll sebagai objektif sampinganpenyelidik maka penggunaan kriterium C" dengan8 = 0.7 menghasilkanD' = x[ 281.4 689.6 1690]- M P 0.23 0.54 0.23Rekabentuk yang terhasil menyarankan supayapengutipan datadilakukan dengan mengambil23%daripadakeseluruhan cerapan masing-masingpadaaras baja 281.4 kg/ha dan 1690.0 kg/ha manakalacerapan selebihnya dikutip pada aras 689.6 kg/ha.Data yang dikutip kemudiannya digunakan umukmenganggar ~ dan juga untuk meramal hasilmaksimum dengan menggunakan kaedah statistikyang sedia ada atau menggunakan pakej statistikseperti SAS (PROC NUN).Sebagai ilustrasi, dua set data yang mengan-dungi 30 bacaan sambutan untuk menganggar pa-rameter model dijana secara simulasi denganmenggunakan kedua-dua jenis rekabentuk iaitu:(i) rekabenruk klasik sama ruang dengan 5bacaan masing-masing dicerap pada aras baja(dalam Kg) 281,553,844,1125, 1405 dan 1690.(ii) rekabentuk yang terhasil dengan meng-9.10]0.4683.120.32ModellFQIP(~" ~" ~,)tetapix [0.95O· -- " - P 0.22rekabentuk di dalam keadaan lain iaitu dengannilai J3 * yang berlainan didapati tidak beruhahuntuk- setiap nilai 8.Sekiranya objektif kedua masih dalam me-ramalkan sambutan tetapi kali ini dengan mengam-bit kim sambutan-sambutan pacta aras faktor yangterkandung di dalam sualu fan tau X, maka dica-dangan kri terium majmuk berikut:HC,..:(GV)"'(IV(x,» ;0,;8';1 (3.3)Sam set rekabentuk I?~l yang dijana dalam rantauX: (0, 12) dapat dilihat dalam Jadual 3. Jelas ba-hawa, tidak seperti kes terdahulu, rekabentuk opti-mum bagi penganggaran parameter (8 ~ I) danperamalan sambutan (8 = 0) masing-masing me-ngandungi tiga titik yang berbeza. Misalnya, dengan8=0155KASSIM HARONgunakan kriterium majrnuk eM dengan e= 0.7 se-perti yang dicadangkan.Proses simulasi dengan 50 lanan dilakukandengan menggunakan nilai parameter ~o - 4.00 xIO-',~, =-4.04 x 10° dan ~, = 8.41 X 10-' yang dian-daikan nHai sebenar parameter-parameler model(Kassim 1989). Ralat ujikaji (andaikan bertaburannormal dengan min sifar dan varians malar) den-gan pekali ubahan, C.v ~ 10% dimasukkan kedalam setiap cerapan yang dijana pada aras-arasbaja tersebut. Analisis statistik untuk mencari nilai-nilai ~'IYtobk' Var ~ dan Var ~.w.~ [Ylol••=(2~+ ~ ) ] kemudiannya dijalankan ke atas set-setda~ tersebut. Hasilnya dimuatkan di dalamjadual4(a) danJaduaI4(b)JADUAL4(A)Nilai-nilai ~ dan YM* yang dihasilkan oleh rekabentukdwitujuan I?M dan rekabemuk klasik sarna ruang I?KRekabentukPenganggarDM D,eo-, -,3.908 x 10 4.02 x 10e, -, -,-4.12 x 10 -4.11 x 10e,-, -,8.50 xlO 8.11 x 10~"" 13,310.77 13,386.05• jADUA!; 4(B)Nilai·niJai Var Pdan Var Y"Uk yang dihasilkan olehrekabentuk dwit~juan J?~ dan rekabentuk klasik sarnaruang 1?RRekabentukVarian !?M !?,eo 0.0784 0.2275e, 0.6475 0.5289e. 0.3565 0.6649~""" • •10.76 x 10 18. 72 x 10Jika diteliti kedua-<iuaJadual4(a) danJadual4(b), jelas bahawa untuk menganggar parameterdan juga meramal hasil maksimum secara serentakmaka rekabentuk yang dicadangkan mernpunyaikelebihan berbanding dengan rekabentuk klasiksarna ruang kerana varians yang sepadan yang di-catat masing-rnasing didapati lebih kedl. Denganitu rekabentukdwitujuan boleh dikatakan suaN re-kabenlUk yang lebih tinggi daIjah kecekapannya.5. KESIMPtlLA>"IHasil daripada kajian yang dijalankanjelas menun-jukkan bahawa kriterium majrnuk berdaya darabmampu berfungsi sebagai alat dalam menerbitkansuatu rekabentuk unik yang dapat dimanfaatkanbagi memenuhi dua objektif penting iaitu meng-anggar parameter dan meramal sambutan. Disamping itu, dapat juga disimpulkan bahawa kri-terium individu yang membentuk kriteriummajmuk memainkan peranan yang penting dalammenentukan kecekapan rekabentuk dwitujuanOM /DM ' yang terjana.RUJUKANATKINSON,A.G. and D.R. Cox. 1974. Planning Experi-ments for Discriminating between Models,journalofthe Royal Statistical Society Ser. B 36: 321-348.Box, G.E.P andJ.S. HUNTER. 1957. Multifactor Experi-mental Designs for Exploring Response Surfaces.Annals ofMathematical Statistics 28: 195-241.Box, G.E.P and H.L. LUCAS. 1959. Design of Experi-ments in Nonlinear Situations. Bio'ffU!tn"ka 46: 77-90.COOK, R.D. and GJ. NACIITSHEIM. 1982. Model Robust,Linear-optimal Designs. Technometrics 24: 49-54.HILL, WJ., w.G. HUNTER and D.W. WICHERN. 1968. AJoint Design Criterion for the Dual Problem ofModel Discrimination and Parameter Estimation.Techncmetrics 10: 145-160.KAssIM, H. 1989. Optimal Design Spacing for theEstimation of the Quadratic Inverse PolynomialParameters. Pertnnika 12(2): 239-244.LAUTER, E. 1974. Experimental Design in a Glass ofModels. Mathematuche operationsforschung undStatistik 5: 379-398.MORRlSON,j, M.Y.JACKSON and P.E. SPARROW. 1980.The Response of Perennial Ryegrass to FertiliserNitrogen in Relation to Climate and Soil. Grass-land Research Institute Technical Report 27.NELDER,jA 1966. Inverse Polynomials, a Useful GroupofMultifactor Response Functions. Biometrics 22:128-141.156 PERTANlKA VOL. 15 NO.2, 1992KRJTERIUM MAJMUK BAGl PENGANGGARAN PARAMETER DAN PERAMAlAN SAMBUTANRAsCH, D. and G. HERRENDORFER. 1980. ExperimentalDesign and Optimal Decision 11. Estimation ofthe Slope ofthe LinearRegression. MathematischeOperalicrnsJamhungundStatislik, Ser. Stat. 11: 125-136.SMmI,].R. 1987. Design and Experiments for the Pre-cise Estimation ofthe Optimum, Economic Opti-mum and Parameters for One FactorInverse Poly-nomial Model. UnpublisludPh.D thesis, UniversityofReading.(Diterima 29Jun 1989)PERTANlKA VOL. 15 NO.2. 1992 157

Kriterium Majmuk bagi Penganggaran Parameter dan Peramalan Sambutan

http://psasir.upm.edu.my/2999/1/Kriterium_Majmuk_bagi_Penganggaran_Parameter_dan_Peramalan_Sambutan.pdf

Kriterium Majmuk bagi Penganggaran Parameter dan Peramalan Sambutan

Abstract

Similar works

Full text

Available Versions

Universiti Putra Malaysia Institutional Repository