spdep (version 0.2-14)

geary: Compute Geary's C

Description

A simple function to compute Geary's C, called by geary.test and geary.mc; $$C=\frac{(n-1)}{2\sum _{i=1}^n\sum _{j=1}^n{w}_{ij}}\frac{\sum _{i=1}^n\sum _{j=1}^n{w}_{ij}(x_i-x_j{)}^2}{\sum _{i=1}^n(x_i-\overline{x}{)}^2}$$ geary.intern is an internal function used to vary the similarity criterion.

Usage

geary(x, listw, n, n1, S0, zero.policy=FALSE)

Arguments

x

a numeric vector the same length as the neighbours list in listw

listw

a listw object created for example by nb2listw

n

number of zones

n1

n - 1

S0

global sum of weights

zero.policy

if TRUE assign zero to the lagged value of zones without neighbours, if FALSE assign NA

Value

• a list with
• C
  Geary's C
• K
  sample kurtosis of x

References

Cliff, A. D., Ord, J. K. 1981 Spatial processes, Pion, p. 17.

See Also

geary.test, geary.mc, sp.mantel.mc

Examples

data(oldcol)
col.W <- nb2listw(COL.nb, style="W")
str(geary(COL.OLD$CRIME, col.W, length(COL.nb), length(COL.nb)-1,
 Szero(col.W)))