# NOT RUN {
# ndvilist <- list(ndvi_30, ndvi_40, ndvi_50)
# su <- c(30,40,50) ## sizes of spatial units
# gdlist <- list() ## list of all geographical detectors results
### set optional parameters of optimal discretization
### optional methods: equal, natural, quantile, geometric, sd and manual
# discmethod <- c("equal","natural","quantile")
# discitv <- c(4:6)
### "gdm" function
# for (i in 1:length(su)){
# ndvidata <- ndvilist[[i]]
# gdlist[[i]] <- gdm(y = ndvidata[,1], xcategorical = ndvidata[,2:3],
# xcontinuous = ndvidata[,c(4:7)],
# discmethod = discmethod, discitv = discitv)
# }
# sesu(gdlist, su) ## size effects of spatial units
# }
Run the code above in your browser using DataLab