馬麗麗,李晉如

(齊齊哈爾大學(xué) 理學(xué)院,黑龍江 齊齊哈爾 161006)

0 引言

隨著社會(huì)網(wǎng)絡(luò)理論研究取得突破性進(jìn)展,廣大學(xué)者發(fā)現(xiàn)國(guó)際貿(mào)易與社會(huì)網(wǎng)絡(luò)高度契合,屬于典型的社會(huì)網(wǎng)絡(luò)關(guān)系,運(yùn)用社會(huì)網(wǎng)絡(luò)分析方法研究國(guó)際貿(mào)易領(lǐng)域的相關(guān)現(xiàn)象成為新熱點(diǎn).A.Wilhite[1]發(fā)現(xiàn)國(guó)家之間的雙邊貿(mào)易交流最終會(huì)演變成具有小世界特征的國(guó)際貿(mào)易網(wǎng)絡(luò);M.E.J.Newman[2]總結(jié)了能夠表達(dá)網(wǎng)絡(luò)相關(guān)拓?fù)涮匦缘娜舾蓽y(cè)度指標(biāo),包括網(wǎng)絡(luò)平均最短距離、節(jié)點(diǎn)度分布、度相關(guān)等;王純杰等[3-4]運(yùn)用神經(jīng)網(wǎng)絡(luò)分段指數(shù)模型,研究了具有不同協(xié)變量的受試者危險(xiǎn)率變化情況,通過Logistic回歸模型和神經(jīng)網(wǎng)絡(luò)模型兩種方法對(duì)咽喉反流疾病進(jìn)行預(yù)測(cè);邢艷春等[5]將動(dòng)態(tài)回歸模型ARIMAX和BP神經(jīng)網(wǎng)絡(luò)算法進(jìn)行結(jié)合,設(shè)計(jì)了一種組合預(yù)測(cè)模型,研究人口結(jié)構(gòu)指標(biāo)對(duì)我國(guó)居民消費(fèi)結(jié)構(gòu)的影響發(fā)展趨勢(shì).目前,社會(huì)網(wǎng)絡(luò)分析方法在國(guó)際貿(mào)易領(lǐng)域的應(yīng)用已相對(duì)成熟,學(xué)者們主要以某一區(qū)域或特定產(chǎn)業(yè)為切入點(diǎn),從以下兩個(gè)方面展開研究.一方面,探究貿(mào)易網(wǎng)絡(luò)的結(jié)構(gòu)特征,測(cè)算具有代表性的網(wǎng)絡(luò)拓?fù)渲笜?biāo),并且劃分網(wǎng)絡(luò)的社團(tuán)結(jié)構(gòu),分析不同國(guó)家在貿(mào)易網(wǎng)絡(luò)中的地位和作用.例如,殷偉等[6]對(duì)全球茶葉貿(mào)易網(wǎng)絡(luò)進(jìn)行動(dòng)態(tài)分析,重點(diǎn)探究網(wǎng)絡(luò)的整體結(jié)構(gòu)特征、重要節(jié)點(diǎn)和貿(mào)易線變化以及我國(guó)茶葉出口市場(chǎng)的格局演變.另一方面,實(shí)證分析影響貿(mào)易網(wǎng)絡(luò)演化的因素和機(jī)制,所采用的模型有引力模型、QAP模型和指數(shù)隨機(jī)圖模型等.例如,王雁斌和程寶棟[7]從內(nèi)生自組織結(jié)構(gòu)、國(guó)家屬性和其他關(guān)系網(wǎng)絡(luò)等角度對(duì)原木貿(mào)易網(wǎng)絡(luò)的影響因素進(jìn)行實(shí)證檢驗(yàn);吳紅梅等[8]構(gòu)建了全球林產(chǎn)品貿(mào)易網(wǎng)絡(luò),采用指數(shù)隨機(jī)圖模型實(shí)證分析了影響該網(wǎng)絡(luò)形成的主要因素;段文奇等[9]發(fā)現(xiàn)國(guó)際貿(mào)易網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的異質(zhì)性下降;G.Fagiolo等[10]基于1981—2000年期間共計(jì)159個(gè)國(guó)家的貿(mào)易數(shù)據(jù)構(gòu)建加權(quán)網(wǎng)絡(luò),得出大部分國(guó)家貿(mào)易聯(lián)系稀疏,只有少數(shù)國(guó)家貿(mào)易互動(dòng)十分密集的結(jié)論;陳銀飛[11]利用社會(huì)網(wǎng)絡(luò)分析方法研究了國(guó)際貿(mào)易格局及其在次貸危機(jī)前后的變化.

