# NOT RUN {
library(GD)
data <- sim[, 4:6]
data.disc <- apply(data, 2, FUN = function(x) disc(x, 4, "quantile"))
layers <- do.call(cbind, lapply(1:ncol(data), function(x)
data.frame(cut(data[, x], data.disc[[x]]$itv, include.lowest = TRUE))))
names(layers) <- names(data)
fo <- fuzzyoverlay(y = sim[,1], layers = layers, method = "fuzzyAND")
qo <- qsoverlay(x = data, xoverlay = fo$fuzzylayer,
location = sim[, c("lo","la")])
# }
Run the code above in your browser using DataLab