空間計量經濟學創造性地處理了經典計量方法在面對空間數據時的缺陷,考察了數據在地理觀測值之間的關聯。近年來在人文社會科學空間轉向的大背景下,空間計量已成為空間綜合人文學和社會科學研究的基礎理論與方法,尤其在區域經濟、房地產、環境、人口、旅遊、地理、政治等領域,空間計量成為開展定量研究的必備技能。
1、數據介紹:
boston.c數據有506行和20列。它包含Harrison和Rubinfeld(1978)的數據,修正了一些小錯誤,並增加了觀測的經度和緯度。Gilley和Pace還指出,MEDV被審查,在$ 50,000或以上的中值被設置為$ 50,000。沒有修改的原始數據集也包括在包mlbench作為BostonHousing。
數據描述:
2、查看數據
> str(boston.c)'data.frame': 506 obs. of 20 variables: $ TOWN : Factor w/ 92 levels "Arlington","Ashland",..: 54 77 77 46 46 46 69 69 69 69 ... $ TOWNNO : int 0 1 1 2 2 2 3 3 3 3 ... $ TRACT : int 2011 2021 2022 2031 2032 2033 2041 2042 2043 2044 ... $ LON : num -71 -71 -70.9 -70.9 -70.9 ... $ LAT : num 42.3 42.3 42.3 42.3 42.3 ... $ MEDV : num 24 21.6 34.7 33.4 36.2 28.7 22.9 27.1 16.5 18.9 ... $ CMEDV : num 24 21.6 34.7 33.4 36.2 28.7 22.9 22.1 16.5 18.9 ... $ CRIM : num 0.00632 0.02731 0.02729 0.03237 0.06905 ... $ ZN : num 18 0 0 0 0 0 12.5 12.5 12.5 12.5 ... $ INDUS : num 2.31 7.07 7.07 2.18 2.18 2.18 7.87 7.87 7.87 7.87 ... $ CHAS : Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1 1 1 ... $ NOX : num 0.538 0.469 0.469 0.458 0.458 0.458 0.524 0.524 0.524 0.524 ... $ RM : num 6.58 6.42 7.18 7 7.15 ... $ AGE : num 65.2 78.9 61.1 45.8 54.2 58.7 66.6 96.1 100 85.9 ... $ DIS : num 4.09 4.97 4.97 6.06 6.06 ... $ RAD : int 1 2 2 3 3 3 5 5 5 5 ... $ TAX : int 296 242 242 222 222 222 311 311 311 311 ... $ PTRATIO: num 15.3 17.8 17.8 18.7 18.7 18.7 15.2 15.2 15.2 15.2 ... $ B : num 397 397 393 395 397 ... $ LSTAT : num 4.98 9.14 4.03 2.94 5.33 ...3、此數據框包含以下列:
4、OLS回歸分析
data(boston)
hr0 <- lm(log(MEDV) ~ CRIM + ZN + INDUS + CHAS + I(NOX^2) + I(RM^2) + AGE + log(DIS) + log(RAD) + TAX + PTRATIO + B + log(LSTAT), data = boston.c)summary(hr0)logLik(hr0)gp0 <- lm(log(CMEDV) ~ CRIM + ZN + INDUS + CHAS + I(NOX^2) + I(RM^2) + AGE + log(DIS) + log(RAD) + TAX + PTRATIO + B + log(LSTAT), data = boston.c)summary(gp0)logLik(gp0)lm.morantest(hr0, nb2listw(boston.soi))5、空間回歸模型
5.1 空間誤差回歸模型
gp1 <- errorsarlm(log(CMEDV) ~ CRIM + ZN + INDUS + CHAS + I(NOX^2) + I(RM^2) + AGE + log(DIS) + log(RAD) + TAX + PTRATIO + B + log(LSTAT), data=boston.c, nb2listw(boston.soi), method="Matrix", control=list(tol.opt = .Machine$double.eps^(1/4)))結果為:
5.2 空間滯後回歸模型
gp1 <- errorsarlm(log(CMEDV) ~ CRIM + ZN + INDUS + CHAS + I(NOX^2) + I(RM^2) + AGE + log(DIS) + log(RAD) + TAX + PTRATIO + B + log(LSTAT), data=boston.c, nb2listw(boston.soi), method="Matrix", control=list(tol.opt = .Machine$double.eps^(1/4)))結果為:
1、數據介紹:
Produc
相關數據來源於:Munnell A (1990). 「Why Has Productivity Growth Declined? Productivity and Public Investment.」 New England Economic Review, 3–22.
2、面板數據分析
3、空間面板數據分析
3.1隨機效應空間面板分析
spml.re <-spml(fm, data = Produc, index =c("state","year"),listw=mat2listw(usaww),model ="random",effect ="individual",lag=T,spatial.error="b")summary(spml.re)結果為:
3.2固定效應空間面板分析
spml.fe <-spml(fm, data = Produc, index =c("state","year"),listw=mat2listw(usaww),model ="within",effect ="individual",lag=T,spatial.error="b")summary(spml.fe)結果為:
3、參考資料:
空間計量八部曲:首套空間計量體系課程
基礎篇:《空間計量經濟學入門:空間計量及Geoda應用》
第一部請點擊:《零基礎|輕鬆搞定空間計量:空間計量及GeoDa、Stata應用》
第二部請點擊:《空間計量及Matlab課程》
第三部請點擊:《第三期空間計量及stata應用精品課程》
第四部請點擊:《空間計量及ArcGis應用課程》
第五期:《空間計量經濟學》(註:第五期更新內容現免費贈送給第一期學員等,可以在第一期課程主頁目錄免費觀看學習)
第六期:《空間計量經濟學及R應用》
第七期:《空間計量經濟學及Python應用》
第八期:《高級空間計量經濟學》