matlab软件moran值,用matlab求moran并作图的具体步骤
现在有从论坛上下载的moran 的代码,如下:求问这个代码如何修正,运行总是有问题。这个代码只能做单变量的moran 吗?
function[I,I_standard,Z_I,Z_alpha,result]=moran(X,W,alpha)
%%X为列向量,W为权重矩阵
%alpha为显著性水平
n = length(X);
Xmean = mean(X);
XminusX_mean = X - Xmean;
S0 = sum(sum(W));
%一、计算,moran'I值
%1、W未标准化,为0-1矩阵
I =XminusX_mean'*W*XminusX_mean*n/(XminusX_mean'*XminusX_mean*S0);
%2、W标准化,为(0,1)矩阵
W_rowSum = sum(W');
W_rowSum_temp = W_rowSum'*ones(1,n);
W_standard = W./W_rowSum_temp;
I_standard =XminusX_mean'*W_standard*XminusX_mean/(XminusX_mean'*XminusX_mean);
%二、moran'I值的显著性检验
Z_I = moran_test(I_standard,W_standard);
Z_alpha = norminv(1-alpha,0,1);
%三、画moran散点图
%y = ax + b;线性拟合
%1、X标准化,W标准化
Xstd = std(X,1);%X的标准差
X_std = (X-Xmean)/Xstd;%X标准化
W_standard_X_standard = W_standard*X_std;
%result = myls(X_std,WX_standard);
result = moranScatterPlot(X_std, W_standard_X_standard );
title('X标准化,W标准化');
% %2、X标准化,W未标准化
% Xstd = std(X,1);%X的标准差
% X_std = (X-Xmean)/Xstd;%X标准化,X_std代表横坐标
% W_X_standard = W*X_std;%WX_standard代表纵坐标
% % result = myls(X_std,WX_standard);
% result = moranScatterPlot(X_std,W_X_standard);
% title('X标准化,W未标准化');
%3、X为未标准化,W标准化
% W_standard_X = W_standard*X;%WX_standard代表纵坐标
% % result = myls(X_std,WX_standard);
% result = moranScatterPlot(X,W_standard_X)
% title('X为未标准化,W标准化');
% % 4、X未标准化,W未标准化
% WX = W*X;%WX_standard代表纵坐标
% result = moranScatterPlot(X,WX);
% title('X未标准化,W未标准化');
Moran I 检验
function [Z_I,Z_alpha] =moran_test(I,W,alpha)
%%I为计算出来的moran'I值
%W为计算相应moran'I值的权重矩阵
n = size(W,1);
E_I = -1/(n-1);
S0 = sum(sum(W));
S1 = sum( sum( (W+W').^2 ) )/2;
S2 =sum( sum((W+W')').^2 );
Var_I =(n*n*S1-n*S2+3*S0*S0)/((n*n-1)*S0*S0)-E_I*E_I;
Z_I = (I-E_I)/sqrt(Var_I);
Z_alpha = norminv(1-alpha,0,1);
Moran I plot
function result = moranScatterPlot(X,WX)
%%X为列向量,W为权重矩阵
result = myls(X,WX);
X_lift = min(X) - 0.5*abs(min(X));
X_right = max(X) + 0.5*max(X);
xx = linspace(X_lift,X_right,1000);
yy = result(1)*xx + result(2);
plot(X,WX,'ro');
hold on
plot(xx,yy)
grid on
axis equal
end
matlab软件moran值,用matlab求moran并作图的具体步骤相关推荐
- Matlab软件有哪些版本,matlab软件功能简要介绍及版本描述
<matlab软件功能简要介绍及版本描述>由会员分享,可在线阅读,更多相关<matlab软件功能简要介绍及版本描述(2页珍藏版)>请在人人文库网上搜索. 1.MATLAB功能简 ...
- matlab显示灰度值的概率,求Matlab统计灰度值的代码
满意答案 ppwfhbrf 2013.07.06 采纳率:47% 等级:12 已帮助:11643人 理,直接显示直方图,这样它统计的就是灰度值的出现次数.例如: x=imread('D:\ebo ...
- 基于matlab测幅值,基于MATLAB的分析测试系统
基于MATLAB的分析测试系统 来源:华强电子网 作者:华仔 浏览:663 时间:2016-08-10 14:18 标签: 摘要: 摘 要:介绍基于matlab的分析测试系统的结构,详细阐述了系统的工 ...
- matlab软件与科学计算,MATLAB与科学计算教程
目 录 第1章 绪论11.1 MATLAB简介11.1.1 21世纪的科学计算语言11.1.2 MATLAB的发展历史31.1.3 MATLAB的应用和网上资源31.2 MATLAB ...
- 天津大学matlab软件许可,天津大学《MATLAB基础和应用》课程教学大纲.PDF
天津大学<MATLAB 基础与应用>课程教学大纲 MATLAB 基础与 课程编号: 2160283 课程名称: 应用 学 时: 48 学 分: 2.5 学时分配: 授课: 28 上机: 2 ...
- matlab软件有著作权,软件著作权matlab样例
软件著作权matlab样例是马景涛创作的一部小说,软件著作权matlab样例小说内容简介: ↙[软萌睿智女vs桀骜忠犬男,爽文互宠,双强双处,一生一世一双人]王忆锦,针灸世家传人,素手纤纤,弹指间风云 ...
- matlab 反向二值化,MATLAB:图像二值化、互补图(反运算)(im2bw,imcomplement函数)...
图像二值化.反运算过程涉及到im2bw,imcomplement函数,反运算可以这么理解:原本黑的区域变为白的区域,白的区域变为黑的区域. 实现过程如下: close all;%关闭当前所有图形窗口, ...
- 用matlab开发软件开发,Matlab软件应用与开发new
<Matlab软件应用与开发new>由会员分享,可在线阅读,更多相关<Matlab软件应用与开发new(13页珍藏版)>请在人人文库网上搜索. 1.Matlab软件应用与开发实 ...
- matlab使用_重磅!哈工大、哈工程无法使用 MATLAB 软件
近日,有消息称:哈工大.哈工程的师生们已经无法使用 MATLAB 软件了. 随后,有哈工大的学生表示收到了正版MATLAB 软件取消激活的通知,学校在与 MATLAB 软件的开发公司 MathWork ...
最新文章
- 数据库备份与还原处理
- Deformable ConvNets--Part2: Spatial Transfomer Networks(STN)
- 黑马程序员java笔记之二-----多线程
- Linux 浏览网址汇集
- Java07-day07【类和对象、对象内存图、成员变量与局部变量、private、this、封装、构造方法】
- AS工程下gradle各文件(gradle.propertie等)与gradle配置(缓存等)讲述
- mysql relay log时间_如何得到Slave应用relay-log的时间
- 4比较三个数的大小输出最大的值并从小到大排序输出
- 爬虫之Xpath的使用
- n个物品,Vi的价值,Ti有n个,求所有的可能的价值:
- vue+springboot实现登录验证码(前后端分离)
- 快压下载|快压软件官方下载
- CentOS 7下载及安装教程
- SolidWorks二次开发-工程图-预定义视图
- 禁止每次打开edge自动在开始菜单创建Microsoft Edge 的 Edge 栏
- kali linux怎么安装无线网卡驱动,Kali Linux 安装BCM43142网卡驱动
- 服务器文档链接电脑,服务器怎么链接电脑
- arp攻击----arpspoof
- DateADD日期Sql
- 单例模式的五种实现形式(懒汉式,饿汉式,双空判断,内部类,枚举)
热门文章
- 通过一个简单的例子学习Angular Injection Token工作原理
- 对具有外部依赖的Angular服务类(service class)进行单元测试的几种方式
- React应用里Invalid hook call错误消息的处理
- 使用nodejs对Marketing Cloud的contact主数据进行修改操作
- SAP CRM ABAP Product search authorization check - 产品搜索的权限检查机制
- SAP CRM WebClient UI on new focus工作原理
- SAP Gateway OData服务的语言决定逻辑 - language determination
- 简述大数据的数据管理方式_智能销售服务商荣登TOP100榜,大数据线上获客方式再得认可...
- 拼接 结果集_MetaQuast:评估宏基因组拼接
- fiddler抓包工具简介