基于多元線性回歸分析方法的物流流通水平研究

潘新宇，薛 洋

(1.唐山學(xué)院 經(jīng)濟管理系，河北 唐山 063000；2.天津大學(xué) 管理與經(jīng)濟學(xué)部，天津 300072)

?

基于多元線性回歸分析方法的物流流通水平研究

潘新宇1，薛 洋2

(1.唐山學(xué)院 經(jīng)濟管理系，河北 唐山 063000；2.天津大學(xué) 管理與經(jīng)濟學(xué)部，天津 300072)

根據(jù)2013年全國31個省市的物流水平相關(guān)指標(biāo)，分別以經(jīng)營單位所在地分的進出口總額、按境內(nèi)目的地和貨源地分的商品進出口總額作為因變量對物流流通水平進行線性回歸分析，得到了多元線性回歸模型，由此找出了流通水平的影響因素，并有針對性地提出了對策建議。

多元線性回歸；物流；流通水平

0 引言

流通業(yè)的物流構(gòu)成消費品供應(yīng)鏈的終端環(huán)節(jié)，也是形成商品使用價值的最重要環(huán)節(jié)，涵蓋從商品采購到將商品銷售給消費者的整個物流全過程，具體包括商品采購、進貨運輸、商品儲存、加工配送、門店物流以及為各環(huán)節(jié)提供配套服務(wù)的信息系統(tǒng)等。中國流通產(chǎn)業(yè)競爭力較低一直是比較突出的問題。產(chǎn)業(yè)規(guī)模通常被用作衡量流通產(chǎn)業(yè)競爭力的基礎(chǔ)性指標(biāo)，它既反映目前流通產(chǎn)業(yè)的總量規(guī)模，又為未來的產(chǎn)業(yè)擴張?zhí)峁┢脚_基礎(chǔ)[1]?！吨袊B鎖零售企業(yè)經(jīng)營狀況分析報告》采用流通產(chǎn)業(yè)增加值、平均資產(chǎn)規(guī)模和社會消費品零售總額等指標(biāo)，綜合評價了中國各省區(qū)的流通產(chǎn)業(yè)規(guī)模競爭力。

物流是一個國家實現(xiàn)流通現(xiàn)代化的重要途徑，也是衡量一個國家和地區(qū)流通現(xiàn)代化程度的重要標(biāo)志。發(fā)展現(xiàn)代物流，實現(xiàn)流通現(xiàn)代化，是提升我國經(jīng)濟素質(zhì)和效益、降低成本、服務(wù)民生的迫切需要。物流在流通業(yè)的表現(xiàn)決定了流通業(yè)服務(wù)水平的高低和商品的最終價格。本文在宏觀上通過研究全國各地區(qū)2013年31個省市的物流水平相關(guān)指標(biāo)，對物流流通水平進行了線性回歸分析，其中能反映物流流通水平的相關(guān)因素，包括國內(nèi)生產(chǎn)總值、3種進出口總額(按經(jīng)營單位所在地分、按境內(nèi)目的地分和貨源分的外商投資企業(yè)貨物)、連鎖零售企業(yè)通體配送商品購進額、億元以上商品交易市場數(shù)量、貨運量、運輸業(yè)人員數(shù)量(鐵路、道路、水上、航空、裝卸搬運和運輸代理、郵政業(yè)、倉儲業(yè))、居民消費水平和各地區(qū)的人口數(shù)。

1 現(xiàn)有研究成果及相關(guān)理論

1.1 物流流通的概念

