Q.map.test: Compute the QE and QI tests of Equivalence and Independence...

In rominsal/spqdata: Testing for spatial structure of qualitative data in cross section

Description

This function compute the QE and QI tests for maps comparison based on symbolic entropy.

Usage

Q.map.test(formula = formula, data = data, coor = NULL, m = m, r = 1, type = "combinations", control = list())

Arguments

formula	a symbolic description of the two factors.
data	(optional) a data frame or a sf object containing the variables to testing for.
coor	(optional) a 2xN vector with coordinates.
m	length of m-surrounding.
r	maximum overlapping between any two m-surroundings (default = 1).
type	Type of symbols: "permutations" or "combinations". Default "combinations"
control	Optional argument. See Control Argument section.

Details

Aquí Antonio escribe una linda historia ...
Two alternative sets of arguments can be included in this function to compute the spatial runs test:

Alternative 1	A factor (xf) and a list of neighborhood (listw) of the class knn.
Alternative 2	A sf object (data) and formula to specify the factor. A list of neighborhood (listw)

Value

A list with two objects of the class htest. The first one is the QE test of Equivalence between maps and the second one is the QI test of independence between maps.

method	a character string giving description of the method.
data.name	a character string giving the name(s) of the data.
statistic	the value of the statistic QE or/and QI
alternative	a character string describing the alternative hypothesis.
p.value	p-value for QE or QI
gl	free degree of the statistic for QE or QI
mh	m-surrounding of th map.
Tm	number of maps (ONLY 2).
R	number of symbolized observations.
nsk	a matrix Tm x symbols with the frequency of the number of symbols of each map

Control arguments

Several parameters to construct the m-surrounding

dtmaxabs

Delete degenerate surrounding based on the absolute distance between observations.

dtmaxpc

A value between 0 and 1. Delete degenerate surrounding based on the distance. Delete m-surrounding when the maximum distance between observation is upper than k percentage of maximum distance between anywhere observation.

dtmaxknn

A integer value 'k'. Delete degenerate surrounding based on the near neighborhood criteria. Delete m-surrounding is a element of the m-surrounding is not include in the set of k near neighborhood of the first element

seedinit

seed to select the initial element to star the algorithm to compute the m-surroundings.

Author(s)

Fernando López	[e-mail]
Román Mínguez	[e-mail]
Antonio	[e-mail]
Manolo	[e-mail]

@references

• Ruiz M, López FA and A Páez (2011). Comparison of Thematic Maps Using Symbolic Entropy. International Journal of Geographical Information Science, 26, 413-439.

• Ruiz, M., López, F., and Páez, A. (2010). Testing for spatial association of qualitative data using symbolic dynamics. Journal of Geographical Systems, 12(3), 281-309.0.

See Also

dgp.spq, m.surround

Examples

# Case 1:
rm(list = ls())
N <- 200
cx <- runif(N)
cy <- runif(N)
x <- cbind(cx,cy)
listw <- spdep::nb2listw(spdep::knn2nb(knearneigh(cbind(cx,cy), k = 4)))
p <- c(1/6,3/6,2/6)
rho = 0.5
QY1 <- dgp.spq(p = p, listw = listw, rho = rho)
rho = 0.8
QY2 <- dgp.spq(p = p, listw = listw, rho = rho)
dt = data.frame(QY1,QY2)
m = 3
r = 1
formula <- ~ QY1 + QY2
control <- list(dtmaxknn = 10)
qmap <- Q.map.test(formula = formula, data = dt, coor = x, m = m, r = r, type ="combinations", control = control)
print(qmap)
plot(qmap)
plot(qmap, ci=.6)
plot(qmap[[1]]$mh)
summary(qmap[[1]]$mh)

control <- list(dtmaxknn = 20)
qmap <- Q.map.test(formula = formula, data = dt, coor = x, m = m, r = r, type ="permutations", control = control)
print(qmap)
plot(qmap)
plot(qmap[[1]]$mh)
qmap <- Q.map.test(formula = formula, data = dt, coor = x, m = m, r = r, type ="combinations")
print(qmap)
plot(qmap)
control <- list(dtmaxknn = 10)
qmap <- Q.map.test(formula = formula, data = dt, coor = x, m = m, r = r, type ="combinations", control = control)
print(qmap)
plot(qmap)

# Case 2:
rm(list = ls())
data(Spain)
m = 3
r = 1
formula <- ~ Coast + Older65
qmap <- Q.map.test(formula = formula, data = spain.sf, m = m, r = r, type ="combinations")
print(qmap)
plot(qmap)
plot(qmap[[1]]$mh)

control <- list(dtmaxknn = 6)
qmap <- Q.map.test(formula = formula, data = spain.sf, m = m, r = r, type ="combinations", control = control)
print(qmap)
plot(qmap[[1]]$mh)
summary(qmap[[1]]$mh)

# Case 3:
rm(list = ls())
data(Newark)
m = 3
r = 1
formula <- ~ IRISH + GERMAN
control <- list(dtmaxknn = 6)
qmap <- Q.map.test(formula = formula, data = Newark, coor = cbind(Newark$X,Newark$Y), m = m, r = r, type ="combinations", control = control)
print(qmap)
plot(qmap)
plot(qmap[[1]]$mh)

rominsal/spqdata documentation built on Jan. 24, 2021, 11:11 a.m.