吳紅梅，王 蕊，殷沐陽(yáng)

（北京林業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，中國(guó) 北京 100083）

一、引言與文獻(xiàn)綜述

隨著中國(guó)林產(chǎn)品貿(mào)易的快速發(fā)展，林產(chǎn)品貿(mào)易水平不斷提升，2019年中國(guó)林產(chǎn)品貿(mào)易總額超1 106億美元，其中林產(chǎn)品出口貿(mào)易額超591 億美元，均居世界首位，林產(chǎn)品出口貿(mào)易是我國(guó)出口貿(mào)易的重要組成部分。在逆全球化和貿(mào)易保護(hù)主義盛行的背景下，從全球林產(chǎn)品出口貿(mào)易格局及影響因素角度研究中國(guó)林產(chǎn)品出口具有重要的理論和現(xiàn)實(shí)意義。

復(fù)雜網(wǎng)絡(luò)理論起源于圖論，作為一種能夠研究大量錯(cuò)綜復(fù)雜系統(tǒng)的有效方法，近年來(lái)已廣泛應(yīng)用于各種領(lǐng)域，如物理學(xué)［1］、生物學(xué)［2］、交通學(xué)［3］、投資學(xué)［4］和人口學(xué)［5］等。國(guó)際貿(mào)易涉及國(guó)家之間復(fù)雜的關(guān)系，復(fù)雜網(wǎng)絡(luò)理論也被國(guó)內(nèi)外學(xué)者逐漸用于國(guó)際貿(mào)易研究中，涉及多種行業(yè)和產(chǎn)品，主要集中于以下兩方面：一是利用復(fù)雜網(wǎng)絡(luò)分析方法分析貿(mào)易特征。有學(xué)者從整體貿(mào)易角度出發(fā)，如Fagiolo等［6］從網(wǎng)絡(luò)密度、中心性等指標(biāo)分析1981—2000 年世界159 個(gè)樣本國(guó)家貿(mào)易網(wǎng)絡(luò)演變，結(jié)果表明該貿(mào)易網(wǎng)絡(luò)為“核心-邊緣”結(jié)構(gòu)，發(fā)達(dá)國(guó)家貿(mào)易聯(lián)系更加緊密和集中，各年份網(wǎng)絡(luò)屬性穩(wěn)定；Ward等［7］通過(guò)構(gòu)建1990—2008 年世界貿(mào)易網(wǎng)絡(luò)，分析其拓?fù)浣Y(jié)構(gòu)，并探究網(wǎng)絡(luò)中節(jié)點(diǎn)之間的依賴(lài)關(guān)系，研究發(fā)現(xiàn)網(wǎng)絡(luò)中的高度依賴(lài)關(guān)系能夠顯著促進(jìn)貿(mào)易規(guī)模的增長(zhǎng)。還有部分學(xué)者從具體產(chǎn)品貿(mào)易角度出發(fā)，應(yīng)用復(fù)雜網(wǎng)絡(luò)方法分析工程機(jī)械類(lèi)產(chǎn)品［8］、棉花［9］、天然氣［10］、能源［11］等具體產(chǎn)品的貿(mào)易網(wǎng)絡(luò)結(jié)構(gòu)特征。二是建立模型實(shí)證分析貿(mào)易網(wǎng)絡(luò)的影響因素。其中，應(yīng)用比較多的如二次指派程序（QAP）分析方法和指數(shù)隨機(jī)圖模型（ERGM）。QAP 分析方法是一種應(yīng)用于網(wǎng)絡(luò)矩陣之間相關(guān)性與回歸性，可以研究多個(gè)“屬性數(shù)據(jù)”矩陣與一個(gè)“關(guān)系數(shù)據(jù)”矩陣之間的回歸關(guān)系［12］。王介勇等［13］將QAP 分析方法應(yīng)用于小麥、水稻、玉米世界三大主糧構(gòu)建的全球糧食貿(mào)易網(wǎng)絡(luò)，定量評(píng)估經(jīng)濟(jì)社會(huì)、資源稟賦、貿(mào)易政策、距離交通、語(yǔ)言文化、政治管理等對(duì)全球糧食貿(mào)易網(wǎng)絡(luò)的影響；馬遠(yuǎn)等［14］應(yīng)用此方法評(píng)估經(jīng)濟(jì)、能源、陸地相鄰與否、語(yǔ)言、人口、物流績(jī)效水平對(duì)絲綢之路經(jīng)濟(jì)帶沿線(xiàn)國(guó)家石油貿(mào)易網(wǎng)絡(luò)的影響。但QAP 分析方法無(wú)法解釋網(wǎng)絡(luò)中關(guān)系與關(guān)系之間的依賴(lài)性。而ERGM 模型是社會(huì)網(wǎng)絡(luò)統(tǒng)計(jì)模型，可以用于解釋網(wǎng)絡(luò)關(guān)系的出現(xiàn)和不出現(xiàn)，因此可用于分析網(wǎng)絡(luò)結(jié)構(gòu)［15］。不同于現(xiàn)有的經(jīng)濟(jì)計(jì)量模型，ERGM 模型可以對(duì)國(guó)際貿(mào)易網(wǎng)絡(luò)的結(jié)構(gòu)控制和調(diào)整提供依據(jù)［16］。此外，ERGM模型可以將多個(gè)層次的網(wǎng)絡(luò)變量納入一個(gè)模型中，綜合考慮形成機(jī)理，這是其相比于一般統(tǒng)計(jì)模型的最大優(yōu)勢(shì)［17］。在實(shí)際應(yīng)用中，吳鋼［16］基于復(fù)雜網(wǎng)絡(luò)理論構(gòu)建并描述1962—2010 年的國(guó)際貿(mào)易網(wǎng)絡(luò)和國(guó)際人文關(guān)系網(wǎng)絡(luò)，在此基礎(chǔ)上應(yīng)用ERGM 模型考察國(guó)際人文關(guān)系網(wǎng)絡(luò)對(duì)不同等級(jí)維度、不同年份國(guó)際貿(mào)易網(wǎng)絡(luò)的影響效應(yīng)；許和連等［18］將復(fù)雜網(wǎng)絡(luò)分析方法用于高端制造業(yè)，構(gòu)建2012 年“一帶一路”國(guó)家高端制造業(yè)貿(mào)易網(wǎng)絡(luò)，進(jìn)行網(wǎng)絡(luò)特征分析，并基于ERGM 模型分析影響該網(wǎng)絡(luò)形成的因素；牛華等［19］從網(wǎng)絡(luò)視角分析全球高技術(shù)產(chǎn)品出口的動(dòng)態(tài)特征并利用ERGM 模型探究研發(fā)投入和貿(mào)易自由化對(duì)高技術(shù)產(chǎn)品出口的影響。

