樓文高,樓際通,宋雷娟,王浪慶

(1. 上海商學(xué)院 財經(jīng)學(xué)院，上海 200235；2. 北卡羅來納大學(xué)教堂山分校, NC 27514 )

?

中小企業(yè)稅務(wù)稽查投影尋蹤建模與實證分析

樓文高1,樓際通2,宋雷娟1,王浪慶1

(1. 上海商學(xué)院 財經(jīng)學(xué)院，上海 200235；2. 北卡羅來納大學(xué)教堂山分校, NC 27514 )

從上海市某區(qū)386家中小企業(yè)申報的15項稅收指標(biāo)數(shù)據(jù)中篩選出對判定企業(yè)納稅情況具有重要影響的10個評價指標(biāo)，并將全部386個樣本分成性質(zhì)相似的建模樣本和測試樣本(其中測試樣本個數(shù)占45%)，建立了基于投影尋蹤分類(PPC)技術(shù)的稅務(wù)稽查評價模型．與多元線性回歸(MLR)、判別分析(MDA)、Logistic和支持向量機(SVM)模型相比，PPC模型的識別錯誤率最低，建模樣本和測試樣本的平均分類錯誤率低于6%，改進型PPC模型包含的評價指標(biāo)少，兩類錯誤率很接近，非常適用于實際企業(yè)的稅務(wù)稽查評估研究和實踐.對339家待判斷企業(yè)納稅情況的判定結(jié)果研究表明，建立的改進型PPC模型具有很好的泛化能力和魯棒性.

稅務(wù)稽查；投影尋蹤分類技術(shù)；分類錯誤率；樣本分組

1 引 言

中小企業(yè)在國家創(chuàng)新經(jīng)濟發(fā)展模式和解決就業(yè)問題中占有越來越重要的位置，量大面廣，給基層稅務(wù)稽查和納稅評估工作帶來了很大的風(fēng)險.因此，建立實用性強和可靠的稅務(wù)稽查評價模型，既能幫助企業(yè)提高涉稅風(fēng)險的防控能力，又能幫助稅務(wù)部門足額征收稅款，日益受到政府有關(guān)部門(稅務(wù)局等)和學(xué)界的重視[1，2].樓文高等[3]對Tobit模型、層次分析法(AHP)、主成分法(PCA)、判別分析(MDA)、Logistic模型和多元線性回歸(MLR)等傳統(tǒng)統(tǒng)計模型以及新興的多層感知器神經(jīng)網(wǎng)絡(luò)(BPNN)、概率神經(jīng)網(wǎng)絡(luò)(PNN)、支持向量機(SVM)、自組織神經(jīng)網(wǎng)絡(luò)(SOM)等數(shù)據(jù)挖掘技術(shù)與傳統(tǒng)統(tǒng)計模型的組合模型[1-2, 4-8]的優(yōu)缺點、適用情況以及現(xiàn)有文獻存在的問題等進行了詳細的評述，并應(yīng)用廣義回歸神經(jīng)網(wǎng)絡(luò)(GRNN)和多重交叉檢驗法，建立了適用于小樣本情況的稅務(wù)稽查GRNN模型，分類錯誤率10%左右，明顯低于傳統(tǒng)統(tǒng)計模型和SVM模型，取得了較好的效果.但是，由于GRNN建模過程中確定合理的光滑因子值是相當(dāng)繁瑣的，而且GRNN模型是隱性模型[3, 9-10]，無法顯性地直接揭示出企業(yè)納稅情況與各個評價指標(biāo)之間的非線性關(guān)系，給后續(xù)的稅務(wù)稽查工作(判定、研究企業(yè)納稅情況)以及企業(yè)如何制定合理的納稅策略、降低涉稅風(fēng)險帶來不便.

另一方面，投影尋蹤分類(Projection Pursuit Clustering，簡稱PPC)技術(shù)是一種適用于高維、非線性、非正態(tài)分布數(shù)據(jù)處理的新興統(tǒng)計建模方法[11-14]，不僅數(shù)學(xué)意義清晰，而且是顯性模型，便于對樣本和評價指標(biāo)的重要性進行排序和分類研究.本文首次將PPC技術(shù)引入到企業(yè)稅務(wù)稽查研究中，212個建模樣本和174個測試樣本(占45%)的平均分類錯誤率低于6.00%，低于MLR、MDA、Logistic等傳統(tǒng)統(tǒng)計模型和SVM模型，建立了更加簡潔、實用、可靠和有效的稅務(wù)稽查模型，應(yīng)優(yōu)先用于中小企業(yè)的稅務(wù)稽查研究和實踐中.

