In graveja0/hcmarket: Defininng and Analyzing Health Care Markets

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(hcmarket)
library(tidyverse)
options(dplyr.summarise.inform = FALSE)
library(ggdendro)
data(geo_hosp_A)
M_z <- geo_hosp_A %>% detect_markets(1)

M_z$dendro %>% ggdendrogram(rotate = TRUE)

M_j <- geo_hosp_A %>% detect_markets(2)

res <-
    0:M_j$max_height %>%
    purrr::map(~(calculate_firm_hhi(A = geo_hosp_A,markets=M_j$markets[[paste0(.x)]])$market_hhi))

res %>%
    bind_rows(.id = "level")  %>%
    select(level,N,hhi,hhi_km) %>%
    gather(measure,value,-level,-N) %>%
    as_tibble() %>%
    group_by(measure, level) %>%
    summarise(value = weighted.mean(value,w=N)) %>%
    mutate(level = as.numeric(paste0(level))) %>%
    ggplot(aes(x = level, y = value, colour= measure)) +
    geom_line() +
    hrbrthemes::theme_ipsum()