目前，對(duì)林產(chǎn)品貿(mào)易的研究主要集中在單獨(dú)國(guó)家［20-21］、雙方國(guó)家間［22-23］或一定區(qū)域內(nèi)［24-25］，鮮有研究將貿(mào)易對(duì)象確定為全球。探究林產(chǎn)品貿(mào)易影響因素的文章使用的研究方法以引力模型［26-27］為主，方法較單一，缺乏創(chuàng)新性?；诖耍疚膮⒄諒?fù)雜網(wǎng)絡(luò)理論在國(guó)際貿(mào)易中的應(yīng)用，建立全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)，從宏觀和微觀兩個(gè)層面分析其格局；同時(shí)，建立ERGM 模型從內(nèi)生結(jié)構(gòu)變量、節(jié)點(diǎn)屬性和網(wǎng)絡(luò)協(xié)變量三方面探究林產(chǎn)品出口貿(mào)易影響因素。

二、理論分析

（一）要素稟賦理論

要素稟賦指一國(guó)所擁有的生產(chǎn)資源狀況。在新古典經(jīng)濟(jì)中，國(guó)家之間要素稟賦的差異是指土地、勞動(dòng)和資本等生產(chǎn)要素相對(duì)量的不同。要素稟賦理論認(rèn)為這種要素稟賦差異引起國(guó)家之間自由貿(mào)易，在開(kāi)展自由貿(mào)易后，一國(guó)出口密集使用其豐裕要素的產(chǎn)品，進(jìn)口密集使用其稀缺要素的產(chǎn)品［28］。因此，要素稟賦理論建立在國(guó)家之間專(zhuān)業(yè)化分工的基礎(chǔ)上，這種分工經(jīng)濟(jì)使國(guó)家之間依存關(guān)系更加密切，從而表現(xiàn)為國(guó)際貿(mào)易網(wǎng)絡(luò)［16］。根據(jù)要素稟賦理論可知，一國(guó)的森林資源是否豐富，資本積累情況也會(huì)影響一國(guó)在全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中的地位。

（二）關(guān)稅同盟理論

關(guān)稅同盟是指在自由貿(mào)易區(qū)的基礎(chǔ)上，成員國(guó)之間實(shí)行自由貿(mào)易，并且統(tǒng)一對(duì)非成員國(guó)實(shí)行進(jìn)口關(guān)稅或其他貿(mào)易政策措施，即建立共同關(guān)稅壁壘［28］。關(guān)稅同盟的靜態(tài)效應(yīng)包括貿(mào)易創(chuàng)造效應(yīng)、貿(mào)易轉(zhuǎn)移效應(yīng)和福利效應(yīng)等，動(dòng)態(tài)效應(yīng)包括擴(kuò)大出口效應(yīng)、規(guī)模經(jīng)濟(jì)效應(yīng)、促進(jìn)競(jìng)爭(zhēng)效應(yīng)、刺激投資效應(yīng)和資源配置效應(yīng)等［28］。本文從動(dòng)態(tài)效應(yīng)的規(guī)模經(jīng)濟(jì)效應(yīng)和促進(jìn)競(jìng)爭(zhēng)效應(yīng)兩個(gè)方面進(jìn)行分析：第一，規(guī)模經(jīng)濟(jì)效應(yīng)。關(guān)稅同盟的建立為成員國(guó)之間產(chǎn)品的相互出口創(chuàng)造良好的條件，所有成員國(guó)的國(guó)內(nèi)市場(chǎng)組成一個(gè)統(tǒng)一的區(qū)域性市場(chǎng)，這種交易市場(chǎng)范圍的擴(kuò)大使成員國(guó)企業(yè)可以擁有更大的需求市場(chǎng)，市場(chǎng)需求的擴(kuò)大使各企業(yè)可以進(jìn)行林產(chǎn)品貿(mào)易的專(zhuān)業(yè)化分工，提高資源利用效率，實(shí)現(xiàn)企業(yè)規(guī)模經(jīng)濟(jì)，最大化成本優(yōu)勢(shì)。第二，促進(jìn)競(jìng)爭(zhēng)效應(yīng)。建立關(guān)稅同盟后，一方面各成員國(guó)林產(chǎn)品企業(yè)面臨本國(guó)和其他成員國(guó)同類(lèi)企業(yè)的競(jìng)爭(zhēng)，這種良性競(jìng)爭(zhēng)不僅促進(jìn)產(chǎn)品流通，還倒逼企業(yè)提高要素生產(chǎn)率，進(jìn)行創(chuàng)新以降低成本，提高經(jīng)濟(jì)福利；另一方面，提高成員國(guó)林產(chǎn)品市場(chǎng)的透明度，使得資源配置效率改善，從而提高經(jīng)濟(jì)效率、增進(jìn)社會(huì)利益。綜上，通過(guò)簽署區(qū)域貿(mào)易協(xié)定實(shí)現(xiàn)貿(mào)易自由化，影響國(guó)家間林產(chǎn)品貿(mào)易情況，進(jìn)而影響全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)的形成。

