请选择 进入手机版 | 继续访问电脑版
楼主: ruamingxi111
8896 6

空间面板回归——建立权重矩阵的方法-距离倒数的权重矩阵 [推广有奖]

  • 0关注
  • 1粉丝

初中生

28%

还不是VIP/贵宾

-

威望
0
论坛币
4 个
通用积分
7.9729
学术水平
3 点
热心指数
3 点
信用等级
0 点
经验
279 点
帖子
8
精华
0
在线时间
10 小时
注册时间
2016-3-15
最后登录
2020-4-29

ruamingxi111 发表于 2019-3-27 17:02:11 |显示全部楼层 |坛友微信交流群
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币


引言: 权重矩阵,是做空间面板回归、空间界面回归必须的文件



从陈强第二版的ussaww.spmat文件(做空间面板回归的)[url=],可以看到,[/url][x1] 是一个48个城市的行正规化的文件。



我们要做回归,50%的工作是要生成这样一个权重矩阵。


理论上 ,3种的计算模式,有临近、经济距离和距离倒数(我用的这个,其他两个做不出);



我的素材是,选了特定45个城市,希望做空间面板,首先想生成这样一个权重矩阵文件.


需要的数据是45个城市的X\Y坐标(百度 全国各城市经纬度坐标就找到);


需要的软件有:arcgis10.0; stata15; GEODA; EXCEL


视频教程软件: 华东理工大学,空间面板公开课。


http://mooc1.chaoxing.com/course/101569745.html?edit=false&knowledgeId=101614908&module=2#content





我实现的方法如下:



1  生成shp文件;by arcgis—— 必经之路。



adddata,导入经纬度坐标文件。  xx.xls (excel要是2003版本的文件簿文件最好,这样能识别)



显示坐标, x ,y已经自动的导入



添加投影仪坐标 wgs84( 如果不添加,后期添加到stata中,stata不自动识别经纬度坐标)



导出shp文件







这样我们就拿到了shp文件 ——全部放到stata的根目录下




# 此处 第一个目标,形成地理坐标的shp文件生成;




2、用stata



2、1


spshape2dta Export_Output_1,replace


# 把 Export_Output_1整体转换成stata的文件dta( 要能创建 CX,CYcoordinate,这样才能成功,这样才能dta 带有stata自身的地理坐标格式)





2.2  


spset _ID,modify replace


# 目的是,设置成空间数据,  _ID 是 北京—xx城市的序号名词( LONG 数据类型);


#错误解决1: 已经set了; 那就 spset,clear 清除掉,再set


# 错误解决2: 有可能spset _ID 就能set 了。



这是我们要出现的结果: coordinate 那里要有 _CX ,CY (此时默认的话,stata这两列没有数据,——2.3步我们自己补上)



2.3 统一说一下原始数据,stata的数据





我们要把X Y 数据 自己复制到,CXCY里面。




3、生成矩阵




结果如下。






4、反思路径



4.1 第三步,为什么不能用邻居矩阵?




显示没有发现neighbor,我觉得用GEODA软件的那个邻居可能有。


一开始我直接用xls 导入的 GEODA,





在这个 4‘ 的地方, 有个最近的neighbor,用xls 生成的文件我看不懂。


可能用shp是另外的生成方式,。


并且GEODA 导出的文件——还通过某一个步骤,转到stata 就能直接用


。我想突破了这个转换的步骤,by用shp; xls ,  就能生成新的权重矩阵—— 临近的。


#




[x1]查看方法,

先把 usaww.spmat (陈强空间面板中的 权重矩阵)

导入stata

前置的操作是,spmat的命令导入

(help spmat)——找到该命令,然后click here install 再用即可。)

使用文件: spmat use usaww usingusaww.spmat

导出成txt: spmat export usaww usingdaochu.txt











二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝


image015.png
image014.png
image012.png
image010.png
image008.png
image007.png
image006.png
image005.png
image003.png
image001.png
image016.jpg
image013.jpg
image011.jpg
image009.jpg
image004.jpg
image002.jpg
已有 1 人评分论坛币 学术水平 热心指数 收起 理由
20041910120 + 5 + 3 + 3 精彩帖子

总评分: 论坛币 + 5  学术水平 + 3  热心指数 + 3   查看全部评分

ruamingxi111 发表于 2019-3-27 17:02:50 |显示全部楼层 |坛友微信交流群
查看陈强那个权重矩阵的方法,

先把 usaww.spmat (陈强空间面板中的 权重矩阵)
导入stata
前置的操作是,spmat的命令导入
(help spmat)——找到该命令,然后click here install 再用即可。)

使用文件: spmat use usaww using usaww.spmat
导出成txt: spmat export usaww using daochu.txt

使用道具

楼主很认真。点赞,即使我没看懂

使用道具

花芽子 发表于 2019-5-28 15:40:43 |显示全部楼层 |坛友微信交流群
ruamingxi111 发表于 2019-3-27 17:02
查看陈强那个权重矩阵的方法,

先把 usaww.spmat (陈强空间面板中的 权重矩阵)
楼主,想问一下,如果是面板的空间权重矩阵,是不是也是一个n*n矩阵,我现在要做一个社会经济特征空间矩阵,首先我先计算出地理距离权重矩阵,后又算出经济对角阵,将两个矩阵相乘得到excel格式的社会经济特征矩阵是一个n*n的矩阵,接下来是不是把这个矩阵的数据放到stata中,以dta形式保存,然后用spamt命令生成spmat格式,这个就是spmat的面板空间权重矩阵?  不知道这样做对不对,困扰了我很久了
命令:
spmat dta w11 var1-var31
save w11.spmat

使用道具

花芽子 发表于 2019-5-28 15:40:46 |显示全部楼层 |坛友微信交流群
ruamingxi111 发表于 2019-3-27 17:02
查看陈强那个权重矩阵的方法,

先把 usaww.spmat (陈强空间面板中的 权重矩阵)
楼主,想问一下,如果是面板的空间权重矩阵,是不是也是一个n*n矩阵,我现在要做一个社会经济特征空间矩阵如图 ((ODSGPQN7%[}IZ%U_PV[AY.png
首先我先计算出地理距离权重矩阵,后又算出经济对角阵,将两个矩阵相乘得到excel格式的社会经济特征矩阵是一个n*n的矩阵,接下来是不是把这个矩阵的数据放到stata中,以dta形式保存,然后用spamt命令生成spmat格式,这个就是spmat的面板空间权重矩阵?  不知道这样做对不对,困扰了我很久了
命令:
spmat dta w11 var1-var31
save w11.spmat

使用道具

Yals_ 发表于 2020-2-22 11:06:13 |显示全部楼层 |坛友微信交流群
我在arcgus生成shp文件后导入stata,显示无法导入r603,为什么啊
求大神解答!!!非常急!!

使用道具

tianwk 发表于 2020-3-3 17:13:30 |显示全部楼层 |坛友微信交流群
thanks for sharing

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-28 16:36