tests/testthat/data/get_scb_distr_data.R

In bjornerstedt/scbapi: Import data from SCB (Statistics Sweden) using a JSON qurey

# NOTE: scripts works from the current directory as working directory
# Currently not used, as SCB does not have district data in JSON format!!!
# devtools::install_github("bjornerstedt/scbapi")
library(scbapi)
library(tidyr)
library(dplyr)
library(stringr)
library(jsonlite)
library(httr)
# Set to TRUE to download and save data as json files
# If FALSE, it uses the saved data files instead
get_data <- TRUE

# Population
pop <- get_scb("pop_distr", get_data) %>%
  group_by(Tid, Region, Kon) %>%
  summarise(value = sum(value)) %>%
  spread(Kon, value) %>%
  select(year=Tid, district = Region, popwomen = kvinnor, popmen = män) %>%
  ungroup() %>%
  mutate(
    year = as.numeric(year),
    district = str_sub(district, end = 2)
    )

# CPI
cpi <- get_scb( "cpi1", get_data) %>%
    bind_rows(get_scb("cpi2", get_data) ) %>%
    rename(cpi = indata) %>%
    separate(time, c("year", "month"), sep = "M", convert = TRUE)

# Real GDP
gdp <- get_scb("gdpreal", get_data) %>%
    rename(GDPreal = indata) %>%
    separate(time, c("year", "quarter"), sep = "K", convert = TRUE) %>%
    na.omit

# Sick leave, pre 2005
sick1 <- get_scb("sick1", get_data, colnames = c("sm", "sw", "time"))

# sick leave from 2005
sick2 <- get_scb("sick2", get_data, colnames = c("sm", "sw", "time")) %>%
    mutate(
        sm = 10 * sm,
        sw = 10 * sw
    )
sick <- bind_rows(sick1, sick2) %>%
    separate(time, c("year", "month"), sep = "M", convert = TRUE)

hexp <- get_scb("health_expend_qu") %>%
    rename(healthexp = indata) %>%
    separate(time, c("year", "quarter"), sep = "K", convert = TRUE)

################## Combine Data ####################

scbdata <- cpi %>%
    mutate(quarter = (month - 1) %/% 3 + 1) %>%
    left_join(gdp, by=c("year","quarter")) %>%
    left_join(pop, by="year") %>%
    left_join(sick, by=c("year","month")) %>%
    left_join(hexp, by=c("year","quarter")) %>%
    select(-quarter) %>%
    filter( year>= 1995)

save(scbdata, file = "../scbdata.Rdata")
saveRDS(scbdata, file = "../scbdata.Rds")