對于流通，文獻[2]指出，流通銜接生產(chǎn)和消費，包括兩個過程：人的轉(zhuǎn)移和商品的實質(zhì)轉(zhuǎn)移。田島義博[3]亦持類似的觀點，他認為商品流通就是商品的轉(zhuǎn)移過程，通常包含人格性轉(zhuǎn)移和物理性轉(zhuǎn)移。鈴木武[4]則宏觀地將流通界定為“制造商與消費者之間的架橋活動”。宋則[5]從微觀層面界定了流通，即“在實體經(jīng)濟范疇內(nèi)，由商品流通直接引起或與商品流通直接有關(guān)的商流、物流、信息流和資金流的總和”。對于流通產(chǎn)業(yè)，郭冬樂、宋則[6]指出，流通產(chǎn)業(yè)是指與商品流通和商業(yè)直接相關(guān)聯(lián)，或為商品流通和商業(yè)提供必要條件的各種投資領(lǐng)域的總稱，它主要包括農(nóng)產(chǎn)品、工業(yè)消費品和工業(yè)生產(chǎn)資料等商品的購銷體系、商業(yè)設(shè)施、倉儲業(yè)、運輸業(yè)、包裝裝卸、流通加工、流通技術(shù)以及與此相關(guān)的信息產(chǎn)業(yè)、服務(wù)業(yè)等。夏春玉[7]認為，流通就是專門以商品流通為經(jīng)營內(nèi)容的贏利性事業(yè)。

1.2 物流流通水平的影響因素

流通產(chǎn)業(yè)規(guī)模：產(chǎn)業(yè)規(guī)模是衡量流通產(chǎn)業(yè)發(fā)展水平高低的基石。

流通結(jié)構(gòu)：主要反映流通產(chǎn)業(yè)的要素內(nèi)部結(jié)構(gòu)狀況。包括人員組成、零售額比重，等等。

流通產(chǎn)業(yè)效率：反映流通產(chǎn)業(yè)運行效率和發(fā)展?fàn)顩r的核心指標(biāo)。包括交易市場等信息。

流通產(chǎn)業(yè)貢獻：改革開放以來，流通產(chǎn)業(yè)逐漸成為支撐國民經(jīng)濟或區(qū)域經(jīng)濟運行的基礎(chǔ)部門，它決定著其他產(chǎn)業(yè)部門經(jīng)濟活動的發(fā)展水平。包括進出口額等。

2 數(shù)據(jù)收集和解析

2.1 宏觀統(tǒng)計數(shù)據(jù)收集和解析

通過《中國統(tǒng)計年鑒》中的數(shù)據(jù)，搜集到了2013年31個省市的生產(chǎn)總值、貨物進出口總額、外商投資企業(yè)貨物進出口總額、連鎖零售企業(yè)基本情況、億元以上商品交易市場基本情況、貨運量、交通運輸業(yè)人員數(shù)、倉儲業(yè)人員數(shù)、郵政業(yè)人員數(shù)、各地區(qū)居民消費水平、各地區(qū)人口數(shù)等數(shù)據(jù)。

將所有數(shù)據(jù)按順序整理到一個大表格中，用于之后的數(shù)據(jù)分析。數(shù)據(jù)解析如下：

生產(chǎn)總值：各地區(qū)經(jīng)濟情況的度量，是指在一定時期內(nèi)(一個季度或一年)，一個地區(qū)的經(jīng)濟中所生產(chǎn)出的全部最終產(chǎn)品和勞務(wù)的價值。

貨物進出口總額：各地區(qū)實際進出我國國境的貨物總金額。

外商投資企業(yè)貨物進出口總額：各地區(qū)外商投資企業(yè)實際進出我國國境的貨物總金額。

連鎖零售企業(yè)基本情況：各地區(qū)連鎖零售企業(yè)數(shù)量、營業(yè)額等。

億元以上商品交易市場基本情況：各地區(qū)消費品綜合市場、農(nóng)產(chǎn)品市場、工業(yè)消費品市場、生產(chǎn)資料綜合市場、工業(yè)生產(chǎn)資料市場、農(nóng)業(yè)生產(chǎn)資料市場的總數(shù)量。

貨運量：各地區(qū)運輸企業(yè)在一定的時期內(nèi)實際運送的貨物數(shù)量。

