求助！！怎么检验全局Moran指数呀？？急~~~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

貌似没有办法直接得出p检验值吧？如果用radomization做模拟的话，每次进行“run”都会得出不同的结果。我自己之前曾经做过一次空间分析的文章，当时是自己根据相应原理自创了一个得出p值的办法，但是得出的结果与用radomization模拟得出的结果不太一样，所以想请教各位高人帮忙判断下是不是可行？

     就是用变量x的空间滞后变量做因变量，x的标准值做自变量做了一个回归模型，然后在那里得出了自变量的系数（就是全局Moran指数）及其相应的p检验值，大家觉得这样做是不是对的呢？？或者说还有没有更直接的办法呀？？现在在做毕业论文特别着急，所以希望大家帮帮忙呀！！万分感激！！！！~~~~

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Moran指数 moran RAN 希望大家帮帮忙 ATION 毕业论文 因变量 自变量 文章 空间

自己先顶一个~~希望早日得到高人指点呀~~这个问题解决不了心理总是觉得很别扭~~纠结呀……

you can try R package "spdep"
Spatial dependence: weighting schemes, statistics and models
###############
library(spdep)
data(oldcol)
test1=moran.test(COL.OLD$CRIME, nb2listw(COL.nb, style="W"))
>test1
        Moran's I test under randomisation

data:  COL.OLD$CRIME  
weights: nb2listw(COL.nb, style = "W")  

Moran I statistic standard deviate = 5.6341, p-value = 8.797e-09
alternative hypothesis: greater
sample estimates:
Moran I statistic       Expectation          Variance
      0.510951264      -0.020833333       0.008908762

######## moran's I - exact test
test2=moran.test(COL.OLD$CRIME, nb2listw(COL.nb, style="W"),
+ randomisation=FALSE)
>test2
        Moran's I test under normality

data:  COL.OLD$CRIME  
weights: nb2listw(COL.nb, style = "W")  

Moran I statistic standard deviate = 5.6754, p-value = 6.92e-09
alternative hypothesis: greater
sample estimates:
Moran I statistic       Expectation          Variance
      0.510951264      -0.020833333       0.008779831

####自行运算如下
statistic standard deviate =(0.5109513-( -0.020833333 ))/sqrt( 0.008779831)=5.6754

alternative == "greater"
p=pnorm(5.6754,lower.tail=FALSE)=6.918248e-09

请问这个是要在stata里面做的吗？？有些复杂啊……

没接触过R,看似复杂其实简单.
约十五分钟一切OK.
http://cran.r-project.org/
Download and Install R
1.Windows
2.base
3.Download R 2.10.1 for Windows (32 megabytes)

安装完后,开启R,
在命令窗口键入install.packages("spdep")
就会自动安装packages("spdep")

package安装完后,就可开始运行,
library(spdep)
data(oldcol)
coords.OLD <- cbind(COL.OLD$X, COL.OLD$Y)
a=moran.test(COL.OLD$CRIME, nb2listw(COL.nb, style="W"),
  randomisation=FALSE)
a

详细function用法,参考spdep.pdf

我也想知道，stata中如何对moran 指数进行显著性检验

同问啊。如果用metlab.又怎么求

全局的莫兰指数很简单，要是求局部的莫兰指数及其相应的p值请看Anselin的1995年的一篇文章

全局的莫兰指数很简单，要是求局部的莫兰指数及其相应的p值请看Anselin的1995年的一篇文章

有谁帮忙做一下用matlab做莫兰指数的操作？必有重谢