- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 28138 个
- 通用积分
- 11.1500
- 学术水平
- 0 点
- 热心指数
- 1 点
- 信用等级
- 1 点
- 经验
- 283 点
- 帖子
- 24
- 精华
- 0
- 在线时间
- 196 小时
- 注册时间
- 2009-7-8
- 最后登录
- 2024-11-28
已卖:17份资源
本科生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 28138 个
- 通用积分
- 11.1500
- 学术水平
- 0 点
- 热心指数
- 1 点
- 信用等级
- 1 点
- 经验
- 283 点
- 帖子
- 24
- 精华
- 0
- 在线时间
- 196 小时
- 注册时间
- 2009-7-8
- 最后登录
- 2024-11-28
|
20论坛币
|
最近在学习用R做空间计量经济分析,采取的步骤是先用GeoDa生成gal格式文件,关键字选择用province,生成结果为附件中的“province_mailand_rook.txt”。然后,再用R将其转化为空间权重矩阵,代码是:
- library(readxl)
- p0<-read_excel("F:/空间计量练习/省区练习数据.xlsx",sheet=1,na="NA")
- library(spdep)
- polgal<-read.gal("F:/GIS基础数据/省区空间权重/province_mailand_rook.gal",override.id=TRUE)#读取gal文件
- swmat1_pr<-nb2mat(polgal,zero.policy=TRUE)View(swmat1_pr)
复制代码 在生成的swmat1_pr中,行的名字是Anhui、Henan、Hubei、Jiangsu等省份英文名,但列的名字是V1、V2、V3等,不是省份英文名(见附件中的swmat1_pr.xls)。是不是我在转换的时候,漏掉了重要参数,或者是其他原因?能否把列名也生成为省份英文名?还望各位前辈不吝赐教,多谢!
swmat1_pr.xls
(34 KB)
province_mailand_rook.txt
(1.42 KB)
|
|