（三）交易成本理論

交易成本泛指所有為促成交易發(fā)生而形成的成本，包括搜尋成本、信息成本、議價(jià)成本、監(jiān)督成本、違約成本等。因此，國(guó)家之間地理距離的遠(yuǎn)近，使用的語(yǔ)言和貨幣是否一致，都會(huì)導(dǎo)致不同國(guó)家間進(jìn)行林產(chǎn)品貿(mào)易的交易成本的不同，這也會(huì)進(jìn)一步影響全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)的形成。

三、模型選擇、變量選取與數(shù)據(jù)來(lái)源

（一）模型選擇

指數(shù)隨機(jī)圖模型（Exponential Random Graph Models，ERGM）是利用網(wǎng)絡(luò)統(tǒng)計(jì)量解釋網(wǎng)絡(luò)形成概率的統(tǒng)計(jì)模型，該模型基于指數(shù)形式來(lái)研究網(wǎng)絡(luò)結(jié)構(gòu)，解釋網(wǎng)絡(luò)中連接關(guān)系形成的原因。ERGM模型與傳統(tǒng)線(xiàn)性模型不同，其假定一條關(guān)系出現(xiàn)的概率取決于其他關(guān)系是否出現(xiàn)，強(qiáng)調(diào)網(wǎng)絡(luò)中關(guān)系間的依賴(lài)性。觀察到的網(wǎng)絡(luò)圖只是所有隨機(jī)網(wǎng)絡(luò)圖中的一張，其出現(xiàn)的概率受到網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)、節(jié)點(diǎn)屬性和網(wǎng)絡(luò)協(xié)變量因素的影響［29］。而林產(chǎn)品貿(mào)易網(wǎng)絡(luò)正是多個(gè)經(jīng)濟(jì)體貿(mào)易微觀行為相互作用形成，因此利用ERGM模型研究全球林產(chǎn)品貿(mào)易形成機(jī)制具有合理性。

根據(jù)Robins 等［30］、吳鋼［16］和許和連等［18］的研究，用Pr(Y=y|θ)表示在θ條件下網(wǎng)絡(luò)圖y出現(xiàn)的概率，則ERGM模型的一般形式為：

式（1）中，gα(y)表示可能影響網(wǎng)絡(luò)中關(guān)系形成的純網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量，gβ(y，x)表示與節(jié)點(diǎn)屬性相關(guān)的網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量，gγ(y，z)表示與其他外部網(wǎng)絡(luò)相關(guān)的網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量，θα、θβ和θγ是相應(yīng)的估計(jì)參數(shù)，κ是用于確保模型具有適當(dāng)概率分布的標(biāo)準(zhǔn)化常量。

ERGM模型主要通過(guò)估計(jì)、診斷、仿真、比較和改進(jìn)等步驟，檢驗(yàn)?zāi)男┮蛩仫@著影響網(wǎng)絡(luò)的形成。目前，主要采用馬爾可夫鏈蒙特卡羅極大似然估計(jì)法（MCMC MLE）對(duì)模型進(jìn)行估計(jì)檢驗(yàn)，通過(guò)不斷對(duì)模型進(jìn)行模擬和參數(shù)修正，使得模型的參數(shù)估計(jì)趨于穩(wěn)定。在模型評(píng)價(jià)方面，采用t統(tǒng)計(jì)量檢驗(yàn)參數(shù)的顯著性，一般要求p值小于0.1，并采用赤池信息準(zhǔn)則（AIC）和貝葉斯信息準(zhǔn)則（BIC）評(píng)估模型的擬合優(yōu)度。本文采用R語(yǔ)言中ergm包構(gòu)建模型并完成ERGM估計(jì)。

（二）變量選取

基于本文研究的是林產(chǎn)品貿(mào)易網(wǎng)絡(luò)，ERGM 模型的被解釋變量Pr(Y=y|θ)是林產(chǎn)品貿(mào)易網(wǎng)絡(luò)形成的概率，解釋變量是一系列林產(chǎn)品貿(mào)易網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量，包括內(nèi)生結(jié)構(gòu)變量gα(y)、節(jié)點(diǎn)屬性變量gβ(y，x)和網(wǎng)絡(luò)協(xié)變量gγ(y，z)，具體如表1所示。

