全部版块 我的主页
› 论坛 › 计量经济学与统计论坛 五区 › 计量经济学与统计软件 › 经管代码库
MATLAB实现：moran'I值、moran'I值检验、moran散点图
天外来客2013
61884 147
2015-02-04
一、用MATLAB计算moran'I值


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

函数头：
function moranI=moranI1(X,W)
%%X为列向量，W为权重矩阵
相应MATLAB代码：
moranI1.rar
大小:(297 Bytes)

只需: 5 个论坛币  马上下载

本附件包括：

  • moranI1.m



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

相应MATLAB代码：
moranI2.rar
大小:(311 Bytes)

只需: 5 个论坛币  马上下载

本附件包括：

  • moranI2.m



二、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









三、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






四、总文件和例子演示


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

    moran总文件程序：
moran.rar
大小:(298.86 KB)

只需: 25 个论坛币  马上下载

本附件包括：

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



    实例演示截图:

实例演示1


实例演示2



实例演示3
                          moran散点图


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

附件列表

moran.rar

大小:1.29 KB

只需: 20 个论坛币  马上下载

本附件包括：

  • moran.m
  • moranScatterPlot.m
  • moran_test.m

二维码

扫码加我 拉你入群

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

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

全部回复
changjiang11
2015-2-6 20:51:54
elhorst的网站有源代码 你还这么贵！！！！
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

天外来客2013
2015-3-4 16:29:06
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;
二维码

扫码加我 拉你入群

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

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

sunnycelia
2015-3-16 17:13:56
楼主运行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:21:21
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.
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群