基于Logit-ISM模型的山東省馬鈴薯種植戶倉儲行為分析

任麗，王士海

（山東農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，山東 泰安 271018）

2015 年農(nóng)業(yè)部提出我國將啟動馬鈴薯主糧化戰(zhàn)略，使馬鈴薯成為繼水稻、小麥、玉米之后的另一種主糧，倉儲馬鈴薯對保障我國的糧食安全有重要作用；此外倉儲馬鈴薯具有一定的經(jīng)濟(jì)效益，對增加農(nóng)戶收入具有重要作用。因此研究農(nóng)戶倉儲馬鈴薯的行為無論從農(nóng)戶微觀還是國家宏觀來說都具有重要意義。

目前對馬鈴薯種植戶倉儲行為的研究相對較少。因此本文從農(nóng)戶儲糧行為影響因素進(jìn)行文獻(xiàn)概述。關(guān)于農(nóng)戶儲糧行為影響因素，學(xué)者普遍認(rèn)為農(nóng)戶個體特征及家庭特征、生產(chǎn)經(jīng)營特征、區(qū)域環(huán)境因素等均會影響農(nóng)戶的儲糧行為［1］。戶主年齡、收入水平、農(nóng)業(yè)勞動力占比、家庭常住人口等特征對農(nóng)戶糧食儲藏數(shù)量具有顯著影響［2，3］，其中收入對農(nóng)戶儲糧行為的影響具有爭議，有學(xué)者認(rèn)為收入對農(nóng)戶小麥儲備數(shù)量有顯著的負(fù)向影響，若家庭收入增加，則小麥儲備數(shù)量會減少［4］，另有觀點認(rèn)為農(nóng)民收入對農(nóng)戶的糧食儲備有不同的影響，對低收入階層是正相關(guān)，對高收入階層是負(fù)相關(guān)［5］。農(nóng)戶生產(chǎn)經(jīng)營特征會影響農(nóng)戶的儲糧行為，糧食產(chǎn)量是影響農(nóng)戶儲糧最重要的因素，糧食復(fù)種指數(shù)、農(nóng)戶家庭兼業(yè)類型都對農(nóng)戶糧食儲備產(chǎn)生影響［6-8］。農(nóng)戶糧食產(chǎn)出規(guī)模對農(nóng)戶儲糧行為具有顯著影響；糧食商品化率則與農(nóng)戶儲糧行為呈負(fù)相關(guān)［9］。農(nóng)戶的儲糧行為還會受到區(qū)域環(huán)境特征的影響，其中地勢類型、與最近縣城的距離［3］、氣候和耕作制度影響農(nóng)戶的儲糧決策［8，9］。對于山區(qū)農(nóng)戶，交通便利程度、儲糧設(shè)施優(yōu)劣對儲藏數(shù)量具有顯著正向影響［2］。糧食儲備效果、地區(qū)因素是影響農(nóng)戶糧食儲備行為的最關(guān)鍵因素［10］。

社會環(huán)境各種因素及農(nóng)戶個人的性格特征等心理因素也會影響農(nóng)戶的售糧儲糧行為［11］。名義利率、通貨膨脹率、糧食價格等因素對農(nóng)戶儲糧行為有顯著影響［12］。在農(nóng)戶收入水平較低、農(nóng)村經(jīng)濟(jì)欠發(fā)達(dá)的情況下，糧食市場預(yù)期供給、當(dāng)前市場價格均會對農(nóng)戶的儲糧行為產(chǎn)生影響［13，14］。農(nóng)戶的儲糧決策和行為存在明顯的模仿性和趨同性［15，16］，間接說明了不確定市場環(huán)境對農(nóng)戶儲糧行為產(chǎn)生了影響。

總體來看，學(xué)者對農(nóng)戶儲糧行為影響因素研究較為豐富，但對各因素間的關(guān)聯(lián)層次結(jié)構(gòu)研究較少。本文在已有研究成果的基礎(chǔ)上，基于山東省馬鈴薯種植戶的問卷調(diào)查資料，運用Logit、ISM模型分析馬鈴薯種植戶倉儲行為的影響因素。

