我正在绘制阿拉斯加州各县的肥胖值,但有些县是灰色的,不代表该县的值。我无法找出灰色县的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")
1条答案
按热度按时间6yt4nkrj1#
在找到缺失县的名称(Chugach人口普查区[02063]和Copper River人口普查区[02066])后,我发现这些县创建于2019年。Valdez-Cordova人口普查区县在2019年被拆分为这两个新县,我的数据创建于2019年之前,所以它仍然有Valdez-Cordova人口普查区县。我最终在数据中创建了Chugach人口普查区和Copper River人口普查区,与Valdez-Cordova人口普查区具有相同的值。