楊 宏

（牡丹江師范學(xué)院 經(jīng)濟(jì)與管理學(xué)院，黑龍江 牡丹江 157000）

0 引言

隨著社會(huì)經(jīng)濟(jì)水平的發(fā)展，越來越多的人享受外出旅行，旅游市場也隨之蓬勃發(fā)展。生態(tài)環(huán)境既是安居樂業(yè)的重要保證，同時(shí)也對(duì)旅游產(chǎn)業(yè)的發(fā)展有著重要影響，分析生態(tài)環(huán)境與旅游產(chǎn)業(yè)發(fā)展間的相互影響對(duì)改善生態(tài)環(huán)境、協(xié)調(diào)旅游行業(yè)健康發(fā)展具有重要意義，不少學(xué)者對(duì)此進(jìn)行了研究，獲得了豐富的研究成果[1-6]。

通過對(duì)目前生態(tài)環(huán)境與旅游產(chǎn)業(yè)發(fā)展相關(guān)的研究發(fā)現(xiàn)，關(guān)注重點(diǎn)主要集中在旅游開發(fā)與環(huán)境保護(hù)的關(guān)系、各區(qū)域旅游產(chǎn)業(yè)發(fā)展與生態(tài)環(huán)境的耦合關(guān)系等方面，關(guān)于生態(tài)環(huán)境對(duì)旅游產(chǎn)業(yè)的影響和影響程度的相關(guān)研究則比較少。本文則以此為切入點(diǎn)，基于門檻模型實(shí)證分析生態(tài)環(huán)境與旅游經(jīng)濟(jì)產(chǎn)業(yè)發(fā)展的門檻效應(yīng)，為生態(tài)環(huán)境保護(hù)和旅游規(guī)劃提供進(jìn)一步的理論指導(dǎo)。

1 分析方法選擇

1.1 門檻效應(yīng)初步分析

不同的地區(qū)生態(tài)環(huán)境水平具有一定的差異，而生態(tài)環(huán)境水平與旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展二者之間的相互影響可能不是嚴(yán)格的線性關(guān)系。當(dāng)前，旅游市場蓬勃發(fā)展，生態(tài)環(huán)境綜合水平較高的地方旅游產(chǎn)業(yè)經(jīng)濟(jì)可能呈現(xiàn)良好的發(fā)展態(tài)勢，而生態(tài)環(huán)境綜合水平中等或偏低的地區(qū)旅游產(chǎn)業(yè)經(jīng)濟(jì)也可能發(fā)展不錯(cuò)。具體而言，生態(tài)環(huán)境對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展產(chǎn)生的很可能不是線性影響，而是表現(xiàn)出相應(yīng)的門檻效應(yīng)；當(dāng)生態(tài)環(huán)境綜合水平超過一定的門檻值時(shí)，其對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展造成的影響會(huì)出現(xiàn)減弱或逆轉(zhuǎn)。

門檻效應(yīng)通常的定義如下：在定義的門檻變量區(qū)間內(nèi)，解釋變量對(duì)被解釋變量的影響會(huì)隨著區(qū)間的不同呈現(xiàn)不同的作用程度和作用方向[7]。在旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的因素當(dāng)中，生態(tài)環(huán)境是十分重要的一個(gè)方面，因此，本文將從門檻效用角度出發(fā)，研究生態(tài)環(huán)境綜合水平在不同的門檻區(qū)間內(nèi)對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展產(chǎn)生影響的程度和大小。需要指出的是，本文做了如下假設(shè)：生態(tài)環(huán)境綜合水平對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的影響不是線性的，而是呈現(xiàn)出門檻效應(yīng)，門檻值的個(gè)數(shù)可能是一個(gè)、兩個(gè)或三個(gè)[8]。為了對(duì)上述假設(shè)進(jìn)行驗(yàn)證，本文將計(jì)算我國30個(gè)省份的生態(tài)環(huán)境綜合水平，并分析其對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的門檻效應(yīng)。如果假設(shè)成立，將確定門檻值的個(gè)數(shù)和門檻值的大小，并揭示生態(tài)環(huán)境綜合水平對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展在不同門檻區(qū)間內(nèi)影響的大小和方向。

1.2 模型選擇