1 數(shù)據(jù)來源與變量選擇

1.1 數(shù)據(jù)來源

山東是馬鈴薯種植大省，其農(nóng)業(yè)氣候條件優(yōu)越，馬鈴薯品質(zhì)和產(chǎn)量一直在全國處于領(lǐng)先地位。其中，滕州市馬鈴薯種植面積約占全省面積的1／4，產(chǎn)量占1／3，有“中國馬鈴薯之鄉(xiāng)”的稱號，此外膠州、肥城、鄒城、樂陵也是山東省馬鈴薯的主要產(chǎn)區(qū)。本研究所用數(shù)據(jù)由山東農(nóng)業(yè)大學(xué)學(xué)生于2019年6—9月通過實地調(diào)查獲得。采取隨機(jī)抽樣的方法，被調(diào)查對象為山東省五個馬鈴薯主產(chǎn)區(qū)種植戶，受訪者為戶主。本次調(diào)研共發(fā)放問卷230份，回收有效問卷200份，問卷有效率87%。其中滕州75份，占總體樣本的37.5%，肥城65份，占總體樣本的32.5%，膠州39份，占總體樣本的19.5%，樂陵和鄒城兩個地區(qū)馬鈴薯種植戶較少，共回收問卷21份，占總體樣本的10.5%。

1.2 變量選擇

基于已有研究結(jié)果，影響農(nóng)戶倉儲行為的因素可以大致總結(jié)為被調(diào)查農(nóng)戶的個體特征及家庭特征、生產(chǎn)經(jīng)營特征、區(qū)域環(huán)境特征、被調(diào)查農(nóng)戶素質(zhì)與態(tài)度、市場特征五個方面，本文從這五個方面選取15個變量來研究馬鈴薯種植戶倉儲行為的影響因素，變量含義及描述統(tǒng)計如表1所示。

被調(diào)查農(nóng)戶的個體與家庭特征主要包括農(nóng)戶的年齡、上學(xué)年限、收入來源、是否有借款、家有老人數(shù)量。出于投機(jī)獲利的想法，年齡越小的農(nóng)戶，就越會進(jìn)行倉儲，因為年齡大的農(nóng)戶思想趨于保守，與年齡小的農(nóng)戶相比，對待風(fēng)險的態(tài)度更為消極。但也有學(xué)者提出，由于年齡較大的農(nóng)戶經(jīng)歷過糧食短缺的特殊時期，為了保障糧食安全，一直都有存糧的習(xí)慣。農(nóng)戶的上學(xué)年限越高，其接受和學(xué)習(xí)新知識的能力越強(qiáng)，對信息收集和分析越充分，在是否選擇倉儲馬鈴薯時也更為理性。主要收入來源為務(wù)農(nóng)的農(nóng)戶經(jīng)濟(jì)收入單一，更可能選擇倉儲馬鈴薯，通過倉儲再賣出以期獲得更高的收入，且他們相較非務(wù)農(nóng)群體，更為重視經(jīng)營馬鈴薯所獲得的收益，不論是在種植、倉儲還是銷售方面都更為“專業(yè)”。有借款的農(nóng)戶，在馬鈴薯收獲后會更加迫切地將實物轉(zhuǎn)化為現(xiàn)金，他們對于風(fēng)險的態(tài)度接近“零容忍”，選擇倉儲馬鈴薯的可能性小。家里老人數(shù)量這一變量測量了農(nóng)戶家庭贍養(yǎng)壓力。通常家庭面對的贍養(yǎng)壓力越大，越會迫使經(jīng)營決策者采用多方式的經(jīng)營策略，去謀求更高的收益，因此，家里老人數(shù)量越多的農(nóng)戶越傾向于倉儲馬鈴薯。