表1 ERGM 變量說(shuō)明

1.內(nèi)生結(jié)構(gòu)變量 內(nèi)生結(jié)構(gòu)變量是指由網(wǎng)絡(luò)結(jié)構(gòu)內(nèi)部關(guān)系來(lái)表示的局部網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量，是用于檢驗(yàn)和控制網(wǎng)絡(luò)內(nèi)生的自我組織特征的變量，包括邊數(shù)、互惠性等。借鑒許和連等［18］的研究，選取邊數(shù)和互惠性?xún)煞N變量作為網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)變量。其中，邊數(shù)類(lèi)似線(xiàn)性回歸模型的常數(shù)項(xiàng)，是控制其他條件情況下關(guān)系發(fā)生的基本傾向或發(fā)生概率。互惠性表示林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中國(guó)家間是否存在互相出口的關(guān)系，參數(shù)估計(jì)值為正，則表明林產(chǎn)品出口網(wǎng)絡(luò)中的國(guó)家間傾向于存在互相出口的關(guān)系。

2.節(jié)點(diǎn)屬性變量 在網(wǎng)絡(luò)中，節(jié)點(diǎn)自身的一些屬性對(duì)網(wǎng)絡(luò)關(guān)系的形成也具有重要影響。在探討的林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中，節(jié)點(diǎn)屬性變量用于測(cè)量國(guó)家屬性的變量，國(guó)家自身的經(jīng)濟(jì)發(fā)展水平、資源稟賦等屬性也會(huì)影響林產(chǎn)品貿(mào)易網(wǎng)絡(luò)關(guān)系的形成。因此，本文所選取的節(jié)點(diǎn)屬性變量為主效應(yīng)、同配性和發(fā)送者效應(yīng)。

主效應(yīng)用于檢驗(yàn)網(wǎng)絡(luò)中是否某種屬性越強(qiáng)，越容易出口?？紤]國(guó)內(nèi)生產(chǎn)總值（GDP）可以反映一國(guó)經(jīng)濟(jì)規(guī)模，人均國(guó)內(nèi)生產(chǎn)總值可以反映一國(guó)的經(jīng)濟(jì)發(fā)展水平，因此選取GDP 和人均GDP 作為主效應(yīng)，分析是否一國(guó)經(jīng)濟(jì)規(guī)模越大，經(jīng)濟(jì)發(fā)展水平越高越容易出口。同配性檢驗(yàn)屬性相近的國(guó)家之間發(fā)生貿(mào)易關(guān)系的傾向，即是否存在“同性相吸”的情形。參考吳鋼［16］的處理方法，選取森林面積作為同配性，且將森林面積分為高、中、低三類(lèi)。具體地，將72 個(gè)國(guó)家中森林面積排在前三分之一的歸為森林面積高，前三分之一至三分之二歸為森林面積中，后三分之一歸為森林面積低?？疾旆謱偃?lèi)內(nèi)部的國(guó)家之間是否有發(fā)生貿(mào)易關(guān)系的傾向。發(fā)送者效應(yīng)則是檢驗(yàn)具有某種屬性的國(guó)家相比其他國(guó)家是否發(fā)送出更多關(guān)系，即具有更高的出口傾向，出口的國(guó)家數(shù)量更多。在發(fā)送者效應(yīng)方面，為考察森林面積是否會(huì)影響一國(guó)的林產(chǎn)品出口，分別對(duì)發(fā)送者效應(yīng)（森林面積中）和發(fā)送者效應(yīng)（森林面積高）進(jìn)行估計(jì)?？疾焯幱谏置娣e高和森林面積中這兩類(lèi)的國(guó)家相比森林面積低這類(lèi)的國(guó)家是否具有更高的出口傾向。

3.網(wǎng)絡(luò)協(xié)變量 外部環(huán)境對(duì)網(wǎng)絡(luò)形成的影響通常用網(wǎng)絡(luò)協(xié)變量進(jìn)行測(cè)算。網(wǎng)絡(luò)協(xié)變量就是檢驗(yàn)其他網(wǎng)絡(luò)關(guān)系是否會(huì)促進(jìn)貿(mào)易關(guān)系形成的變量。本文選擇貿(mào)易自由化網(wǎng)絡(luò)、距離網(wǎng)絡(luò)、共同語(yǔ)言網(wǎng)絡(luò)和共同貨幣網(wǎng)絡(luò)，分別檢驗(yàn)貿(mào)易自由化、距離、語(yǔ)言及貨幣對(duì)林產(chǎn)品出口網(wǎng)絡(luò)的影響。

（三）數(shù)據(jù)來(lái)源