本文基于社會(huì)網(wǎng)絡(luò)視角,利用進(jìn)出口貿(mào)易數(shù)據(jù)構(gòu)建全球化妝品貿(mào)易網(wǎng)絡(luò),并進(jìn)一步運(yùn)用指數(shù)隨機(jī)圖模型探討影響世界化妝品貿(mào)易網(wǎng)絡(luò)演化的因素,以期對(duì)中國(guó)化妝品貿(mào)易提出一些有益的建議.

1 研究方法和數(shù)據(jù)說明

1.1 研究方法

1.1.1 社會(huì)網(wǎng)絡(luò)分析

本文將化妝品貿(mào)易網(wǎng)絡(luò)中的國(guó)家抽象為節(jié)點(diǎn)Vi(i=1,2,…,n),使用鄰接矩陣A(t)=[aij](t)來表示各個(gè)節(jié)點(diǎn)之間的貿(mào)易關(guān)系(i,j=1,2,…,n;t表示時(shí)間),若存在從國(guó)家Vi出口到國(guó)家Vj的貿(mào)易流動(dòng)時(shí),aij=1,否則為0.

互惠性指標(biāo)即雙向比,是兩兩節(jié)點(diǎn)之間雙向關(guān)系的總數(shù)占網(wǎng)絡(luò)所有關(guān)系的比例,數(shù)值越大代表網(wǎng)絡(luò)中各成員國(guó)的貿(mào)易互通性越強(qiáng).假設(shè)有向圖中共有M條邊,其中有MB條是雙向邊,雙向比的計(jì)算公式為

(1)

度是刻畫單個(gè)節(jié)點(diǎn)屬性的指標(biāo)中最簡(jiǎn)單而又最重要的概念之一.在沒有自環(huán)和多重邊的無向無權(quán)貿(mào)易網(wǎng)絡(luò)中,某一節(jié)點(diǎn)的度數(shù)ki是與該節(jié)點(diǎn)有直接貿(mào)易關(guān)系的國(guó)家或地區(qū)的數(shù)目.

1.1.2 指數(shù)隨機(jī)圖模型

指數(shù)隨機(jī)圖模型(exponential random graph model,ERGM)旨在分析出影響網(wǎng)絡(luò)實(shí)現(xiàn)概率的特殊結(jié)構(gòu),找出較為重要的網(wǎng)絡(luò)局部關(guān)系構(gòu)建過程.給定一個(gè)有向網(wǎng)絡(luò)G(V,E),其中V為網(wǎng)絡(luò)中節(jié)點(diǎn)的集合,E為網(wǎng)絡(luò)中關(guān)系的集合.以集合J表示該網(wǎng)絡(luò)節(jié)點(diǎn)集V中所有節(jié)點(diǎn)對(duì)可能出現(xiàn)的關(guān)系,即J={(i,j)|i,j∈V,i≠j}.定義隨機(jī)變量Yij為關(guān)系變量,(i,j)∈E時(shí),Yij=1,(i,j)?E時(shí),Yij=0.由此構(gòu)成隨機(jī)的鄰接矩陣Y=[Yij],記y=[yij]是Y的一個(gè)特定實(shí)現(xiàn).ERGM的基本形式如下:

(2)

本文運(yùn)用指數(shù)隨機(jī)圖模型估計(jì)的步驟是首先選取合適的網(wǎng)絡(luò)統(tǒng)計(jì)量,然后采用馬爾科夫鏈蒙特卡羅極大似然估計(jì)法(MCMC MLE)進(jìn)行模型估計(jì),最后對(duì)模型展開有效性檢驗(yàn),直至模型的最終仿真網(wǎng)絡(luò)結(jié)果的各種結(jié)構(gòu)特征能夠與觀測(cè)網(wǎng)絡(luò)保持很好的一致性.

2 世界化妝品貿(mào)易網(wǎng)絡(luò)影響因素分析

2.1 變量選取與數(shù)據(jù)來源