被調(diào)查農(nóng)戶的家庭生產(chǎn)特征主要包括種植經(jīng)驗、馬鈴薯產(chǎn)量、受災(zāi)程度和之前是否倉儲。種植經(jīng)驗豐富的農(nóng)戶一般能根據(jù)作物的長勢對馬鈴薯產(chǎn)量以及市場價格的波動進(jìn)行一定的判斷，其豐富的經(jīng)驗可以為倉儲決策提供較為準(zhǔn)確的指導(dǎo)。馬鈴薯產(chǎn)量大的農(nóng)戶選擇性更多，可能會選擇先賣出部分馬鈴薯，來滿足下一期生產(chǎn)需要和日常開支，將剩下的馬鈴薯倉儲起來，待價而沽；而對于那些產(chǎn)量大銷售渠道有限的農(nóng)戶來說，為了減少馬鈴薯自然損耗，不得不將馬鈴薯倉儲起來逐漸賣出。若農(nóng)戶種植馬鈴薯的受災(zāi)程度嚴(yán)重，市場上馬鈴薯供不應(yīng)求，價格上漲，農(nóng)戶對市場價格滿意就會選擇將馬鈴薯直接賣出；若預(yù)期市場價格還會繼續(xù)上漲，農(nóng)戶也可能選擇將馬鈴薯倉儲起來，以期獲取更多收益。此外出于糧食安全考慮，受災(zāi)程度嚴(yán)重時農(nóng)戶可能會選擇將馬鈴薯倉儲起來，供自己家庭食用。之前倉儲過的農(nóng)戶為了繼續(xù)獲取較高的收益，很有可能選擇倉儲馬鈴薯，但同時之前倉儲損失較大或倉儲后價格下跌的農(nóng)戶可能不會選擇倉儲馬鈴薯；此外之前倉儲過馬鈴薯的農(nóng)戶，也可能存在行為慣性，會選擇再將馬鈴薯倉儲起來。

被調(diào)查農(nóng)戶的區(qū)域環(huán)境特征主要包括村莊地形。對于地形為山地或是丘陵的村莊，交通不便利，交易成本高，出于多種考慮，農(nóng)戶一般會選擇將馬鈴薯直接入庫，之后由專業(yè)戶或倉儲企業(yè)統(tǒng)一出售。

農(nóng)戶素質(zhì)與態(tài)度主要包括農(nóng)戶對倉儲技術(shù)的了解、當(dāng)季價格滿意度與家庭成員態(tài)度。農(nóng)戶對倉儲技術(shù)了解程度會影響其對于倉儲技術(shù)有用性和易用性的認(rèn)知，進(jìn)而影響倉儲行為。家庭成員的態(tài)度會影響農(nóng)戶的主觀規(guī)范，若家庭成員支持倉儲馬鈴薯，戶主就有很大的可能會去倉儲。農(nóng)戶對馬鈴薯當(dāng)季價格較為滿意，出于規(guī)避風(fēng)險的心理，害怕未來價格下跌，會選擇將馬鈴薯全部賣出去，而不倉儲馬鈴薯；另一方面農(nóng)戶出于投機(jī)心理，可能會覺得馬鈴薯價格會繼續(xù)上漲，而選擇將馬鈴薯倉儲起來，等到價格更高時再賣出。

市場特征主要為銷售渠道。若農(nóng)戶銷售馬鈴薯的渠道較少，為了減少自然損耗，只能選擇將馬鈴薯倉儲起來。

表1 變量含義及其統(tǒng)計描述

1.3 樣本描述分析

樣本基本特征描述如表2，樣本農(nóng)戶年齡主要集中于41～50歲，占比46.0%；被調(diào)查農(nóng)戶的受教育程度主要是6～9年，占到總體樣本的43.0%，說明農(nóng)戶的受教育程度普遍不高，但具有一定的學(xué)習(xí)能力，能夠處理信息并作出決策。從家庭勞動力人數(shù)來看，92.5%的家庭有2～4個勞動力。從收入來源來看，農(nóng)戶的家庭收入大都來源于務(wù)農(nóng)，且接近一半的農(nóng)戶家庭收入集中于5萬～10萬元。從種植經(jīng)驗來看，馬鈴薯種植戶種植年限多集中于10年內(nèi)，說明農(nóng)戶對馬鈴薯的種植已有一定程度的了解，種植經(jīng)驗較為豐富，對馬鈴薯的生產(chǎn)狀況較為熟悉。