不同模型具有各自的適用范圍，簡單將存在結(jié)構(gòu)突變的非線性模型做簡化處理（如假設(shè)結(jié)構(gòu)突變不存在）進(jìn)行線性分析，可能會(huì)造成模型不能完全呈現(xiàn)各變量間的相互作用，甚至造成模型不適用。對(duì)此，有不少研究人員對(duì)模型結(jié)構(gòu)變化的問題進(jìn)行了分析，改進(jìn)的方法主要有以下幾種：采取分組回歸的形式[9]、增加多項(xiàng)式項(xiàng)數(shù)強(qiáng)化解釋變量的說明[10]、設(shè)置交乘項(xiàng)和虛擬變量[11]等。盡管上述幾種方法有所改進(jìn)，但還是存在著不同程度的不足，如分組回歸或交乘項(xiàng)設(shè)置難以確定結(jié)構(gòu)突變點(diǎn)的數(shù)值大小，進(jìn)而造成最終的分析結(jié)果出現(xiàn)差異。而增加多項(xiàng)式的項(xiàng)數(shù)則會(huì)導(dǎo)致分析模型出現(xiàn)多重共線性問題[12]。

通過上文的分析可知，生態(tài)環(huán)境綜合水平與旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展之間的相互作用很有可能是非線性的，表現(xiàn)出明顯的區(qū)間效應(yīng)[13]。為了消除分組回歸方法中人為劃分效應(yīng)區(qū)間所產(chǎn)生的偏差，本文將利用一種數(shù)據(jù)特征劃分效應(yīng)區(qū)間的方法——面板門檻回歸模型，確定生態(tài)環(huán)境綜合水平對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的區(qū)間，并分析其作用方式。

1.3 門檻值估計(jì)及檢驗(yàn)

確定分析模型的種類之后，需要對(duì)門檻值的大小進(jìn)行確定，即門檻值預(yù)估。下面將以單一門檻模型為例，對(duì)其門檻值的預(yù)估和檢驗(yàn)方法進(jìn)行介紹，雙重門檻及三重門檻的門檻值確定方法類似。

門檻值的預(yù)估方法遵從“數(shù)據(jù)對(duì)話”的原則，對(duì)單一門檻而言，首先需要在不同的層次上設(shè)定結(jié)構(gòu)變化點(diǎn)，這一過程會(huì)產(chǎn)生如0或1的虛擬變量，這樣模型就會(huì)變成常見的線性形式，只是模型中增加了一些交乘項(xiàng)；接著就可以使用普通的分析方法（如最小二乘法）進(jìn)行估計(jì)，根據(jù)估算產(chǎn)生的RSS（殘差平方和）及R2（擬合優(yōu)度）結(jié)果，擬合程度和誤差最小的那個(gè)點(diǎn)就是模型中結(jié)構(gòu)突變點(diǎn)[14]。而R2在突變點(diǎn)附近變化不顯著，因此，模型分析中更多選擇RSS作為擬合結(jié)果的依據(jù)，RSS最小值對(duì)應(yīng)的變量數(shù)值即為門檻變量值。整個(gè)分析過程中采用網(wǎng)格搜索進(jìn)行。

上述方法對(duì)于單門檻非常適用，但對(duì)于雙重門檻或者三門檻而言，門檻值采用網(wǎng)格搜索的方式則會(huì)由于嵌套循環(huán)的存在使得計(jì)算量急劇增加，造成估值困難；而對(duì)于大數(shù)據(jù)集的分析，采用常規(guī)的網(wǎng)格搜索方法幾乎不可能。對(duì)于這種情況，Hansen于1999年提出了一種改進(jìn)的搜索方法，此方法計(jì)算處理量幾何級(jí)的增大，因此能夠很好地解決傳統(tǒng)網(wǎng)格搜索存在的不足。

通過上述方式獲得粗略的門檻值后需要進(jìn)行兩方面的檢驗(yàn)。第一個(gè)檢驗(yàn)是效應(yīng)的顯著性檢驗(yàn)，即確定模型是否具有非線性特征；第二個(gè)檢驗(yàn)是真實(shí)性檢驗(yàn)，即檢查預(yù)估的門檻值與真實(shí)值存在的偏差。

