guery: Mutivariable and univariable Geary's C statistic

In amvallone/estdaR: Exploratory Spatio-Temporal Data Analysis in R

Mutivariable and univariable Geary's C statistic

Description

Compute the univariable Anselin (1995) and multivarible Anselin (2017) Geary's C statistic

Usage

geary(x, W, nsim = NULL, type = "uni", nbcom = length(W$neighbours), ...)

Arguments

x	a ventor, matrix or data frame containing the variables and spatial units
W	a listw object
nsim	number of random permutations used to the compute the pseudo p value. By default it is NULL
type	a character indicating the Geary's C to be compute. If it is set as "uni" the standart Geary's C statistic is calulated. Set type as "multi" to cumpute the multivariable Geary's C.
nbcom	number of comparisons use in the Bonferroni multiple comparisons correction. By default it is set aas the number of spatial unit to use.
...	other argument to quad function. See quad for more information.

Details

later...

Value

a data frame (more explanation in Details later)

References

Anselin, L. (1995). Local Indicators of Spatial Association—LISA. Geographical Analysis, 27(2), 93–115. https://doi.org/10.1111/j.1538-4632.1995.tb00338.x Anselin, L (2017). A Local Indicator of Multivariate Spatial Association: Extending Geary’s c. Center for Spatial Data Science working papers. University of Chicago.

Examples

data(Guerry)
w1queen<-nb2listw(poly2nb(Guerry))
pcc <- cbind(Guerry$pc1,Guerry$pc2)
## Not run: e1 <- geary(Guerry$pc1,w1queen,nsim=999)
e2 <- geary(pcc,w1queen,nsim=999)
e3 <- geary(pcc,w1queen,type="multi",nsim=999)
## End(Not run)

amvallone/estdaR documentation built on March 30, 2024, 9:38 p.m.