test/test_read_acs1year.R

In GL-Li/totalcensus: Extract Decennial Census and American Community Survey Data

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")