请在此处找到我需要可视化的数据。我想揭示是否有任何变量与参考变量显着不同。我绘制条形图,列出将相互比较的变量,进行统计并绘制显着性线,如下所示。到目前为止一切都很好。
library(tidyverse)
library(here)
library(readxl)
library(ggforce)
library(ggplot2)
library(ggsignif)
library(ggpubr)
library(rstatix)
data <- read.csv(here("data", "data.csv"), sep = ";")
# List the variables to be compared with each other in pairwise comparison
my_comparisons <- list( c("%0", "250 µg/ml"), c("%0", "100 µg/ml"), c("%0", "50 µg/ml"), c("%0", "25 µg/ml"))
#Statistics
stat.test <- data %>%
t_test(OD ~ Material_Concentration, p.adjust.method = "bonferroni", comparisons = my_comparisons) %>%
add_significance("p.adj")
stat.test
# add significance
stat.test <- stat.test %>%
add_xy_position(x = "Material_Concentration")
#plot the data
p <- ggplot(data, aes(Material_Concentration, OD, colour = Material_Concentration))+
stat_summary(fun = mean, geom = "col", aes(fill = Material_Concentration), show.legend = FALSE)+
stat_summary(fun.data = mean_se, geom = "errorbar", width = 0.2, show.legend = FALSE)+
xlab("Treatment Concentrations")+
ylab("Average OD")+
theme(panel.background = element_blank(), plot.title = element_text(color="red", size=40, face="bold.italic", hjust = 0.5),
axis.text.x = element_text(size = 20, angle = 45, hjust = 1),
axis.text.y = element_text(size = 15),
axis.title.y = element_text(size = 15),
axis.title.x = element_text(size = 15))+
stat_pvalue_manual(stat.test, label = "p.adj.signif")
p
但是,当我重新排列变量在 x 轴上的位置时,通过使用scale-x-discrete
有效线的 x 位置不会改变。
library(tidyverse)
library(here)
library(readxl)
library(ggforce)
library(ggplot2)
library(ggsignif)
library(ggpubr)
library(rstatix)
data <- read.csv(here("data", "data.csv"), sep = ";")
# List the variables to be compared with each other in pairwise comparison
my_comparisons <- list( c("%0", "250 µg/ml"), c("%0", "100 µg/ml"), c("%0", "50 µg/ml"), c("%0", "25 µg/ml"))
#Statistics
stat.test <- data %>%
t_test(OD ~ Material_Concentration, p.adjust.method = "bonferroni", comparisons = my_comparisons) %>%
add_significance("p.adj")
stat.test
# add significance
stat.test <- stat.test %>%
add_xy_position(x = "Material_Concentration")
#plot the data
p <- ggplot(data, aes(Material_Concentration, OD, colour = Material_Concentration))+
stat_summary(fun = mean, geom = "col", aes(fill = Material_Concentration), show.legend = FALSE)+
stat_summary(fun.data = mean_se, geom = "errorbar", width = 0.2, show.legend = FALSE)+
xlab("Treatment Concentrations")+
ylab("Average OD")+
theme(panel.background = element_blank(), plot.title = element_text(color="red", size=40, face="bold.italic", hjust = 0.5),
axis.text.x = element_text(size = 20, angle = 45, hjust = 1),
axis.text.y = element_text(size = 15),
axis.title.y = element_text(size = 15),
axis.title.x = element_text(size = 15))+
scale_x_discrete(limits = c("250 µg/ml", "100 µg/ml", "50 µg/ml", "25 µg/ml", "%0", "%10"))+
stat_pvalue_manual(stat.test, label = "p.adj.signif")
p
我怎么解决这个问题?如果您能揭露我的错误,我将非常高兴。
谢谢