plot_credible_intervals <- readRDS(here("data/3_prime_batch_1/fast_pipeline_results/differential_composition/plot_credible_intervals.rds"))
plot_credible_intervals[[1]]$credible_intervals_1D
plot_credible_intervals[[1]]$boxplot[[1]] +
guides(fill=guide_legend(nrow=4,byrow=TRUE))
## Warning: Removed 1512 rows containing non-finite values (stat_summary).
Testing evolution of blood composition during time
COVID Atlas - https://atlas.fredhutch.org/fredhutch/covid/
atlas_composition = readRDS(here("data/covid_atlas/plot_composition_analysis.rds"))
knitr::include_graphics(here("data/covid_atlas/plot_atlas.pdf"))