R语言地理探测器--因子探测器实现

上一篇gdm包实现整个地理探测器分析流程，但是整个包直接计算不同数据会产生很多问题，本人现阶段只需要因子探测，因此单独拿出来因子探测实现脚本，上代码！

install.packages("readxl")
library(readxl)

# 读取excel文件
beibuwan <- read_excel("F:/beibuwan.xlsx")

# 查看数据前三行
head(beibuwan)[1:3,]
# 数据转化
beibuwan$X1=as.numeric(beibuwan$X1)
beibuwan$X2=as.numeric(beibuwan$X2)
beibuwan$X3=as.numeric(beibuwan$X3)
beibuwan$X4=as.numeric(beibuwan$X4)
beibuwan$X5=as.numeric(beibuwan$X5)
beibuwan$X6=as.numeric(beibuwan$X6)
beibuwan$X7=as.numeric(beibuwan$X7)
beibuwan$X8=as.numeric(beibuwan$X8)
beibuwan$X9=as.numeric(beibuwan$X9)
beibuwan$X10=as.numeric(beibuwan$X10)
beibuwan$X11=as.numeric(beibuwan$X11)
beibuwan$X12=as.numeric(beibuwan$X12)
beibuwan$X13=as.numeric(beibuwan$X13)
beibuwan$X14=as.numeric(beibuwan$X14)
beibuwan$X15=as.numeric(beibuwan$X15)
beibuwan$X16=as.numeric(beibuwan$X16)
# 下载并导入GD包
# install.packages("GD")
library(GD)
# 数