plot.sprunstest: Plot the empirical distribution of runs

In spqdep: Testing for Spatial Independence of Qualitative Data in Cross Section

Plot the empirical distribution of runs

Description

Plot the empirical distribution of runs

Usage

## S3 method for class 'sprunstest'
plot(x, ...)

Arguments

x	A object of class sprunstest.
...	further arguments passed to or from other methods.

Details

Plot the histogram with the empirical distribution of the runs

Value

No return value, called for side effects

Author(s)

Fernando López	[e-mail]
Román Mínguez	[e-mail]
Antonio Páez	[e-mail]
Manuel Ruiz	[e-mail]

See Also

sp.runs.test.

Examples

# Example 1: Fastfood example. sf (points)
data("FastFood.sf")
x <- sf::st_coordinates(sf::st_centroid(FastFood.sf))
listw <- spdep::knearneigh(x, k = 2)
formula <- ~ Type
srq <- sp.runs.test(formula = formula, data = FastFood.sf, listw = listw, nsim = 299)
plot(srq)

# Example 2: Spain example (poligons with 0 neinghbourhood)
data("provinces_spain")
sf::sf_use_s2(FALSE)
listw <- spdep::poly2nb(as(provinces_spain,"Spatial"), queen = FALSE)
provinces_spain$Older <- cut(provinces_spain$Older, breaks = c(-Inf,19,22.5,Inf))
levels(provinces_spain$Older) = c("low","middle","high")
formula <- ~ Older
srq <- sp.runs.test(formula = formula, data = provinces_spain, listw = listw, nsim = 299)
plot(srq)
provinces_spain$Male2Female <- factor(provinces_spain$Male2Female > 100)
levels(provinces_spain$Male2Female) = c("men","woman")
formula <- ~ Male2Female
srq <- sp.runs.test(formula = formula, data = provinces_spain, listw = listw, nsim = 299)
plot(srq)

spqdep documentation built on March 28, 2022, 5:06 p.m.