2 建立中小企業(yè)稅務(wù)稽查評價指標(biāo)體系

本課題組選取上海某區(qū)稅務(wù)局管轄的木制家具制造中小企業(yè)增值稅納稅情況進行研究，共有2 000多家企業(yè)，考慮到數(shù)據(jù)的真實可靠性，選取一般納稅人作為研究對象.在2010～2012年度有較完整財務(wù)數(shù)據(jù)和增值稅納稅申報數(shù)據(jù)的正常經(jīng)營企業(yè)有725家(次)，成立時間短則1年多，長則16年多.筆者查閱了40多篇有關(guān)稅務(wù)稽查(納稅評估等)的文獻和文件，先后被采用過的評價指標(biāo)有180多個.在比較研究的基礎(chǔ)上，筆者根據(jù)現(xiàn)有企業(yè)財務(wù)數(shù)據(jù)和增值稅納稅申報表可獲得數(shù)據(jù)兩個方面，建立了由如下15個指標(biāo)組成的企業(yè)增值稅稅務(wù)稽查評價指標(biāo)體系：全部收入稅負(fù)率(x1)、成本費用稅負(fù)率(x2)、銷售毛利率(x3)、成本費用利潤率(x4)、銷售利潤率(x5)、總資產(chǎn)收益率(x6)、銷售費用率(x7)、管理費用率(x8)、財務(wù)費用率(x9)、流動資產(chǎn)周轉(zhuǎn)率(x10)、總資產(chǎn)周轉(zhuǎn)率(x11)、流動比率(x12)、速動比率(x13)、現(xiàn)金比率(x14)和資產(chǎn)負(fù)債率(x15).經(jīng)區(qū)稅務(wù)局多位經(jīng)驗豐富的同志和課題組稅務(wù)學(xué)方面專家共同研究和認(rèn)真、仔細的甄別，在上述725家企業(yè)中，判定211家企業(yè)進行了誠實納稅，175家企業(yè)沒有誠實納稅，對其他339家企業(yè)的納稅情況進行初步判定，希望借助本課題組建立的模型進行輔助稽查判斷.

3 投影尋蹤分類(PPC)建模原理

Friedman等[11]于1974年提出了一種適用于非線性、高維和非正態(tài)分布數(shù)據(jù)處理的新興統(tǒng)計建模方法——投影尋蹤分類(PPC)模型，即把高維數(shù)據(jù)投影到低維(1～3維)子空間上，通過分析低維子空間上的數(shù)據(jù)規(guī)律以達到揭示高維數(shù)據(jù)特征的目的.由于在常用的一維PPC模型建模過程中須確定合理的窗口半徑R值，而如何確定合理、有效的R值，迄今還缺乏理論依據(jù)和指導(dǎo)[11-14]，并且不同的R值，往往得到不同的建模結(jié)果.因此，為避免出現(xiàn)PPC模型結(jié)果的不確定性和非唯一性，消除R值對建模結(jié)果的不利影響，本文采用另一種投影尋蹤目標(biāo)函數(shù)型式.

根據(jù)稅務(wù)專家的先驗知識和經(jīng)驗，已經(jīng)事先判定了部分企業(yè)的納稅情況y(i)(即誠實納稅或非誠實納稅).因此，根據(jù)要使所有樣本點盡可能形成若干個團(類)(誠實納稅類或非誠實納稅類)，類與類之間盡可能分散和類內(nèi)樣本點盡可能密集的要求，可以構(gòu)建如下投影尋蹤目標(biāo)函數(shù)[12]：

(1)

4 建立中小企業(yè)稅務(wù)稽查PPC模型

4．1 輸入數(shù)據(jù)歸一化預(yù)處理和模型輸出期望值的設(shè)定

