基于R语言GD包的Risk Map制作（批量生成）

短短的代码…踩了多少坑
过去完全没有接触过R语言，很多本该知道的基础知识都不清楚，遇到了盲点就要花很长的时间解决。

一个完成的流程：

1)载入GD包

.libPaths('D:/Rpackages')


library('GD')

2)导入数据，按照年份对所有变量分组

data<-read.table(file="D:/GWRR/final.csv",header=TRUE,sep=',')

factors2011<-data[,c(4,7,9:39)]
factors2016<-data[,c(5,8,9:39)]


discmethod <- c("equal","natural","quantile")
discitv <- c(4:6)


name2011<-names(factors2011)
name2016<-names(factors2016)

3. 对2011年和2016年的数据，分别运行GD包，生成Riskmap

ndvigdm <- gdm(income2011 ~ .,continuous_variable = name2011[-1],
                data = factors2011,
                discmethod = discmethod, discitv = discitv) 
### Risk values of each units

name2011<-name2011[-1]

for(i in 1:length(name2011)){
   
  
  b <- ndvigdm$Risk.mean[[name2011