- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 1672 个
- 通用积分
- 19.1517
- 学术水平
- 1 点
- 热心指数
- 1 点
- 信用等级
- 0 点
- 经验
- 261 点
- 帖子
- 11
- 精华
- 0
- 在线时间
- 160 小时
- 注册时间
- 2016-9-3
- 最后登录
- 2020-9-8
本科生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
- 1672 个
- 通用积分
- 19.1517
- 学术水平
- 1 点
- 热心指数
- 1 点
- 信用等级
- 0 点
- 经验
- 261 点
- 帖子
- 11
- 精华
- 0
- 在线时间
- 160 小时
- 注册时间
- 2016-9-3
- 最后登录
- 2020-9-8
| 开心 2017-9-3 17:40:34 |
---|
签到天数: 223 天 连续签到: 1 天 [LV.7]常住居民III
|
30论坛币
我自己写了一个程序如下:
x=read.csv("X.csv",sep = ';' , header = F) #数据读入
w=read.csv("w.csv",header = F) #矩阵读入
library(spdep)
w<-as.matrix(w) #转换为matrix
listw<-mat2listw(w) #转换为listw
a<-lagsarlm(V1~V2+V3+V4+V5+V6+V7+V8,data = x,listw = listw)
summary(a)
b<-errorsarlm(V1~V2+V3+V4+V5+V6+V7+V8,data = x,listw = listw)
summary(b)
测试数据见附件(不能传csv,我转化为xls传,大神们转换一下)
有几个问题
1.首先想问一下,这样跑SLM和SEM可以么?如果不行应该怎么设置?
2.那个出来的Warning message:
In errorsarlm(V1 ~ V2 + V3 + V4 + V5 + V6 + V7 + V8, data = x, listw = listw) :
inversion of asymptotic covariance matrix failed for tol.solve = 1e-10
倒条件数=1.46095e-14 - using numerical Hessian.
这个警告是什么意思?
3.如果跑可以的话,还想问问,这样summary里出来的东西,都是什么东西?(我只能看懂系数,标准差,Z值)如何对应下面这个表?(主要是R方,SC这些哪里找?)
4.最后还想问一下,哪些指标可以看空间相关性?选用SLM或者SEM是看哪个指标的?
谢谢大家的解答!!
|
|