售后

编辑本页后请点击“保存”。请参阅 syntax 了解维基语法。只有在您能改进该页面的前提下才编辑它。如果您想尝试一些东西，请先到 playground 热身。

====== 富集结果桑基图展示 ======
<code>
suppressMessages({
library(tidyverse)
library(ggsankey)
library(cols4all)
#BiocManager::install("dittoSeq")
library(dittoSeq)
library(argparser)})

argv <- arg_parser('')
argv <- add_argument(argv,"--stat", help="the stat file")
argv <- add_argument(argv,"--type", help="the plot type")
argv <- add_argument(argv,"--prefix", help="the prefix of outfile")
argv <- parse_args(argv)

stat <- argv$stat
type <- argv$type
prefix <- argv$prefix

stat_frame <- read.delim(stat,head=TRUE,sep='\t',quote='')
if (type=='GO') {
  data <- dplyr::select(stat_frame, geneName, Description)%>%
    separate_rows(geneName, sep = '/')
  colnames(data) <- c('gene', 'pathway')
}

if (type=='KEGG') {
  data <- dplyr::select(stat_frame, geneName, Description)%>%
    separate_rows(geneName, sep = '/')
  colnames(data) <- c('gene', 'pathway')
}

data2 <- make_long(data, gene, pathway )

p <- ggplot(data2, aes(x = x,
                      next_x= next_x,
                      node= node,
                      next_node= next_node,
                      fill= node,
                      label= node)) +
  geom_sankey(flow.fill="#DFDFDF",
              flow.color="grey60",
              node.fill=dittoColors()[1:length(unique(data2$node))],
              width=0.15) +
  geom_sankey_text(size = 4,
                   color= 'black',
                   hjust=1) +
  theme_void()

pdf(file=paste(prefix,'.pdf',sep=''))
p
dev.off()

svg(filename=paste(prefix,'.svg',sep=''))
p
dev.off()
ggsave(file=paste(prefix,'.png',sep=''),type='cairo-png',plot=p)
</code>
===== 展示图 =====

右图为文献中的图，可以将左图与结果中气泡图结合在一起，复现右图

编辑摘要

当您选择开始编辑本页，即寓示你同意将你贡献的内容按下列许可协议发布： CC Attribution-Share Alike 4.0 International

售后

用户工具

站点工具

页面工具