表2 樣本基本特征描述

2 模型構(gòu)建

2.1 Logit模型

受訪農(nóng)戶是否倉儲馬鈴薯只有“是”、“否”兩種答案，這是典型的二值選擇問題，適合選擇Probit模型或 Logit模型。Probit模型和 Logit模型在本質(zhì)上無太大區(qū)別，二者差異之處在于它們趨于0或1的速率不同，Logit模型分布函數(shù)是正態(tài)分布函數(shù)，相比Probit模型會以更快的速率趨于0或1，所得出的結(jié)果也更符合實際情況。除此之外，由于Logit模型的累積分布函數(shù)有解析表達(dá)式，Probit模型沒有，故計算 Logit通常比Probit更為方便。因此，本研究采用 Logit模型來分析馬鈴薯種植戶倉儲行為的影響因素。農(nóng)戶倉儲馬鈴薯情況賦值為1，農(nóng)戶不倉儲馬鈴薯情況賦值為0。模型設(shè)定如下：

式中，p表示因變量為1的概率，即農(nóng)戶倉儲馬鈴薯的概率；1-p表示因變量為0的概率，即農(nóng)戶不倉儲馬鈴薯的概率；β0為常數(shù)項；i表示解釋變量的個數(shù)；βi為解釋變量的回歸系數(shù)；xi表示解釋變量；μ為隨機(jī)誤差項。

2.2 ISM模型

ISM是美國沃菲爾德教授開發(fā)的一種方法，目的是為分析復(fù)雜社會經(jīng)濟(jì)系統(tǒng)的結(jié)構(gòu)問題。其原理是通過確定影響系統(tǒng)的各種因素及其相互關(guān)系，利用圖論中的關(guān)聯(lián)矩陣原理和計算機(jī)技術(shù)，對不同因素及其相互關(guān)系的信息進(jìn)行處理，以明確因素間的關(guān)聯(lián)性和層次性，從而發(fā)現(xiàn)主要因素及其內(nèi)在聯(lián)系。馬鈴薯種植戶倉儲行為影響因素之間相互獨立又相互關(guān)聯(lián)。通過ISM模型分析農(nóng)戶倉儲行為影響因素之間的關(guān)系，分析其深層影響因素。本文用S0表示馬鈴薯種植戶倉儲行為，Si（i＝1，2，…，k）表示馬鈴薯種植戶倉儲行為的影響因素。鄰接矩陣R由公式（2）定義。

通過（3）式計算可達(dá)矩陣：

（3）式中 I為單位矩陣，2≤λ≤k，以上矩陣運算遵循的是布爾運算法則。最高層級因素確定如下：

（4）式中，P（Si）為可達(dá)集，包含可達(dá)矩陣中從 Si出發(fā)可以直接或間接到達(dá)的全部要素；Q（Si）為先行集，包含可達(dá)矩陣中可以直接或間接到達(dá)Si的全部要素。P（Si）和 Q（Si）的表達(dá)式如下：

（5）式中 mij和 mji都是可達(dá)矩陣的元素。通過（4）、（5）式得到頂層要素集后，在可達(dá)矩陣中刪掉L1層中元素，得到矩陣M′，再對矩陣M′進(jìn)行上述的運算，可得位于第二層要素的集合L2。依照此種方法，再得到位于第三層直至最后一層的要素集合。依據(jù)層次關(guān)系，將各層要素用箭頭連接起來，即可得到馬鈴薯種植戶倉儲行為影響因素間的層次結(jié)構(gòu)。

3 結(jié)果與分析

3.1 模型估計結(jié)果

