Hi, George. Tacking this on to Matt H's reply, `tidycensus::load_variables()` will give you only the combined labels, which can be difficult to sift through to find what you're looking for. For example, S0501_C01_088 has the label "Estimate!!Total!!EARNINGS IN THE PAST 12 MONTHS (IN 2019 INFLATION-ADJUSTED DOLLARS) FOR FULL-TIME, YEAR-ROUND WORKERS!!Population 16 years and over with earnings!!Median earnings (dollars) for full-time, year-round workers:!!Female".
I have an easier time seeing the hierarchy in the tables when I split up the label into separate parts:
load_variables(year = 2019, dataset = "acs5/subject") %>%
separate(label, into = paste0("label", 1:9), sep = "!!", fill = "right", remove = FALSE)
That will break apart the label components and give you new variables: label1 would be "Estimate", label2 would be "Total", ..., and label6 would be "Female". Then it's easier to filter on different dimensions of each table.
Something else that may be helpful: A given variable name in one year may represent something completely different in another year (unlike in the detailed tables, which are often replaced by tables with a different name when categories change substantially). And even when the variable names are the same from year to year, the format of the labels sometimes does. I do like the subject tables for getting more accurate margins of error for certain fields (rather than using the formulas for calculating them ourselves), but it can be hard to grab them over time.
Good luck!
--Matt