區(qū)域玉米生產(chǎn)成本影響因素的實(shí)證分析*

盧德成

(中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)信息研究所,北京 100081)

0 引言

在中國(guó)，玉米不僅是“三大主糧”之一，同時(shí)也是畜牧業(yè)主要飼料來(lái)源和重要的工業(yè)原料。2017年，中國(guó)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革進(jìn)入深化之年，玉米供給側(cè)改革取得顯著成效，階段性供給過(guò)剩的局面得到抑制，但同時(shí)也意味著傳統(tǒng)玉米優(yōu)勢(shì)產(chǎn)區(qū)會(huì)因此而承擔(dān)更多的糧食供給風(fēng)險(xiǎn)，而當(dāng)前主產(chǎn)省玉米“高成本”的局面不僅不利于玉米生產(chǎn)的穩(wěn)定發(fā)展，而且還會(huì)對(duì)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革的進(jìn)程造成一定影響。通過(guò)研究發(fā)現(xiàn)有效降低玉米生產(chǎn)成本的方法具有重要的現(xiàn)實(shí)意義。研究以玉米為研究對(duì)象，旨在通過(guò)全國(guó)各區(qū)域玉米生產(chǎn)成本影響因素的實(shí)證分析，發(fā)現(xiàn)影響各區(qū)域玉米生產(chǎn)成本變化的主要因素，并據(jù)此為優(yōu)化區(qū)域玉米生產(chǎn)成本提供參考建議。

1 玉米區(qū)域劃分及成本測(cè)算

1.1 玉米生產(chǎn)區(qū)域劃分

國(guó)內(nèi)專家、學(xué)者在玉米生產(chǎn)區(qū)域劃分方面進(jìn)行了大量且深入的研究。根據(jù)不同省份產(chǎn)量，閆麗珍等[1]將14個(gè)主產(chǎn)省劃為北段、中段和南段3個(gè)區(qū)域； 根據(jù)不同產(chǎn)區(qū)的資源、地理、氣候等條件，岳德榮[2]、郭慶海[3]將種植區(qū)劃分為北方春播玉米區(qū)、黃淮海夏播玉米區(qū)、西南山地玉米區(qū)、南方丘陵玉米區(qū)、西北灌溉玉米區(qū)和青藏高原玉米區(qū)6個(gè)區(qū)域； 柴斌鋒[4]據(jù)此將13個(gè)玉米主產(chǎn)省劃分為東北春玉米、黃淮海夏玉米區(qū)和西南山地丘陵玉米區(qū)3個(gè)區(qū)域。

結(jié)合產(chǎn)量和資源條件兩方面，選取玉米產(chǎn)量較高的20個(gè)省份，按資源、地理和氣候條件的不同劃分為5個(gè)區(qū)域：東北春玉米區(qū)(簡(jiǎn)稱東北區(qū))，包含黑、吉、遼、蒙4省； 黃淮海夏玉米區(qū)(簡(jiǎn)稱黃淮海區(qū))：包含晉、魯、豫、冀4省； 西北灌溉玉米區(qū)(簡(jiǎn)稱西北區(qū))：包含陜、甘、寧、新4省； 南方丘陵玉米區(qū)(簡(jiǎn)稱南方區(qū))：包含皖、江、鄂、桂4省； 西南山地玉米區(qū)(簡(jiǎn)稱西南區(qū))，包含云、黔、川、渝4省。

1.2 區(qū)域成本測(cè)算方法

由于該研究所分區(qū)域內(nèi)的省份數(shù)均為4個(gè)，因此設(shè)某產(chǎn)區(qū)的每667m2生產(chǎn)成本為C，區(qū)域生產(chǎn)總投入為C0，各省份的每667m2生產(chǎn)成本分別為C1、C2、C3、C4； 且設(shè)區(qū)域總播種面積為A0，對(duì)應(yīng)的各省播種面積為A1、A2、A3、A4。則從生產(chǎn)成本的定義出發(fā)，該區(qū)域的平均生產(chǎn)成本等于總投入與總播種面積的比值。

(1)

由式(1)可見(jiàn)，區(qū)域生產(chǎn)成本及各構(gòu)成項(xiàng)目的數(shù)據(jù)能通過(guò)各省、市、自治區(qū)的玉米播種面積加權(quán)后計(jì)算得出。

