for (i in seq(5,62,2)) {
  value1 <- names(phe)[i]
  value2 <- names(phe)[i+1]
  phe %>% ggplot(aes(!!sym(value1),!!sym(value2)))+geom_point()+theme_bw()+stat_cor(method = "spearman")+
    theme(axis.text.x = element_text(size = 24,colour = "black",face = "bold",family = "serif"))+
    theme(axis.text.y = element_text(size = 24,colour = "black",face = "bold",family = "serif"))+
    theme(axis.title.x = element_text(size = 24,face = "bold",family = "serif"))+
    theme(axis.title.y = element_text(size = 24,face = "bold",family = "serif"))
  ggsave(filename = paste0(value1,".png"))
}