本人準大四生，之前為了完成發展經濟學的論文接觸了一下空間計量的內容，使用的軟件是stata，折騰了好久算是有點心得，想和和我一樣的純新手分享一下。其實人大論壇這邊還是有一些關于利用stata做空間計量的內容，本人也參考了不少，以下有參考的都盡可能列出來。
首先講一講如何利用stata制作分布地圖，本人的課程論文是關于老齡化的區域分布，所以上一個2012年的老齡人口分布圖（省域）：
方法參考自
首先要有地圖的數據，在那個帖子中介紹了一個命令： ；
大致的意思即下載了中國的地圖數據，類似的好像還有一個可以找到世界地圖的數據，不過記不清了，輸入這個命令的過程有點類似于安裝命令，安裝完之后會在stata的安裝文件夾中出現一個.zip的壓縮包，解壓后會出現兩組數據，.dta和.dta，都可以打開看一下，前者有省份的標識，后者并沒有，然后當時本人使用的是第一個文件，我們可以看一下里面的數據：
name一欄可能是中文無法顯示的原因，不過我們可以根據英文名來確定不同的省份。然后這里的原理本人也并不是非常清楚，只介紹本人是怎么做的。我們可以發現一個省份可能對應若干行，所以在這里做的事就是將我們的數據對應填進這個stata數據里，這會花掉一點時間，不過在excel里整理起來還是比較快的。在這里如果我們使用的是面板數據的話，可以把我們的變量一行一行整理，而不是一個個，像我就一下子輸入了幾年的數據，到時候作圖換個變量就可以實現。
給大家截一個excel整理的圖：（請忽視變量）
數據到手，我們還需要一個命令，相信很多人都知道的一個spmap的命令，下載安裝：ssc spmap；
在這里用2012年的數據做一個范例（變量名我加了個a，直接數字似乎不可行）：
輸入命令：spmap a2012 using ".dta",id(id)
然后我們就可以看到彈出來的圖了：
我們在這里可以進行簡單的編輯，然后在這里還有一個值得關注的部分，那就是顏色劃分的范圍，我們可以看到左下角的劃分是自動生成的，如果我們想要自定義怎么辦呢？只要我們在之前的命令上加上后綴即可： spmap ? a2012 using ".dta",id(id) clm(c) clb(0 7 8 9 10 11 16)
參考自：
相信大家能看出這個后綴是什么意思，以及我們應該怎么去定義范圍：
最后再提供一個比較好的帖子，也是關于作圖的，當然要比我的解說專業得多：
；
感覺上畫這個圖非常ZB，做出的論文老師不置可否，然后說：你這圖畫的不錯，我二十年前就開始畫了/(ㄒoㄒ)/~~
畫完了圖，我們接下來就開始正式進入空間計量部分，一般情況下，我們都會先接觸一個Moran's I指數，或者叫莫蘭指數，這是一個衡量空間自相關的變量，如果沒有空間上的相關性，那么我們做空間計量也就沒有意義了。而為了計算莫蘭指數，乃至為了之后的空間計量模型部分，我們都需要一個空間權重矩陣W；
關于空間權重矩陣部分，主要有兩種計算的方式：一、以相鄰為1，不相鄰為0；二、以距離小于固定值的為1，距離大于固定值的為0；由于本人使用的數據是省域，于是恬不知恥地直接找了一個論壇里分享的矩陣完事，計算方法是第一種，附件中會分享；而如果是縣域或者其他空間分布的矩陣的話，本人并沒有嘗試過，但是有一個和我一起嘗試做空間計量的同學倒是折騰過，折騰了很久，最后還是以失敗告終，所以在這里就不誤人子弟了，僅提供一個相對詳細的帖子作參考：；
于是我們有了矩陣（說得好像真的有了一樣(#‵′)凸），關于矩陣的形式stata定義面板數據不行，以我的數據為例，就是31×31的一個矩陣，我們需要將其變為dta數據，其實貼進stata，然后另存為就可以，例如我們取名為W.dta。
然后我們就需要空間計量的命令了，在這里推薦一組spat開頭的系列命令：
（截圖來自陳強的高級計量經濟學及stata應用第二版，這本書有一章就在講空間計量，后期摸索我主要就靠這本書，由于論壇好像已經有這本書的分享，而且是收費的，我也不好全部免費發上來，僅免費提供空間計量一章）
在調用命令之前我們需要將W.dta放到我們的安裝文件夾里，然后輸入命令： using W.dta,name(W);
正常情況下，我們的空間權重矩陣就存在了，不過這里可能會出現問題，而出現問題的還是我那個悲催的同學，我的電腦是，然后她是mac，結果她做到這一步時發現無法調用W.dta，也就是說無法生成矩陣，雖然的確把W.dta放到安裝文件里了。不過后來我們嘗試了一下，發現一種解決辦法，這一點等一會兒再講。
對于像我一樣的新手有一點必須要說明的就是權重矩陣和我們的數據應該是兩個文件，而在我們導入矩陣之后，并不會占用我們stata儲存數據的那一塊區域，所以我們要跑回歸的數據甚至是可以直接復制粘貼進去的，其實本人就是這么干的，我們可以看一下矩陣和數據的形式：
（注意不要把第一列復制進去）
本人的數據是多年份的，所以是這樣的，在這里，我們并不需要一個ID列，stata會自動默認，即我的矩陣有31列，這里的數據也應當有31列，于是就會自動對應；
之前我們已經引入了權重矩陣，因而接下來直接輸moran指數的命令：
a2012,(W) moran//全局莫蘭指數，a2012是我的變量，即2012年的各省老齡化水平，別忘記變量右面有個逗號，不然會出錯，此命令默認單尾；
a2012,(W) moran//局部莫蘭指數
a2013 ,(W) moran graph (moran) (n)//莫蘭散點圖
關于莫蘭散點圖有一點需要說明的就是當時做出的散點圖是黑底的，直接貼到word上非常難看，但是在stata之中卻提醒無法編輯，搞了很久沒找到解決辦法，后來機智的我把黑底圖貼進word之后，利用word自帶的對圖片的編輯功能設置了黑色為透明色，然后略加調整，散點圖就已經有點像模像樣了。不過只能說奇技淫巧，遇到這個問題又實在解決不了的可以考慮下。
關于莫蘭指數的部分就到此為止了，問題并不大，除了我那個悲催的同學。我那個同學出現的問題是無法調用數據生成權重矩陣，而我們后來的解決辦法是把權重矩陣的數據粘貼進stata之后stata定義面板數據不行，直接生成矩陣： ,name(W)；剛才試了一下，這個辦法在我的電腦上并不行，不知道出了什么問題，但是我們的確是按照直接生成矩陣的思路來的，生成矩陣之后再把矩陣數據清空，粘貼進我們要使用的數據即可；這一塊可能有待大神補充；
整體上這塊內容在我搞到陳強的那本書之后很快就出來了，主要是那本書上有個例子，在附件中我也會將例子的數據附上，大家可以按照那里的思路學習；
接下來的部分是空間模型部分，由于我做的是面板的數據，所以這一部分并沒有完全接觸過，當然兩者之間在操作上沒有本質的差異，這一塊內容大家可以參考陳強的教材，摸索一下，很快就能掌握；
好了，最后就是我的空間面板部分了，事實上，這部分內容在陳強的教材占也有講到，但是但是我照著他的教材操作，仍然出現了一點問題（事實證明，那是我的問題，而不是教材的/(ㄒoㄒ)/~~，所以作為經驗貼我覺得有必要講一下，哪些地方可能會出現問題）
首先下一個非官方命令xsmle:net xsmle;
然后第一步仍然是導入矩陣，這一點可以參考之前的操作，陳強的教材中又提供了一條新命令，但就不贅述了，本質是一樣的： using W.dta,name(W);
然后大家可以照著陳強的內容開始操作，我以我自己的命令舉例：
using .dta,name(W)//打開文件后導入矩陣
xtset id 年份//設定面板
xsmle 老年化率 醫療機構數 人均GDP 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sdm) nolog
（本人是，所以支持中文變量名，大家看到中文變量名不要感到奇怪）
然后你就會喜聞樂見地發現出錯了，stata會提醒： not ；
哪里出了問題呢？其實是我們在我們進行這一輪操作的時候，我們導入的矩陣必須是標準化的！不要問我為什么，我也只是查來的/(ㄒoㄒ)/~~
using .dta,name(W) //打開文件后導入矩陣，標準化以防出現 not 錯誤
xtset id 年份//設定面板
xsmle 老年化率 醫療機構數 人均GDP 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sdm) nolog
在標準化之后，就不會再出現錯誤提醒了；
具體操作大家都可以參考陳強的教材；
然后大家應該知道空間計量有多個可選模型，SDM，SAR和SEM，對于如何選擇模型我不想多說，因為我也完全不懂啊，但是在三個模型的命令上其實存在一個非常微小的差別，我們先來看一看正確的命令：
xsmle 老年化率 醫療機構數 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sdm) nolog
xsmle 老年化率 醫療機構數 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sar) nolog
xsmle 老年化率 醫療機構數 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,emat(W)model(sem) nolog
不知道大家看出差別沒有，關鍵在于SEM模型在調用權重矩陣的命令是：emat，而其他兩個模型為wmat，關于這一點，我們可以help xsmle查看幫助
如果還是和我一樣直接將括號里的sdm,sar,sem互換一下的話，就會出現： () with SEM ；
最后再提醒一個關于霍斯曼檢驗的問題，其實和空間計量問題不大，不過我當時遇到了就提醒一下；
我們之前在跑回歸的時候都是默認隨機模型，而且都是穩健結果（），但是在霍斯曼檢驗時，就需要將（）去掉：
xsmle 老年化率 醫療機構數 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sdm) ( 普通高中學校數 ) nolog
est sto re
xsmle 老年化率 醫療機構數 普通高中學校數 在職職工參加養老保險人數萬人 人口自然增長率,wmat(W)model(sdm) ( 普通高中學校數 ) nolog fe
est sto fe
fe re
如果穩健（）的話，就會出現： be used with vce(), vce( cvar), or p- data；
說實話，本人是為了完成那一篇課程論文才接觸的空間計量，一開始完全處于一臉懵逼的狀態，摸索了許久才算有點成果，不過也就僅限于這個帖子里所說的啦，所以要是提問什么高深的問題我也是完全不會的啊。而這個帖子呢，也不算什么技術貼，只能說是一個經驗貼，相信還是能夠為那些和我一樣剛接觸空間計量，啥都不懂的小白一點幫助。
最后說明一下附件的情況。
“我的數據”里包括三部分，第一部分是一個31×31的空間權重矩陣（以相鄰為1，不相鄰為0，然后海南省不與任何一個省份相鄰）；第二部分是1997-2013年的人口老齡化數據（65歲以上老齡人口比重）；第三部分是自變量和因變量；
“陳強 空間計量部分”是我從陳強的《高級計量經濟學及stata應用》的第二版中截取的空間計量的部分內容；
然后附上的數據是陳強教材中使用到的四個數據，有兩個不讓我傳，于是打成壓縮包了
壓縮包.zip(25.59 KB)
我的數據.xlsx(24.58 KB)
陳強 空間計量部分.pdf(10.03 MB)