2 區(qū)域玉米生產(chǎn)成本的演變特征

根據(jù)測(cè)算方法，計(jì)算得出1978～2016年各年度、各區(qū)域的玉米生產(chǎn)成本(表1)。

表1 1978～2016年區(qū)域玉米生產(chǎn)成本數(shù)量及增長(zhǎng)率變化情況 元/667m2、%

區(qū)域玉米生產(chǎn)成本總體呈上升趨勢(shì)，并表現(xiàn)出東部低于西部、北方低于南方的特點(diǎn)。各區(qū)域每667m2玉米生產(chǎn)成本由1978年的41.84元、47.36元、40.71元、46.86元和57.13元增長(zhǎng)至2016年的662.09元、784.96元、824.68元、1 189.84元和1 251.12元，除個(gè)別年份外，成本數(shù)量基本符合西南區(qū)>西北區(qū)>全國(guó)平均>南方區(qū)>黃淮海區(qū)>東北區(qū)的態(tài)勢(shì)。成本增長(zhǎng)呈現(xiàn)西南區(qū)>南方區(qū)>西北區(qū)>全國(guó)平均>黃淮海區(qū)>東北區(qū)的態(tài)勢(shì)， 1978～2016年年均增幅分別為8.5%、8.2%、8.2%、7.7%和7.5%。

近年來(lái)，不同區(qū)域玉米生產(chǎn)成本的變化趨勢(shì)發(fā)生新的變化。以2004年為節(jié)點(diǎn)， 2004～2016年，全國(guó)玉米生產(chǎn)成本年均增長(zhǎng)率達(dá)到8.4%，高于1978～2016年的7.7%，不同區(qū)域玉米生產(chǎn)成本年均增長(zhǎng)率有所不同。其中，(1)黃淮海區(qū)成本年均增長(zhǎng)率為9.3%，不僅增幅為各區(qū)域中的最高，而且增速也明顯快于1978～2016年的7.67%； (2)西北區(qū)和西南區(qū)的成本年均增速有所上升，年均增速較1978～2016年分別上升11.9%、9.8%； (3)東北區(qū)和南方區(qū)的成本年均增長(zhǎng)率略低于1978～2016年的成本年均增速。

近幾年，中國(guó)經(jīng)濟(jì)和社會(huì)發(fā)展逐漸進(jìn)入新常態(tài)，隨著農(nóng)業(yè)生產(chǎn)方式的逐漸改變，得益于結(jié)構(gòu)性改革成效， 2013年以來(lái)，全國(guó)各區(qū)域玉米的生產(chǎn)成本增速迅速放緩。2013～2016年，東北區(qū)和南方區(qū)玉米生產(chǎn)成本實(shí)現(xiàn)負(fù)增長(zhǎng)，分別為-0.9%、-0.8%； 西南區(qū)、黃淮海區(qū)和西北區(qū)雖保持增長(zhǎng)，但增速已經(jīng)大幅降低，黃淮海區(qū)2016年也實(shí)現(xiàn)成本負(fù)增長(zhǎng)，當(dāng)年降幅3.7%，西北區(qū)和西南區(qū)則仍處于微增狀態(tài)。

3 區(qū)域玉米生產(chǎn)成本的影響因素分析

3.1 變量選取及數(shù)據(jù)來(lái)源

影響區(qū)域間玉米生產(chǎn)成本變動(dòng)的因素分為內(nèi)部因素和外部因素兩個(gè)方面。其中，內(nèi)部因素是指直接作用于玉米生產(chǎn)過(guò)程的物質(zhì)資料投入，因使用或租用勞動(dòng)工具進(jìn)行作業(yè)而產(chǎn)生的投入(簡(jiǎn)稱租賃作業(yè)投入)和勞動(dòng)力投入。物質(zhì)資料投入使用最新農(nóng)產(chǎn)品成本核算體系中的種子費(fèi)、化肥費(fèi)、農(nóng)藥費(fèi)和農(nóng)家肥費(fèi)加總來(lái)代指； 生產(chǎn)工具投入使用同體系中的租賃作業(yè)費(fèi)來(lái)代指； 勞動(dòng)力投入則是使用人工成本來(lái)代指。外部因素則主要是影響玉米綜合生產(chǎn)能力水平的規(guī)模化、生產(chǎn)技術(shù)水平和資源型產(chǎn)品等因素。生產(chǎn)規(guī)模為其中的主要因素，以玉米播種面積來(lái)代指。

