I have data such as this:
dat_in <- data.frame(Happy=c(1,1,0,NA,0,0,1), Exercise=c(1,3,4,NA,1,1,3)) %>% mutate_all(as.factor)%>% mutate(across(matches(\