交通運輸：各地區(qū)海陸空運輸業(yè)、倉儲和郵政業(yè)人員數(shù)量。

居民消費水平：各地區(qū)居民在物質(zhì)產(chǎn)品和勞務(wù)的消費過程中，對滿足人們生存、發(fā)展和享受需要方面所達到的程度。

人口數(shù)：各地區(qū)人口數(shù)量。

2.2 各地區(qū)物流流通水平多元線性回歸分析

2.2.1 模型的建立

進出口是體現(xiàn)流通水平的一個重要方面，所以在模型中，分別將按經(jīng)營單位所在地分的進出口總額、按境內(nèi)目的地和貨源地分的進出口總額，作為因變量進行線性回歸分析。

Y=β0+β1X1+β2X2+β3X3+β4X4+β5X5+β6X6+β7X7+β8X8+β9X9+β10X10+β11X11+β12X12+β13X13+β14X14。

變量的定義如下：

Y1：按經(jīng)營單位所在地分的進出口總額(萬美元)；

Y2：按境內(nèi)目的地和貨源地分的進出口總額(萬美元)；

X1：外商投資企業(yè)貨物進出口總額(萬美元)；

X2：連鎖零售企業(yè)通體配送商品購進額(億元)；

X3：億元以上商品交易市場數(shù)量(個)；

X4：貨運量(萬噸)；

X5：鐵路運輸業(yè)人員數(shù)量(人)；

X6：道路運輸業(yè)人員數(shù)量(人)；

X7：水上運輸業(yè)人員數(shù)量(人)；

X8：航空運輸業(yè)人員數(shù)量(人)；

X9：裝卸搬運和運輸代理業(yè)人員數(shù)量(人)；

X10：郵政業(yè)人員數(shù)量(人)；

X11：倉儲業(yè)人員數(shù)量(人)；

X12：國內(nèi)生產(chǎn)總值(億元)；

X13：居民消費水平(元)；

X14：人口數(shù)(人)。

2.2.2Y1的多元線性回歸分析

2.2.2.1 參數(shù)估計

采用強行進入法進行參數(shù)的最小二乘估計，求得回歸系數(shù)如表1所示。

表1 回歸系數(shù)

因變量：按經(jīng)營單位所在地分進出口總額(萬美元)。

初步得到本問題的多元線性回歸模型為：

Y1=9 087 629.563+0.891X1-177.301X2+23 513.380X3+15.816X4-30.797X5+30.555X6-183.836X7+183.065X8+270.562X9+89.168X10+300.212X11+620.642X12-789.551X13-526.122X14。2.2.2.2 三種檢驗

檢驗結(jié)果如表2所示。因文章篇幅關(guān)系，相關(guān)檢驗圖表省略。

回歸方程的擬合度檢驗：復(fù)相關(guān)系數(shù)R、判定系數(shù)、調(diào)整判定系數(shù)、估計值誤差和DW，這些值反映了自變量解釋因變量的能力。由表2可知R的平方是0.989，調(diào)整后的R的平方等于0.979，即判定系數(shù)和調(diào)整后的判定系數(shù)都很接近于1，故認定其通過了擬合度檢驗，且模型擬合度很好。

回歸關(guān)系的整體性檢驗(F檢驗)：檢驗因變量的方差來源、方差平方和、自由度、均方、F檢驗統(tǒng)計量的觀測值和顯著性水平。得出F=92.605，且顯著性水平是0.000。因為本模型的F統(tǒng)計量服從F(14，16)分布，由此查表得臨界值F(14，16)≈3.4，而本模型的F值92.605大于此臨界值，且顯著性水平0.000小于0.05，故拒絕原假設(shè)，則表明此回歸模型的整體回歸模式顯著，具有統(tǒng)計上的意義。

