550 likes | 832 Views
Part1. 繪製地圖 安裝繪製地圖所需程式. ssc install spmap ssc install shp2dta ssc install mif2dta. 讀取 CHN_adm1 地圖檔。創出 chinaprovince 以及 coord 兩個檔案. shp2dta using CHN_adm1, database (chinaprovince) coordinates(coord) genid(id) gencentroids(c). *打開 chinaprovince 檔案,刪除 Paracel Islands.
E N D
Part1.繪製地圖安裝繪製地圖所需程式 ssc install spmap ssc install shp2dta ssc install mif2dta
讀取CHN_adm1地圖檔。創出chinaprovince以及coord兩個檔案讀取CHN_adm1地圖檔。創出chinaprovince以及coord兩個檔案 shp2dta using CHN_adm1, database (chinaprovince) coordinates(coord) genid(id) gencentroids(c)
繪出2007年中國FDI分布圖 spmap FDI_2007 using coord, id(id) clnumber(5)
Part2.建立空間計量矩陣 findit spmat
生成距離矩陣 spmat idistance idistance_province longitude latitude, id(id) dfunction(euclidean) normalize(row)
生成stata可讀之檔案idistance_province.spmat spmat save idistance_province using idistance_province.spmat (此時stata資料夾出現idistance_province.spmat)
將剛剛的idistance_province.spmat檔案轉為我們可閱讀的txt檔將剛剛的idistance_province.spmat檔案轉為我們可閱讀的txt檔 spmat export idistance_province using idistance_province.txt
生成相鄰矩陣 spmat contiguity contiguity_province using coord, id(id) normalize(row)
生成stata可讀之檔案 spmat save contiguity_province using contiguity_province.spmat
將剛剛的contiguity_province.spmat檔案轉為我們可閱讀的txt檔將剛剛的contiguity_province.spmat檔案轉為我們可閱讀的txt檔 spmat export contiguity_province using contiguity_province.txt
Part.3 跑空間計量 findit xsmle
spmat use W using "idistance_province.spmat " (使用距離矩陣跑空間計量) spmat use W using "contiguity_province.spmat" (使用相鄰矩陣跑空間計量)
告知stata要跑panel data,才可以跑空間計量 xtset id year
開啟資料data,公布年度_上課用 xsmle FDI FT GDP WAGE OPEN CFDI , wmat(W) mode(sdm) re (其中第一個變數stata會自動設定為Y,之後都是X,re為隨機效果模型,fe為固定效果模型。mode可選sar與sdm)
Sem模型 xsmle FDI FT GDP WAGE OPEN CFDI , ematrix(W) mode(sem) re
將W改為相鄰矩陣 spmat use W using "contiguity_province.spmat" ,replace
重複剛剛的步驟 xsmle FDI FT GDP WAGE OPEN CFDI , wmat(W) mode(sdm) re
檢定使用sdm模型或是sar模型如果拒絕H0,則使用sdm檢定使用sdm模型或是sar模型如果拒絕H0,則使用sdm test [Wx]FT=[Wx]GDP=[Wx]WAGE=[Wx]OPEN=[Wx]CFDI=0
檢定使用sdm模型或是sem模型如果拒絕H0,則使用sdm檢定使用sdm模型或是sem模型如果拒絕H0,則使用sdm testnl ([Wx]FT = -[Spatial]rho*[Main]FT) ([Wx]GDP = -[Spatial]rho*[Main]GDP) ([Wx]WAGE = -[Spatial]rho*[Main]WAGE) ([Wx]OPEN = -[Spatial]rho*[Main]OPEN) ([Wx]CFDI = -[Spatial]rho*[Main]CFDI)
使用hausmam test檢定固定效果模型或隨機效果模型
先跑隨機效果模型,並儲存結果 xsmle FDI FT GDP WAGE OPEN CFDI , wmat(W) mode(sdm) re estimates store re
再先跑固定效果模型,並且儲存結果 xsmle FDI FT GDP WAGE OPEN CFDI , wmat(W) mode(sdm) fe estimates store fe
檢定使用固定效果模型或隨機效果模型 hausman re fe
Part4. Moran’s I findit spatwmat
讀取chinaprovince spatwmat, name(Chinaweights) xcoord(longitude) ycoord(latitude) band(0 8)
下載繪製Moran’s I工具 findit splagvar