地理探测器的R语言实现

最新推荐文章于 2024-10-09 18:35:41 发布

四九马赫

最新推荐文章于 2024-10-09 18:35:41 发布

阅读量1.2k

点赞数 16

文章标签： r语言

本文链接：https://blog.youkuaiyun.com/qq_64891272/article/details/136546826

版权

文章讲述了作者在使用R语言的GD包进行数据离散化分析时遇到错误，但最终仍能得到结果的过程，重点在于问题解决和分析过程的描述。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

记录一下自己的学习记录，有一部分报错了，但是还是出来结果了，初学还是有点不懂。

library(GD)

# 设置工作目录
setwd("G:/CITY_TONGJI3/西宁")

# 读取数据
data <- read.csv("2010-2020西宁.CSV", as.is = TRUE)

# 定义离散化方法和区间数
discmethod <- c("equal", "natural", "quantile", "geometric", "sd")
discitv <- c(3:10)

# 存储结果
results <- list()

# 使用 gdm() 函数进行分析
for (i in discitv) {
result <- gdm(Y ~ X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 + X11 + X12 + X13 + X14 + X15 + X16,
continuous_variable = c("X1", "X2", "X3", "X4", "X5", "X6", "X7", "X8", "X9", "X10", "X11", "X12", "X13", "X14", "X15", "X16"),
data = data,
discmethod = discmethod,
discitv = i)

results[[as.character(i)]] <- result
}

# 打印结果
for (i in discitv) {
print(paste("Result for discitv =", i))
print(results[[as.character(i)]])
}

# 绘制图表
for (i in discitv) {
plot(results[[as.character(i)]])
}

有这个报错，不知道为什么，但是还是跑出结果了