R数据分析|地理集中指数

地理集中指数主要用于衡量地理要素的空间分布的集中程度，是地理学研究中一个常见指标。本文主要介绍如何利用R语言计算地理集中指数。

1、公式介绍

公式和数据主要参考冯亚芬等（2017）
在地理科学发表的《广东省传统村落空间分布特征及影响因素研究》 一文。演示数据来自文中表1数据。数据也可以在后台回复【20220527】获取。具体如下：

公式数据

2、指标计算

主要利用dplyr包中的mutate() 函数与summarise() 函数进行计算。

setwd("C:\\Users\\Acer\\Desktop\\R学习\\常用数据")
data <- readxl::read_xlsx("不均衡指数.xlsx")
library(tidyverse)
head(data)
# A tibble: 6 x 2
#  region value
#  <chr>  <dbl>
#1 梅州      40
#2 清远      13
#3 广州      12
#4 湛江      12
#5 肇庆       9
#6 惠州       7

# 计算
data %>% mutate(prop = value/sum(value), #实际占比
                mean_prop = mean(value)/sum(value) #均匀分布占比
                ) %>% 
  summarise(index = sqrt(sum(prop^2)) *100, #地理集中指数
            mean_index = sqrt(sum(mean_prop^2)) *100 #均匀分布地理集中指数
            ) %>% 
  as.data.frame()

#    index mean_index
#   38.0456   21.82179

计算结果显示，地理集中指数为38.0456，均匀分布时的地理集中指数为21.82179，与原文结果一致。

原文结果

3、其他

• 关于空间基尼系数的计算可参考R数据分析|空间基尼系数。
• 关于不均衡指数计算可参考R数据分析|不均衡指数。
• 更多内容可阅读公众号其他文章。

如有帮助请多多点赞哦！