为什么在geoda 内生成生成空间权重矩阵阵时,总是出现duplicate values

苹果/安卓/wp
积分 776, 距离下一级还需 24 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 269 天连续签到: 2 天[LV.8]以坛为家I
& && &&&新人第一次发帖!!& && &&&我相信计算Moran指数是每一位做空间计量的同学都需要掌握的技能,其中空间计量专门化软件Geoda为Moran指数的计算提供了相当便捷的方法。虽然lesage的工具包也可以做,但是毕竟没有Geoda专业,也无法做出想要的空间图。
& && &&&然而一开始我便遇到了一些问题,对于一个shp文件,我们需要把自己的变量和数据导进去,但事实上Geoda并未提供此功能。我在网上找了很多帖子都是建议用ArcGIS编辑shp文件然后导入数据,再返回到Geoda中。
& && &&&其实Geoda很方便的就可以导入数据,我以前不知道都是看教程用ArcGis处理的现在想想真的蠢哭了
& && &&&首先,你需要一个shp文件,和一个已经输入好你的数据的csv文件,注意:截面里每一个区域的数据编号和数据要对应好
& && &&&然后打开数据表格,右击选择merge table data,选择你的csv文件然后merge by record order就可以啦,当然你要是定义了key也可以用键来导入。
& && &&&这样就可以了有木有!!别告诉我只有我一个人现在才知道!
& && &&&楼主大三学生一枚,对计量尤其是空间计量十分感兴趣,大部分属于自学,遇到困难也只能求助于论坛和文献,实在不易。所以希望这个帖子能帮到同我一样对空间计量感兴趣的同学~也希望可以多多交流~~共勉
支持楼主:、
购买后,论坛将把您花费的资金全部奖励给楼主,以表示您对TA发好贴的支持
载入中......
很好,解决了一个难题!
鼓励积极发帖讨论
总评分:&论坛币 + 25&
学术水平 + 2&
热心指数 + 2&
信用等级 + 1&
本帖被以下文库推荐
& |主题: 12303, 订阅: 29
shgy88 发表于
只有分析数据,没有合适地图,只做空间自相关分析,不在地图上表现。按论坛经验:http://bbs.pinggu.org/th ...那我建议你用Matlab工具箱做就可以了,很方便的,如果你不需要在地图上表现的话。
制作shp文件得用ArcGIS,好像是要一个坐标数据一个变量数据,具体我也不会。自建权重矩阵Geoda有shp文件就可以自己生成,如果没有shp文件你可以建立简单的0-1权重矩阵,就是相邻为1不相邻为0,在标准化,如果你想更深刻一点,可以在权重矩阵里加入经济意义,参考林光平的《我国地区经济收敛的空间计量实证分析》,这就是所谓的经济空间权重矩阵。
shgy88 发表于
只有分析数据,没有合适地图,只做空间自相关分析,不在地图上表现。按论坛经验:http://bbs.pinggu.org/th ...噢噢你说的是自建权重文件。。。是.gal格式的对吧。。那我也不会额,可以请教论坛里高人
alvinyin 发表于
楼主,我正在使用merge table data导入数据,您说的程序我也懂,但是我在导入的时候出现了问题,我做的是中 ...你把shp和dbf文件里面的地级市用数字编号,作为key试一试
alvinyin 发表于
谢谢楼主的提点,我的问题解决了,么么哒嘿嘿有问题多多交流哦~
jyt8866989 发表于
楼主,您说的需要一个shp文件,和一个已经输入好你的数据的csv文件。shp文件可以下载的到,可是怎样找csv文 ...csv文件的数据根据你要研究的内容自己收集呀,比如你要想研究经济发展空间聚集性就是吧各个省的GDP数据导入呀
能源经济小白鼠 发表于
楼主你好,我也同是大三学生一枚~最近在学geoda做图,因为写论文需要。不过我按照楼主的方法保存csv文件后, ...你看看这样行不行,
1.csv文件名最好是数字或者英文
2.文件路径不能有中文,或者你直接保存在桌面上试试
3.确定一下你merge data选择的文件类型是csv数据吗,不是的话你就选所有数据类型然后选中你的文件
不知道这样能不能解决
jyt8866989 发表于
谢谢,我明白您的意思,可是具体怎么生成csv文件呢求指教?额,你只需要在excel里另存为csv就可以了。。。。
能源经济小白鼠 发表于
我在将Excel的数据保存为CSV的时候总是会出现问题。具体情况是每次在保存的适合,系统会提示说我原本的 ...谢谢你~~~嘿嘿
不过说实话Geoda做空间计量的局限性太大了,stata作图可以,具体分析如果靠两个埃及狂魔编的程序也不行,matlab无法作图,分析的话LeSage并没有考虑内生性问题,面板的2SLS、GMM等估计也没有提供。强烈建议用R,因为你研究的深了肯定会需要用到R来做自己独一无二的分析的,如果以后反正要用,还不如现在就开始学,虽然学习曲线要比Geoda陡峭的多,但物有所值。
<font color="#13133 发表于
亲 看了你的帖子 有三个问题 (问题有点白,刚学)
1、shp文件怎么来?论坛下么?还是自己怎么生成呢?
2 ...shp文件自己去下载,现在很容易下的,你要精确到地级市的现在也有很多,但是数据匹配起来有点复杂。
200多个市的话我研究的不深,不太敢说,但是你要知道现在中国有300个左右地级市,你要是做这个的话将会面临相当大的缺失数据问题,这个解决方法Lesage(1999)的书里提到过,但这个也是现在空间计量被诟病的地方,具体你自己去看吧~你要是用城市间的邻接矩阵,还有一个问题就是会得到一个稀疏矩阵,对这个的处理又是一门学问,所以相当复杂。
stata是可以做Moran散点图的,你findit spatlsa 和 spmap,自己摸索去吧
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师GeoDa软件关于空间权重矩阵
核心提示:GeoDa软件关于空间权重矩阵主要分为两类:基于邻接关系和基于距离关系的空间权重矩阵。一、基于邻接(Contiguity)关系的空间权重1.rook...
GeoDa软件关于空间权重矩阵主要分为两类:基于邻接关系和基于距离关系的空间权重矩阵。一、基于邻接(Contiguity)关系的空间权重1.rook:表示共边为邻接;【与arcgis中的Polygoncontiguity(edgesonly)对应】2.queen:表示共边或共点为邻接;【与arcgis中的Polygon contiguity(edgesandcorners)对应】3.bishop(软件中没有):表示共点为邻接。二、基于距离(Distance)关系的空间权重1.threshold:表示既定距离下的相关,一般软件有默认一个最小值,但可视实际情况调整(一般应大于最小值);【arcgis中的inverse distance/fixed distance都要设置threshold距离】2.K-nearest:表示指定某个多边形周围的多边形个数(K=3,4,5,6...)。【与arcgis中的K-nearestneighbors对应】PS:1.在建立空间权重矩阵过程中,对于关键字变量【俗称的ID变量】的选择十分重要,这是两个软件同时提到的,因为ID是用于判别每个变量位置的唯一性,所以应该是每个变量都不一样的,大多情况下以数字表示(N=0,1,2,3.......):基于邻接关系和基于距离关系的空间权重矩阵。一、基于邻接(Contiguity)关系的空间权重1.rook:表示共边为邻接;【与arcgis中的Polygoncontiguity(edgesonly)对应】2.queen:表示共边或共点为邻接;【与arcgis中的Polygon contiguity(edgesandcorners)对应】3.bishop(软件中没有):表示共点为邻接。二、基于距离(Distance)关系的空间权重1.threshold:表示既定距离下的相关,一般软件有默认一个最小值,但可视实际情况调整(一般应大于最小值);【arcgis中的inverse distance/fixed distance都要设置threshold距离】2.K-nearest:表示指定某个多边形周围的多边形个数(K=3,4,5,6...)。【与arcgis中的K-nearestneighbors对应】PS:1.在建立空间权重矩阵过程中,对于关键字变量【俗称的ID变量】的选择十分重要,这是两个软件同时提到的,因为ID是用于判别每个变量位置的唯一性,所以应该是每个变量都不一样的,大多情况下以数字表示(N=0,1,2,3.......)苹果/安卓/wp
积分 917, 距离下一级还需 458 积分
权限: 自定义头衔, 签名中使用图片, 隐身
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 设置帖子权限道具: 提升卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 81 天连续签到: 1 天[LV.6]常住居民II
第一次接触空间计量。请问一下,对于面板数据,空间权重矩阵怎么生成?已经有了全国省份相邻与否的excel表,但只适用一个截面,面板数据怎么办。。。求解答
载入中......
But some things are worth fighting for
同求。有人说要看这篇文章QML Estimation of Spatial Dynamic Panel Data Models with Time Varying Spatial Weights Matrices。网址:/doi/pdf/10.2.#.U-bPUNL59m0。我没有下载权限,你试试。
热心帮助其他会员
总评分:&热心指数 + 1&
解决了吗楼主,同求
菜园一块田 发表于
解决了吗楼主,同求我当时实在没招 就手动拉了一个对角矩阵
之后也没有去钻研呢。。
要是解决了的话,记得share一下呀
<font color="#5102137 发表于
同求。有人说要看这篇文章QML Estimation of Spatial Dynamic Panel Data Models with Time Varying Spatia ...资料我下了,但是还没看。。
18:05:59 上传
But some things are worth fighting for
kocoma 发表于
资料我下了,但是还没看。。这篇文章太专业,看不懂
有人会了吗&&好惆怅啊
请问楼主解决这个问题了吗? 我只有30个地区 30*30的相邻矩阵,不知道面板的矩阵应该怎样得到。好像要用spmat命令,但具体怎么用那?请指教谢谢
http://econpapers.repec.org/software/bocbocode/s457325.htm
听说这个可以做, SPWEIGHT。仅供参考
o5104369 发表于
http://econpapers.repec.org/software/bocbocode/s457325.htm
听说这个可以做, SPWEIGHT。仅供参考不好意思,看出了,我这个是R中做截面权重矩阵的
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师苹果/安卓/wp
积分 2127, 距离下一级还需 98 积分
权限: 自定义头衔, 签名中使用图片, 隐身, 设置帖子权限
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡下一级可获得
权限: 设置回复可见道具: 沉默卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 205 天连续签到: 1 天[LV.7]常住居民III
本帖最后由 wanghaidong918 于
17:20 编辑
大家好,最近在做空间面板回归,但是生成空间权重后做回归,老是显示“Weight Matrix is not Square”求问各位达人有没有解决办法?因为原来做geoda时感觉geoda生成空间权重比较简单,是否有办法把geoda生成的权重文件.gal引入到stata。非常谢谢
载入中......
Geoda可以生成rook、queen权重矩阵,但转换到stata或matlab里面需要处理下,根据Geoda生成的gal文件手工形成邻接矩阵太麻烦,也容易产生误差,可以java的gal2WM的小程序转换。
本帖被以下文库推荐
& |主题: 12303, 订阅: 29
Geoda可以生成rook、queen权重矩阵,但转换到stata或matlab里面需要处理下,根据Geoda生成的gal文件手工形成邻接矩阵太麻烦,也容易产生误差,可以java的gal2WM的小程序转换。
10:08:05 上传
热心帮助其他会员
总评分:&经验 + 10&
论坛币 + 10&
我下载了,试试看,谢谢分享
下载一下,学一下
点了转换也没动静,楼主帮忙看一看?
11:06:46 上传
yhw1234 发表于
Geoda可以生成rook、queen权重矩阵,但转换到stata或matlab里面需要处理下,根据Geoda生成的gal文件手工形成 ...最近在做空间计量分析,但是空间权重矩阵还是不知道如何从gal文件导入stata中,麻烦可以说明下如何使用gal2WM程序吗……
需要安装Java运行环境程序
spwmatrix creates distance-based spatial weights and imports adjacency-based weights (gal files saved in sparse format)
lishujuan 发表于
点了转换也没动静,楼主帮忙看一看?我也是点了转换没动静,请问你后来解决问题了么,可以告诉下不
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师苹果/安卓/wp
积分 1277, 距离下一级还需 98 积分
权限: 自定义头衔, 签名中使用图片, 隐身
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 设置帖子权限道具: 提升卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
开心签到天数: 19 天连续签到: 1 天[LV.4]偶尔看看III
本帖最后由 fei355 于
18:54 编辑
以下都用中国的例子
(一)得到数据库文件和坐标文件
1.之所以这一步是由于在实际分析中会删除个别省份,然后基于剩余省份制作shp文件后生成空间权重文件。找到中国地图省级、省会城市、城市、县的dbf和shp文件,注意:各级地图的dbf和shp文件的名称必须相同且是英文,如province.dbf和province.shp,provcap.dbf和provcap.shp等2.在stata中安装spmap、shp2dta和mif2dta命令是ssc install spmap\& &shp2dta\mif2dta
3.在指定目录下,运行命令shp2dta using ProCap, data(&ProCap_data&)&&coor(&ProCap_coordinates&),其中下划线部分分别是使用的数据库和坐标名称,转换后的数据库和坐标文件名称,然后用命令 merge ID using procap_coordinates,将文件复制到txt文档中编辑成符合geoda格式(第一行写上截面单位数和变量数,第二行是变量名称,用英文下双引号引起,用逗号隔开,剩下行是观测值也要用逗号隔开,当然可以先在excel中编辑再拷过来)然后打开geoda~tools~shape~points from ascii,主要设置位置,然后保存为shp文件,最后基于此文件生成空间权重矩阵,当然质点文件只能生成knearest和threshold distance距离权重。
(二)如果用geoda生成了权重数据,那么可以通过以下方法转换成stata文件
1.用geoda生成不同类型的权重文件,比如province_knear.gwt,打开另存为province_knear.txt
2.用记事本打开后复制出来或者在stata中输入命令insheet using province_knear.txt,删除关键词,保留位置变量和邻近数量变量,并重命名为v1和v2,并另存为spweight3.dta
3.在stata中安装spweightxt命令
<font color="#.输入命令& &
& && & clear all
& && &&&sysuse spweight3.dta, clear
& && &&&list v1 v2
& && &&&spweightxt v1 v2 , panel(34)&&time(8) matrix(W) eigw table
就可以生成中国34个省份观测单位8年的权重指标了。
注意:其中panel()项可以设置你要生成的权重观测单位数一般与原始文件一致,如果想减少几个省份,则先在原始文件减少,time()设置生成时间的维数,一般空间面板数据用的上。然后生成的权重数据分为横截面和面板类型分别保存在 Cross Section Weight Matrix File (Wcs.dta)和 Panel Weight Matrix File& & (Wxt.dta)中要查看只需输入 use Wcs.dta,clear 或者use Wxt.dta ,clear即可
(三)空间面板回归
主要说gmm估计
命令是 spgmmxt y&&x1 x2 , nc(7) wmfile(SPWxt) gmm(1) mfx(lin) test
暂时就说这么多,有空了再写详细点,主要想把在STATA中如何进行空间面板的操作写一下!
载入中......
对论坛有贡献
热心帮助其他会员
鼓励积极发帖讨论
总评分:&经验 + 180&
论坛币 + 90&
学术水平 + 15&
热心指数 + 14&
信用等级 + 14&
本帖被以下文库推荐
& |主题: 753, 订阅: 199
& |主题: 151, 订阅: 78
& |主题: 106, 订阅: 32
& |主题: 105, 订阅: 43
非常感谢楼主!期待您写出运用stata分析空间面板数据模型的详细步骤,如原始数据的导入、几类空间面板权重的生成、moran\G指数的计算,真是大功德一件了!
纵一苇之所如,凌万顷之茫然。
首先非常感谢楼主提供的宝贵经验!其次我仔细揣摩你提供的方法,有两点疑问,希望你能够帮忙解答。其一,你在第一步转换得到的数据库和坐标文件的作用是什么,或者说它与第二步的关系是什么。其二,我们知道SPWEIGHT命令的V1和V2只能是0或者1,那么请问若是构建地理距离权重矩阵或者经济距离权重矩阵,V1和V2又该如何设定呢,若是不能设定,那么如何得出面板数据的空间权重矩阵呢,这些问题一直困扰我,所以希望你能帮忙解答一下。总结来说,就是我该如何通过STATA构建地理/经济距离权重矩阵,详细步骤是什么,麻烦你了。
鼓励积极发帖讨论
总评分:&论坛币 + 10&
热心指数 + 1&
zhichi&&ge&&xiexie
寻梦?撑一支长篙,向青草更青处漫溯;满载一船星辉,在星辉斑斓里放歌。
it is goooooooooooooood
不错,期待LZ接下来的详解
楼主 我在stata中输入spweightxt v1 v2 , panel(34)&&time(5) matrix(W) ptable,得不出结果。
显示matrix __00000K11 would have missing values
楼主帮我看看。geoda生成的权重文件在附件中。谢谢!
16:57:43 上传
geoda生成的权重文件(按v1和v2整理的文件)。
17:03:51 上传
热心帮助其他会员
总评分:&论坛币 + 20&
学术水平 + 3&
谢谢好经验~!!!!!!!!!!
非常感谢您的分享,但是我想请问, spweightxt v1 v2 , panel(34)&&time(8) matrix(W) eigw table中的“ eigw table”是什么意思呢?
鼓励积极发帖讨论
总评分:&经验 + 10&
热心指数 + 1&
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师

我要回帖

更多关于 如何生成空间权重矩阵 的文章

 

随机推荐