依據(jù)ERGM的理論原理,本文以1996—2021年有向無權(quán)的世界化妝品貿(mào)易網(wǎng)絡(luò)鄰接矩陣中各個(gè)國(guó)家之間是否存在貿(mào)易關(guān)系為模型的被解釋變量,網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)、節(jié)點(diǎn)屬性特征和外部環(huán)境因素三大方面入手選取解釋變量.

2.1.1 網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)

互惠性(reciprocity)即雙向比,若其參數(shù)估計(jì)值顯著為正,則表明有向網(wǎng)絡(luò)中的節(jié)點(diǎn)間傾向于存在互惠性關(guān)系.幾何加權(quán)邊共享伙伴(gwesp)是指網(wǎng)絡(luò)中形似三角形狀的特殊結(jié)構(gòu),可以檢驗(yàn)國(guó)家之間的化妝品貿(mào)易是否存在集聚性和傳遞性.

2.1.2 節(jié)點(diǎn)屬性特征

本文以各國(guó)化妝品貿(mào)易合作伙伴數(shù)量(節(jié)點(diǎn)度數(shù))作為節(jié)點(diǎn)屬性變量,考察是否具有同配效應(yīng).節(jié)點(diǎn)度數(shù)可以反映一國(guó)在化妝品貿(mào)易網(wǎng)絡(luò)中的活躍度,因此可以探究活躍度相似的國(guó)家之間是否更易建立貿(mào)易關(guān)系.

2.1.3 外部環(huán)境因素

通過外部協(xié)變量網(wǎng)絡(luò)可以檢驗(yàn)網(wǎng)絡(luò)節(jié)點(diǎn)之間發(fā)生的關(guān)系是否建立在其他關(guān)系的基礎(chǔ)上.以下從地理、政治、經(jīng)濟(jì)、文化、貿(mào)易政策五大方面選取外部網(wǎng)絡(luò)指標(biāo),其具體含義和數(shù)據(jù)來源見表1,并截選部分?jǐn)?shù)據(jù)加以展示(表2).

表1 解釋變量的數(shù)據(jù)說明

表2 部分?jǐn)?shù)據(jù)示意

2.2 模型構(gòu)建

ERGM的基準(zhǔn)模型如公式(3)所示.它是簡(jiǎn)單隨機(jī)圖模型,僅包含弧數(shù)(Edges)統(tǒng)計(jì)項(xiàng),即網(wǎng)絡(luò)中節(jié)點(diǎn)之間有向邊的數(shù)量,其效應(yīng)如同線性回歸中的截距,可以為后續(xù)復(fù)雜模型有效性的檢驗(yàn)提供參照.

P(Y=y|θ)=exp (θ0Edges)/k.

(3)

經(jīng)過上文所述的變量選取和數(shù)據(jù)收集,本文最終建立的綜合模型如公式(4)所示.其中,各個(gè)參數(shù)所代表變量的含義見表1.

P(Y=y|θ)=exp (θ0Edges+θ1Mutual+θ2Gwesp+θaDegree+

θb1Dis+θb2Ctg+θb3Ins+θb4Eco+θb5Fta+θb6Lan+θb7Rlg)/k.

(4)

2.3 實(shí)證結(jié)果

本文基于1996年、2000年、2006年、2011年、2016年和2021年的數(shù)據(jù),借助R語言中的ergm包進(jìn)行ERGM擬合,估計(jì)結(jié)果見表3.

網(wǎng)絡(luò)內(nèi)生結(jié)構(gòu)變量中,互惠性指標(biāo)均在0.1%的水平上顯著為正,且估計(jì)系數(shù)偏高,說明化妝品貿(mào)易網(wǎng)絡(luò)中的國(guó)家之間傾向于存在雙邊貿(mào)易關(guān)系,許多國(guó)家互有進(jìn)出口貿(mào)易往來.而幾何加權(quán)邊共享伙伴對(duì)化妝品貿(mào)易網(wǎng)絡(luò)的作用方向和強(qiáng)度存在波動(dòng).

節(jié)點(diǎn)屬性特征中,從節(jié)點(diǎn)度數(shù)的同配效應(yīng)來看,其估計(jì)值均在0.1%的水平上顯著為正,說明度數(shù)相似的國(guó)家建立化妝品貿(mào)易關(guān)系的概率不大,即度數(shù)差距大的國(guó)家進(jìn)行化妝品貿(mào)易互動(dòng)的概率更大.

