MATLAB实现：moran'I值、moran'I值检验、moran散点图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

一、用MATLAB计算moran'I值

---

1、权重矩阵(W)未标准化后的moran'I值
   计算公式：

函数头：
function moranI=moranI1(X,W)
%%X为列向量，W为权重矩阵
相应MATLAB代码： moranI1.rar (297 Bytes, 需要: 5 个论坛币) 本附件包括：

• moranI1.m

2015-2-4 18:07:34 上传

权重矩阵未标准化后的moran'I值程序
需要: 5 个论坛币  [购买]

2、权重矩阵(W)标准化后的moran'I值
a、权重矩阵标准化方法：W每行的值都除以该行的和即可。
b、计算公式：
     注：权重标准化后，moran'I可确保在[-1,1]之间，否则不敢保证。
函数头:
function moranI_standard = moranI2(X,W)
%%X为列向量，W为权重矩阵

相应MATLAB代码： moranI2.rar (311 Bytes, 需要: 5 个论坛币) 本附件包括：

• moranI2.m

2015-2-4 18:07:38 上传

权重矩阵标准化后的moran'I值程序
需要: 5 个论坛币  [购买]

二、moran'I值的检验

---

函数头：
function [Z_I,Z_alpha] = moran_test(I,W,alpha)
%%I为计算出来的moran'I值
%W为计算相应moran'I值的权重矩阵

%alpha为显著性水平
相应MATLAB代码： moran_test.rar (356 Bytes, 需要: 5 个论坛币) 本附件包括：

• moran_test.m

2015-2-4 18:19:48 上传

moran'I值的检验程序
需要: 5 个论坛币  [购买]

三、moran散点图

---

函数头：
function  result = moranScatterPlot(X,WX)
    %%X为列向量，W为权重矩阵
    %WX是W与X的乘积
    %X,W根据个人需要选择性的标准化
    %X标准化方法:(X-X_mean)/X_std

    %result为回归结果
相应MATLAB代码： moranScatterPlot.rar (319 Bytes, 需要: 5 个论坛币) 本附件包括：

• moranScatterPlot.m

2015-2-4 18:28:47 上传

moran散点图程序
需要: 5 个论坛币  [购买]

四、总文件和例子演示

---

    moran文件中包含了上述所有的程序和一篇理论性文档指导,文档中有详细的例子说明，在此，仅将例子作为演示，演示如下。

    moran总文件程序： moran.rar (298.86 KB, 需要: 25 个论坛币) 本附件包括：

• www.uni-kassel.de_~rkosfeld_lehre_spatial_SpatialEconometrics2.pdf
• moran.m
• moranScatterPlot.m
• moran_test.m

2015-2-4 18:31:43 上传

包含上述所有程序以及详细的理论说明文档和例子
需要: 25 个论坛币  [购买]

    实例演示截图:

                          moran散点图


    说明:上述MATLAB程序中的理论性知识都是根据moran文件中包含的pdf文档。除了代码中的bug之外，其它情况本人概不负责。

扫码加我 拉你入群

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

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

分享0 收藏76 回帖

关键词：matlab实现 moran散点图 MATLAB moran matla MATLAB moran散点图 moran值 moran值检验

elhorst的网站有源代码 你还这么贵！！！！

你好 moran散点图里面的  result = myls(X,WX)  ，myls函数是什么

joan宝宝 发表于 2015-3-4 13:47
你好 moran散点图里面的  result = myls(X,WX)  ，myls函数是什么

这个是用最小二乘求回归系数的。
源码是这样的：
function result = myls(X,Y)
%X,Y均为列向量
%y = ax + b
n = length(X);
A = [X,ones(n,1)];
result = inv(A'*A)*A'*Y;

楼主运行Z_I = moran_test(I_standard,W_standard);
Z_alpha = norminv(1-alpha,0,1);
出现错误Undefined command/function 'moran_test'.
为什么啊？

sunnycelia 发表于 2015-3-16 17:13
楼主运行Z_I = moran_test(I_standard,W_standard);
Z_alpha = norminv(1-alpha,0,1);
出现错误Undefined ...

之前我忘记放在工作夹里了，但放进去后就出现这样错误了,求解答
Z_I = moran_test(I_standard,W_standard)
??? Input argument "alpha" is undefined.

Error in ==> moran_test at 11
Z_alpha = norminv(1-alpha,0,1);
>> Z_alpha = norminv(1-alpha,0,1)
??? Error using ==> alpha
Too many output arguments.

sunnycelia 发表于 2015-3-16 17:21
之前我忘记放在工作夹里了，但放进去后就出现这样错误了,求解答
Z_I = moran_test(I_standard,W_standar ...

明天早上给你解答啊，亲～

天外来客2013 发表于 2015-3-16 23:51
明天早上给你解答啊，亲～

Z_I = moran_test(I_standard,W_standard)
??? Input argument "alpha" is undefined.
这是因为moran_test函数中少了alpha这个参数，如下所示改就行了：
把你moran.m文件中的：
Z_I = moran_test(I_standard,W_standard);
改为：
Z_I = moran_test(I_standard,W_standard,alpha);

大家好：
       针对上述朋友问的两个问题，在此做个总结：
       问题1：
      “Z_I = moran_test(I_standard,W_standard)
       ??? Input argument "alpha" is undefined.”
       改正：
              把moran.m文件中的：
              Z_I = moran_test(I_standard,W_standard);（大概是第19行左右）
              改为：
              Z_I = moran_test(I_standard,W_standard,alpha);（大概是第19行左右）
       问题2：
        “??? Undefined function or method 'myls' for input
        arguments of type 'double'.”
        改正：
        新建myls.m文档，把如下代码粘贴进去：
        function result = myls(X,Y)
        %X,Y均为列向量
        %y = ax + b
        n = length(X);
        A = [X,ones(n,1)];
        result = inv(A'*A)*A'*Y;

天外来客2013 发表于 2015-3-17 08:36
Z_I = moran_test(I_standard,W_standard)
??? Input argument "alpha" is undefined.
这是因为moran_t ...

那个alpha是假设检验时的显著性水平，一般为0.05。