您好,请问导入的数据是什么数据啊
rm(list=ls())
setwd("your_dir")
# read the file
mydata <- read.csv("your_diff_expr_file",row.names=1)
head(mydata)
# if (!requireNamespace("BiocManager", quietly = TRUE))
# install.packages("BiocManager")
#
# BiocManager::install("msigdbr")
library(clusterProfiler)
library(enrichplot)
library(msigdbr) # Molecular Signatures Database
library(org.Hs.eg.db) #人类GO注释数据