本文借鑒伍海泉等［31］、張洪瑞等［32］對(duì)林產(chǎn)品的界定方法，綜合參考聯(lián)合國(guó)糧食及農(nóng)業(yè)組織（FAO）分類(lèi)、《中國(guó)林業(yè)和草原統(tǒng)計(jì)年鑒》以及《中國(guó)海關(guān)統(tǒng)計(jì)年鑒》，選用HS1992 分類(lèi)標(biāo)準(zhǔn)，數(shù)據(jù)來(lái)源于UN Comtrade 數(shù)據(jù)庫(kù)。本文選取2019 年為研究年份，選取72 個(gè)國(guó)家，其中包括亞洲國(guó)家29 個(gè)、歐洲國(guó)家28個(gè)、非洲國(guó)家5 個(gè)、大洋洲國(guó)家2 個(gè)、美洲國(guó)家8 個(gè)，所選國(guó)家涵蓋全球大多數(shù)國(guó)際組織或聯(lián)盟成員國(guó)，2019 年林產(chǎn)品出口貿(mào)易額占全球比例為96%，林產(chǎn)品進(jìn)口貿(mào)易額占全球比例為92%，能夠較好地代表全球林產(chǎn)品貿(mào)易水平。為避免過(guò)小的貿(mào)易量不足以反映國(guó)家間的貿(mào)易模式，清晰反映林產(chǎn)品貿(mào)易關(guān)系的親密程度，參照許和連等［18］研究，只保留林產(chǎn)品出口額大于100 萬(wàn)美元的貿(mào)易邊，提取主要林產(chǎn)品貿(mào)易關(guān)系網(wǎng)絡(luò)進(jìn)行分析。

基于上述ERGM 模型的變量選取，GDP（美元）、人均GDP（美元）及森林面積（平方千米）均可從世界銀行網(wǎng)站獲得。為得到更好的估計(jì)結(jié)果，對(duì)GDP 和人均GDP 進(jìn)行對(duì)數(shù)處理。貿(mào)易自由化網(wǎng)絡(luò)、距離網(wǎng)絡(luò)、共同語(yǔ)言網(wǎng)絡(luò)和共同貨幣網(wǎng)絡(luò)數(shù)據(jù)來(lái)源于CEPⅡ數(shù)據(jù)庫(kù)，其中貿(mào)易自由化網(wǎng)絡(luò)根據(jù)貿(mào)易雙方國(guó)家是否均為WTO成員國(guó)測(cè)定。

四、結(jié)果與分析

（一）全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)分析結(jié)果

在經(jīng)濟(jì)全球化背景下，全球各國(guó)和地區(qū)間貿(mào)易活動(dòng)形成相互作用和相互依賴(lài)的復(fù)雜貿(mào)易網(wǎng)絡(luò)。林產(chǎn)品貿(mào)易網(wǎng)絡(luò)將各個(gè)國(guó)家和地區(qū)視為點(diǎn)，將雙邊貿(mào)易關(guān)系視為線(xiàn)，進(jìn)而用G=(Vi，Vj，A，W)表示全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)。其中，向量Vi=[vi](i= 1，2，…，n)表示出口國(guó)，向量Vj=[vj](j= 1，2，…，n)表示進(jìn)口國(guó)，鄰接矩陣A=[aij](i，j= 1，2，…，n)表示兩國(guó)之間的貿(mào)易關(guān)系，當(dāng)i國(guó)對(duì)j國(guó)出口大于0 時(shí)，aij= 1，否則aij= 0；權(quán)重矩陣W=[wij](i，j= 1，2，…，n)表示Vi國(guó)對(duì)Vj國(guó)出口的貿(mào)易值。本文從宏觀和微觀兩個(gè)層面對(duì)林產(chǎn)品貿(mào)易網(wǎng)絡(luò)進(jìn)行分析。

1.宏觀層面 在宏觀層面，選取網(wǎng)絡(luò)密度和基于距離的凝聚指數(shù)（緊湊度）兩個(gè)指標(biāo)。網(wǎng)絡(luò)密度可以反映國(guó)家間貿(mào)易關(guān)系的疏密程度，通過(guò)網(wǎng)絡(luò)實(shí)際貿(mào)易關(guān)系數(shù)除以理論最大貿(mào)易關(guān)系數(shù)來(lái)測(cè)算。網(wǎng)絡(luò)密度越大，表示國(guó)家之間貿(mào)易聯(lián)系越緊密。設(shè)網(wǎng)絡(luò)中包含的實(shí)際貿(mào)易關(guān)系數(shù)為m，網(wǎng)絡(luò)中的國(guó)家個(gè)數(shù)為n，則包含的貿(mào)易關(guān)系總數(shù)在理論上的最大可能值為n（n-1），則網(wǎng)絡(luò)密度的計(jì)算公式為：

根據(jù)式（2）計(jì)算出2019年全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)密度為0.501。各國(guó)之間林產(chǎn)品貿(mào)易聯(lián)系較緊密，貿(mào)易關(guān)聯(lián)性較強(qiáng)，但仍存在較大提升空間。

基于距離的凝聚指數(shù)（緊湊度）反映國(guó)家間貿(mào)易往來(lái)的通達(dá)性，是平均距離的倒數(shù)，其中，兩國(guó)的距離指兩國(guó)至少可通過(guò)多少條邊關(guān)聯(lián)在一起。網(wǎng)絡(luò)中包含n（n-1）對(duì)國(guó)家，每對(duì)國(guó)家的距離之和是∑ijdij，因此網(wǎng)絡(luò)中兩兩國(guó)家之間的平均距離計(jì)算公式為：