回歸系數(shù)的顯著性檢驗(T檢驗)：本模型中樣本數(shù)是31個，自變量數(shù)是14個，則本模型各回歸系數(shù)的T統(tǒng)計量應(yīng)服從T(16)的分布，由檢驗結(jié)果得到其臨界值約為2.12，而本模型的15個回歸系數(shù)的T值有的大于此臨界值，有的則小于臨界值，所以沒有全部通過T檢驗。

2.2.2.3 多重共線性分析

通過自變量的方差膨脹因子分析的方法判斷自變量間是否存在共線性。

容許度在0～1之間變化，其值越接近于0，說明自變量間的共線性越強；其值越接近于1，說明自變量間的共線性越弱。所以，若方差膨脹因子VIF接近于1，則表明自變量間的共線性較弱；VIF大于10，則表明自變量之間存在多重共線性。

由檢驗結(jié)果可知，其方差膨脹因子值有的大于10，有的小于10，故本模型中的自變量之間存在多重共線性。

2.2.2.4 殘差的回歸假設(shè)檢驗

(1)殘差的自相關(guān)分析。表2為Y1的模型摘要。由表2可知，此回歸模型的DW值為1.962。而其樣本個數(shù)是31，自變量個數(shù)是14，顯著水平為0.05。由于解釋變量數(shù)目過多，查不到相關(guān)數(shù)據(jù)。但是可見DW值接近2，因此可以認為沒有自相關(guān)性。

(2)殘差值統(tǒng)計量。表3為殘差統(tǒng)計。從表3中可以得到此回歸模型殘差值的描述性統(tǒng)計量，包括預(yù)測值、殘差、標(biāo)準(zhǔn)預(yù)測值和標(biāo)準(zhǔn)殘差等。從中可以看出，最小和最大殘差為-5.491 445 500 000 0E6和4.623 251 500 000 0E6，且標(biāo)準(zhǔn)化后，其均在3倍的標(biāo)準(zhǔn)差以內(nèi)，故此模型殘差值的正態(tài)分布性較好。

表2 Y1的模型摘要

表3 殘差統(tǒng)計

因變量：按經(jīng)營單位所在地分進出口總額(萬美元)。

(3)標(biāo)準(zhǔn)化殘差直方圖。本模型的標(biāo)準(zhǔn)化殘差直方圖如圖1所示，這是為了觀測模型的殘差是否符合正態(tài)分布。由圖1可知，本模型標(biāo)準(zhǔn)化殘差的分布結(jié)果基本符合正態(tài)分布的假設(shè)，且殘差值都在3個標(biāo)準(zhǔn)差內(nèi)。

圖1 標(biāo)準(zhǔn)化殘差直方圖

2.2.2.5 逐步修正法

因為原模型并不是所有系數(shù)都通過了T檢驗，所以對模型進行逐步回歸，得到表4。

得到3個回歸模型：

Y1=2 638 292.88+1.753X1；

Y1=-2 009 306.170+1.381X1+194.345X10；

Y1=-1.072E7+1.14X1+218.808X10+574.69X13。

對上述3個模型進行檢驗，結(jié)果如表5，表6所示，從中可以看出，模型1，模型2，模型3都通過了T檢驗和F檢驗。

從表6擬合優(yōu)度R平方來看，模型3的擬合優(yōu)度最高，其次是模型2，模型1。

表4 修正后的模型系數(shù)

因變量：按經(jīng)營單位所在地分進出口總額(萬美元)。

表5 Y1的ANOVA模型

表6 修正后的Y1模型摘要

2.2.3 Y2的多元線性回歸分析

本節(jié)根據(jù)“2.2.2 Y1的多元線性回歸分析”同理可得4個回歸模型：

Y2=1 136 849.844+1.986X1；

Y2=-1 572 064.174+1.858X1+19 697.532X3；

Y2=-4 026 313.219+1.673X1+17 061.435X3+109.158X13；

Y2=-1 503 281.847+1.726X1+20 772.420X3+200.505X6-52.332X9。