3.1.1 農(nóng)戶馬鈴薯倉儲行為的影響因素 本文運用Stata 14軟件對山東省馬鈴薯種植戶倉儲行為進(jìn)行Logit模型處理，實證結(jié)果如表3所示?？芍匦?、馬鈴薯產(chǎn)量、受災(zāi)程度、之前是否倉儲、倉儲技術(shù)了解、銷售渠道、家庭成員態(tài)度均對農(nóng)戶的倉儲行為有顯著的影響。

表3 馬鈴薯種植戶倉儲行為影響因素實證結(jié)果

3.1.2 農(nóng)戶倉儲行為影響因素的解釋性結(jié)構(gòu)基于 Logit模型的實證結(jié)果，用 S1、S2、S3、S4、S5、S6、S7表示受災(zāi)程度、地形、馬鈴薯產(chǎn)量、之前是否倉儲、倉儲技術(shù)了解、銷售渠道、家庭成員態(tài)度，“A”表示列因素對行因素有影響，“V”表示行因素對列因素有影響，各影響因素間的邏輯關(guān)系如圖1所示。

圖1 影響因素間的邏輯關(guān)系

根據(jù)（2）式，可得鄰接矩陣 R（略），再根據(jù)（3）式并利用Matlab軟件得出可達(dá)矩陣M，如（6）式。根據(jù)可達(dá)矩陣，通過（4）式和（5）式，首先得到頂層要素集，即L1＝（S0）。其次依次確定其他層因素，依次得到第二層、第三層、第四層、第五層要素，分別是 L2＝（S4，S5，S7），L3＝（S3），L4＝（S1，S6），L5＝（S2），其中，第二層作為表層因素，第三、四層作為間接因素，第五層作為最深層次因素。根據(jù) L1，L2，L3，L4，L5將可達(dá)矩陣 M的行與列重新排序得到矩陣B，如（7）式。

根據(jù)各因素間的邏輯關(guān)系以及上文得出的層次結(jié)構(gòu)，得出馬鈴薯種植戶倉儲行為影響因素間關(guān)聯(lián)層次結(jié)構(gòu)圖，如圖2所示。

圖2 馬鈴薯種植戶倉儲行為影響因素間關(guān)聯(lián)層次結(jié)構(gòu)示意圖

3.2 結(jié)果分析

由圖2可見，在農(nóng)戶馬鈴薯倉儲行為影響因素中，農(nóng)戶之前是否倉儲、農(nóng)戶對倉儲技術(shù)的了解、農(nóng)戶家庭成員的態(tài)度這三個因素是影響農(nóng)戶倉儲行為的表層因素，受災(zāi)程度、農(nóng)戶銷售渠道、馬鈴薯產(chǎn)量是影響馬鈴薯種植戶倉儲行為的間接因素，地形是影響農(nóng)戶倉儲行為的最深層次因素。

農(nóng)戶對倉儲技術(shù)的了解通過1%的顯著性水平檢驗，且方向為正，這與之前的假設(shè)一致。農(nóng)戶對倉儲技術(shù)了解程度越高，越會影響其對于倉儲技術(shù)有用性和易用性的認(rèn)知，最終影響其倉儲行為意向。因此農(nóng)戶對倉儲技術(shù)越了解越可能做出倉儲行為。

之前是否倉儲也通過1%的顯著性水平檢驗，且方向為正?？赡艿慕忉屖?，之前倉儲的農(nóng)戶存在行為慣性，會選擇繼續(xù)倉儲馬鈴薯。此外之前倉儲過馬鈴薯的農(nóng)戶對馬鈴薯的倉儲過程更加了解，可以降低倉儲損失程度，會選擇繼續(xù)倉儲馬鈴薯。

銷售渠道通過5%的顯著性水平檢驗，且符號為負(fù)。說明若農(nóng)戶銷售馬鈴薯的渠道很少，馬鈴薯賣不出，為減少馬鈴薯的自然損耗，只能選擇將馬鈴薯倉儲起來逐漸賣出。

農(nóng)戶家庭成員的態(tài)度通過5%的顯著性水平檢驗，且符號為正，說明家庭成員的態(tài)度會影響農(nóng)戶的主觀規(guī)范，若家庭成員對倉儲行為是支持的態(tài)度，農(nóng)戶很有可能會倉儲馬鈴薯。