當(dāng)網(wǎng)絡(luò)中的所有國(guó)家兩兩之間存在直接貿(mào)易關(guān)系時(shí)，緊湊度等于1；不存在任何關(guān)系時(shí)，緊湊度等于0。緊湊度越接近1，表明網(wǎng)絡(luò)越具有凝聚力，貿(mào)易網(wǎng)絡(luò)任意兩國(guó)間普遍存在直接貿(mào)易關(guān)系，貿(mào)易往來(lái)的通達(dá)性越強(qiáng)。

根據(jù)式（3）和式（4）分別計(jì)算出2019 年全球林產(chǎn)品出口網(wǎng)絡(luò)的平均距離為1.481，緊湊度為0.696，較為接近1，因此該網(wǎng)絡(luò)任意兩國(guó)基本普遍存在直接貿(mào)易關(guān)系，網(wǎng)絡(luò)凝聚力較強(qiáng)，各國(guó)林產(chǎn)品貿(mào)易往來(lái)通達(dá)性較強(qiáng)。

2.微觀層面 選取出度中心度和入度中心度、中間中心度兩個(gè)指標(biāo)。一國(guó)的出度中心度和入度中心度分別反映一國(guó)發(fā)出和接收貿(mào)易關(guān)系的能力，即一國(guó)出口和進(jìn)口的能力。i國(guó)的點(diǎn)入度是出口到i國(guó)的其他國(guó)的個(gè)數(shù)，即i國(guó)得到的直接關(guān)系數(shù)；i國(guó)的點(diǎn)出度是i國(guó)直接發(fā)出的關(guān)系數(shù)。為便于比較，參照劉軍［33］的研究，利用直接關(guān)系數(shù)除以最大可能度數(shù)（n-1）進(jìn)行計(jì)算。數(shù)值越大，表示該國(guó)的貿(mào)易對(duì)象越廣泛。設(shè)i國(guó)直接發(fā)出的關(guān)系數(shù)∑jaij，i國(guó)直接接收的關(guān)系數(shù)為∑jaji，則i國(guó)出度中心度（ODCi）和入度中心度（IDCi）的計(jì)算公式分別為：

一國(guó)的中間中心度反映一國(guó)控制其他國(guó)家貿(mào)易往來(lái)的能力。這里定義捷徑為一國(guó)出口至另一國(guó)的最短路徑，假設(shè)t國(guó)出口到s國(guó)存在多條捷徑，其中有一條可能經(jīng)過(guò)i國(guó)。其中間中心度可以用t國(guó)向s國(guó)出口并且經(jīng)過(guò)i國(guó)的捷徑數(shù)與t國(guó)出口到s國(guó)的捷徑總數(shù)之比來(lái)計(jì)算。

為便于比較，參照羅家德［34］的研究，利用上述比值除以最大可能度數(shù)（n-1）（n-2）進(jìn)行計(jì)算。數(shù)值越大，表示該國(guó)影響其他國(guó)家之間貿(mào)易往來(lái)的能力越強(qiáng)。設(shè)dtis為t國(guó)出口到s國(guó)并且經(jīng)過(guò)i國(guó)的捷徑數(shù)，dts是t國(guó)出口到s國(guó)的捷徑總數(shù)，則i國(guó)的中間中心度的計(jì)算公式為：

表2所示為2019年全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)出度中心度和入度中心度、中間中心度指標(biāo)排名前20的國(guó)家。表2顯示，中國(guó)、德國(guó)和意大利各中心性指標(biāo)均排在前10名，表明這三國(guó)出口和進(jìn)口貿(mào)易對(duì)象都最為廣泛，影響其他國(guó)家貿(mào)易往來(lái)的能力最強(qiáng)，處于網(wǎng)絡(luò)核心地位，發(fā)揮著重要的“樞紐”作用。西班牙、法國(guó)、荷蘭、土耳其、比利時(shí)、英國(guó)、美國(guó)和韓國(guó)的各項(xiàng)中心性指標(biāo)均位列前20 名，說(shuō)明這些國(guó)家在此網(wǎng)絡(luò)中具有一定的影響力，但并沒(méi)有位于核心位置。中亞國(guó)家所占比例相較其他國(guó)家比例較小，處于林產(chǎn)品貿(mào)易網(wǎng)絡(luò)邊緣。

表2 2019年全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中心性指標(biāo)前20名國(guó)家

（二）ERGM 估計(jì)結(jié)果

