p = ggplot(to_plot_fr) + geom_col(mapping = aes(x = year, y = patient_num, fill = miss_or_nonmiss)) + facet_wrap(~province)