該研究所使用的物質(zhì)資料投入、生產(chǎn)工具投入及勞動(dòng)力投入等原始成本數(shù)據(jù)均來(lái)自歷年《中國(guó)農(nóng)產(chǎn)品成本收益資料匯編》； 播種面積數(shù)據(jù)來(lái)源自歷年《中國(guó)統(tǒng)計(jì)年鑒》和《中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒》。

3.2 模型構(gòu)建

設(shè)玉米生產(chǎn)成本為物質(zhì)資料類投入、勞動(dòng)工具類投入、勞動(dòng)力投入、規(guī)模化程度相關(guān)的函數(shù)，則建立雙對(duì)數(shù)線性回歸模型：

lnY=β1lnXM+β2lnXT+β3lnXL+β4lnXA+μ

(2)

式(2)中，Y為“生產(chǎn)成本”；XM為“物質(zhì)資料類投入”；XT為“租賃作業(yè)類投入”；XL為“勞動(dòng)力投入”；XA為“規(guī)模化程度”。

對(duì)數(shù)線性模型的優(yōu)勢(shì)在于可以在一定程度上消除原有數(shù)據(jù)中可能存在的多重共線性和自相關(guān)性，使分析結(jié)果更準(zhǔn)確。并且在使用取自然對(duì)數(shù)的方式對(duì)原始數(shù)據(jù)進(jìn)行相應(yīng)處理后，該模型所反映的數(shù)學(xué)概念即可變?yōu)槊繂挝蛔宰兞康淖儎?dòng)可以引起的因變量的變動(dòng)幅度，其變量系數(shù)β1、β2、β3和β4即為彈性系數(shù)。因此，可將對(duì)數(shù)后的指標(biāo)lnY、lnXM、lnXT、lnXL和lnXA直接作為因變量和自變量，通過(guò)對(duì)其彈性系數(shù)的分析，反映因各影響因素的變化而產(chǎn)生的生產(chǎn)成本變化情況。

3.3 模型檢驗(yàn)

由于原始面板數(shù)據(jù)T>N，可能產(chǎn)生偽回歸。為避免分析誤差，須先對(duì)數(shù)據(jù)做ADF單位根檢驗(yàn)。使用Eviews8.0對(duì)原始數(shù)據(jù)取自然對(duì)數(shù)后進(jìn)行單位根檢驗(yàn)，得出各區(qū)域的變量數(shù)據(jù)在水平(Level)上多數(shù)為非平穩(wěn)序列，但所有變量數(shù)據(jù)的一階差分(1st Difference)序列均為平穩(wěn)序列。一階差分序列的ADF檢驗(yàn)結(jié)果如表2。

表2 原始數(shù)據(jù)一階差分ADF檢驗(yàn)分析結(jié)果

變量區(qū)域ADF統(tǒng)計(jì)值1%臨界參考值5%臨界參考值10%臨界參考值結(jié)論Y東北區(qū)-35422-36329-29484-261285%平穩(wěn)XM-47443-36210-29434-261031%平穩(wěn)XT-39373-36793-29678-262301%平穩(wěn)XL-59651-36210-29434-261031%平穩(wěn)XA-48606-36210-29434-261031%平穩(wěn)Y黃淮海區(qū)-36665-36210-29434-261031%平穩(wěn)XM-61972-36210-29434-261031%平穩(wěn)XT-62762-36268-29458-261151%平穩(wěn)XL-46200-42733-35578-321241%平穩(wěn)XA-42335-42268-35366-320031%平穩(wěn)Y西北區(qū)-43173-42268-35366-320031%平穩(wěn)XM-78075-36210-29434-261031%平穩(wěn)XT-79962-36210-29434-261031%平穩(wěn)XL-64485-36210-29434-261031%平穩(wěn)XA-56776-36210-29434-261031%平穩(wěn)Y南方區(qū)-41482-36210-29434-261031%平穩(wěn)XM-91897-36210-29434-261031%平穩(wěn)XT-71339-36210-29434-261031%平穩(wěn)XL-65931-36210-29434-261031%平穩(wěn)XA-43563-36210-29434-261031%平穩(wěn)Y西南區(qū)-42874-36210-29434-261031%平穩(wěn)XM-68408-42268-35366-320031%平穩(wěn)XT-70696-36210-29434-261031%平穩(wěn)XL-66883-36210-29434-261031%平穩(wěn)XA-49026-36210-29434-261031%平穩(wěn)

