library(totalcensus)
# read_acs1year_filesegment ===================================================
aaa <- totalcensus:::read_acs1year_filesegment_(2005, "RI", "0003")
# read_acs1year_1_file_tablecontents ==========================================
aaa <- totalcensus:::read_acs1year_1_file_tablecontents_(
2005, "RI", "0003", c("B02005_042", "C02003_001"), "m"
)
# read_acs1year_tablecontents_ ==================================================
aaa <- totalcensus:::read_acs1year_tablecontents_(
2005, "RI", c("B02005_042", "C02003_001", "B01001_002", "B01001_026")
)
# read_acs1year_geoheader_file_ ================================================
aaa <- totalcensus:::read_acs1year_geoheader_file_(2005, "RI")
# read_acs1year_geo_() =========================================================
aaa <- totalcensus:::read_acs1year_geo_(
year = 2005,
state = "RI"
)
cols <- c("GEOID", "STUSAB", "NAME", "LOGRECNO", "SUMLEV", "GEOCOMP")
stopifnot(names(aaa) == cols)
aaa = totalcensus:::read_acs1year_geo_(
2010, "RI", geo_headers = c("STUSAB", "STATE", "CBSA")
)
stopifnot(dim(aaa) == c(38, 8))
aaa <- totalcensus:::read_acs1year_geo_(
year = 2005,
state = "RI",
geo_headers = dict_acs_geoheader_2005_1year$reference
)
cols <- c("GEOID", "STUSAB", "NAME", "LOGRECNO", "SUMLEV", "GEOCOMP",
"FILEID", "US", "REGION", "DIVISION", "STATE", "COUNTY", "MCD",
"PL", "AINDN", "CBSA", "METDIV", "CSA", "CNECTA", "NECTA", "NECTADIV",
"UA", "CD2000", "BST", "PUMA5", "SDELM", "SDSEC", "SDUNI", "UR",
"MEMI", "MEMIR", "PCI", "CD1990", "COUSUB", "PLACE")
stopifnot(names(aaa) == cols)
aaa <- totalcensus:::read_acs1year_geo_(
year = 2017,
state = "DC",
geo_headers = dict_acs_geoheader_2011_now$reference
)
stopifnot(aaa[7, "PLACE"] == "50000")
aaa <- totalcensus:::read_acs1year_geo_(
year = 2009,
state = "US",
geo_headers = dict_acs_geoheader_2009_1year$reference
)
stopifnot(dim(aaa) == c(1275, 45))
# read_acs1_geoheaders_() ======================================================
aaa <- totalcensus:::read_acs1year_geoheaders_(
year = 2010,
states = "RI",
geo_headers = c("STATE", "CBSA")
)
cols <- c("GEOID", "NAME", "STATE", "CBSA", "GEOCOMP", "SUMLEV", "state",
"STUSAB", "lon", "lat")
stopifnot(names(aaa) == cols)
aaa <- totalcensus:::read_acs1year_geoheaders_(
year = 2017,
states = "US",
table_contents = c("B01001_002", "B01001_026"), # reference only
geo_headers = c("STATE", "COUNTY")
)
stopifnot(aaa[7, 6] == 550273)
# read_acs1year_areas_() =======================================================
aaa <- totalcensus:::read_acs1year_areas_(
year = 2016,
states = c("UT", "RI"),
table_contents = c("B01001_002", "B01001_026"),
areas = c("Salt Lake City city, UT",
"Providence city, RI",
"PLACE = RI19180"),
summary_level = "place",
with_margin = TRUE
)
stopifnot(aaa[3, 5] == 40454)
# read_acs1year() ==============================================================
# provide geo_headers
aaa <- read_acs1year(
year = 2015,
states = c("UT", "RI"),
table_contents = c("male = B01001_002", "female = B01001_026"),
geo_headers = c("COUNTY", "CBSA"),
summary_level = "county",
with_margin = TRUE,
dec_fill = "dec2010"
)
stopifnot(dim(aaa) == c(10, 16))
# provided areas
aaa <- read_acs1year(
year = 2016,
states = c("UT", "RI"),
table_contents = c("male = B01001_002", "female = B01001_026"),
areas = c("Salt Lake City city, UT",
"Providence city, RI",
"PLACE = RI19180"),
summary_level = "place",
with_margin = TRUE
)
cols <- c("area", "GEOID", "NAME", "PLACE", "population", "population_margin",
"male", "male_margin", "female", "female_margin", "GEOCOMP",
"SUMLEV", "state", "STUSAB", "lon", "lat")
stopifnot(names(aaa) == cols)
aaa <- read_acs1year(
year = 2017,
states = "CA",
table_contents = c("male = B01001_002", "female = B01001_026"),
summary_level = "place"
)
stopifnot(dim(aaa) == c(142, 11))
# _________________________ ===================================================
# check all file segments =====================================================
read_segment <- function(year, state, fs, est_marg){
aaa <- totalcensus:::read_acs1year_filesegment_(
year = year,
state = state,
file_seg = fs,
est_marg = est_marg
)
}
read_all_segment <- function(year, state = "RI", first_seg = 1, est_marg = "e"){
look <- get(paste0("lookup_acs1year_", year))
n_fs <- max(as.integer(look$file_segment))
for (i in first_seg:n_fs){
if (i < 10){
fs <- paste0("000", i)
} else if (i < 100){
fs <- paste0("00", i)
} else {
fs <- paste0("0", i)
}
#print(fs)
read_segment(year, state, fs, est_marg)
}
}
read_all_segment(2018, "RI", 1, "e")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.