本文采用R 語(yǔ)言中ergm 包構(gòu)建模型并完成ERGM 估計(jì)，估計(jì)結(jié)果表明，互惠性、GDP 和人均GDP在1%的水平上顯著為正，互惠性估計(jì)系數(shù)為0.864，說(shuō)明全球林產(chǎn)品貿(mào)易具有較明顯的互惠性，許多國(guó)家之間互有貿(mào)易進(jìn)口和出口關(guān)系；GDP 估計(jì)系數(shù)為正，表明一國(guó)經(jīng)濟(jì)規(guī)模越大，越容易展開(kāi)林產(chǎn)品出口貿(mào)易；人均GDP 估計(jì)系數(shù)為正，但略小于GDP 的估計(jì)系數(shù)，表明經(jīng)濟(jì)發(fā)展水平越高，越容易展開(kāi)林產(chǎn)品出口貿(mào)易，但影響強(qiáng)度不如GDP。模型（1）考察包含網(wǎng)絡(luò)的內(nèi)生結(jié)構(gòu)變量和節(jié)點(diǎn)屬性中的主效應(yīng)，在模型（1）的基礎(chǔ)上引入同配性考察森林面積相近的國(guó)家之間是否更傾向于發(fā)生貿(mào)易關(guān)系。同配性（森林面積高）和同配性（森林面積中）在1%的水平上顯著為正，且前者系數(shù)大于后者，而同配性（森林面積低）在5%的水平上顯著為負(fù)。這表明森林面積在高等和中等水平的國(guó)家之間更傾向于發(fā)生林產(chǎn)品貿(mào)易，且森林面積越大，這種傾向就越強(qiáng)，而在低等森林面積水平的國(guó)家之間開(kāi)展的傾向較小。造成這種現(xiàn)象的一種可能原因是，低等森林面積水平的國(guó)家林產(chǎn)品種類(lèi)和數(shù)量普遍偏低，因此不傾向發(fā)生林產(chǎn)品貿(mào)易。

表3 ERGM 回歸結(jié)果

在模型（2）的基礎(chǔ)上，模型（3）～模型（6）分別加入貿(mào)易自由化網(wǎng)絡(luò)、距離網(wǎng)絡(luò)、共同語(yǔ)言網(wǎng)絡(luò)和共同貨幣網(wǎng)絡(luò)，考察4 種網(wǎng)絡(luò)對(duì)林產(chǎn)品貿(mào)易網(wǎng)絡(luò)的影響，即貿(mào)易自由化、距離、語(yǔ)言和貨幣在林產(chǎn)品貿(mào)易關(guān)系中的作用。估計(jì)結(jié)果顯示貿(mào)易自由化網(wǎng)絡(luò)、共同語(yǔ)言網(wǎng)絡(luò)和共同貨幣網(wǎng)絡(luò)在1%水平上顯著為正，距離網(wǎng)絡(luò)在1%水平上顯著為負(fù)。這表明林產(chǎn)品貿(mào)易網(wǎng)絡(luò)與貿(mào)易自由化網(wǎng)絡(luò)和共同貨幣網(wǎng)絡(luò)存在明顯重疊，對(duì)林產(chǎn)品貿(mào)易網(wǎng)絡(luò)具有重要的正向影響，貿(mào)易自由化網(wǎng)絡(luò)的影響更強(qiáng)；共同語(yǔ)言網(wǎng)絡(luò)與林產(chǎn)品貿(mào)易網(wǎng)絡(luò)間存在重疊，并且共同語(yǔ)言網(wǎng)絡(luò)對(duì)于林產(chǎn)品貿(mào)易網(wǎng)絡(luò)具有正向影響，但不如前兩個(gè)網(wǎng)絡(luò)影響顯著；距離網(wǎng)絡(luò)的影響效應(yīng)相對(duì)較小。模型（3）的估計(jì)結(jié)果顯示，貿(mào)易自由化估計(jì)系數(shù)高達(dá)1.382，表明貿(mào)易自由化對(duì)全球林產(chǎn)品貿(mào)易發(fā)展產(chǎn)生極大促進(jìn)作用。其可能原因?yàn)椋嘿Q(mào)易自由化使關(guān)稅降低，擴(kuò)大貿(mào)易。根據(jù)模型（4），距離對(duì)林產(chǎn)品貿(mào)易影響較小。林產(chǎn)品貿(mào)易主要依賴(lài)的運(yùn)輸方式為陸運(yùn)和海運(yùn)，距離對(duì)林產(chǎn)品貿(mào)易的貿(mào)易發(fā)展的限制作用仍然存在，但因?yàn)檫\(yùn)輸成本并不高，所以其影響較小。共同語(yǔ)言和共同貨幣代表林產(chǎn)品貿(mào)易中的交易成本，模型（5）和模型（6）的估計(jì)結(jié)果表明，具有共同語(yǔ)言和共同貨幣更有利于國(guó)家間開(kāi)展林產(chǎn)品貿(mào)易。這是因?yàn)槭褂孟嗤Z(yǔ)言有利于貿(mào)易雙方溝通交流，易于減少和解決貿(mào)易過(guò)程中出現(xiàn)的問(wèn)題；使用相同貨幣可以規(guī)避匯率風(fēng)險(xiǎn)，有利于貿(mào)易關(guān)系的形成。

模型（7）綜合上述內(nèi)生結(jié)構(gòu)變量、主效應(yīng)、同配性和網(wǎng)絡(luò)協(xié)變量，變量估計(jì)結(jié)果均在1%的水平上顯著，且符號(hào)方向均未改變。相比于模型（1）～模型（6），模型（7）中AIC和BIC數(shù)值減小，擬合優(yōu)度進(jìn)一步提高。上述模型可以說(shuō)明森林面積在高等和中等水平的國(guó)家間容易展開(kāi)林產(chǎn)品貿(mào)易，森林面積在低等水平的國(guó)家間不易展開(kāi)林產(chǎn)品貿(mào)易，但并未說(shuō)明森林面積處于較高水平的國(guó)家是否具有更高的出口傾向。為此，模型（8）在模型（7）的基礎(chǔ)上進(jìn)一步加入發(fā)送者效應(yīng)（森林面積中）和發(fā)送者效應(yīng)（森林面積高），估計(jì)結(jié)果顯示發(fā)送者效應(yīng)（森林面積中）和發(fā)送者效應(yīng)（森林面積高）均在1%的水平上顯著為正且后者大于前者，這表明高等和中等森林面積水平的國(guó)家相對(duì)于低等森林面積水平的國(guó)家具有更高的林產(chǎn)品貿(mào)易出口傾向，且森林面積越大的國(guó)家林產(chǎn)品出口傾向越高。

