运用软件:Rstudio
运用的软件包:tidyverse、clusterProfiler(4.6.1以上)
library(readr)
library(tidyverse)
library(clusterProfiler)
emapper <- read_delim("E:/转录/workspace/cds1.emapper.annotations",
delim = "\t", escape_double = FALSE,
trim_ws = TRUE)
emapper <- dplyr::select(emapper,
GID = 'query',
SYMBOL = "Preferred_name",
GO = "GOs",
KO = "KEGG_ko",
Pathway = "KEGG_Pathway",
Gene_Name = "PFAMs")
gene_info <- dplyr::select(emapper,GID,Gene_Name,SYMBOL) %>%
dplyr::filter(!is.na(Gene_Name)) %>%
dplyr::filter(!is.na(SYMBOL))
gene_info$GID <- str_spli