為了消除各指標(biāo)量綱不同對建模結(jié)果的不利影響，應(yīng)對各評價指標(biāo)的輸入數(shù)據(jù)進行線性歸一化預(yù)處理，轉(zhuǎn)化為 [0, 1]區(qū)間內(nèi)的值.考慮到分類要求的對稱性，筆者設(shè)定誠實納稅和非誠實納稅樣本的模型理論期望輸出值分別為-0.5和0.5.

4.2 篩選出對稅務(wù)稽查具有重要影響的評價指標(biāo)

大量研究表明，如果模型中包含不重要的(或者稱為不相關(guān)的)指標(biāo)，不僅不能提高模型的精度(性能)，反而會降低模型的精度，并無畏增加建模和收集樣本數(shù)據(jù)的成本，可以說是有害無益的[1, 3，4,6 , 15，16].本文采用簡單、實用和有效的靈敏度分析法來篩選重要評價指標(biāo)[6, 15，16]，借助于Statsoft公司出品的商品化軟件STATISTICA Neural Network (以下簡稱SNN)軟件GRNN模型的靈敏度分析功能[6, 9]，刪除了x13、x8、x9、x12和x14共五個指標(biāo)，即從15個指標(biāo)中篩選出了對稅務(wù)稽查具有重要影響的10個評價指標(biāo).

4.3 把樣本數(shù)據(jù)進行合理分組提高模型的泛化能力和魯棒性

建立的模型是否可靠和有效，必須用非建模樣本(或稱測試樣本)的性能(精度和有效性)來衡量，因此，必須把樣本數(shù)據(jù)分成建模樣本和測試樣本.如果樣本數(shù)量足夠多，無論采用什么方法，隨機抽取就能確保各組樣本具有相似(同)的性質(zhì)(即均值和方差基本相等)，但如果樣本數(shù)量不是很多，隨機分組就可能出現(xiàn)偏差，這樣建立的模型很可能就是無效和不可靠的.因此，國外文獻早在1960年代就開始研究如何分組才能確保各組樣本具有相似的性質(zhì)，其中SOM方法分組效果較好，也便于實現(xiàn)[6,15，16]．為此，本課題組選用SNN的SOM方法進行樣本分組.一般地，模型的內(nèi)插性能總是好于外插性能，即要求建模樣本必須包含每個指標(biāo)的最大值和最小值(包括誠實納稅和非誠實納稅兩種情況)．對于本例，包含10個誠實納稅樣本和10個非誠實納稅樣本.據(jù)此，將上述386個樣本分成了具有相似性質(zhì)的建模樣本和測試樣本(174個，約占45%)，數(shù)據(jù)特征值如表1所示.

表1 用SOM方法進行分組的兩類樣本的均值和標(biāo)準(zhǔn)差

*注：“1”和“2”分別表示誠實納稅和非誠實納稅樣本，Tr和Te分別指建模樣本和測試樣本；誠實納稅樣本211個，其中Tr和Te分別為117和94個，非誠實納稅樣本175個，其中Tr和Te分別為95和80個.

4．4 建立中小企業(yè)稅務(wù)稽查PPC模型

(2)

再次調(diào)用筆者編制的最優(yōu)化程序，求得全局最優(yōu)解，目標(biāo)函數(shù)值Q(b)=16.378 6，b1=1.233 1，b2=-6.229 3.

5 結(jié)果與分析

5.1 PPC建模結(jié)果及其樣本屬性的判定(誠實納稅或非誠實納稅)

PPC模型每個樣本的輸出值是一個實數(shù)，因此，必須采用某種規(guī)則來判定每個樣本究竟屬于誠實納稅還是非誠實納稅.對于結(jié)果基本服從(右半支)正態(tài)分布N(μ1,σ1)(誠實納稅樣本)，和(左半支)正態(tài)分布N(μ2,σ2)(非誠實納稅樣本)的兩類樣本來說，區(qū)分他們的分界值μF可用下式確定[12]：

(3)

因此，如果模型輸出值小于分界值μF，就可判定該樣本為誠實納稅，離分界值越遠，誠實納稅的概率就越大，反之為非誠實納稅.根據(jù)上述判斷準(zhǔn)則，可以判定本例212個建模樣本和174個測試樣本屬性，其I類、II類和平均分類錯誤率如表2所示.

表2 不同建模方法得到的兩類樣本的分類錯誤率