馬鈴薯產(chǎn)量通過5%的顯著性水平檢驗，符號也為正，說明產(chǎn)量大的農(nóng)戶更傾向于倉儲馬鈴薯。產(chǎn)量大的農(nóng)戶有更多的選擇，家庭決策后更有可能出于理性考慮先賣出部分馬鈴薯，來滿足下一期生產(chǎn)需要和日常開支，將剩下的馬鈴薯倉儲起來，待價而沽，以獲得更多收益；此外對于那些產(chǎn)量大、銷售渠道有限的農(nóng)戶來說，馬鈴薯賣不出去時，為了減少馬鈴薯自然損耗，也只能將馬鈴薯倉儲起來逐漸賣出。

受災(zāi)程度通過10%的顯著性水平檢驗，且符號為正，與前面的假設(shè)一致。若農(nóng)戶種植馬鈴薯的受災(zāi)程度嚴(yán)重，市場上馬鈴薯價格多會上漲，農(nóng)戶與其家庭成員出于對市場價格的預(yù)期，會選擇將馬鈴薯倉儲起來，以期獲取更多收益。且若受災(zāi)程度嚴(yán)重，農(nóng)戶可能會選擇將馬鈴薯倉儲起來，供自己家庭食用。

地形通過1%的顯著性水平檢驗，且符號為正，這與前面的假設(shè)一致。說明在山區(qū)的農(nóng)戶更傾向于倉儲馬鈴薯?？赡艿慕忉屖巧絽^(qū)的農(nóng)戶，信息閉塞，馬鈴薯銷售渠道較少，此外山區(qū)相較于平原地區(qū)受災(zāi)程度可能更嚴(yán)重，因此農(nóng)戶更有可能做出倉儲行為。

綜上所述，以上7個影響因素既獨立發(fā)揮作用，又相互關(guān)聯(lián)，形成完整的農(nóng)戶倉儲行為影響因素系統(tǒng)。作為深層影響因素，地形直接影響農(nóng)戶銷售渠道、受災(zāi)程度、馬鈴薯產(chǎn)量，沿著“地形”→“農(nóng)戶銷售渠道、受災(zāi)程度、馬鈴薯產(chǎn)量”→“農(nóng)戶之前是否倉儲、農(nóng)戶對倉儲技術(shù)的了解、農(nóng)戶家庭成員的態(tài)度”→“農(nóng)戶倉儲行為”的正向傳導(dǎo)關(guān)系，影響馬鈴薯種植戶的倉儲行為。

4 結(jié)論與建議

研究結(jié)果表明：地形、馬鈴薯產(chǎn)量、受災(zāi)程度、之前是否倉儲、倉儲技術(shù)了解、銷售渠道、家庭成員態(tài)度均顯著影響馬鈴薯種植戶倉儲行為，除銷售渠道外，其他因素影響均為正向。此外，通過ISM模型分析，得出地形是影響農(nóng)戶倉儲行為的最深層次因素且正向影響顯著，表明從地形出發(fā)，基于平原、山區(qū)的不同情況來有效增加農(nóng)戶的倉儲行為是提高我國糧食安全的關(guān)鍵。

因此山區(qū)應(yīng)加快倉儲點建設(shè)，改善倉儲條件，提高農(nóng)戶倉儲的質(zhì)量；對于條件便利、農(nóng)戶倉儲行為少的平原地區(qū)，應(yīng)加強(qiáng)宣傳教育，增強(qiáng)農(nóng)戶糧食安全意識，促進(jìn)農(nóng)戶儲糧行為。農(nóng)戶對倉儲技術(shù)的了解對其倉儲行為有顯著的正向影響，應(yīng)積極開展馬鈴薯倉儲技術(shù)的相關(guān)培訓(xùn)，推廣科學(xué)倉儲技術(shù)和設(shè)備，加強(qiáng)對農(nóng)戶的技術(shù)指導(dǎo)，增強(qiáng)農(nóng)戶倉儲積極性。