Question: Gene set enrichment analysis (GSEA) of some lncRNAs(known for Gencode and newly assembled)
gravatar for Ashley
2.0 years ago by
Ashley50 wrote:

Hello all, I have the FPKM values from NORMAL vs TUMOR.I am doing GSEA analysis using GSEA guy.

Below is my input(lncRNA_up.gct):


130 100 NAME DESCRIPTION SRR934732.sorted.bam SRR934781.sorted.bam RP3-340N1.5 na 0.405482833 0.709856841 AC116614.1 na 0.583072119 0.777716411 XLOC_000564 na 523.3443013 53.65888594 …… I have a gmt format file from GSEA(c2.cp.kegg.v6.0.symbols.gmt).When I try to Run GSEA,I get the following error: 4891 [INFO ] Begun importing: Dataset from: lncRNA_up.gct

of elements = 1

lncRNA_up.cls 5032 [INFO ] Begun importing: Template from: lncRNA_up.cls 5049 [INFO ] No dataset collapsing was done .. using original as is to parse>c2.cp.kegg.v6.0.symbols.gmt< got: [c2.cp.kegg.v6.0.symbols.gmt] 5049 [INFO ] Begun importing: GeneSetMatrix from: c2.cp.kegg.v6.0.symbols.gmt 5083 [INFO ] Got gsets: 186 now preprocessing them ... min: 15 max: 500 5101 [INFO ] Done preproc for smaller than: 15 xtools.api.param.BadParamException: After pruning, none of the gene sets passed size thresholds. at xtools.api.param.ParamFactory.checkAndBarfIfZeroSets( at xtools.gsea.Gsea.execute( at xtools.api.AbstractTool.tool_main( at xtools.gsea.Gsea.main( 5111 [INFO ] Renaming rpt dir on error to: /mnt/pub/work/liuzhe/projs-cas-ibp/xuelab/GSCs_lncRNA/GSEA_analysis/Up-reagulated/apr07/error_my_analysis.Gsea.1491526934954

My parameters are: java -Xmx1024m -cp /mnt/pub/work/liuzhe/tools/gsea2-2.2.4.jar xtools.gsea.Gsea -res lncRNA_up.gct -cls lncRNA_up.cls -gmx c2.cp.kegg.v6.0.symbols.gmt -collapse false

Has anyone encountered similar Error? Thank you very much!

ADD COMMENTlink written 2.0 years ago by Ashley50

I've used it previously for microarray data(form GSEA example dataset P53) and it worked fine.

ADD REPLYlink written 2.0 years ago by Ashley50

It seems that the filters to select the gene sets are too stringent, and you end up without any gene set to test with. Try increasing the limits of genes in gene sets. 5 in the lower band and 5000 in the upper band

ADD REPLYlink written 2.0 years ago by Lluís R.810

Thank you! From your comment, I known add a parameter min size and/or max size can overcome the question. Thanks for your kindly help!

ADD REPLYlink written 2.0 years ago by Ashley50
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1859 users visited in the last hour