library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(coriverse)
library(DBI)

con <- connect_to_db('sch_source')
latent_inno <- cori.db::read_db(con, "goetz_latent_innovation")
dbDisconnect(con)

readr::write_csv(latent_inno, '../export/data/goetz_latent_innovation.csv')

Rural Latent Innovation

source('../scripts/utils-data.R')
cbsa <- getCBSA()


rural_latent_inno <- dplyr::left_join(
    latent_inno,
    cbsa,
    by = c("fips" = "geoid")
  ) %>%
  dplyr::filter(
    rural_cbsa_flag == 1
  )

readr::write_csv(rural_latent_inno, '../export/data/goetz_latent_innovation_rural.csv')
sum_stats <- latent_inno %>%
  dplyr::left_join(
    cbsa,
    by = c("fips" = "geoid")
  ) %>%
  dplyr::group_by(metropolitan_designation) %>%
  dplyr::summarise(
    mean_innovation = mean(innovation, na.rm = TRUE),
    sd_innovation = sd(innovation, na.rm = TRUE),
    sum_innovation = sum(innovation, na.rm = TRUE)
  )

knitr::kable(sum_stats)
metropolitan_designation mean_innovation sd_innovation sum_innovation
Metro 0.4124808 0.9376038 486.72734
Micro -0.0993061 0.8819528 -65.54204
Non-CBSA -0.3237543 0.9825001 -418.93804
NA -0.3158961 0.7329809 -1.57948