“第八届空间计量研讨会”由南通祐禟计量有限公司主办,ASEF承办,将于2025年1月16~17日采用网络课堂的形式在江苏省南通市举办。诚邀您届时参会!过大年,选年货!32讲空间计量课程就是最好的知识年货哦!第二阶段缴费,依然有一定的优惠幅度。详情参见“第八届空间计量研讨会”第一号公告。
一、生成空间权重矩阵
Create a row-standardized binary spatial weights matrix assuming spherical coordinates and a distance cut-off of 10 miles(生成行标准化的二元空间权重矩阵)
. use 31_provinces_longitude&latitude,clear
.spwmatrix gecon latitude longitude, wn(wbin) wtype(bin) db(0 550) r(3958.761) row
.clear
.svmat wbin
.save wbin
Generate an inverse distance squared spatial weights matrix using projected latitudes and longitudes(生成反距离空间权重矩阵)
. use 31_provinces_longitude&latitude,clear
. spwmatrix gecon latitude longitude, wname(winvsq) wtype(inv) alpha(1) dband(0 100) cart
. clear
. svmat winvsq
. save winvsq
(一)计算全域Moran's I
1.使用反距离空间权重矩阵
. use winv.dta,clear
. spatwmat using winv.dta,name(winv) standardize
. use 31_provinces_core_variable_chn.dta,clear
. spatgsa pcrgdp_2002 pcrgdp_2003 pcrgdp_2004 pcrgdp_2005 pcrgdp_2006 pcrgdp_2007 pcrgdp_2008 pcrgdp_2009 pcrgdp_2010 pcrgdp_2011 pcrgdp_2012 pcrgdp_2013 pcrgdp_2014 pcrgdp_2015 pcrgdp_2016 ,weights(winv) moran(见图1)
图1
2.使用二元空间权重矩阵
. use wbin.dta,clear
. spatwmat using wbin.dta,name(wbin) standardize
. use 31_provinces_core_variable_chn.dta,clear
. spatgsa pcrgdp_2002 pcrgdp_2003 pcrgdp_2004 pcrgdp_2005 pcrgdp_2006 pcrgdp_2007 pcrgdp_2008 pcrgdp_2009 pcrgdp_2010 pcrgdp_2011 pcrgdp_2012 pcrgdp_2013 pcrgdp_2014 pcrgdp_2015 pcrgdp_2016 ,weights(wbin) moran(见图2)
图2
(二)计算局域Moran's I
. use winv.dta,clear
. spatwmat using wbin.dta,name(winv) standardize
. use 31_provinces_core_variable_chn.dta,clear
. spatlsa pcrgdp_2002, w(winv) moran twotail(见图3)
图3
(一)使用二元空间权重矩阵
. use wbin.dta,clear
. spatwmat using wbin.dta,name(wbin) standardize
. use 31_provinces_core_variable_chn.dta,clear
. spatlsa pcrgdp_2002,weights(wbin) moran graph(moran) symbol(id) id(prvn)(见图4)
. spatlsa pcrgdp_2002,weights(wbin) moran graph(moran) symbol(id) id(prvn) twotail
图4
. use winv.dta,clear
. spatwmat using winv.dta,name(winv) standardize
. use 31_provinces_core_variable_chn.dta,clear
. spatlsa pcrgdp_2002,weights(winv) moran graph(moran) symbol(id) id(prvn)(见图5)
. spatlsa pcrgdp_2002,weights(winv) moran graph(moran) symbol(id) id(prvn) twotail
图5