weizhouUMICH / SAIGE

GNU Lesser General Public License v3.0
187 stars 72 forks source link

SAIGE memory error #358

Closed lanjblue closed 2 years ago

lanjblue commented 3 years ago

I was running SAIGE on multiple chromosomes. All went well except chr4. Here's the error. Do you know what could be the cause of this type of error and how to fix it? Thank you so much for your suggestions!

     user    system   elapsed
 588492.4 1280165.4  183101.4
[1] 2250000
numPassMarker:  2070290

 *** caught bus error ***
address 0x7f1bbf4b3030, cause 'non-existent physical address'

Traceback:
 1: SPAGMMATtest(vcfFile = opt$vcfFile, vcfFileIndex = opt$vcfFileIndex,     vcfField = opt$vcfField, bgenFile = opt$bgenFile, bgenFileIndex = opt$bgenFileIndex,     savFile = opt$savFile, savFileIndex = opt$savFileIndex, idstoExcludeFile = opt$idstoExcludeFile,     idstoIncludeFile = opt$idstoIncludeFile, rangestoExcludeFile = opt$rangestoExcludeFile,     rangestoIncludeFile = opt$rangestoIncludeFile, chrom = opt$chrom,     start = opt$start, end = opt$end, IsDropMissingDosages = opt$IsDropMissingDosages,     sampleFile = opt$sampleFile, GMMATmodelFile = opt$GMMATmodelFile,     varianceRatioFile = opt$varianceRatioFile, SAIGEOutputFile = opt$SAIGEOutputFile,     minMAF = opt$minMAF, minMAC = opt$minMAC, numLinesOutput = opt$numLinesOutput,     IsOutputAFinCaseCtrl = opt$IsOutputAFinCaseCtrl, IsOutputNinCaseCtrl = opt$IsOutputNinCaseCtrl,     condition = opt$condition, maxMAFforGroupTest = opt$maxMAFforGroupTest,     groupFile = opt$groupFile, sparseSigmaFile = opt$sparseSigmaFile,     minInfo = opt$minInfo, cateVarRatioMinMACVecExclude = cateVarRatioMinMACVecExclude,     cateVarRatioMaxMACVecInclude = cateVarRatioMaxMACVecInclude,     IsSingleVarinGroupTest = opt$IsSingleVarinGroupTest, dosageZerodCutoff = opt$dosageZerodCutoff,     IsOutputPvalueNAinGroupTestforBinary = opt$IsOutputPvalueNAinGroupTestforBinary,     IsAccountforCasecontrolImbalanceinGroupTest = opt$IsAccountforCasecontrolImbalanceinGroupTest,     method = opt$method, kernel = opt$kernel, weights.beta.rare = weights.beta.rare,     weights.beta.common = weights.beta.common, weightMAFcutoff = opt$weightMAFcutoff,     r.corr = opt$r.corr, weightsIncludeinGroupFile = opt$weightsIncludeinGroupFile,     weights_for_G2_cond = weights_for_G2_cond, IsOutputBETASEinBurdenTest = opt$IsOutputBETASEinBurdenTest,     SPAcutoff = opt$SPAcutoff)
An irrecoverable exception occurred. R is aborting now ...
Bus error
weizhouUMICH commented 2 years ago

Sorry for the late reply! We have just released a new version 1.0.0. It has substantial computational efficiency improvements for both Step 1 and Step 2 for single-variant and set-based tests and clearer log output. We have created a new program github page https://github.com/saigegit/SAIGE with the documentation provided https://saigegit.github.io/SAIGE-doc/ The program will be maintained by multiple SAIGE developers there. The docker image has been updated. Please feel free to try the version 1.0.0 and report issues if any.

Thanks! Wei