Analýza projektů s nehierarchickou lokalizací

prjs_geostatus <- read_parquet(here::here("data-processed",
                                          "projects-geo-check-groups.parquet"))
places_geostatus <- read_parquet(here::here("data-processed",
                                            "places-geo-check.parquet"))
ops <- read_parquet(here::here("data-processed", 
                               "op-codes.parquet"))
prj_meta <- read_parquet(here::here("data-processed", 
                                    "prj-esif-meta.parquet"))

head(prjs_geostatus$prj_id) # memo jak vypadá číslo projektu

prjs_groupable <- prjs_geostatus %>% 
  separate(prj_id, sep = "/", into = c("prog", "obj", "what", "vyzva", "proj"), 
           remove = F) %>% 
  separate(vyzva, sep = "_", into = c("vyzva_rok", "vyzva_cislo"), remove = F)

prj_ids_nonhierarchical <- prjs_geostatus$prj_id[prjs_geostatus$geostatus == "více míst nehierarchicky"]

library(ggupset)

gst_for_upset <- dtl %>% 
  filter(typ == "id" & !is.na(value),
         prj_id %in% prj_ids_nonhierarchical) %>% 
  select(op_id, prj_id, level) %>%
  distinct() %>% 
  mutate(level = as.character(level)) %>% 
  group_by(prj_id) %>% 
  summarise(geos = list(level)) %>%
  left_join(dt %>% filter(prj_radek == 1) %>% select(prj_id, op_id)) %>% 
  # mutate(nr = map_int(geos, nrow)) %>%
  # arrange(desc(nr)) %>% 
  filter(TRUE)

ggplot2::update_geom_defaults("line", list(colour = "black"))
ggplot(gst_for_upset) +
  geom_bar(mapping = aes(x = geos, fill = op_id), position = "stack") +
  scale_x_upset(order_by = "degree", sets = geolevels) +
  # scale_y_continuous(limits = c(0, 6e4)) +
  # scale_y_log10(limits = c(1, 1e5), n.breaks = 10,
  #               labels = scales::label_number(1)) +
  # scale_y_continuous(trans = scales::log_trans(), breaks = scales::pretty_breaks()) +
  ptrr::theme_ptrr() +
  scale_fill_brewer(type = "qual")

ggplot2::update_geom_defaults("line", list(colour = "darkblue"))

ggplot2::update_geom_defaults("line", list(colour = "black"))

ggplot(gst_for_upset) +
  geom_bar(mapping = aes(x = geos)) +
  scale_x_upset(order_by = "degree", sets = geolevels) +
  # scale_y_log10(n.breaks = 6, labels = scales::label_number(1)) +
  # scale_y_continuous(trans = scales::log_trans(), breaks = scales::pretty_breaks()) +
  ptrr::theme_ptrr(multiplot = T) +
  facet_wrap(~op_id, scales = "fixed")

ggplot2::update_geom_defaults("line", list(colour = "darkblue"))

dtl_with_geostatus <- dtl %>% 
  left_join(prjs_geostatus %>% 
              select(prj_id, geostatus)) %>% 
  select(op_id, prj_id, level, typ, geostatus, value) %>% 
  filter(typ %in% c("nazev", "id"), geostatus == "více míst nehierarchicky") %>% 
  drop_na(value)

dtl_with_geostatus_places <- dtl %>% 
  filter(typ == "id") %>% 
  drop_na(value) %>% 
  left_join(places_geostatus %>% 
              select(prj_id, value, level, value_has_valid_parent)) %>% 
  left_join(prjs_geostatus %>% 
              select(prj_id, geostatus)) %>% 
  select(op_id, prj_id, level, typ, value_has_valid_parent, value, 
         prj_geostatus = geostatus)

data_geocheck_file <- here::here("data-processed", "dt-geohierarchy-check.parquet")

data_geocheck <- read_parquet(data_geocheck_file)

# data_geocheck %>% 
#   filter(prj_id %in% prj_ids_nonhierarchical) %>% 
#   arrange(prj_id)

data_geocheck %>% 
  group_by(prj_id, level, value) %>% 
  filter(all(!levels_ok)) %>% 
  # filter(str_detect(prj_id, "CZ.04")) %>% 
  ungroup() %>% 
  count(prj_id, sort = T)