0

请在此处找到我需要可视化的数据。我想揭示是否有任何变量与参考变量显着不同。我绘制条形图,列出将相互比较的变量,进行统计并绘制显着性线,如下所示。到目前为止一切都很好。

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

此图像中重要线的定位不正确

我怎么解决这个问题?如果您能揭露我的错误,我将非常高兴。

谢谢

4

0 回答 0