Last chance! 50% off unlimited learning
Sale ends in
Public use microdata areas (PUMAs) are decennial census areas that have been defined for the tabulation and dissemination of Public Use Microdata Sample (PUMS) data, American Community Survey (ACS) data, and ACS period estimates. For the 2010 Census, the State Data Centers (SDCs) in each state, the District of Columbia, and the Commonwealth of Puerto Rico were given the opportunity to delineate PUMAs within their state or statistically equivalent entity. All PUMAs must nest within states and have a minimum population threshold of 100,000 persons. 2010 PUMAs were built on census tracts and cover the entirety of the United States, Puerto Rico, Guam, and the U.S. Virgin Islands. Because they do not meet the minimum population requirement, the Commonwealth of the Northern Mariana Islands and American Samoa do not contain any 2010 PUMAs.
pumas(state = NULL, cb = FALSE, year = NULL, ...)The two-digit FIPS code (string) of the state you want. Can also
be state name or state abbreviation. When NULL and combined with
cb = TRUE, a national dataset of PUMAs will be returned when
year = 2019 only.
If cb is set to TRUE, download a generalized (1:500k) states file. Defaults to FALSE (the most detailed TIGER/Line file)
the data year (defaults to 2020).
arguments to be passed to the underlying `load_tiger` function, which is not exported.
Options include class, which can be set to "sf" (the default) or "sp" to
request sf or sp class objects, and refresh, which specifies whether or
not to re-download shapefiles (defaults to FALSE).
https://www.census.gov/programs-surveys/geography/guidance/geo-areas/pumas.html
Other general area functions:
block_groups(),
blocks(),
counties(),
county_subdivisions(),
places(),
school_districts(),
states(),
tracts(),
zctas()
# NOT RUN {
library(tigris)
us_states <- unique(fips_codes$state)[1:51]
continental_states <- us_states[!us_states %in% c("AK", "HI")]
pumas_list <- lapply(continental_states, function(x) {
pumas(state = x, cb = TRUE, year = 2017)
})
us_pumas <- rbind_tigris(pumas_list)
plot(us_pumas$geometry)
# }
Run the code above in your browser using DataLab