使用tidycensus的get_decennial()获取2010年街区组数据报错

看起来你遇到的问题是tidycensus在处理2010年block group层级的人口普查数据时的兼容性或参数配置问题，结合你的代码和报错信息，我整理了几个针对性的解决方案：

1. 显式指定普查汇总文件（sumfile参数） #

2010年人口普查分为SF1（全样本数据）和SF3（抽样数据）两个核心汇总文件。虽然get_decennial通常会自动匹配地理层级对应的文件，但针对block group这类细分层级，显式指定sumfile = "sf1"可能会解决"变量不可用"的报错：

bg_10 <- get_decennial(
  geography = "block group",
  variables = "P001001",
  year = 2010,
  state = "WY",
  county = "Teton",
  geometry = FALSE,
  sumfile = "sf1"
)

2. 更新tidycensus及依赖包到最新版本 #

你遇到的gather_函数报错大概率是因为tidycensus与dplyr等依赖包的版本不兼容。旧版本的tidycensus在处理2010年block group数据时存在已知bug，更新到最新版本可以直接修复这类问题：

# 更新tidycensus包
install.packages("tidycensus")
# 同步更新依赖的dplyr包（避免版本冲突）
install.packages("dplyr")

3. 验证变量在目标地理层级的可用性 #

如果上述方法仍无效，可以先确认P001001变量是否确实存在于2010年block group层级的SF1数据中：

# 加载2010年SF1的变量列表
vars_2010 <- load_variables(year = 2010, dataset = "sf1")
# 搜索目标变量
dplyr::filter(vars_2010, name == "P001001")

根据人口普查官方文档，P001001是2010年SF1中block group层级的总人口变量，理论上应该存在；如果搜索结果为空，你需要使用对应的替代变量名。

4. 先获取州级数据再筛选目标县 #

有时候指定county参数可能触发边界处理的小bug，你可以先获取整个怀俄明州的2010年block group数据，再筛选Teton县：

# 获取怀俄明州所有block group数据
bg_10_wy <- get_decennial(
  geography = "block group",
  variables = "P001001",
  year = 2010,
  state = "WY",
  geometry = FALSE,
  sumfile = "sf1"
)

# 筛选Teton县（FIPS代码为021）的数据
bg_10 <- dplyr::filter(bg_10_wy, COUNTY == "021")

优先尝试更新包和指定sumfile参数，这两个是解决该问题最常见的途径。

内容的提问来源于stack exchange，提问作者Pukitta