Segfault encountered in using read.maf to read mutation and Gistic results at the same simultaneously:
>maf=read.maf(fn,gisticAllLesionsFile='all_lesions.conf_99.txt',gisticAmpGenesFile='amp_genes.conf_99.txt',gisticDelGenesFile='del_genes.conf_99.txt')
reading maf..
Mutation_Status not found. Assuming all variants are Somatic and validated.
Excluding 2831 silent variants.
ID N
1: Samples 102
2: 3'UTR 337
3: 5'Flank 232
4: 5'UTR 262
5: RNA 310
6: Silent 1588
7: lincRNA 102
Processing Gistic files..
*** caught segfault ***
address 0xf1000013, cause 'memory not mapped'
Traceback:
1: rbindlist(l, use.names, fill, idcol)
2: data.table::.rbind.data.table(...)
3: rbind(deparse.level, ...)
4: rbind(ampGenes, delGenes)
5: readGistic(gisticAllLesionsFile = gisticAllLesionsFile, gisticAmpGenesFile = gisticAmpGenesFile, gisticDelGenesFile = gisticDelGenesFile, isTCGA = isTCGA)
6: read.maf(fn, gisticAllLesionsFile = "all_lesions.conf_99.txt", gisticAmpGenesFile = "amp_genes.conf_99.txt", gisticDelGenesFile = "del_genes.conf_99.txt")
The above error occurred when both amp_genes.conf_XX.txt and del_genes.conf_XX.txt are provided simultaneously. If provided only amp_genes.conf_XX.txt or del_genes.conf_XX.txt, it works.
Anyone here encountered the same issue?
Best regards,
Xiangchun