吳紅梅，王 蕊，殷沐陽

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

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

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

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

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

二、理論分析

（一）要素稟賦理論

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

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

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

（三）交易成本理論

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

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

（一）模型選擇

指數(shù)隨機(jī)圖模型（Exponential Random Graph Models，ERGM）是利用網(wǎng)絡(luò)統(tǒng)計量解釋網(wǎng)絡(luò)形成概率的統(tǒng)計模型，該模型基于指數(shù)形式來研究網(wǎng)絡(luò)結(jié)構(gòu)，解釋網(wǎng)絡(luò)中連接關(guān)系形成的原因。ERGM模型與傳統(tǒng)線性模型不同，其假定一條關(guān)系出現(xiàn)的概率取決于其他關(guān)系是否出現(xiàn)，強(qiáng)調(diào)網(wǎng)絡(luò)中關(guān)系間的依賴性。觀察到的網(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ò)正是多個經(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)計量，gβ(y，x)表示與節(jié)點(diǎn)屬性相關(guān)的網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計量，gγ(y，z)表示與其他外部網(wǎng)絡(luò)相關(guān)的網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計量，θα、θβ和θγ是相應(yīng)的估計參數(shù)，κ是用于確保模型具有適當(dāng)概率分布的標(biāo)準(zhǔn)化常量。

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

（二）變量選取

基于本文研究的是林產(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)計量，包括內(nèi)生結(jié)構(gòu)變量gα(y)、節(jié)點(diǎn)屬性變量gβ(y，x)和網(wǎng)絡(luò)協(xié)變量gγ(y，z)，具體如表1所示。

表1 ERGM 變量說明

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

2.節(jié)點(diǎn)屬性變量 在網(wǎng)絡(luò)中，節(jié)點(diǎn)自身的一些屬性對網(wǎng)絡(luò)關(guān)系的形成也具有重要影響。在探討的林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中，節(jié)點(diǎn)屬性變量用于測量國家屬性的變量，國家自身的經(jīng)濟(jì)發(fā)展水平、資源稟賦等屬性也會影響林產(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)，越容易出口?？紤]國內(nèi)生產(chǎn)總值（GDP）可以反映一國經(jīng)濟(jì)規(guī)模，人均國內(nèi)生產(chǎn)總值可以反映一國的經(jīng)濟(jì)發(fā)展水平，因此選取GDP 和人均GDP 作為主效應(yīng)，分析是否一國經(jīng)濟(jì)規(guī)模越大，經(jīng)濟(jì)發(fā)展水平越高越容易出口。同配性檢驗(yàn)屬性相近的國家之間發(fā)生貿(mào)易關(guān)系的傾向，即是否存在“同性相吸”的情形。參考吳鋼［16］的處理方法，選取森林面積作為同配性，且將森林面積分為高、中、低三類。具體地，將72 個國家中森林面積排在前三分之一的歸為森林面積高，前三分之一至三分之二歸為森林面積中，后三分之一歸為森林面積低。考察分屬三類內(nèi)部的國家之間是否有發(fā)生貿(mào)易關(guān)系的傾向。發(fā)送者效應(yīng)則是檢驗(yàn)具有某種屬性的國家相比其他國家是否發(fā)送出更多關(guān)系，即具有更高的出口傾向，出口的國家數(shù)量更多。在發(fā)送者效應(yīng)方面，為考察森林面積是否會影響一國的林產(chǎn)品出口，分別對發(fā)送者效應(yīng)（森林面積中）和發(fā)送者效應(yīng)（森林面積高）進(jìn)行估計。考察處于森林面積高和森林面積中這兩類的國家相比森林面積低這類的國家是否具有更高的出口傾向。

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

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

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

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

四、結(jié)果與分析

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

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

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

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

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

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

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

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

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

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

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

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

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

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

表3 ERGM 回歸結(jié)果

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

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

整體來看，模型（8）綜合所有變量，與前面的估計結(jié)果相比，AIC和BIC均為最低值，模型擬合優(yōu)度最高。模型（1）～模型（8）變量估計結(jié)果符號一致，驗(yàn)證模型的穩(wěn)健性。節(jié)點(diǎn)屬性變量中發(fā)送者效應(yīng)對網(wǎng)絡(luò)促進(jìn)作用最強(qiáng)，協(xié)變量網(wǎng)絡(luò)中貿(mào)易自由化網(wǎng)絡(luò)對林產(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，國家之間林產(chǎn)品貿(mào)易聯(lián)系緊密，貿(mào)易往來的通達(dá)性強(qiáng)；從網(wǎng)絡(luò)微觀層面分析，中國、德國和意大利的網(wǎng)絡(luò)中心性指標(biāo)位居世界前列，進(jìn)出口貿(mào)易對象最為廣泛，影響其他國家貿(mào)易往來的能力最強(qiáng)，在林產(chǎn)品貿(mào)易網(wǎng)絡(luò)中發(fā)揮著重要的“樞紐”作用，眾多中亞國家處于林產(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)易具有較明顯的互惠性，許多國家之間互有貿(mào)易進(jìn)出口關(guān)系；節(jié)點(diǎn)屬性變量中，一國GDP和人均GDP越高，越容易展開林產(chǎn)品出口貿(mào)易，人均GDP較GDP影響強(qiáng)度??；高等和中等水平森林面積的國家之間更傾向于發(fā)生林產(chǎn)品貿(mào)易，且相對于低等森林面積水平的國家具有更高的林產(chǎn)品貿(mào)易出口傾向，森林面積越大的國家這兩種傾向就越強(qiáng)；網(wǎng)絡(luò)協(xié)變量中，貿(mào)易自由化網(wǎng)絡(luò)、共同貨幣網(wǎng)絡(luò)和共同語言網(wǎng)絡(luò)對林產(chǎn)品貿(mào)易網(wǎng)絡(luò)存在重疊和正向影響，影響強(qiáng)度依次減弱，距離網(wǎng)絡(luò)的影響較小。

（二）政策建議

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