R语言处理CCK8或MTT数据，一步绘制生长曲线

CCK8数据主要是统计OD450的吸收值，一般的数据形式如下，保存为csv格式

cck.csv

cck8 <- read.csv("~/Desktop/cck8.csv") 
## 宽数据转换为长数据
library(reshape2) 
cck8_long<-melt(cck8,
                id.vars = c('Time'),
                variable.name='Sample',
                value.name='OD450')
## 新增一个分组
cck8_long$Group=rep(c('control', 'treat'), each = 15)  ##分为两组，每组15个数据
### 使用ggpubr画图，可以使用prism主题
library(ggpubr)
library(ggprism)
ggline(cck8_long, x = "Time", y = "OD450", 
      color = "Group",shape='Group', 
      add = "mean_sd",ylab = "OD450 value",
      ggtheme = theme_prism(base_size = 12))

一步出图

加个P值，调整坐标，适当美化

• 首先预统计看下结果

compare_means(OD450~Group, data=cck8_long, group.by = "Time",method = 'anova')
####################################################################################
# A tibble: 5 x 7
#  Time  .y.           p   p.adj p.format p.signif method
#                     
# 1 day1  OD450 0.0823    0.082   0.08228  ns       Anova 
# 2 day2  OD450 0.00701   0.014   0.00701  **       Anova 
# 3 day3  OD450 0.000371  0.0015  0.00037  ***      Anova 
# 4 day4  OD450 0.0000675 0.00034 6.7e-05  ****     Anova 
# 5 day5  OD450 0.00309   0.0093  0.00309  **       Anova 
####################################################################################

结果还不错，添加显著标志

ggline(cck8_long, x = "Time", y = "OD450", 
       color = "Group",shape='Group', title = 'xxx cell',
       add = "mean_sd", palette = 'lancet',ylab = "OD450 value",legend = c(0.1,0.9),
       ggtheme = theme_prism(base_size = 12))+stat_compare_means(label = "p.signif",
                                                                 aes(group=Group),
                     method = "anova",  
                     hide.ns = TRUE)

加标志，加配色

细节地方AI处理一下,当然如果你安装了esquisse，你也可以直接导出成PPT

esquisse:::ggplot_to_ppt()