用R做空间自相关分析-Moran'I遇到问题求大神解答！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

最近在学习用R做空间统计分析，按照[基于R语言的空间统计分析研究与应用-杨中庆]的论文来做的。

然后前面要下载的包我都下了：
  > search()
[1] ".GlobalEnv"        "package:spdep"     "package:Matrix"   
[4] "package:sp"        "package:stats"     "package:graphics"
[7] "package:grDevices" "package:utils"     "package:datasets"
[10] "package:methods"   "Autoloads"         "package:base"   

然后进行到这一步：看多边形化后的哥伦布市地图和各区的ID值时，总是出现这个错误：
> plot(polys.OLD)
错误于plot(polys.OLD) :
  在为'plot'函数选择方法时评估'x'参数出了错: 错误: 找不到对象'polys.OLD'   
请问这是什么意思，要怎么解决啊？第一次发帖，请各位大神帮帮忙！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：空间自相关分析 空间自相关 moran 相关分析 自相关 methods package search Matrix 空间

求不沉！自己顶！

R中没有怎么处理过，Matlab有现成的：

function result = moran(y,x,W);
% PURPOSE: computes Moran's I-statistic for spatial correlation
%          in the residuals of a regression model
% ---------------------------------------------------
%  USAGE: result = moran(y,x,W)
%  where: y = dependent variable vector
%         x = independent variables matrix
%         W = contiguity matrix (standardized or unstandardized)
% ---------------------------------------------------
%  RETURNS: a  structure variable
%         result.meth   = 'moran';
%         result.morani = e'*W*e/e'*e (I-statistic)
%         result.istat  = [i - E(i)]/std(i), standardized version
%         result.imean  = E(i),   expectation
%         result.ivar   = var(i), variance
%         result.prob   = std normal marginal probability
%         result.nobs   = # of observations
%         result.nvar   = # of variables in x-matrix
% ---------------------------------------------------
% NOTES: istat > 1.96, => small prob,
%                     => reject HO: of no spatial correlation
% This function standardizes the weight matrix and produces
% results based on a standardized weights matrix
% ---------------------------------------------------
% See also: lmerrs, walds, lratios
% ---------------------------------------------------

% written by:
% James P. LeSage, Dept of Economics
% University of Toledo
% 2801 W. Bancroft St,
% Toledo, OH 43606
% jlesage@spatial-econometrics.com


if nargin ~= 3
error('Wrong # of arguments to moran');
end;

[n k] = size(x);

% standardize the weight matrix
W = normw(W);

% do ols to get residuals
b = inv(x'*x)*x'*y;
e = y - x*b;
epe = e'*e;
mi = (e'*W*e)/epe;
M = eye(n) - x*(inv(x'*x))*x';
tmw = trace(M*W);

meani = tmw/(n-k);
vari =  trace((M*W)*(M*W')) + trace((M*W)*(M*W)) + tmw*tmw;
vari = vari/((n-k)*(n-k+2));
vari = vari - meani*meani;
mis = (mi-meani)/sqrt(vari);
prob = norm_prb(mis);

result.meth = 'moran';
result.nobs = n;
result.nvar = k;
result.morani = mi;
result.istat  = mis;
result.imean  = meani;
result.ivar   = vari;
result.prob   = prob;

wwwujl 发表于 2014-12-4 22:31
R中没有怎么处理过，Matlab有现成的：

function result = moran(y,x,W);

谢谢你啊，可是我老师要我必须用R  >_<

某个包没有load吧

qoiqpwqr 发表于 2014-12-5 09:29
某个包没有load吧

我也遇到同样的问题，请问是要下那个包呢？

易老师喵了个咪 发表于 2015-11-12 19:14
我也遇到同样的问题，请问是要下那个包呢？

楼上说的是可能没有LOAD某包，就是说你所主要用的那个包所要求必须提前先加载好的包，看包的DOC，其实一般用用还是R STUDIO方便

tzj110 发表于 2015-11-13 18:34
楼上说的是可能没有LOAD某包，就是说你所主要用的那个包所要求必须提前先加载好的包，看包的DOC，其实一般 ...

好的，谢谢！

wwwujl 发表于 2014-12-4 22:31
R中没有怎么处理过，Matlab有现成的：

function result = moran(y,x,W);

在MATLAB里边用这个程序算出来的不是因变量的MORAN值？

wwwujl 发表于 2014-12-4 22:31
R中没有怎么处理过，Matlab有现成的：

function result = moran(y,x,W);

这个求出来的是因变量的moran值吗？