外部網(wǎng)絡(luò)協(xié)變量中,代表國(guó)家是否相鄰的變量Ctg和地理距離Dis的系數(shù)在0.1%的水平上分別顯著為正和顯著為負(fù),說明地理相鄰有助于促進(jìn)兩國(guó)成為化妝品貿(mào)易伙伴,這與傳統(tǒng)貿(mào)易理論的結(jié)論一致,各國(guó)更傾向與土地接壤或者距離較近的國(guó)家進(jìn)行化妝品貿(mào)易.

政治環(huán)境變量Ins的系數(shù)估計(jì)值一直在0.1%的水平上顯著為正,說明不同國(guó)家政治制度的差異對(duì)于兩國(guó)化妝品貿(mào)易聯(lián)系是否建立不構(gòu)成阻礙.

經(jīng)濟(jì)規(guī)模差異Eco基本均在0.1%的水平上顯著為負(fù),且其系數(shù)估計(jì)值在各個(gè)研究年份都最高,說明經(jīng)濟(jì)水平差距對(duì)化妝品貿(mào)易往來具有顯著的抑制作用,經(jīng)濟(jì)水平距離越大的國(guó)家越不易開展化妝品貿(mào)易互動(dòng),發(fā)展水平越相近的國(guó)家更易進(jìn)行化妝品貿(mào)易往來.

表示雙方是否簽訂貿(mào)易協(xié)定的變量Fta一直在0.1%的水平上顯著為正,而且其估計(jì)系數(shù)呈現(xiàn)下降趨勢(shì),說明簽訂了同一區(qū)域貿(mào)易協(xié)定的國(guó)家之間發(fā)生化妝品貿(mào)易的可能性更大.

語言Lan和宗教Rlg也總是顯著為正,說明使用同一官方語言、宗教背景越相近的國(guó)家進(jìn)行化妝品貿(mào)易的可能性更大.但是二者的顯著性水平在逐漸下降,體現(xiàn)出在全球融合的大趨勢(shì)下,語言和宗教這類文化因素對(duì)國(guó)際化妝品貿(mào)易的影響力逐漸減弱.

2.4 模型有效性檢驗(yàn)

為檢驗(yàn)?zāi)Ｐ偷挠行?以下采用兩種方法加以驗(yàn)證.

一方面,通過赤池信息準(zhǔn)則(AIC)和貝葉斯信息準(zhǔn)則(BIC)衡量模型擬合數(shù)據(jù)的優(yōu)良性.BIC略大于AIC,兩個(gè)數(shù)值越小,說明模型的擬合優(yōu)度越高.表4展示了基準(zhǔn)模型和綜合模型回歸結(jié)果的赤池信息量和貝葉斯信息量,通過對(duì)比可以看出,與基準(zhǔn)模型相比,綜合模型AIC和BIC的值均明顯降低,說明擬合效果優(yōu)良.

表4 模型擬合優(yōu)良性評(píng)價(jià)結(jié)果

另一方面,模型診斷能夠幫助判斷估計(jì)算法是否已經(jīng)收斂還是存在近似退化的問題,模型是否還需要進(jìn)行相應(yīng)的修改調(diào)整.以2021年為例,模型中的網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量在最后一次迭代期間所表現(xiàn)出的情況見圖1,可以看出模型中每一個(gè)網(wǎng)絡(luò)結(jié)構(gòu)統(tǒng)計(jì)量的圖表現(xiàn)出以0為中心隨機(jī)變化狀態(tài),說明模型達(dá)到了穩(wěn)定收斂的結(jié)果.

3 結(jié)語

本文基于社會(huì)網(wǎng)絡(luò)分析方法,對(duì)1996—2021年國(guó)際化妝品貿(mào)易的影響因素展開分析,得到世界化妝品貿(mào)易網(wǎng)絡(luò)受多種因素的共同影響.例如,互惠性指標(biāo)起顯著正向驅(qū)動(dòng)作用,兩國(guó)之間更傾向于建立雙邊貿(mào)易關(guān)系.世界化妝品貿(mào)易網(wǎng)絡(luò)呈現(xiàn)負(fù)向匹配特征,擁有較多貿(mào)易伙伴的國(guó)家傾向于與貿(mào)易伙伴數(shù)量少的國(guó)家建立化妝品貿(mào)易聯(lián)系.經(jīng)濟(jì)水平差距更有利于化妝品貿(mào)易關(guān)系的構(gòu)建等.