Deleted:Customize biological processes shown in KEGG analysis dotplot
0
0
Entering edit mode
@3b25d11d
Last seen 3.2 years ago
United States

I am semi-new to R and bioinformatics, but I want to be able to customize the pathways that are shown in my dotplot after gseKEGG pathway analysis. Currently, only 15 or so pathways are differentially expressed according to this KEGG pathway analysis. However, I still want to show pathways that are not differentially enriched. I have attached what my current dotplot looks like, as well as the code for it, and a dot plot of what I want my data to look like. Any help would be greatly appreciated.

KEGG Pathway Analysis

#place ENSEMBL into ENTREZ id format
ids<-bitr(names(original_gene_list), fromType = "ENSEMBL", toType = "ENTREZID", OrgDb=organism)
dedup_ids = ids[!duplicated(ids[c("ENSEMBL")]),]
df2 = g[g$gene_id %in% dedup_ids$ENSEMBL,]
df2$Y = dedup_ids$ENTREZID
# Create a vector of the gene unuiverse
kegg_gene_list <- df2$log2FoldChange
#Name vector with ENTREZ ids
names(kegg_gene_list) <- df2$Y
head(df2$Y)
list(df2)
# omit any NA values 
kegg_gene_list<-na.omit(kegg_gene_list)

# sort the list in decreasing order (required for clusterProfiler)
kegg_gene_list = sort(kegg_gene_list, decreasing = TRUE)
kegg_organism = "mmu"
kk2 <- gseKEGG(geneList     = kegg_gene_list,
               organism     = kegg_organism,
               nPerm        = 10000,
               minGSSize    = 5,
               maxGSSize    = 800,
               pvalueCutoff = 0.05,
               pAdjustMethod = "none",
               keyType       = "ncbi-geneid")
dotplot(kk2, showCategory = 10, title = "Enriched Pathways" , split=".sign") + facet_grid(.~.sign) + scale_color_gradient(low = "#56B1F7", high = "#132B43")
ridgeplot(kk2, showCategory = 30) + labs(title = "KEGG_Pathway Analysis", x = "enrichment distribution") + scale_fill_gradient(low = "#56B1F7", high = "#132B43")

This is what I currently have: enter image description here

This is what I want it to look like (I will be looking at non-ciliary pathways, but hopefully you get the idea): enter image description here

RNASeqR GO clusterProfiler rnaseq KEGG • 3.2k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 417 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6