注：*I類錯誤率是指把誠實納稅的樣本判定為非誠實納稅的百分率，反之是II類錯誤率；**采用MDA、MLR和Logistic建模時，刪除了共線性變量x1、x3、x7、x10和x15；&為改進型PPC模型的結(jié)果.

從PPC模型建模樣本和測試樣本的I類和II類錯誤率可以看出，本例的樣本分組是合理和有效的.如果分組不合理，就會導(dǎo)致建模樣本和測試樣本的分類錯誤率相差很大，從而導(dǎo)致模型沒有泛化能力和實用價值.

5.2 與MLR、MDA、Logistic和SVM等模型計算結(jié)果的對比

針對上述兩類樣本數(shù)據(jù)，筆者也建立了MDA、MLR、Logistic和SVM模型，其I類、II類和平均分類錯誤率也列于表2中.

由表2可知， PPC模型的分類錯誤率最低，MDA、MLR和SVM模型的分類錯誤率都出現(xiàn)了大于10%的情況，尤其是MDA和SVM模型的II類錯誤率較高，危害更大.PPC模型的分類錯誤率與Logistic模型基本相當(dāng).但從原理上講，MDA、MLR、Logistic和SVM模型只有對服從正態(tài)分布規(guī)律數(shù)據(jù)的建模結(jié)果才具有較好的可靠性，否則可靠性難以保證.顯然，本例的各個評價指標(biāo)數(shù)據(jù)都不服從正態(tài)分布，因此，這些傳統(tǒng)統(tǒng)計模型是否能用于企業(yè)實際納稅情況的判定，還有待進一步研究分析.

犯II類錯誤的危害性(即把非誠實納稅企業(yè)判定為誠實納稅，將導(dǎo)致稅收流失)遠遠大于犯I類錯誤.因此，II類錯誤率低于I類錯誤率的模型較為合理.

5.3 其他339家木制家具制造企業(yè)納稅情況的判定

根據(jù)前述，另有339家企業(yè)的納稅情況還有待判定.339家企業(yè)各個指標(biāo)的均值和方差如表3所示，作為對比，上述212個建模樣本及其誠實納稅樣本與非誠實納稅樣本的均值和方差也列于表3中.顯然，這339家企業(yè)的均值與建模樣本中誠實納稅樣本和非誠實納稅樣本的均值都相差較大，說明這339家企業(yè)中肯定既有誠實納稅的，也有非誠實納稅的.為此，把上述339家企業(yè)10個評價指標(biāo)的歸一化數(shù)據(jù)導(dǎo)入上述建立的PPC模型，得到了PPC模型的輸出值，其中190家企業(yè)的輸出值小于分界值μF(-0.096 2)，這些企業(yè)被判定為誠實納稅，另外149家企業(yè)被判定為非誠實納稅.為了進一步驗證PPC模型的有效性、泛化能力和魯棒性，表3也列出了339家企業(yè)中被判定為誠實納稅和非誠實納稅企業(yè)的各個指標(biāo)的均值和標(biāo)準(zhǔn)差.從表3可以看出，對于指標(biāo)x1、x2、x7和x15來講，339家企業(yè)與212個建模樣本的性質(zhì)是矛盾的(如339家企業(yè)中，誠實納稅樣本x1的均值顯著大于非誠實納稅樣本，而212個建模樣本的情況，則正好相反).

5.4 建立稅務(wù)稽查改進型PPC模型

表3 建模樣本和339家企業(yè)及其被PPC模型判定為誠實納稅、非誠實納稅樣本各個評價指標(biāo)的均值和標(biāo)準(zhǔn)差

表4 建模樣本和339家企業(yè)及其被改進型PPC模型判定為誠實納稅、非誠實納稅樣本指標(biāo)的均值和標(biāo)準(zhǔn)差

5.5 各個評價指標(biāo)的重要性排序和分類

