calotriton_resp1 %>%
filter(time_point!="0") %>%
ggplot(aes(x = time_point, y = o2_normalized, color=time_point, fill=time_point)) +
geom_boxplot()+
geom_jitter() +
scale_color_manual(
name = "Time point",
breaks = c("1", "2"),
values = c('#ffa500', "#b37400")
) +
scale_fill_manual(
name = "Time point",
breaks = c("1", "2"),
values = c('#ffa50050', "#b3740050")
)+
theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
facet_wrap(~ factor(temperature)+elevation)+
stat_compare_means(size=3)+
theme(legend.position="right")+
labs(x = "Time_point", y= "Normalized o2 consumption (ml/s/g)")
calotriton_resp1 %>%
filter(time_point!="0") %>%
ggplot(aes(x = time_point, y = o2_normalized, color=time_point, fill=time_point)) +
geom_boxplot()+
geom_jitter() +
scale_color_manual(
name = "Time point",
breaks = c("1", "2"),
values = c('#ffa500', "#b37400")
) +
scale_fill_manual(
name = "Time point",
breaks = c("1", "2"),
values = c('#ffa50050', "#b3740050")
)+
theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
facet_grid(temperature~factor(elevation)+population)+
stat_compare_means(size=3)+
theme(legend.position="right")+
labs(x = "Time_point", y= "Normalized o2 consumption (ml/s/g)")