以下是我在近三年做各類計(jì)量和統(tǒng)計(jì)分析過程中感受最深的東西，或能對(duì)大家有所幫助。當(dāng)然，它不是ABC的教程，也不是細(xì)致的數(shù)據(jù)分析方法介紹，它只 是“總結(jié)”和“體會(huì)”。由于我所學(xué)所做均甚雜，我也不是學(xué)統(tǒng)計(jì)、數(shù)學(xué)出身的，故本文沒有主線，只有碎片，且文中內(nèi)容僅為個(gè)人觀點(diǎn)，許多論斷沒有數(shù)學(xué)證明， 望統(tǒng)計(jì)、計(jì)量大牛輕拍。
關(guān)于軟件
于我個(gè)人而言，所用的數(shù)據(jù)分析軟件包括EXCEL、SPSS、STATA、EVIEWS。在分析前期可以使用EXCEL進(jìn)行數(shù)據(jù)清洗、數(shù)據(jù)結(jié)構(gòu)調(diào) 整、復(fù)雜的新變量計(jì)算（包括邏輯計(jì)算）；在后期呈現(xiàn)美觀的圖表時(shí)，它的制圖制表功能更是無可取代的利器；但需要說明的是，EXCEL畢竟只是辦公軟件，它 的作用大多局限在對(duì)數(shù)據(jù)本身進(jìn)行的操作，而非復(fù)雜的統(tǒng)計(jì)和計(jì)量分析，而且，當(dāng)樣本量達(dá)到“萬”以上級(jí)別時(shí)，EXCEL的運(yùn)行速度有時(shí)會(huì)讓人抓狂。
SPSS是擅長(zhǎng)于處理截面數(shù)據(jù)的傻瓜統(tǒng)計(jì)軟件。首先，它是專業(yè)的統(tǒng)計(jì)軟件，對(duì)“萬”甚至“十萬”樣本量級(jí)別的數(shù)據(jù)集都能應(yīng)付自如；其次，它是統(tǒng)計(jì)軟 件而非專業(yè)的計(jì)量軟件，因此它的強(qiáng)項(xiàng)在于數(shù)據(jù)清洗、描述統(tǒng)計(jì)、假設(shè)檢驗(yàn)（T、F、卡方、方差齊性、正態(tài)性、信效度等檢驗(yàn)）、多元統(tǒng)計(jì)分析（因子、聚類、判 別、偏相關(guān)等）和一些常用的計(jì)量分析（初、中級(jí)計(jì)量教科書里提到的計(jì)量分析基本都能實(shí)現(xiàn)），對(duì)于復(fù)雜的、前沿的計(jì)量分析無能為力；第三，SPSS主要用于 分析截面數(shù)據(jù)，在時(shí)序和面板數(shù)據(jù)處理方面功能了了；最后，SPSS兼容菜單化和編程化操作，是名副其實(shí)的傻瓜軟件。
STATA與EVIEWS都是我偏好的計(jì)量軟件。前者完全編程化操作，后者兼容菜單化和編程化操作；雖然兩款軟件都能做簡(jiǎn)單的描述統(tǒng)計(jì)，但是較之 SPSS差了許多；STATA與EVIEWS都是計(jì)量軟件，高級(jí)的計(jì)量分析能夠在這兩個(gè)軟件里得到實(shí)現(xiàn)；STATA的擴(kuò)展性較好，我們可以上網(wǎng)找自己需要 的命令文件（.ado文件），不斷擴(kuò)展其應(yīng)用，但EVIEWS就只能等著軟件升級(jí)了；另外，對(duì)于時(shí)序數(shù)據(jù)的處理，EVIEWS較強(qiáng)。
綜上，各款軟件有自己的強(qiáng)項(xiàng)和弱項(xiàng)，用什么軟件取決于數(shù)據(jù)本身的屬性及分析方法。EXCEL適用于處理小樣本數(shù)據(jù)，SPSS、 STATA、EVIEWS可以處理較大的樣本；EXCEL、SPSS適合做數(shù)據(jù)清洗、新變量計(jì)算等分析前準(zhǔn)備性工作，而STATA、EVIEWS在這方面 較差；制圖制表用EXCEL；對(duì)截面數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析用SPSS，簡(jiǎn)單的計(jì)量分析SPSS、STATA、EVIEWS可以實(shí)現(xiàn)，高級(jí)的計(jì)量分析用 STATA、EVIEWS，時(shí)序分析用EVIEWS。
關(guān)于因果性
做統(tǒng)計(jì)或計(jì)量，我認(rèn)為最難也最頭疼的就是進(jìn)行因果性判斷。假如你有A、B兩個(gè)變量的數(shù)據(jù)，你怎么知道哪個(gè)變量是因（自變量），哪個(gè)變量是果（因變量）？
早期，人們通過觀察原因和結(jié)果之間的表面聯(lián)系進(jìn)行因果推論，比如恒常會(huì)合、時(shí)間順序。但是，人們漸漸認(rèn)識(shí)到多次的共同出現(xiàn)和共同缺失可能是因果關(guān) 系，也可能是由共同的原因或其他因素造成的。從歸納法的角度來說，如果在有A的情形下出現(xiàn)B，沒有A的情形下就沒有B，那么A很可能是B的原因，但也可能 是其他未能預(yù)料到的因素在起作用，所以，在進(jìn)行因果判斷時(shí)應(yīng)對(duì)大量的事例進(jìn)行比較，以便提高判斷的可靠性。
有兩種解決因果問題的方案：統(tǒng)計(jì)的解決方案和科學(xué)的解決方案。統(tǒng)計(jì)的解決方案主要指運(yùn)用統(tǒng)計(jì)和計(jì)量回歸的方法對(duì)微觀數(shù)據(jù)進(jìn)行分析，比較受干預(yù)樣本與 未接受干預(yù)樣本在效果指標(biāo)（因變量）上的差異。需要強(qiáng)調(diào)的是，利用截面數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，不論是進(jìn)行均值比較、頻數(shù)分析，還是方差分析、相關(guān)分析，其結(jié)果 只是干預(yù)與影響效果之間因果關(guān)系成立的必要條件而非充分條件。類似的，利用截面數(shù)據(jù)進(jìn)行計(jì)量回歸，所能得到的最多也只是變量間的數(shù)量關(guān)系；計(jì)量模型中哪個(gè) 變量為因變量哪個(gè)變量為自變量，完全出于分析者根據(jù)其他考慮進(jìn)行的預(yù)設(shè)，與計(jì)量分析結(jié)果沒有關(guān)系?？傊?，回歸并不意味著因果關(guān)系的成立，因果關(guān)系的判定或 推斷必須依據(jù)經(jīng)過實(shí)踐檢驗(yàn)的相關(guān)理論。雖然利用截面數(shù)據(jù)進(jìn)行因果判斷顯得勉強(qiáng)，但如果研究者掌握了時(shí)間序列數(shù)據(jù)，因果判斷仍有可為，其中最經(jīng)典的方法就是 進(jìn)行“格蘭杰因果關(guān)系檢驗(yàn)”。但格蘭杰因果關(guān)系檢驗(yàn)的結(jié)論也只是統(tǒng)計(jì)意義上的因果性，而不一定是真正的因果關(guān)系，況且格蘭杰因果關(guān)系檢驗(yàn)對(duì)數(shù)據(jù)的要求較高 （多期時(shí)序數(shù)據(jù)），因此該方法對(duì)截面數(shù)據(jù)無能為力。綜上所述，統(tǒng)計(jì)、計(jì)量分析的結(jié)果可以作為真正的因果關(guān)系的一種支持，但不能作為肯定或否定因果關(guān)系的最 終根據(jù)。
科學(xué)的解決方案主要指實(shí)驗(yàn)法，包括隨機(jī)分組實(shí)驗(yàn)和準(zhǔn)實(shí)驗(yàn)。以實(shí)驗(yàn)的方法對(duì)干預(yù)的效果進(jìn)行評(píng)估，可以對(duì)除干預(yù)外的其他影響因素加以控制，從而將干預(yù)實(shí)施后的效果歸因?yàn)楦深A(yù)本身，這就解決了因果性的確認(rèn)問題。
關(guān)于實(shí)驗(yàn)
在隨機(jī)實(shí)驗(yàn)中，樣本被隨機(jī)分成兩組，一組經(jīng)歷處理?xiàng)l件（進(jìn)入干預(yù)組），另一組接受控制條件（進(jìn)入對(duì)照組），然后比較兩組樣本的效果指標(biāo)均值是否有差 異。隨機(jī)分組使得兩組樣本“同質(zhì)”，即“分組”、“干預(yù)”與樣本的所有自身屬性相互獨(dú)立，從而可以通過干預(yù)結(jié)束時(shí)兩個(gè)群體在效果指標(biāo)上的差異來考察實(shí)驗(yàn)處 理的凈效應(yīng)。隨機(jī)實(shí)驗(yàn)設(shè)計(jì)方法能夠在最大程度上保證干預(yù)組與對(duì)照組的相似性，得出的研究結(jié)論更具可靠性，更具說服力。但是這種方法也是備受爭(zhēng)議的，一是因 為它實(shí)施難度較大、成本較高；二是因?yàn)樵诟深A(yù)的影響評(píng)估中，接受干預(yù)與否通常并不是隨機(jī)發(fā)生的；第三，在社會(huì)科學(xué)研究領(lǐng)域，完全隨機(jī)分配實(shí)驗(yàn)對(duì)象的做法會(huì) 涉及到研究倫理和道德問題。鑒于上述原因，利用非隨機(jī)數(shù)據(jù)進(jìn)行的準(zhǔn)實(shí)驗(yàn)設(shè)計(jì)是一個(gè)可供選擇的替代方法。準(zhǔn)實(shí)驗(yàn)與隨機(jī)實(shí)驗(yàn)區(qū)分的標(biāo)準(zhǔn)是前者沒有隨機(jī)分配樣 本。
通過準(zhǔn)實(shí)驗(yàn)對(duì)干預(yù)的影響效果進(jìn)行評(píng)估，由于樣本接受干預(yù)與否并不是隨機(jī)發(fā)生的，而是人為選擇的，因此對(duì)于非隨機(jī)數(shù)據(jù)，不能簡(jiǎn)單的認(rèn)為效果指標(biāo)的差異 來源于干預(yù)。在剔除干預(yù)因素后，干預(yù)組和對(duì)照組的本身還可能存在著一些影響效果指標(biāo)的因素，這些因素對(duì)效果指標(biāo)的作用有可能同干預(yù)對(duì)效果指標(biāo)的作用相混 淆。為了解決這個(gè)問題，可以運(yùn)用統(tǒng)計(jì)或計(jì)量的方法對(duì)除干預(yù)因素外的其他可能的影響因素進(jìn)行控制，或運(yùn)用匹配的方法調(diào)整樣本屬性的不平衡性——在對(duì)照組中尋 找一個(gè)除了干預(yù)因素不同之外，其他因素與干預(yù)組樣本相同的對(duì)照樣本與之配對(duì)——這可以保證這些影響因素和分組安排獨(dú)立。
隨機(jī)實(shí)驗(yàn)需要至少兩期的面板數(shù)據(jù)，并且要求樣本在干預(yù)組和對(duì)照組隨機(jī)分布，分析方法就是DID（倍差法，或曰雙重差分法）；準(zhǔn)實(shí)驗(yàn)分析用截面數(shù)據(jù)就 能做，不要求樣本在干預(yù)組和對(duì)照組隨機(jī)分布，分析方法包括DID（需兩期的面板數(shù)據(jù)）、PSM（傾向性得分匹配法，需一期的截面數(shù)據(jù)）和PSM- DID（需兩期的面板數(shù)據(jù)）。從準(zhǔn)確度角度來說，隨機(jī)實(shí)驗(yàn)的準(zhǔn)確度高于準(zhǔn)實(shí)驗(yàn)和非實(shí)驗(yàn)分析。
關(guān)于分析工具的選擇
如果根據(jù)理論或邏輯已經(jīng)預(yù)設(shè)了變量間的因果關(guān)系，那么就無需使用實(shí)驗(yàn)方法。我對(duì)非實(shí)驗(yàn)數(shù)據(jù)分析工具的選擇原則如下。
隨著統(tǒng)計(jì)和計(jì)量經(jīng)濟(jì)學(xué)的發(fā)展，各種前沿分析工具層出不窮，但我認(rèn)為最靠譜的分析工具不外乎以下四種：DID（針對(duì)隨機(jī)實(shí)驗(yàn)），多元線性回歸，固定效 應(yīng)變截距模型（FE，針對(duì)面板數(shù)據(jù)），Logit模型或Probit模型（針對(duì)分類因變量數(shù)據(jù)）。其他方法或適用條件苛刻，或分析過程折騰，或方法本身不 可靠（尤其是聚類分析、判別分析，超級(jí)不靠譜），因此能用以上四種方法分析問題時(shí)，不必為“炫方法”而瞎折騰。
關(guān)于擬合優(yōu)度、變量選擇原則及估計(jì)值絕對(duì)大小的意義
在人人的“數(shù)據(jù)分析”小站中，某同學(xué)提出這樣一個(gè)問題：“多元回歸分析中，怎么選擇自變量和因變量，可以使R方達(dá)到80%以上？”
很顯然，問這個(gè)問題的同學(xué)要么沒學(xué)好計(jì)量，要么就是犯了功利主義的錯(cuò)誤，或者二者皆有。擬合優(yōu)度的大小很大程度上取決于數(shù)據(jù)本身的性質(zhì)。如果數(shù)據(jù)是 時(shí)序數(shù)據(jù)，只要拿有點(diǎn)相關(guān)關(guān)系的變量進(jìn)行回歸就能使擬合優(yōu)度達(dá)到80%以上，但這樣的高R方根本說明不了什么，很可能使分析者陷入偽回歸的陷阱，嚴(yán)謹(jǐn)?shù)淖?法當(dāng)然是做平穩(wěn)性檢驗(yàn)和協(xié)整檢驗(yàn)；如果是截面數(shù)據(jù)，根本沒必要追求R方到80%的程度，一般來說，有個(gè)20%、30%就非常大了。
如果一定要增大R方，那么最應(yīng)該做的的確是對(duì)納入模型的變量進(jìn)行選擇。選擇納入模型的原則我認(rèn)為有三條。第一，從理論和邏輯出發(fā)，將可能影響因變量 的變量作為自變量納入模型，即理論上或邏輯上能影響因變量的自變量必須納入模型，即使該自變量的回歸系數(shù)不顯著。第二，奧姆剃刀原則——如無必要，勿增實(shí) 體，即理論上或邏輯上不能影響因變量的自變量不能納入模型，即使該自變量的回歸系數(shù)顯著。第三，防止納入具有多重共線性的自變量。
前面說了，對(duì)截面數(shù)據(jù)進(jìn)行計(jì)量分析，R方能達(dá)到20%、30%是非常了不起的事情。但是，如果擬合優(yōu)度（或類似擬合優(yōu)度的指標(biāo)）在20%、30%或 更低時(shí)，回歸系數(shù)只具有定性或定序上的意義，強(qiáng)調(diào)其絕對(duì)數(shù)值的大小沒什么意義。譬如lnY=alnA+blnB+…+zlnZ+c回歸的R方為20%，a 為0.375，b為0.224，且二者的T檢驗(yàn)顯著，那么我們可以說，A、B對(duì)Y有影響，也可以說一百分點(diǎn)的A變化對(duì)Y的影響大于一百分點(diǎn)的B變化對(duì)Y的 影響（控制其他因素的情況下），但說一百分點(diǎn)的A變化對(duì)Y的影響較一百分點(diǎn)的B變化對(duì)Y的影響大0.151%，就沒什么意義了。
其他一些建議或忠告
用心思考變量間的因果關(guān)系：是A影響了B還是B影響了A？A、B之間是否真的有因果關(guān)系？是否存在C，使C既影響A又影響B(tài)，而A、B本身無直接關(guān)系？
仔細(xì)選擇自變量，不要遺漏重要變量，否則會(huì)造成內(nèi)生性問題。如果遇上了內(nèi)生性問題，先不要忙著尋找工具變量或使用2SLS，尋找被遺漏的變量才是最 重要的事情。如果被遺漏的變量即使找到卻囿于各種困難無法納入分析，而你又忽然想到了一個(gè)絕佳的工具變量，那么恭喜你，你可以在核心期刊發(fā)文章了！
一定要控制其他可能對(duì)因變量產(chǎn)生影響的因素，并認(rèn)識(shí)到對(duì)回歸系數(shù)和偏相關(guān)分析結(jié)果的解釋都是建立在“其他條件不變”的情況之下。
看到R方很大時(shí)不要忙著高興，如果F檢驗(yàn)顯著而T檢驗(yàn)不顯著，很可能存在多重共線性?？吹絫值很大時(shí)，也不要忙著高興，因?yàn)檫@很可能是偽回歸的產(chǎn)物；如果此時(shí)DW值很?。ㄐ∮?.5），那么偽回歸的可能性進(jìn)一步變大。
均值比較雖然簡(jiǎn)單卻考驗(yàn)分析者的嚴(yán)謹(jǐn)性。兩個(gè)看似不同的平均數(shù)、中位數(shù)或比率是否意味著高下有別？樣本取自獨(dú)立總體還是相關(guān)總體？方差“齊”或“不齊”？比較的是平均數(shù)、中位數(shù)還是比率差異？
樣本量限制了所能做的分析，小樣本時(shí)請(qǐng)珍惜自由度；不要用小于30個(gè)樣本的數(shù)據(jù)進(jìn)行計(jì)量分析（尤其是時(shí)序分析）和復(fù)雜的統(tǒng)計(jì)分析；不要以為能從小于或等于5期的數(shù)據(jù)中看出什么“發(fā)展趨勢(shì)”；不要沒有依據(jù)的使用復(fù)雜的模型和分析方法；不要將一目了然的簡(jiǎn)單問題故意復(fù)雜化。
最重要的，不要造假！不對(duì)數(shù)據(jù)本身造假，也不對(duì)分析結(jié)果造假！數(shù)據(jù)分析前可以進(jìn)行一定的清洗，將奇異值去掉，也可以嘗試對(duì)未預(yù)料到的分析結(jié)果進(jìn)行探討和解釋，但如果去改數(shù)據(jù)改分析結(jié)果，那還有什么必要進(jìn)行數(shù)據(jù)分析呢？直接編文章編報(bào)告不就得了？某些“詭異的”、不合常理的數(shù)據(jù)分析結(jié)果，很可能就是研究最重要的所得。
以上，如有錯(cuò)誤，敬請(qǐng)指正；如有補(bǔ)充，歡迎留言，我會(huì)加進(jìn)文中。
后記：過完年就要去工作了，每想及此就頗為傷感。在北大兩年多，除了讓自己更加理想主義外，除了愛上燕園的學(xué)術(shù)氛圍和結(jié)識(shí)到一些好朋友好師長(zhǎng)外，我 學(xué)到并且目前還記得的知識(shí)并不多，且這些知識(shí)大多停留在“術(shù)”的層面。當(dāng)然，“術(shù)”之道亦博大精深，我所掌握的不足萬一。之所以還敢寫下上面的文字貽笑大 方，是想為需要的人提供幫助，也是以此形式悼念自己逝去的大學(xué)時(shí)光。感謝丁延慶老師，感謝邵宜航老師，感謝所有給我以指導(dǎo)和幫助的師長(zhǎng)與朋友。
聯(lián)系客服
