先从excel中粘贴权重数据(不含地区名)到stata,另存为W.DTA(stata文件),存入stata安装文件中
读取权重文件:spatwmatusing'C:\Users\张燕\Downloads\计量软件\Stata14\DataSets-Stata-1111531048_364888\W.DTA', name(w)standardize。标黄部分是W.DTA文件存放位置。权重矩阵进行了行标准化。
读取变量文件:import excel'C:\Users\张燕\Desktop\workingpaper\“荷兰病”空间效应\荷兰病\论文写作\Moran'sI检验数据\1994年.xlsx',sheet('Sheet1') firstrow clear。标黄部分是变量(1994年)表格存放位置。表格中不含地区名,第一行是变量名。Moran'sI检验只能单个变量分年检验,所以有几年数据就要分开做几次。
Moran'sI检验:spatgsaRA,weights(w) moran。RA是变量,只需替换变量名就可以检验其他的变量。结果如下:
Moran's I
--------------------------------------------------------------
Variables | I E(I) sd(I) z p-value*
--------------------+-----------------------------------------
RA | 0.285 -0.034 0.108 2.959 0.002标红字体就是M值和P值。
--------------------------------------------------------------
*1-tailtest(单侧检验)
5. 局部Moran'sI散点图:spatlsaRA ,weight(w) moran graph (moran) symbol(n)
Symbol(n)表示每个点用序号表示(变量表格中没有地区序号,采用系统默认序号),可以换成“Symbol(id),id(state)”,此时变量表格中第一列要有“state”数据,即各省的名称。
Graph可以进行编辑,把西文字体调成新罗马。复制到Word中自动的变成白底黑字。
联系客服