Biostar Beta. Not for public use.
clusterprofiler: Pathview and dotplot
0
Entering edit mode
16 months ago
nsmaan • 10

Hello Experts. I have processed my RNA-seq data in a well-formatted text file (atleast i thought so..) that looks like this (~100+ genes) where i have Ensembl, entrezgene, gene(ID), group (up/down regulated), keggID, logFC, and padj, etc

> head(ids)
Ensembl entrezgene  gene    group   kegg    logFC   padj
ENSRNOG00000019500  24296   Cyp1a1  upregulated rno24296    2.441405    0.002797
ENSRNOG00000003280  24411   Grin2c  upregulated rno24411    6.245858    0.000071
ENSRNOG00000049590  24988   RT1-M2  downregulated   rno24988    -2.451328   0.000000
ENSRNOG00000014018  25280   St8sia1 downregulated   rno25280    -3.139369   0.000017
ENSRNOG00000001527  25408   Cd80    upregulated rno25408    2.669957    0.000050
ENSRNOG00000008157  29179   Syn2    upregulated rno29179    2.737939    0.001181
ENSRNOG00000011815  29517   Sgk1    downregulated   rno29517    -2.073134   0.000152
-----------------------------

But no matter i much i try i can't get the Pathview and dotplot working on this file based on clusterprofiler documentation (which actually has excellent description, but i am limited by my basic R skills).

I would appreciate anyone's help in plotting Pathview and dotplot on a file like mine...

Thank you, -N

RNA-Seq • 240 views
ADD COMMENTlink
0
Entering edit mode

If you show your code we can see where it goes wrong. Without any code it is impossible to diagnose.

ADD REPLYlink
0
Entering edit mode

My input looks like this:

data<-read.delim("test.txt") # File is tab-delimited
#file contains columns of ENSEMBL, entrez, GENENAME, logFC, padj
genes<-data$entrez[abs(data$logFC)>1.5 & data$padj<0.01 ]
ego<- enrichGO(genes,OrgDb= org.Rn.eg.db,ont= "CC",  pAdjustMethod = "BH", pvalueCutoff  = 0.01,qvalueCutoff  = 0.05, readable= TRUE)
cnetplot(ego, categorySize="pvalue", foldChange=genes)
#barplot(ego, showCategory=8)

I can see cnetplot output but do not see the lines connecting cluster centers with genes?

Also i am not sure how to generate 'compareCluster' data to select up and downregulated groups (based on logFC and padj) for dotplot(formula_res) example provided in ....doc/clusterProfiler.html

I would appreciate any help in this... Thank you

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1