整體來(lái)看，模型（8）綜合所有變量，與前面的估計(jì)結(jié)果相比，AIC和BIC均為最低值，模型擬合優(yōu)度最高。模型（1）～模型（8）變量估計(jì)結(jié)果符號(hào)一致，驗(yàn)證模型的穩(wěn)健性。節(jié)點(diǎn)屬性變量中發(fā)送者效應(yīng)對(duì)網(wǎng)絡(luò)促進(jìn)作用最強(qiáng)，協(xié)變量網(wǎng)絡(luò)中貿(mào)易自由化網(wǎng)絡(luò)對(duì)林產(chǎn)品出口貿(mào)易網(wǎng)絡(luò)的正向影響最強(qiáng)。

五、主要結(jié)論與政策建議

（一）主要結(jié)論

基于構(gòu)建的2019 年全球林產(chǎn)品出口貿(mào)易網(wǎng)絡(luò)，運(yùn)用網(wǎng)絡(luò)分析方法分析其格局，并建立ERGM 模型探究影響該網(wǎng)絡(luò)形成的因素。得出如下主要結(jié)論：第一，從網(wǎng)絡(luò)宏觀層面分析，全球網(wǎng)絡(luò)密度和緊湊度均接近1，國(guó)家之間林產(chǎn)品貿(mào)易聯(lián)系緊密，貿(mào)易往來(lái)的通達(dá)性強(qiáng)；從網(wǎng)絡(luò)微觀層面分析，中國(guó)、德國(guó)和意大利的網(wǎng)絡(luò)中心性指標(biāo)位居世界前列，進(jìn)出口貿(mào)易對(duì)象最為廣泛，影響其他國(guó)家貿(mào)易往來(lái)的能力最強(qiáng)，在林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中發(fā)揮著重要的“樞紐”作用，眾多中亞國(guó)家處于林產(chǎn)品貿(mào)易網(wǎng)絡(luò)的邊緣。第二，全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò)的形成是由網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)變量、節(jié)點(diǎn)屬性變量和網(wǎng)絡(luò)協(xié)變量共同影響的。分析發(fā)現(xiàn)：內(nèi)生結(jié)構(gòu)變量中，全球林產(chǎn)品貿(mào)易具有較明顯的互惠性，許多國(guó)家之間互有貿(mào)易進(jìn)出口關(guān)系；節(jié)點(diǎn)屬性變量中，一國(guó)GDP和人均GDP越高，越容易展開(kāi)林產(chǎn)品出口貿(mào)易，人均GDP較GDP影響強(qiáng)度小；高等和中等水平森林面積的國(guó)家之間更傾向于發(fā)生林產(chǎn)品貿(mào)易，且相對(duì)于低等森林面積水平的國(guó)家具有更高的林產(chǎn)品貿(mào)易出口傾向，森林面積越大的國(guó)家這兩種傾向就越強(qiáng)；網(wǎng)絡(luò)協(xié)變量中，貿(mào)易自由化網(wǎng)絡(luò)、共同貨幣網(wǎng)絡(luò)和共同語(yǔ)言網(wǎng)絡(luò)對(duì)林產(chǎn)品貿(mào)易網(wǎng)絡(luò)存在重疊和正向影響，影響強(qiáng)度依次減弱，距離網(wǎng)絡(luò)的影響較小。

（二）政策建議

基于以上結(jié)論，得出如下政策建議：第一，中國(guó)應(yīng)繼續(xù)發(fā)揮在林產(chǎn)品貿(mào)易中的“樞紐”作用，一方面通過(guò)政策溝通擴(kuò)大與其他核心國(guó)家的林產(chǎn)品貿(mào)易范圍，深化經(jīng)貿(mào)關(guān)系；另一方面加快建立與林產(chǎn)品貿(mào)易網(wǎng)絡(luò)邊緣國(guó)家的林產(chǎn)品貿(mào)易聯(lián)系，構(gòu)建多層次貿(mào)易網(wǎng)絡(luò)，實(shí)現(xiàn)合作共贏。第二，保持國(guó)家經(jīng)濟(jì)狀況穩(wěn)步提升，經(jīng)濟(jì)平穩(wěn)發(fā)展，積極推進(jìn)貿(mào)易自由化進(jìn)程，堅(jiān)持對(duì)外開(kāi)放，這對(duì)林產(chǎn)品貿(mào)易具有積極的促進(jìn)作用。第三，中國(guó)應(yīng)繼續(xù)在進(jìn)行林產(chǎn)品出口貿(mào)易的過(guò)程中注重森林資源的保護(hù)問(wèn)題，積極踐行“綠水青山就是金山銀山”理念，以此保持森林面積處于世界高等水平地位，穩(wěn)固林產(chǎn)品貿(mào)易“大國(guó)”地位。