书接上回:为了满足要求,需要对数据进行t.test,并标注到图上。
最终作图效果是这样的啦:
为了能够将显著性标记换成P值,并以合适的计数方法输出,运用了t.test()函数。同时,运用一个自定义函数(根据数字的大小选择是否使用科学计数法),感觉应该有更简洁的实现方式,还需要继续学习鸭吼吼吼吼。
只能说chatgpt真是厉害,可以帮忙写代码。
具体代码如下:
library(gridExtra)
library(ggpubr)
library(ggsci)
library(ggsignif)
##install.packages(”ggsci“)
##没有package的话用上一行的命令,下载包。
#带显著性比较,作图,单因素双水平,t.test
#选择L,P行
list3 <- list()
for (i in 1:27){
list3[[i]] <-
list2[[i]] %>%
filter(!!sym(y[i]) %in% c("L", "P"))
}
##太大或太小的数字,选择合适计数法输出
format_number <- function(number, decimal_places &#