qs: Power of spatial determinant (PSD).

In IDSA: An Interactive Detector for Spatial Associations

Description

Power of spatial determinant (PSD).

Usage

qs(y, xh, location)

Arguments

y	A numeric vector of a response variable
xh	A character variable, a data frame or a matrix of explanatory variables
location	A matrix of spatial locations

Value

A power of spatial determinant (PSD) value.

Examples

# an explanatory variable
library(GD)
data.disc <- disc(sim$xa, 4, "quantile")
xh <- cut(sim$xa, data.disc$itv, include.lowest = TRUE)
qs(sim$y, xh, location = sim[, c("lo","la")])
# multiple explanatory variables
data <- sim[,4:6]
data.disc <- apply(data, 2, FUN = function(x) disc(x, 4, "quantile"))
xh <- do.call(cbind, lapply(1:ncol(data), function(x)
  data.frame(cut(data[, x], data.disc[[x]]$itv, include.lowest = TRUE))))
names(xh) <- names(data)
qs(sim$y, xh, location = sim[, c("lo","la")])

IDSA documentation built on April 18, 2021, 5:07 p.m.