已知有13个地区,它们的人均GDP为A=[26998;30840;34892;41099;45444;50467;58204;64491;66940;73856;81658;87475;19141;]
临界矩阵为
W=[0,1,0,0,1,0,1,0,1;1,0,0,1,1,0,0,1,0;0,0,0,0,0,0,1,0,0;0,1,0,0,0,1,0,0,0;1,1,0,0,0,0,1,1,0;0,0,0,1,0,0,0,0,0;1,0,1,0,1,0,0,1,1;0,1,0,0,1,0,1,0,0;1,0,0,0,0,0,1,0,0;1,0,0,1,0,1,0,0,0;0,0,0,0,0,0,0,0,0;0,0,1,0,0,0,0,0,0;0,0,1,0,0,0,1,1,0;]
然后我找到了Moran_I指数的函数
HELP morans_I
PURPOSE: calculate global Moran's I for an input grid (matrix) by calculating all
local Moran's I for a given moving windows size using a weight matrix.
-------------------------------------------------------------------
USAGE: M = moransI(grid, W, s);
where: [grid] is the matrix to analyse
[W] is the normalized weight matrix of the size the local Moran's
I will be calculated for (uneven sized!)
[s] is an optional flag to use zscores of input values for
calculation. Set to 'true' if zscores of local grid should be
calculated. Leave blank if not desired or input values are already
standardized.
-------------------------------------------------------------------------
OUTPUTS:
[M] matrix of all local Moran's I
-------------------------------------------------------------------
NOTES: Weight matrix needs to be 'moving window' style, not contiguity
matrix: Moran's I is calculated and weighted for neighbours to center cell.
Matrix needs to be normalized (weights sum to 1) and center cell weight
will be set to 0 if not already. Uses localmoran.m
-> Use nanmean(M(:)) to get the average global Moran's I.
See Anselin (1995, 'LISA.', Geogr. Analysis 27(2),p.93f) for details on
standardized variables in calculation of local Moran's I.
EXAMPLE: M = moransI(rand(20,20),ones(5,5),'true')
Felix Hebeler, Geography Dept.,de University Zurich, March 2006.
请问使用M=Morans_I(A,W,'ture')计算吗?为什么算出来都是NaN?
求大神帮助