步驟一：導(dǎo)入數(shù)據(jù)原始表如下，數(shù)據(jù)請以時間（1998，1999，2000，2001??）為橫軸，樣本名（北京，天津，河北??）為縱軸將中文地名替換為數(shù)字。注意：表中不能有中文字符，否則會出現(xiàn)錯誤。面板數(shù)據(jù)中不能有空值。去除年份的一行，將其余部分復(fù)制到stata的dataeditor中，或保存為csv格式。打開stata，調(diào)用數(shù)據(jù)。方法一：直接復(fù)制到dataeditor中。方法二：使用口令：insheetusing

文件路徑調(diào)用例如：insheetusing

C:\STUDY\paper\taxi.csv其中csv格式可用excel的“另存為”導(dǎo)出步驟二：調(diào)整格式首先請將代表樣本的var1重命名

口令：renamevar1

樣本名

例如：renamevar1province

也可直接在var1處雙擊，在彈出的窗口中修改:接下來將數(shù)據(jù)轉(zhuǎn)化為面板數(shù)據(jù)的格式

口令：reshapelongvar,i(樣本名)

例如：reshapelongvar,i(province)其中var代表的是所有的年份（var2,var3,var4??）轉(zhuǎn)化成功后繼續(xù)重命名，其中_j這里代表原始表中的年份，var代表該變量的名稱口令例如：rename_jyearrenamevartaxi也可直接在需要修改的名稱處雙擊，在彈出的窗口中修改步驟三：排序

口令：sort

變量名

例如：sortprovinceyear意思為將province按升序排列，然后再根據(jù)排好的province數(shù)列排year這一列最后，保存。至此，一個變量的前期數(shù)據(jù)處理就完成了，請如法炮制的處理所有的變量，也就是說每個變量都做一個dta文件。在處理新變量前請使用口令：clear將stata重置步驟四：合并數(shù)據(jù)任意打開一個處理過的變量的dta文件作為基礎(chǔ)表（推薦使用因變量的dta文件，這里使用so2作為因變量）

口令：

merge

樣本名時間

using

文件路徑例如：mergeprovinceyearusingC:\STUDY\paper\taxi.dta

意思是將taxi的數(shù)據(jù)添加到so2的數(shù)據(jù)表中然后使用口令：tab_merge然后使用口令：drop_merge

將數(shù)據(jù)表中的_merge一列去掉，接著重新使用口令：sort

樣本名時間例如：sortprovinceyear為新生成的表排序。

如法炮制，將所有的變量都添加到基礎(chǔ)表中，最終步驟：回歸首先，使用口令：xtset

樣本名時間定義面板數(shù)據(jù)例如：xtsetprovinceyear然后使用：口令：xtreg

因變量自變量進(jìn)行回歸分析例如：xtregso2taxibusloaddriversroadlength至此，使用stata進(jìn)行面板數(shù)據(jù)回歸分析完成面板模型分為混合回歸模型、固定效應(yīng)模型、隨機效應(yīng)模型固定效應(yīng)分為個體/時點固定效應(yīng)，個體時點雙固定效應(yīng)隨機效應(yīng)分為個體/時點隨機效應(yīng)，個體時點雙隨機效應(yīng)描述性統(tǒng)計：sum標(biāo)準(zhǔn)化：sum(x-均值）/標(biāo)準(zhǔn)差產(chǎn)生新變量：genpol=(pol-均值）/標(biāo)準(zhǔn)差普通回歸命令：regyx1x2一般p<0.05檢驗多重共線性：estatvifvif為方差膨脹因子，vif<10,否則要消除多重共線性相關(guān)系數(shù)矩陣corryx1x2區(qū)分固定效應(yīng)還是隨機效應(yīng)：xtregyx1x2,feeststorefe這一步結(jié)束看結(jié)果最后一行F檢驗p<0.05,排除混合回歸xtregyx1x2,reeststorerehausmanfere,constantsigmamorehausman檢驗P>0.05接受原假設(shè)