library(tidyverse)
library(ggplot2)
library(labelled)
library(forcats)
library(haven)
# Read data or import data
an_data <- read_dta("data/ESPA_2021/ETAN81FLSP.DTA")
sc_data <- read_dta("data/ESPA_2021/ETSC81FLSP.DTA")
fc_data <- read_dta("data/ESPA_2021/ETFC81FLSP.DTA")
fc_data_fct <- to_factor(fc_data)
ggplot(fc_data_fct, aes(x = fct_rev(fct_infreq(factype)))) +
geom_bar(fill = "#2E86AB", width = 0.5) +
coord_flip()+
geom_text(stat = 'count', aes(label = after_stat(count)),
hjust = -0.1, size = 4) +
labs(title = "Healthcare Facility Type Distribution",
subtitle = "ESPA 2021 Facility Data",
x = "Facility Type",
y = "Count") +
theme_classic() +
theme(axis.line.x = element_blank() )