F統(tǒng)計(jì)量是門檻顯著性檢驗(yàn)常用的方法，以單門檻檢驗(yàn)為例進(jìn)行說明，雙重門檻和三門檻模型的檢驗(yàn)方法類似。假設(shè)β1=β2，對(duì)應(yīng)的計(jì)算公式為：

其中，S0表示的時(shí)沒有門檻值時(shí)的RSS（殘差平方和），Si(γ)表示的是單一門檻情況下的RSS，δ2表示的是殘差值與自由度的比值，即δ2=RSS/n（t-1）

考慮到F統(tǒng)計(jì)量分布的不確定性，本文中檢驗(yàn)方法為自抽樣檢驗(yàn)，即對(duì)分析樣本進(jìn)行數(shù)次重復(fù)性的抽樣分析以得到F統(tǒng)計(jì)量，最后獲得P值的大小。檢驗(yàn)過程中，如果假設(shè)不成立，則意味著門檻效應(yīng)明顯，表示生態(tài)環(huán)境對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展呈現(xiàn)明顯的門檻作用。

門檻值真實(shí)性檢驗(yàn)就是確定預(yù)估的與實(shí)際的門檻值差別，檢驗(yàn)常用的方法的是LRtest（似然比統(tǒng)計(jì)量），初始的假設(shè)是，對(duì)應(yīng)的LRtest計(jì)算公式為：

其中，Si(γ，)表示的是其他門檻值計(jì)算得到的殘差平方和數(shù)值，Si(γ)實(shí)際門檻值得到的RSS數(shù)值，δ2與F統(tǒng)計(jì)量中的定義類似，即殘差值與自由度的比值。

同樣的，考慮到似然比統(tǒng)計(jì)量的非標(biāo)準(zhǔn)性，通常利用的是Hansen創(chuàng)建的非拒絕域模型來進(jìn)行檢驗(yàn)，定義為：

其中，α表示的是顯著性水平，如5%水平上相應(yīng)的臨界值大小是7.35。如果LR的計(jì)算值在各自水平所對(duì)應(yīng)的臨界值內(nèi)，則初始假設(shè)成立，意味著預(yù)估門檻值是真實(shí)準(zhǔn)確的。

2 門檻效應(yīng)分析

2.1 數(shù)據(jù)來源及處理

本文中30個(gè)省份（不含西藏和港、澳、臺(tái)）的數(shù)據(jù)整理自各權(quán)威門戶、統(tǒng)計(jì)年鑒等。網(wǎng)站包括國家文物局官網(wǎng)、國家旅游局官網(wǎng)，年鑒數(shù)據(jù)包括各省市區(qū)《國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)（2005—2015）》《中國旅游統(tǒng)計(jì)年鑒（2005—2015）》《中國統(tǒng)計(jì)年鑒（2005—2015）》，對(duì)于缺失的小部分采取均值代替或者是插值處理。

生態(tài)環(huán)境是一個(gè)綜合性的影響因素，涵蓋不同的指標(biāo)，本文中采取的是主成分分析法對(duì)生態(tài)環(huán)境進(jìn)行表征，首先基于SPSS平臺(tái)對(duì)整理得到的30個(gè)省份相關(guān)數(shù)據(jù)進(jìn)行主成分分析，獲取能夠表征生態(tài)環(huán)境的有效成分，然后計(jì)算得到各主成分的綜合值大小來衡量生態(tài)環(huán)境綜合水平，具體的計(jì)算公式為：

其中，左邊ecοli表示的是綜合得分；右邊式子中，i表示某省份，j表示主成分，F(xiàn)ij是j主因素上對(duì)應(yīng)的省市得分，Cij為對(duì)應(yīng)的權(quán)重值大小。

2.2 直觀分析

為了使生態(tài)環(huán)境的解釋變量更真實(shí)準(zhǔn)確，本文在參考大量相關(guān)文獻(xiàn)的基礎(chǔ)上主要從生態(tài)環(huán)境壓力、建設(shè)、及狀態(tài)三個(gè)方面表征生態(tài)環(huán)境水平，進(jìn)而構(gòu)建30個(gè)省份的生態(tài)環(huán)境評(píng)價(jià)體系，具體評(píng)價(jià)體系如表1所示。

表1 生態(tài)環(huán)境評(píng)價(jià)體系

