ndvilist <- list(ndvi_30, ndvi_40, ndvi_50)
su <- c(30, 40, 50) ## sizes of spatial units
## "gdm" function
gdlist <- lapply(ndvilist, function(x){
gdm(NDVIchange ~ Climatezone + Mining, data = x)
})
sesu(gdlist, su) ## size effects of spatial units
Run the code above in your browser using DataLab