group,phylum) %>%
summarise(value_mean=mean(value)) %>%
ungroup() %>%
as.data.frame() %>% drop_na...,aov_data %>% select(1,2,4),by=c("phylum","group")) %>%
ggplot(aes(group,phylum,fill=value))+
geom_tile...rev(RColorBrewer::brewer.pal(3, "RdBu")))+
scale_x_discrete(expand = c(0,0),position = "top")+
scale_y_discrete