需要說明的是，研究過程中對(duì)性質(zhì)為負(fù)的指標(biāo)做倒數(shù)處理，并由SPSS軟件進(jìn)行相應(yīng)的數(shù)據(jù)處理。

通過對(duì)表1的生態(tài)環(huán)境各個(gè)具體因素進(jìn)行主成分分析，可以獲得2005—2015年30個(gè)省份生態(tài)環(huán)境綜合水平初步評(píng)價(jià)得分；同時(shí)，還對(duì)初步處理得到的數(shù)據(jù)進(jìn)行了基于3σ原則的坐標(biāo)移動(dòng)處理，以滿足后面門檻模型分析過程中的需要。本文中以2005年及2015年的生態(tài)環(huán)境水平為例進(jìn)行說明，相關(guān)數(shù)據(jù)結(jié)果匯總?cè)绫?所示。

表2 2005年和2015年國內(nèi)30個(gè)省份的生態(tài)環(huán)境水平結(jié)果

由表2的生態(tài)環(huán)境水平綜合得分可以看出，2005年30個(gè)省份生態(tài)環(huán)境前五依次是海南、青海、北京、江蘇及浙江；同樣，2015年30個(gè)省份生態(tài)環(huán)境前五依次是海南、山東、江蘇、北京、重慶。從整體情況來看，不少省份的生態(tài)環(huán)境水平出現(xiàn)了下降；按照中部、東部、西部三個(gè)區(qū)域劃分來看，中部大部分省份的生態(tài)環(huán)境水平得到了提高，僅有黑龍江、吉林兩個(gè)省份出現(xiàn)了不同程度的下滑，而中部地區(qū)各個(gè)省份的生態(tài)環(huán)境均以上升為主，東部各省份則出現(xiàn)了較多地區(qū)的生態(tài)環(huán)境下降，值得注意。

2.3 面板門檻模型分析

單門檻模型的計(jì)算通式為：

其中，qit表示的是門檻變量，xit表示的是解釋變量，yit表示的是被解釋變量，ui表示每個(gè)因素所產(chǎn)生的影響，εit表征的是誤差項(xiàng)或干擾值；I(.)是一個(gè)條件函數(shù)，當(dāng)門檻條件達(dá)到時(shí)為1，沒有達(dá)到時(shí)則為0；i＝1，2…，N ，表示的是不同因素影響，；t表示的是時(shí)間，取值范圍為1,2，…，T。需要說明的是，雙重門檻或三門檻模型的計(jì)算公式都是基于此公式衍生得到。

本文中門檻變量為生態(tài)環(huán)境綜合水平，而影響旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的因素眾多，本文綜合各因素選取的解釋變量包括區(qū)域旅游資源情況、區(qū)域交通設(shè)施建設(shè)情況、區(qū)域旅游服務(wù)設(shè)施配套以及區(qū)域經(jīng)濟(jì)發(fā)展水平等方面。

基于stata12.1軟件平臺(tái)對(duì)門檻模型分析以確定生態(tài)環(huán)境對(duì)產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的門檻效應(yīng)類型，檢驗(yàn)結(jié)果匯總見下頁表3。可以看出，5%水平下，單門檻模型和雙重門檻模型均通過了顯著性檢驗(yàn)，而三門檻模型沒有通過，所以本文中的門檻效應(yīng)類別確定為雙重門檻。

確定門檻模效應(yīng)的門檻值個(gè)數(shù)之后，接著需要確定兩個(gè)門檻值的大小并對(duì)其真實(shí)性進(jìn)行檢驗(yàn)，具體計(jì)算過程由軟件完成，相關(guān)計(jì)算結(jié)果匯總見下頁表4。

表3 門檻檢驗(yàn)結(jié)果

表4 門檻值預(yù)計(jì)大小

根據(jù)表4中的門檻值計(jì)算結(jié)果和真實(shí)性檢驗(yàn)，可以將生態(tài)環(huán)境綜合水平對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的門檻效應(yīng)劃為3個(gè)區(qū)間進(jìn)行分析，具體的計(jì)算結(jié)果如表5所示。

表5 門檻效應(yīng)分析結(jié)果