從公式(1)和(2)可知，權(quán)重越大的指標(biāo)對納稅情況的影響就越顯著，在這六個指標(biāo)中，x11最重要，其次是x10，然后依次是x4、x3、x5，x6最不重要，這些指標(biāo)的數(shù)值越大，改進型PPC模型的輸出值越大，即企業(yè)存在非誠實納稅的可能性也越高.因此，在判斷研究企業(yè)納稅情況時，首先要重點考察其總資產(chǎn)周轉(zhuǎn)率數(shù)據(jù)是否過高，其次是考察流動資產(chǎn)周轉(zhuǎn)率數(shù)據(jù)是否過高，然后是依次考察成本費用利潤率、銷售成本率、銷售利潤率及其總資產(chǎn)收益率等指標(biāo)數(shù)據(jù)是否過高.根據(jù)各個指標(biāo)的權(quán)重大小，也可以對他們的重要性進行分類，由于x11和x10的權(quán)重大于0.50，其屬性可以歸為“最重要”一類，x4、x3和x5的權(quán)重在[0.3, 0.4]范圍內(nèi)，其屬性可以歸為“重要”一類，x6的權(quán)重小于0.30，其屬性可以歸為“次重要”一類.

6 結(jié)束語

針對上海市某區(qū)386家木制家具制造企業(yè)納稅申報的15項指標(biāo)數(shù)據(jù)，首先采用靈敏度分析方法篩選出對企業(yè)納稅情況具有重要影響的10個指標(biāo)，再采用SOM方法把樣本分成具有相似性質(zhì)的建模樣本和測試樣本(174個，約占45%)，以確保后續(xù)建模的魯棒性和泛化能力.

采用可用于非線性、高維、非正態(tài)分布數(shù)據(jù)處理的投影尋蹤分類(PPC)技術(shù)，針對212個建模樣本，建立了求得真正全局最優(yōu)解的PPC模型，建模樣本和測試樣本的平均分類錯誤率僅為4.72%和4.60%，均低于MLR、MDA、SVM等模型.同時發(fā)現(xiàn)，對339家實際企業(yè)納稅情況的判定表明，有4個指標(biāo)的性質(zhì)是相互矛盾的，為此建立了刪除該四個指標(biāo)后僅包含六個指標(biāo)的改進型PPC模型，建模樣本和測試樣本的平均分類錯誤率分別為6.13%和4.60%，而且與339家實際企業(yè)誠實納稅和非誠實納稅樣本的評價指標(biāo)性質(zhì)完全一致，說明建立的改進型PPC模型具有很高的精度、泛化能力和魯棒性，而且模型包含的評價指標(biāo)更少，實用性更強，可用于實際企業(yè)納稅情況的判定.

采用MLR、MDA和SVM方法建模，不僅模型的平均識別錯誤率高于PPC模型，而且兩類錯誤率差異也較大，說明模型存在一定的偏態(tài)性，再者，從原理上講，在數(shù)據(jù)不服從正態(tài)分布規(guī)律時，無法保證這些模型的可靠性和魯棒性.

PPC模型采用一維連續(xù)實數(shù)輸出方式，便于研究兩類不同樣本輸出值的分布規(guī)律，從而為分析、判定真實樣本屬性提供便利，可同時完成樣本的分類和排序研究.改進型PPC模型建模過程簡捷，意義清晰，無須人為確定某些參數(shù)的合理值，可最大程度規(guī)避人為因素的影響.

[1]RSWU,CSOU,HYLIN,etal. Using data mining technique to enhance tax evasion detection performance [J]. Expert Systems with Applications, 2012, 39(10): 8769-8777．

[2] 何輝, 侯偉. 我國納稅評估存在的問題與路徑選擇[J]. 稅務(wù)研究, 2013，29(2):75-77．

[3] 樓文高, 婁元英, 尹淑平. 基于廣義回歸神經(jīng)網(wǎng)絡(luò)的稅務(wù)稽查選案實證研究[J]. 廣東商學(xué)院學(xué)報, 2013, 28(6):74-80．

[4] J ALM. Measuring, explaining, and controlling tax evasion: lessons from theory, experiments, and field studies [J]. International Tax Public Finance, 2012, 19(1):54-77．

[5] 葉藝勇. 基于支持向量機和領(lǐng)域知識的納稅評估預(yù)警模型[J]. 數(shù)學(xué)的實踐與認(rèn)識, 2014, 44(1):72-77．

[6] P C GONZALEZ, G D VELASQUEZ. Characterization and detection of taxpayers with false invoices using data mining techniques [J]. Expert Systems with Applications, 2013, 40(5):1427-1436．