檢驗結(jié)果如表7，表8所示。從擬合優(yōu)度來看，模型4的擬合優(yōu)度最高，其次是模型3，模型2，模型1。

3 結(jié)論

(1)從分析結(jié)果可以看出，外商投資企業(yè)貨物進出口總額、郵政業(yè)人員數(shù)量、居民消費水平這幾個指標(biāo)在14個指標(biāo)中對按經(jīng)營單位所在地分的進出口總額有較大影響。所以按經(jīng)營單位所在地劃分中，關(guān)于居民日常的貨物流通較為顯著。因此，在進出口時，對于日常類的流通應(yīng)多加管理和監(jiān)控，運用物流管理知識合理規(guī)劃運輸和倉儲，從而提升流通效率。

(2)外商投資企業(yè)貨物進出口總額、億元以上商品交易市場數(shù)量、道路運輸業(yè)人員數(shù)量這幾個指標(biāo)在14個指標(biāo)中對按境內(nèi)目的地分和貨源地分的進出口總額有較大影響。所以按境內(nèi)目的地和貨源地劃分中，關(guān)于工業(yè)、商業(yè)等的大中型非日常的流通較為顯著。因此加強各大企業(yè)和政府部門的管理和監(jiān)控是很有必要的。

表7 Y2的ANOVA模型

表8 模型摘要

(3)在不同的運輸方式中，對進出口的影響并不相同，有的呈現(xiàn)正相關(guān)，有的呈現(xiàn)負相關(guān)。例如，對于按境內(nèi)目的地和貨源地劃分的進出口中，裝卸搬運和運輸代理人員數(shù)量和其呈負相關(guān)，說明裝卸搬運和運輸代理人員數(shù)量并不能提高此類進出口的效率。但在各類運輸方式都是必需的情況下，應(yīng)合理分配運輸人員，以提高流通的整體效率。

[1] 中國流通產(chǎn)業(yè)發(fā)展報告課題組.我國流通產(chǎn)業(yè)的整體競爭力發(fā)展亟待提速[N].經(jīng)濟日報，2005-09-05.

[2] 劉向陽.20世紀(jì)營銷學(xué)的范式演變與特點[J].科技進步與對策，2002(7)：74-76.

[3] 田島義博.流通的活力[M].于淑華，譯.北京：中國商業(yè)出版社，2000.

[4] 鈴木武.現(xiàn)代流通政策和課題[M].王哲，譯.北京：中國商業(yè)出版社，1993.

[5] 宋則.流通現(xiàn)代化及流通業(yè)競爭力研究(上)[J].商業(yè)時代，2006(4)：11-13.

[6] 郭冬樂，宋則.中國商業(yè)理論前沿(二)[M].北京：社會科學(xué)文獻出版社，2001.

[7] 夏春玉.商品流通產(chǎn)業(yè)的歷史發(fā)展及其理論解說[J].財經(jīng)問題研究，1998(2)：71-76.

(責(zé)任編校：夏玉玲)

A Research on Logistical Circulation Efficiency Based on Multivariate Linear Regression Analysis

PAN Xin-yu1， XUE Yang2

(1.Department of Economics and Management， Tangshan University， Tangshan 063000， China； 2.College of Management and Economics， Tianjin University， Tianjin 300072， China)

By studying the logistics-related indexes of the 31 provinces in 2013，the authors of this paper have carried out a linear regression analysis of the logistical distribution in these places，with the total value of imports and exports by business location and the total amount of imports and exports of goods by destination and source location as the dependent variable， obtained the multiple linear regression model， discovered the influencing factors in the circulation efficiency， and put forward the corresponding suggestions.

multivariate linear regression analysis； logistics； circulation efficiency

潘新宇(1982-)，男，河北唐山人，講師，博士，主要從事物流與供應(yīng)鏈管理研究。

F25

A

1672-349X(2016)06-0072-06

10.16160/j.cnki.tsxyxb.2016.06.019