導入要分析的數據，并用切片器選擇因變量和自變量，因變量在左側，并刪除無關數據，包括索引。

testdata <- read.csv(csvfile, header = TRUE, sep = ",") testdata=testdata[3:9]

輸入可供選擇的分類方法

#3 approaches to calculate disgress.discmethod <- c("equal","natural","quantile")

輸入可供選擇的分類基本數量

#classify 4~6 discitv <- c(4:6)

篩選連續的自變量

#continuous variablescontinuous_variable <- colnames(testdata)[-1]

運行gdm函數?

#gdm functiontestgdm <- gdm(WSI ~ .,continuous_variable = continuous_variable ,data = testdata,discmethod = discmethod, discitv = discitv)

分析結果

?完整代碼

Sys.setlocale(category = 'LC_CTYPE', locale = 'C')#This is for using in pycharm. library(GD)filepath <- "E:/BaiduNetdiskDownload/Driving/Yl_naturalBreaks/" setwd(filepath)#Change the workspace temp=list.files(path=filepath,pattern="*.csv")#list csvfilenames in workspace for( csvfile in temp){print(csvfile)testdata <- read.csv(csvfile, header = TRUE, sep = ",")#load datatestdata <- read.csv("2003.csv", header = TRUE, sep = ",")testdata=testdata[3:9]#5 approach to calculate disgress.discmethod <- c("equal","natural","quantile")#select the best one among automatically 3 approaches of discrizetion #classify 4~6 automaticallydiscitv <- c(4:6)#continuous independent variables' name listcontinuous_variable <- colnames(testdata)[-1] #[-1] is for unselect dependent variable #gdm functiontestgdm <- gdm(WSI ~ .,continuous_variable = continuous_variable ,data = testdata,discmethod = discmethod, discitv = discitv)#dependent variable is WSI}print(testgdm)}

?

總結

以上是生活随笔為你收集整理的基于R语言地理探测器包（GD）空间异质性与驱动力分析的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。