[7] C LIN, I LIN, C WU,etal. The application of decision tree and artificial neural network to income tax audit: the examples of profit-seeking enterprise income tax and individual income tax in Taiwan [J]. Journal of the Chinese Institute of Engineers, 2012, 35(4):401-411．

[8] 夏輝, 李仁發(fā). 基于SVM與SOM的稅務(wù)稽查選案模型研究[J]. 科學(xué)技術(shù)與工程, 2009, 9(11):4027-4031．

[9] STATSOFT Inc. Electronic Statistics Textbook [EB/OL]．Tulsa ( http://www.statsoft.com/textbook ) , 2011．

[10]D F SPECHT. A generalized regression neural network [J]. IEEE Transactions on Neural Networks, 1991, 2(6):568-576．

[11]J H FRIEDMAN, J W TUKEY. A projection pursuit algorithm for exploratory data analysis [J]. IEEE Transactions on Computers, 1974, 23(9):881-890.

[12]樓際通, 樓文高, 于秀榮. 商業(yè)銀行個人信用風(fēng)險評價的投影尋蹤建模及其實證研究[J]. 經(jīng)濟數(shù)學(xué), 2013, 30(4):27-33．

[13]H CAUSSINUSS, A RUIZ-GAZEN，G GOVAERT. Exploratory projection Pursuit [C]// G GOVAERT．Data Analysis，London: ISTE Ltd,2009：67-92．

[14]樓文高, 喬龍. 投影尋蹤聚類建模理論的新探索與實證研究[J]. 數(shù)理統(tǒng)計與管理, 2015, 34(1):47-58．

[15]W WU, R J MAY, H R MAIER,etal. A benchmarking approach for comparing data splitting methods for modeling water resources parameters using artificial neural networks [J]. Water Resources Research, 2013, 49(11):7598-7614.

[16]R J MAY, H R MAIER, G C DANDY. Data splitting for artificial neural networks using SOM-based stratified sampling [J]. Neural Networks, 2010, 23(2):283-294．

Tax-Checking Assessment of Small and Medium-Sized Enterprises Applying Projection Pursuit Clustering Technique and Its Positive Research

LOU Wen-gao1，LOU Ji-tong2，SONG Lei-juan1， WANG Lang-qing1

(1.Faculty of Fiscal and Financial, Shanghai Business School, Shanghai 200235,China; 2. The University of North Carolina at Chapel Hill, NC 27514,USA)

Based on the 15 variables’ (indexes’ ) tax-reporting data of 386 wooden-furniture manufacturing small- and medium-sized enterprises (WFMSMEs) located in some districts of Shanghai city, the ten variables mainly influencing the tax-checking situation (tax evasion or compliance) of the 386 WFMSMEs were obtained by applying sensitivity analysis method (SAM) for selecting input variables. The modelling set data and testing set data (about taking up 45%) with similar characteristics - similar mean values and variance-were divided using self-organizing map (SOM) approach. The practical, feasible and effective projection pursuit clustering (PPC) model for tax-checking assessment was thus established. Compared with the multivariate linear regression (MLR), the multivariate discriminant analysis (MDA), Logistic and the support vector machine (SVM), the established PPC model possesses the most accurate and the lowest classification-error percentage (CEP) of the models. The mean CEP of modelling set data and the testing set data is lower than 6%. The improved PPC model including fewer variables is thus suitable to tax-checking assessment and research. The tax-checking situation of the other 339 WFMEs was also assessed and judged, and the results show that the established improved PPC model possesses high generalization and robustness.

tax-checking assessment; projection pursuit clustering (PPC) model; classification-error percentage; samples splitting

2015-07-02

上海高校知識服務(wù)平臺“上海商貿(mào)服務(wù)業(yè)知識服務(wù)中心”建設(shè)子項目“稅收風(fēng)險管理信息系統(tǒng)設(shè)計及開發(fā)”(ZF1226)；上海地方本科院?！笆濉眱?nèi)涵建設(shè)上海商學(xué)院重點學(xué)科專業(yè)建設(shè)“金融學(xué)人才培養(yǎng)模式改革與創(chuàng)新”；2013年國家大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃項目(CXGJ-13-002)資助

樓文高(1964—)， 男,教授，博士

E-mail: wglou@sbs.edu.cn；wlou64@126.com

TV139.1; N945.12

A