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
# <chr> <chr> <dbl> <dbl> <chr> <chr> <chr>
# 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()
