R中moran.test函数怎么使用？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我不知道怎么生成空间权重矩阵
setwd("D:/R工作空间")
mydata<-read.table("分析.txt",stringsAsFactors=FALSE,header=TRUE)
distance<-read.table("省会距离.txt",header=FALSE)
####检查数据是否正确
distance
head(distance)
mydata
head(mydata)
####提取需要进行莫兰检验的数据
morandata<-mydata[,3]
####检验数据是否正确
head(morandata)
#####进行莫兰检验
distanceone<-as.matrix(distance)
library(spdep)
moran.test(morandata,nb2listw（ distanceone）)

请问这个nb2listw中的对象时什么类型的？我这里的distanceone是各省份距离的一个矩阵

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：moran test Est RAN 空间权重矩阵 我不知道 空间

求大神

啊啊啊啊啊啊啊啊

Error in nb2listw(w, style = "w") : Not a neighbours list怎么回事呢？求大神啊

z8938560 发表于 2015-10-22 16:46
Error in nb2listw(w, style = "w") : Not a neighbours list怎么回事呢？求大神啊

我也是这个错误，你解决了吗

求解

Error: unexpected input in "moran.test(morandata,nb2listw?

edsioncjlu 发表于 2017-2-19 08:06
求解

Error: unexpected input in "moran.test(morandata,nb2listw?

1. > x <- runif(20, 0, 10)
   
2. > y <- runif(20, 0, 10)
   
3. > obs <- rnorm(20)
   
4. > coords <- cbind(x, y)
   
5. > my.nb <- dnearneigh(coords, 0, 5)
   
6. > moran.test(obs, nb2listw(my.nb, style = "W"))
   
7. 
   
8.         Moran I test under randomisation
   
9. 
   
10. data:  obs  
    
11. weights: nb2listw(my.nb, style = "W")  
    
12. 
    
13. Moran I statistic standard deviate = -1.62, p-value = 0.9474
    
14. alternative hypothesis: greater
    
15. sample estimates:
    
16. Moran I statistic       Expectation          Variance
    
17.      -0.180124078      -0.052631579       0.006193256

复制代码

moran.test(y1,listw=mat2listw(W))

morantest（x，w）x是向量类型，如果不是用as.vector转化成向量，w空间权重矩阵是listw类型，原来是matric用mat2listw转换，nb用nb2listw转换。我目前是这样理解的，如有错误欢迎指正