從表5可以看出，在旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的幾個(gè)變量之中，經(jīng)濟(jì)發(fā)展情況和旅游資源量對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展具有積極促進(jìn)作用；盡管交通設(shè)施建設(shè)情況及旅游設(shè)施建設(shè)這兩個(gè)方面對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展也具有正面促進(jìn)作用，但在統(tǒng)計(jì)學(xué)上不明顯，即影響不顯著。同時(shí)可以看到，D08的系數(shù)大小為-0.122，表明我國旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展在2008年出現(xiàn)了波動(dòng)，很可能是由汶川地震、經(jīng)濟(jì)危機(jī)等因素造成。另外，在分析過程中對(duì)全部數(shù)據(jù)都進(jìn)行了對(duì)數(shù)化，得到的系數(shù)能夠直接進(jìn)行比較說明?？梢钥吹?，區(qū)域經(jīng)濟(jì)發(fā)展對(duì)應(yīng)的系數(shù)為0.231，而旅游資源量對(duì)應(yīng)的系數(shù)為0.433，意味著豐富地區(qū)旅游資源比發(fā)展地區(qū)經(jīng)濟(jì)更有利于旅游產(chǎn)業(yè)經(jīng)濟(jì)的發(fā)展，在某種程度上解釋了為什么國內(nèi)各省市爭相申遺、加快旅游品牌建設(shè)以豐富旅游資源。

在三個(gè)門檻區(qū)間內(nèi)，生態(tài)環(huán)境對(duì)旅游經(jīng)濟(jì)發(fā)展的回歸系數(shù)分別為0.22、0.331、0.249，說明我國旅游產(chǎn)業(yè)經(jīng)濟(jì)隨著生態(tài)環(huán)境綜合水平的提高先增加后降低，關(guān)系曲線呈現(xiàn)倒“V”。具體來說，當(dāng)生態(tài)環(huán)境綜合水平門檻值大于-1.109時(shí)，生態(tài)環(huán)境綜合水平對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的促進(jìn)作用下降，意味著在該水平上單純的提高生態(tài)環(huán)境不能提高旅游產(chǎn)業(yè)經(jīng)濟(jì)，表明此時(shí)旅游產(chǎn)業(yè)經(jīng)濟(jì)的發(fā)展需要考慮改善或提高其他方面，同時(shí)注重旅游產(chǎn)業(yè)自身質(zhì)量的提高；而當(dāng)生態(tài)環(huán)境綜合水平門檻值小于-1.109時(shí)，提高生態(tài)環(huán)境綜合水平能夠助力旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展。因此，為了更好地引導(dǎo)旅游產(chǎn)業(yè)的經(jīng)濟(jì)發(fā)展，需要提高生態(tài)環(huán)境水平，同時(shí)還需要綜合考慮其他影響旅游產(chǎn)業(yè)發(fā)展的因素。

3 結(jié)論

本文從理論層面分析生態(tài)環(huán)境和旅游經(jīng)濟(jì)產(chǎn)業(yè)發(fā)展相互作用關(guān)系出發(fā)，介紹了面板門檻模型等基本概念，引用我國30個(gè)省份2005年和2015年生態(tài)環(huán)境相關(guān)數(shù)據(jù)，基于主成分分析方法確定了生態(tài)環(huán)境的各主要影響因素，計(jì)算得到了他們的生態(tài)環(huán)境綜合水平，并基于面板門檻模型進(jìn)行了門檻效應(yīng)分析。結(jié)果表明，30個(gè)省份生態(tài)環(huán)境綜合水平呈現(xiàn)明顯的地區(qū)差異，東部及西部區(qū)域的生態(tài)環(huán)境水平要高于中部省份，但隨著中部各省份的生態(tài)環(huán)境水平提高，差距略有縮??；30個(gè)省份生態(tài)環(huán)境對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的作用符合雙重門檻模型，門檻值為-1.152、-1.109，在門檻區(qū)間分析中，隨著生態(tài)環(huán)境水平的提高，其對(duì)旅游產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的作用先增強(qiáng)后減弱，作用曲線為倒“V”。對(duì)于各省份而言，既需要根據(jù)自身?xiàng)l件進(jìn)一步改善生態(tài)環(huán)境水平，同時(shí)還需要從其他方面入手，以保證旅游產(chǎn)業(yè)經(jīng)濟(jì)的持續(xù)健康發(fā)展。