7093_HW7

Read ACS 2019 5 Year data

hw7ddi <- read_ipums_ddi("C:/UTSA//OneDrive - University of Texas at San Antonio//3_W_7093_GIS/GISClasswork/usa_00001.xml")

hw7data <- read_ipums_micro(hw7ddi)

hw7data <- haven::zap_labels(hw7data) 

names(hw7data) <- tolower(names(hw7data))

Get PUMA Geographic data

options(tigris_class = "sf")

ca_pumas <- pumas(state = "CA",
               year = 2019,
               cb = T)

# plot(ca_pumas["GEOID10"],
#     main = "Public Use Microdata Areas in California")

# mapview::mapview(ca_pumas, zcol= "GEOID10")

Recode “migrate1d” variable into “samehouse” binary variable

hw7data$samehouse <- Recode(hw7data$migrate1d, recodes = "10=1; 20:40=0; else=NA")

Survey design and estimate residential stability rate

hw7design <- svydesign(ids = ~cluster,
                       strata = ~strata, 
                       weights = ~perwt, 
                       data = hw7data)

puma_est_samehouse <- svyby (formula = ~samehouse,
                             by = ~puma,
                             design = hw7design, 
                             FUN = svymean,
                             na.rm = TRUE)

Join to PUMA geography data

ca_pumas$puma <- as.numeric(ca_pumas$PUMACE10)

geo7 <- left_join(ca_pumas, puma_est_samehouse, by=c("puma"="puma"))

Map

tmap_mode("view")

geo7 %>%
  tm_shape()+
  tm_polygons("samehouse",
              style = "kmeans",
              n = 8, 
              legend.hist = TRUE)+
  tm_layout(title = "Residential Stability Rate in California, 2015-2019")