回复：A simple function to compute Moran's I, callded by moran.test or moran.mc
        I=(n sum_i sum_j w_ij (x_i-xbar)(x_j-xbar))/(S0 sum_i (x_i-xbar)^2)
        在R中可以用moran (x,listw,n,S0,zero.policy=NULL,NAOK=FALSE)