SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq and GSEA with pre-ranked gene list proteomania Bioinformatics 7 02-17-2015 10:08 AM
GSEA for RNA-seq data sazz Bioinformatics 3 10-20-2013 05:47 AM
Creating Gene Sets - CummeRbund SHeaph Bioinformatics 0 01-17-2013 07:56 AM
From Affy probe sets/gene symbols to genomic coordinates? ETHANol Epigenetics 7 10-25-2010 02:13 AM
Programs comparing gene sets 454andSolid Bioinformatics 0 03-16-2010 02:22 AM

Reply
 
Thread Tools
Old 08-28-2013, 02:59 PM   #1
id0
Senior Member
 
Location: USA

Join Date: Sep 2012
Posts: 130
Default GSEA: None of the gene sets passed the size thresholds

I am trying to use GSEA for some RNA-Seq data. I've used it previously for microarray data and it worked fine. My guess is that using well-defined probe set IDs helped.

For my data, I tried using the provided gene symbols as well as symbols exactly matching what I have. It seems no matter what I do, I end up with "None of the gene sets passed the size thresholds" error.

Using GenePattern (which I assume should be the safest option), I get the following output:
Code:
1286 [INFO ] Begun importing: Chip from: /xchip/gpprod-upload/servers/genepattern/users/uploads/tmp/run4907393233174326312.tmp/chip.platform.file/1/same.chip
1334 [WARN ] Missing chip file: >/xchip/gpprod-upload/servers/genepattern/users/uploads/tmp/run4907393233174326312.tmp/chip.platform.file/1/same.chip<	at edu.mit.broad.vdb.chip.FileInMemoryChip.initHere(?:?)
1502 [INFO ] Parsed from dotchip : 21862
1350 [WARN ] Missing chip file: >/xchip/gpprod-upload/servers/genepattern/users/uploads/tmp/run4907393233174326312.tmp/chip.platform.file/1/same.chip<	at edu.mit.broad.vdb.chip.FileInMemoryChip.initHere(?:?)
1738 [INFO ] Collapsing dataset was done. Original: 21862x2 (ann: 21862,2,same.chip) collapsed: 860x2 (ann: 860,2,GENE_SYMBOL)
to parse>c5.all.v4.0.symbols.gmt< got: [c5.all.v4.0.symbols.gmt]
1763 [INFO ] Begun importing: GeneSetMatrix from: c5.all.v4.0.symbols.gmt
2110 [INFO ] Got gsets: 1454 now preprocessing them ... min: 3 max: 500
Done removeGeneSetsSmallerThan: 3 for: 501 / 1454
Done removeGeneSetsSmallerThan: 3 for: 1001 / 1454
2259 [INFO ] Done preproc for smaller than: 3
2428 [INFO ] Renaming rpt dir on error to: error_.
2276 [WARN ] Could not rename for error to: error_.	at edu.mit.broad.genome.reports.api.ToolReport.setErroredOut(?:?)
Those warnings look questionable, but they are not exactly informative. Why would it say "Missing chip file" when the chip file is obviously present?
id0 is offline   Reply With Quote
Old 10-23-2014, 11:33 AM   #2
snaporaz
Junior Member
 
Location: Italy

Join Date: Oct 2014
Posts: 1
Default

hi!

i just started encountering the same error. ('Renaming rpt dir on error to: error_.')
differently from you i'm still working on some microarray data, and i also tried using collapse option and also symbols that mach exactly my dataset probes.

were you able to resolve the problem? could anyone have some advice?

thanks!

p.s. i'm sure that threshold is well above number of the genes in my geneset.

Last edited by snaporaz; 10-23-2014 at 11:35 AM. Reason: clarification
snaporaz is offline   Reply With Quote
Old 10-23-2014, 06:35 PM   #3
Elsie
Member
 
Location: Australia

Join Date: Mar 2011
Posts: 85
Default

I haven't used the Broads GSEA in a while, but for RNAseq data you could try a Bioconductor package:
http://www.bioconductor.org/packages...l/SeqGSEA.html
Elsie is offline   Reply With Quote
Old 03-08-2017, 07:42 AM   #4
jamesjcai
Junior Member
 
Location: Texas

Join Date: May 2012
Posts: 1
Default

Set 'Collapse dataset to gene symbols' false -- might help.

Last edited by jamesjcai; 03-08-2017 at 07:51 AM.
jamesjcai is offline   Reply With Quote
Reply

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 09:05 AM.


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