鑒于各區(qū)域變量序列在一階差分上平穩(wěn)，采用Engle-Granger兩步法進(jìn)行協(xié)整檢驗(yàn)，以判斷變量間是否具有穩(wěn)定關(guān)系。(1)使用Eviews8.0軟件，運(yùn)用普通最小二乘法(Ordinary Least Square，OLS)對(duì)原始數(shù)據(jù)進(jìn)行線性回歸，分別計(jì)算各區(qū)域的非均衡誤差et1(東北區(qū))、et2(黃淮海區(qū))、et3(西北區(qū))、et4(南方區(qū))和et5(西北區(qū))； (2)分別對(duì)不同區(qū)域的非均衡誤差進(jìn)行ADF單位根檢驗(yàn)，檢驗(yàn)結(jié)果顯示(表3)，除東北區(qū)非均衡誤差項(xiàng)在一階差分上平穩(wěn)外，其他各區(qū)域非均衡誤差項(xiàng)都在水平上平穩(wěn)，通過(guò)協(xié)整檢驗(yàn)，不同變量間存在穩(wěn)定關(guān)系，可以進(jìn)行回歸分析。

表3 非均衡誤差的ADF檢驗(yàn)結(jié)果

條件ADF統(tǒng)計(jì)值1%臨界參考值5%臨界參考值10%臨界參考值結(jié)論et1Level-21930-36268-29458-26115不平穩(wěn)et11stdifference-72784-36268-29458-261151%平穩(wěn)et2Level-46701-36156-29411-260911%平穩(wěn)et3Level-60506-36156-29411-260911%平穩(wěn)et4Level-40209-36156-29411-260911%平穩(wěn)et5Level-44372-36156-29411-260911%平穩(wěn)

3.4 分析結(jié)論

通過(guò)回歸得出不同區(qū)域變量相關(guān)結(jié)果如表4。其中各變量系數(shù)即為彈性系數(shù)，表示該區(qū)域生產(chǎn)成本受變量變化的變化程度。

表4 回歸分析結(jié)果

系數(shù)東北區(qū)黃淮海區(qū)南方區(qū)西北區(qū)西南區(qū)常數(shù)項(xiàng)μ38635??40272??1387129107??19674?lnXMβ103584??03116??03044??01803??02504??lnXTβ202328??01591??00966?02626??00326lnXLβ304526??05427??05702??05512??07012??lnXAβ4-03145??-03312??-00337-02242??-01339AdjustedR20998509987099840999009993F-statistic58617457685629529913386650841197185 注：表4中采用系數(shù)標(biāo)注上標(biāo)的方式表示顯著性水平，其中，??表示1%的顯著性水平；?表示5%的顯著性水平；無(wú)上標(biāo)標(biāo)注的表示不具顯著性

