R语言 usmap未显示阿拉斯加州某些县的数据

mi7gmzs6 于 2023-02-14 发布在其他

关注(0)|答案(1)|浏览(131)

我正在绘制阿拉斯加州各县的肥胖值，但有些县是灰色的，不代表该县的值。我无法找出灰色县的FIPS代码，也无法找出它们为什么是灰色的。我假设我为它们提供的FIPS代码是不正确的，但即使是内部向上Map人口数据也将它们显示为灰色。

# read in data
chd_map = read.csv(file = 
'https://www.countyhealthrankings.org/sites/default/files/media/document/analytic_data2022.csv', 
header = TRUE, sep=",")

# pull out FIPS codes and obesity values
chd_map = chd_map[-c(1,2),c(3,75)]

# make obesity values numeric
chd_map$Adult.obesity.raw.value = as.numeric(chd_map$Adult.obesity.raw.value) * 100

# rename fips column
colnames(chd_map)[1] = 'fips'

# plot data
plot_usmap(
  data = chd_map,
  size = 0.1,
  regions = "counties",
  include = c('AK'),
  values = "Adult.obesity.raw.value",
  color = "black") +
  scale_fill_continuous(type = "viridis", name = "Percentage of Adults with Obesity (2022)", label = scales::comma
  ) +
  labs(title = "Figure X. Distribution of Obesity Prevalence by County") + 
  theme(panel.background = element_rect(color = "white", fill = "white"), legend.position = "bottom")

r

来源：https://stackoverflow.com/questions/75404058/usmap-not-showing-data-for-some-alaska-counties

1条答案

按热度按时间

6yt4nkrj1#

在找到缺失县的名称（Chugach人口普查区[02063]和Copper River人口普查区[02066]）后，我发现这些县创建于2019年。Valdez-Cordova人口普查区县在2019年被拆分为这两个新县，我的数据创建于2019年之前，所以它仍然有Valdez-Cordova人口普查区县。我最终在数据中创建了Chugach人口普查区和Copper River人口普查区，与Valdez-Cordova人口普查区具有相同的值。

valdez_obesity = chd_map$Adult.obesity.raw.value[chd_map$fips == '2261']

chugach = data.frame(2063,valdez_obesity)
names(chugach) = c('fips', 'Adult.obesity.raw.value')

copper = data.frame(2066,valdez_obesity)
names(copper) = c('fips', 'Adult.obesity.raw.value')

obesity_map_data = rbind(obesity_map_data, chugach, copper)

obesity_map_data = obesity_map_data[-which(obesity_map_data$fips == 2261),]

plot_usmap(
  data = chd_map,
  size = 0.1,
  regions = "counties",
  include = c('AK'),
  values = "Adult.obesity.raw.value",
  color = "black") +
  scale_fill_continuous(type = "viridis", name = "Percentage of Adults with Obesity (2022)", label = scales::comma
  ) +
  labs(title = "Figure X. Distribution of Obesity Prevalence by County") + 
  theme(panel.background = element_rect(color = "white", fill = "white"), legend.position = "bottom")

赞(0）回复(0）举报 2023-02-14

我来回答

R语言 usmap未显示阿拉斯加州某些县的数据

1条答案

相关问题

热门标签

最新问答