Question: GISTIC 2.0.23 Error: Floating Point Exception
gravatar for rpaa
4 days ago by
rpaa0 wrote:

Hello everyone,

I have been trying to perform GISTIC analysis in Linux on methylation CNV data of ~650 tumor samples.

I use the following command:

gistic2 -b <directory> -seg gistic.segments -refgene hg19.mat conf -0.95

The head of the gistic.segments file is as follows: ID chrom loc.start loc.end num.mark seg.mean

Sample1 chr1 560684 145225000 862 0.156

Sample1 chr1 145375000 249195311 644 -0.13

Sample1 chr10 105000 107175000 589 -0.02

Sample1 chr10 108425000 135462374 202 -0.273

Sample1 chr11 130000 134873258 881 -0.149

Sample1 chr12 172870 34375000 207 0.121

Sample1 chr12 34425000 34703347 4 -0.424

Sample1 chr12 38178347 133770948 607 -0.107

Sample1 chr13 19110000 115079939 416 -0.046

The tool works perfectly if I run ~430 samples or less. In this case it doesn't matter which samples I select out of the 650. However, if I try to run more samples I run in to an error due to a floating point exception:

Finding peaks on chromosome 1

Finding peaks on chromosome 2

Amplification peak at chr2:204700000-206675000(36452:36650)

q-value = 0.0021198 & residual q-value = 0.0021198

Finding peaks on chromosome 3

Amplification peak at chr3:8500000-9525000(39381:39467)

q-value = 0.005833 & residual q-value = 0.005833

Finding peaks on chromosome 4

Resolving conflicts...

Amplification peak at chr4:54625000-54975000(54479:54514)

q-value = 1.023e-22 & residual q-value = 1.023e-22

Finding peaks on chromosome 5

Finding peaks on chromosome 6

Floating point exception detected at Wed Jan 15 13:52:45 2020


Crash Decoding : Disabled

Current Visual : None

Default Encoding : UTF-8

GNU C Library : 2.27 stable

MATLAB Architecture: glnxa64

MATLAB Root : /usr/share/mcr/v83

MATLAB Version : (R2014a)

Operating System : Linux 4.15.0-51-generic #55-Ubuntu SMP Wed May 15 14:27:21 UTC 2019 x86_64

Processor ID : x86 Family 6 Model 62 Stepping 4, GenuineIntel

Virtual Machine : Java 1.7.0_11-b21 with Oracle Corporation Java HotSpot(TM) 64-Bit Server VM mixed mode

Window System : No active display

Fault Count: 1

Abnormal termination:

Floating point exception

The chromosome at which the tool crashes varies (in this case chromosome 6), but it will crash somewhere if I use more than 430 samples.

Does anybody know what I'm doing incorrectly?

Thank you very much for your help!

gistic • 76 views
ADD COMMENTlink modified 4 days ago • written 4 days ago by rpaa0

It is not an answer at all: what is methylation cnv data? How did you get it?

In my experience all the errors from gistic were caused by faulty samples (one per cohort with one "bad" cnv) or wrong markers file. Can you try to support your run with manually generated markers file, taking into account that all of your cnvs have to contain at least one marker within?

ADD REPLYlink written 4 days ago by German.M.Demidov1.2k
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: 1118 users visited in the last hour