從分析結(jié)果來(lái)看：物質(zhì)資料投入、租賃作業(yè)投入和勞動(dòng)力投入與生產(chǎn)成本呈正相關(guān)關(guān)系，生產(chǎn)規(guī)模則與之呈現(xiàn)負(fù)相關(guān)。其中，(1)勞動(dòng)力投入對(duì)全國(guó)玉米生產(chǎn)成本的影響最大，除東北區(qū)外，區(qū)域成本對(duì)勞動(dòng)力投入的彈性系數(shù)均超過(guò)0.5，西南區(qū)甚至超過(guò)了0.7。(2)物質(zhì)資料投入對(duì)各區(qū)域生產(chǎn)成本均產(chǎn)生顯著影響，且對(duì)東部區(qū)域的影響程度高于西部區(qū)域，東北區(qū)、黃淮海區(qū)和南方區(qū)玉米生產(chǎn)成本對(duì)物質(zhì)資料投入的彈性系數(shù)在0.3以上，西南區(qū)和西北區(qū)則低于這一數(shù)字，分別為0.25和0.18。(3)租賃作業(yè)投入主要對(duì)北方的東北區(qū)、黃淮海區(qū)和西北區(qū)產(chǎn)生影響，其中西北區(qū)生產(chǎn)成本對(duì)租賃作業(yè)投入的彈性系數(shù)最高，為0.26，其次是東北區(qū)，黃淮海區(qū)最低，彈性系數(shù)為0.16； 對(duì)南方的南方區(qū)、西南區(qū)兩區(qū)而言，租賃作業(yè)投入對(duì)生產(chǎn)成本的影響并不顯著，只有南方區(qū)玉米成本與租賃作業(yè)投入在5%的顯著水平下存在相關(guān)性，彈性系數(shù)僅為0.10，影響程度較低，西南區(qū)則不具備顯著的相關(guān)性。(4)生產(chǎn)規(guī)模與租賃作業(yè)投入類似，也是主要對(duì)北方三區(qū)生產(chǎn)成本產(chǎn)生顯著影響，對(duì)南方兩區(qū)則不具有顯著相關(guān)性。其中，黃淮海區(qū)生產(chǎn)成本受生產(chǎn)規(guī)模的影響程度最高，其次是東北區(qū)，彈性系數(shù)的絕對(duì)值均大于0.3，西北區(qū)相對(duì)較低，彈性系數(shù)絕對(duì)值為0.22。

4 結(jié)論及建議

4.1 結(jié)論

該文將全國(guó)20個(gè)玉米主產(chǎn)省劃分為5個(gè)不同區(qū)域，運(yùn)用雙對(duì)數(shù)線性模型對(duì)不同區(qū)域玉米生產(chǎn)成本影響因素進(jìn)行了實(shí)證分析。影響因素分析結(jié)果表明：由于地理?xiàng)l件、生產(chǎn)方式等的不同，不同變量對(duì)不同區(qū)域玉米生產(chǎn)成本的影響程度不一，勞動(dòng)力投入是最主要的影響因素，其次是物質(zhì)資料投入，租賃作業(yè)投入和生產(chǎn)規(guī)模主要對(duì)北方產(chǎn)區(qū)產(chǎn)生影響，對(duì)南方產(chǎn)區(qū)影響并不顯著。不同區(qū)域玉米生產(chǎn)成本對(duì)不同影響因素變動(dòng)的彈性有所不同。

4.2 相關(guān)建議

近年來(lái)，中國(guó)玉米生產(chǎn)出現(xiàn)階段性供過(guò)于求的情況，國(guó)家玉米產(chǎn)業(yè)政策發(fā)生較大轉(zhuǎn)變。2015年11月，農(nóng)業(yè)部發(fā)布《農(nóng)業(yè)部關(guān)于“鐮刀彎”地區(qū)玉米結(jié)構(gòu)調(diào)整的指導(dǎo)意見(jiàn)》，對(duì)玉米種植面積進(jìn)行結(jié)構(gòu)性調(diào)整。2016年起，國(guó)家逐漸退出玉米的臨儲(chǔ)制度，還價(jià)于市，在東北和內(nèi)蒙古率先開(kāi)展“市場(chǎng)化收購(gòu)+補(bǔ)貼”的新政。基于國(guó)家政策的新變化，結(jié)合影響因素分析結(jié)論，建議：

(1)“藏糧于地”，通過(guò)推行適度規(guī)模經(jīng)營(yíng)的生產(chǎn)方式，提高玉米生產(chǎn)集約化、規(guī)模化生產(chǎn)水平，降低人工成本。勞動(dòng)力投入是影響全國(guó)各區(qū)域玉米生產(chǎn)成本變化的最主要因素，提高玉米生產(chǎn)經(jīng)營(yíng)的規(guī)模化水平，可以充分釋放規(guī)模經(jīng)濟(jì)效果，減少玉米生產(chǎn)的單位面積用工數(shù)量，達(dá)到抑制勞動(dòng)力投入增長(zhǎng)的目的。對(duì)于不同區(qū)域，建議在東北區(qū)、黃淮海區(qū)和西北區(qū)內(nèi)有條件的地區(qū)因地制宜推行大規(guī)模種植模式，以充分釋放生產(chǎn)規(guī)模經(jīng)濟(jì)效益，條件尚不足的地方推行適度規(guī)模經(jīng)營(yíng)的生產(chǎn)方式； 而在南方各區(qū)則主要推廣適度規(guī)模經(jīng)營(yíng)的生產(chǎn)方式，以起到抑制勞動(dòng)力投入上漲的作用。

(2)“藏糧于技”，通過(guò)推廣使用良種、良技，減少種子損耗、提高肥料利用率，降低物質(zhì)資料投入品支出。種子、化肥等物質(zhì)資料投入品的使用直接影響玉米生產(chǎn)成本，提高投入品使用效率，抑制投入品數(shù)量和交易成本，可以有效降低成本。建議政府創(chuàng)造條件，主要在東北區(qū)、黃淮海區(qū)、南方區(qū)和西南區(qū)，合理引導(dǎo)種糧主體采取良種栽培、測(cè)土配方施肥等方式投入物質(zhì)資料，以提高種子和化肥的利用率，降低損耗，從而降低種子、化肥等投入品的投入數(shù)量； 同時(shí)，通過(guò)制定出臺(tái)交通優(yōu)惠政策，降低種子和化肥的運(yùn)輸費(fèi)用，盡量減少交易成本，從而達(dá)到減少投入品支出的目的。

(3)在部分“鐮刀彎”玉米生產(chǎn)比較優(yōu)勢(shì)較低的種植區(qū)域，適當(dāng)退出玉米種植，即降低區(qū)域總體成本，也可提高當(dāng)?shù)剞r(nóng)民收益。東北冷涼區(qū)、北方農(nóng)牧交錯(cuò)區(qū)的黑龍江、內(nèi)蒙古，西北風(fēng)沙干旱區(qū)甘肅，西南石漠化區(qū)的云南、貴州的部分地區(qū)氣候和地理?xiàng)l件不利于玉米種植，生產(chǎn)成本反而較高，適當(dāng)調(diào)減退出種植，改種適宜環(huán)境的其他作物，不僅利于降低北方區(qū)整體玉米生產(chǎn)成本，也利于增加當(dāng)?shù)剞r(nóng)民的經(jīng)濟(jì)收益。

[1] 閆麗珍，成升魁，劉愛(ài)民，等.中國(guó)玉米生產(chǎn)成本收益的區(qū)域分布規(guī)律研究.農(nóng)業(yè)技術(shù)經(jīng)濟(jì)， 2003，(6)： 27～34

[2] 岳德榮. 中國(guó)玉米品質(zhì)區(qū)劃及產(chǎn)業(yè)布局.北京：中國(guó)農(nóng)業(yè)出版社， 2004

[3] 郭慶海. 中國(guó)玉米主產(chǎn)區(qū)的演變與發(fā)展.玉米科學(xué)， 2010, 18(1): 139～145

[4] 柴斌鋒， 鄭少鋒，李哲.中國(guó)玉米生產(chǎn)成本地區(qū)差異的實(shí)證分析.商業(yè)研究, 2007，(12)： 110～115

[5] 鄭少鋒， 邵建成.主要糧食作物生產(chǎn)成本影響因素分析.中國(guó)農(nóng)學(xué)通報(bào)， 2003, 19(3): 115～119

[6] 范成方， 史建民.山東省糧食種植成本影響因素的實(shí)證分析——以玉米、小麥為例.中國(guó)農(nóng)業(yè)資源與區(qū)劃， 2014, 35(2): 67～74

[7] 李國(guó)祥. 玉米價(jià)格與生產(chǎn)者收益關(guān)系的研究——基于我國(guó)玉米收儲(chǔ)制度改革背景下的思考.價(jià)格理論與實(shí)踐, 2016，(4)： 53～58

[8] 馬曉河. 中國(guó)農(nóng)業(yè)收益與生產(chǎn)成本變動(dòng)的結(jié)構(gòu)分析.中國(guó)農(nóng)村經(jīng)濟(jì).2011，(5)： 4～11， 56

[9] 王彩峰， 史建民.山東省蘋(píng)果種植收益影響因素的實(shí)證分析.中國(guó)農(nóng)業(yè)資源與區(qū)劃， 2017, 38(10): 215～221

[10]許存興， 魏建中，張芙蓉.基于多元回歸的三種糧食生產(chǎn)成本分析.山東農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2016, 47(5): 779～784