df %>% filter(include == "TRUE") %>% ggplot(aes(var1,fill=re_analysis_dataset))+ geom_histogram(center = 0.5, binwidth = 0.1)