導入需要的包：

> .libPaths('D:/Rpackages') > library('GD')

載入數據

data<-read.table(file="D:/GWRR/統計總表-終版.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)

運行

> ndvigdm <- gdm(income2011 ~ .,continuous_variable = name2011[-1],data = factors2011,discmethod = discmethod, discitv = discitv)

在此研究中，所有的factors都是連續值，需要定義出來。否則GD包會將每個獨立數值理解為要給類別，而不是進行分組計算。（速度會賊慢）

可視化

由于結果通常會有很大的圖，經常會報錯：

Error in plot.new() : figure margins too large

解決方法：
1）將Plots界面盡可能拖到最大
2）使用

par(mar=c(1,1,1,1))

總結

以上是生活随笔為你收集整理的使用R包GD实现地理探测器算法的全部內容，希望文章能夠幫你解決所遇到的問題。

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