Panel A
Acetyle CoA Gene abundance and butyric acid relationship
metabolomics_db <- readRDS("../data/cleaned_data/metabolomics_data.rds")
ggplot(metabolomics_db) +
geom_smooth(aes(x=log(pyruvate), y=Butyric.acid), method=lm, color="#e05f5f") +
geom_point(aes(x=log(pyruvate), y=Butyric.acid),
shape=21, size=3, fill="#e05f5f", alpha=0.8) +
theme_classic()
cor.test(metabolomics_db$Butyric.acid, log(metabolomics_db$pyruvate))
Pearson's product-moment correlation
data: metabolomics_db$Butyric.acid and log(metabolomics_db$pyruvate)
t = 2.3571, df = 62, p-value = 0.02159
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
0.04407955 0.49751891
sample estimates:
cor
0.2867815