基于R语言的地理探测器实现与问题研究

由于本人近期需要用到地理探测器做实验分析，以此记录并学习。话不多说，先上代码。

# 下载并导入包，用以读取excel文件
install.packages("readxl")
library(readxl)
# 读取excel文件
beibuwan <- read_excel("文件路径")
# 展示该文件
View(beibuwan)
# 下载并导入GD包
#install.packages("GD")
library(GD)
# 数据离散方法
discmethod <- c("equal","natural","quantile","geometric")
# 间断点范围
discitv <- c(3:10)
# beibuwandata<-edit(beibuwan)
# 连续型变量选择，这里C的含义是除去导入文件的哪些列
# 值得注意的是，R语言索引从1开始
continuous_variable <- colnames(beibuwan)[-c(1,2,3,4,5,6,7,8,9,10,11,12,13,17)]
continuous_variable
# 直接调用GD中的gdm，其将根据最大q，选取上述discmethod,discitv中最优值
testgdm <- gdm(Y ~ X13+X14+X15,
               continuous_variable = continuous_variable ,
               data = as.data.frame(beibuwandata),
               discmethod = discmethod, discitv = discitv)
testgdm
View(testgdm)

运算一定时间后，4种类型直接在控制台显示。也可以进行视图表达View（testgdm