SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
GSEA: None of the gene sets passed the size thresholds id0 Bioinformatics 3 03-08-2017 07:42 AM
Edge-pro error chariko Bioinformatics 0 03-02-2015 01:25 AM
DESeq and GSEA with pre-ranked gene list proteomania Bioinformatics 7 02-17-2015 10:08 AM
Gene Set Enrichment Analysis uqfgaiti Bioinformatics 1 07-12-2014 03:14 PM

Reply
 
Thread Tools
Old 02-28-2018, 09:52 PM   #1
seqprone
Junior Member
 
Location: USA

Join Date: Oct 2015
Posts: 7
Default GSEA Leading edge analysis fails with gene set not found error

hi all!

I'm trying to run a leading edge analysis on my RNASeq data.

Hi everyone!

I'm trying to run a leading edge analysis on my RNASeq data using GSEA

I first create a ranked gene set file and run GSEA Prerank analysis.

I use the GMT file from http://download.baderlab.org/EM_Gene...018_symbol.gmt

Once I run the GSEA Preranked analysis, I try running the Leading edge analysis and run into this error https://imgur.com/a/zX3YI

It says:


---- Stack Trace ----
# of exceptions: 1
------In GeneSetMatrix: gene_sets.gmt no GeneSet found with name: TRANSMEMBRANE RECEPTOR PROTEIN SERINE\THREONINE KINASE SIGNALING PATHWAY%GOBP%GO:0007178

I open the gene_sets.gmt file in notepad++ and search for the "missing" geneset and I find it to be present:

The line looks like this:

TRANSMEMBRANE RECEPTOR PROTEIN SERINE/THREONINE KINASE SIGNALING PATHWAY%GOBP%GO:0007178 transmembrane receptor protein serine/threonine kinase signaling pathway Ddx5 Sub1 Zyx Bmpr1b Chrd Trps1 Tob1 Scx Inha Nup93 Ptprk Tmem100 Id1 Trf Mir23a Ccl2 Zcchc18 Pml Dand5 Vim Zfyve9 Map3k7 Acvr2b Hfe2 Parp1 Eid2 Acvr1b Itgb5 Lrp4 Map3k1 Bmp6 Acvr1c Bmp5 Mtmr4 Bmp3 Dcp1a Tab1 Akap2 Col1a2 Col3a1 Ppm1l Magi2 Rgmb Twsg1 Src Rbm14 Atoh8 Nceh1 Tgfb3 Tgfbr1 Nlk Creb1 Usp15 Tgfb2 Smad6 Amhr2 Inhba Eng Sptbn1 Fermt2 Tgfbr3 Dusp22 Tgfbr2 Tgfb1 Ecsit Jun Gdf10 Rnf111 Pxn Smad1 Gdf15 Smad2 Btbd11 Smad5 Smad3 Smad4 Pdcd4 Fut8 Cldn5 Egr1 Mapk3 Adam9 Trp53 Acvr1 Fos Tgfb1i1 Hfe Stoml1 4632428N05Rik Bmpr2 Smad7 Sirt1 Lnpep Inhbb Slc33a1 Bmp8a Bmp8b Smurf1 Msx2 Usp9x Bmp4 Cited2 Bmp7 Fam83g Ptk2 Myd88 Bmpr1a Skil Mapk14 Runx2 Akap7 Megf8 Smurf2 Gdf5 Acvrl1 Rgma Meg3 Dlx5 Gdf9 Hpgd Itgb1 Fst Ski Bmp15 Tgif1 Ror2 Tgif2 Nog Zfp128 Acvr2a Ltbp3 Ltbp1 Ltbp4 Gdf11 Stub1 Hipk2 Bambi Fkbp1a Gcnt2 Lefty1 Msx1 Lef1 Cited1 Cdh5 Arrb2 Bmp2 Hexa Smad9

So, if it's present in the gene_sets.gmt, why am I getting this error?

Could someone please help me troubleshoot this?
seqprone is offline   Reply With Quote
Reply

Tags
gsea, leading edge analysis, rnaseq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:33 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO