knitr::opts_chunk$set(echo = TRUE) knitr::opts_chunk$set(message = FALSE) knitr::opts_chunk$set(warning = FALSE) library(ProvidenciaChemo) library(tidyverse) library(patchwork) library(tidybayes) library(modelr) library(emmeans) library(magrittr) library(emmeans) library(rstanarm) theme_set(theme_my)
library(scales) TAsucc <- read_csv('extdata/Figure_2E_3D_TAsucc.csv') %>% filter(food %in% c("OP50", "JUb39")) %>% mutate(food = fct_relevel(food, c("OP50", "JUb39"))) plot <- TAsucc %>% ggplot(aes(x = food, y = AUC)) + #geom_point(aes(colour = food)) + geom_bardots(fillvar = food, dotvar = food) + stat_summary(geom = "errorbar", fun.data = mean_se, width = 0.2) + scale_color_plot("grey-blue", drop = TRUE) + scale_fill_plot("grey-blue", drop = TRUE) + add.n('food', y.pos = 0.5) + labs(x = "food", y = "AUC") + guides(colour = FALSE) + theme(axis.text.x = element_blank()) + facet_grid(.~genotype) + scale_y_continuous(trans = 'log10', breaks = trans_breaks('log10', function(x) 10^x), labels = trans_format('log10', math_format(10^.x))) lin_mod <- lm(TAsucc, formula = AUC ~ genotype * food) lin_mod %>% emmeans::ref_grid() %>% emmeans::contrast(method = "pairwise") plot
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.