US Small-area Life Expectancy Estimates Project (USALEEP)
Estimates of life expectancy at birth e (0) by census tract 2010-2015.
National Center for Health Statistics. 2018. "US Small-Area Life Expectancy Estimates Project (USALEEP): Life Expectancy Estimates File, 2010-2015." Updated September. https://www.cdc.gov/nchs/nvss/usaleep/usaleep.html.
Geography
Polygons are from the USCB TIGER/Line shapefiles.
- GEOIDFQ: Fully-qualified GEOID code composed of a prefix and FIPS code that can be used to unambiguously join with US Census Bureau data from data.census.gov.
- ST (STUSPS): USPS two-letter state code (states, counties, tracts, and ZCTA only)
- Latitude: NAD 1983 latitude of the centroid calculated from the area polygons with st_centroid()
- Longitude: NAD 1983 longitude of the centroid calculated from the area polygons with st_centroid()
- Square_Miles: Land area in square miles calculated by dividing ALAND + AWATER (square meters) by 2,589,988
Data Fields
- Life_Expectancy: Life expectancy at birth e(0)
- Standard_Error: Standard error of the estimate se(e(0))
- Flags Abridged life table flag:
- 1 = Observed age-specific death rates for all age groups
- 2 = Predicted age-specific death rates for all age groups
- 3 = Combination of observed and predicted age-specific death rates
R Processing Code
library(sf)
polygons = st_read("https://michaelminn.net/tutorials/data/2015-2019-acs-tracts.geojson")
polygons = polygons[,c("FactFinder.GEOID", "ST", "Square.Miles", "geometry")]
polygons$Longitude = round(st_coordinates(st_centroid(polygons))[,"X"], 3)
polygons$Latitude = round(st_coordinates(st_centroid(polygons))[,"Y"], 3)
names(polygons) = c("GEOIDFQ", "ST", "Square_Miles", "geometry", "Latitude", "Longitude")
polygons = polygons[,c("GEOIDFQ", "ST", "Square_Miles", "Latitude", "Longitude", "geometry")]
object.size(polygons)
data = read.csv("2015-usaleep.csv")
data$GEOIDFQ = sprintf("1400000US%011.0f", data[,"Tract.ID"])
data = data[,c("GEOIDFQ", "e.0.", "se.e.0..", "Abridged.life.table.flag")]
names(data) = c("GEOIDFQ", "Life_Expectancy", "Standard_Error", "Flags")
polygons = merge(polygons, data, all.x=T)
st_write(polygons, "2015-usaleep.geojson")
# plot(polygons["Life_Expectancy"], border=NA)