saigegit / SAIGE

Development for SAIGE and SAIGE-GENE(+)
GNU General Public License v3.0
65 stars 28 forks source link

Error during conditional analysis #129

Open mohasal12 opened 10 months ago

mohasal12 commented 10 months ago

I am getting the following error

Error: vector::_M_range_check: __n (which is 2538) >= this->size() (which is 2538)
Execution halted

when we reach the chunck number 41, without conditioning the code works fine.

My code

/usr/bin/time -v Rscript ${saigepath}step2_SPAtests.R \
  --vcfFile=${vcfpath}CHR6_PID112437_colon.vcf.gz \
  --vcfFileIndex=${vcfpath}CHR6_PID112437_colon.vcf.gz.csi \
  --vcfField=DS \
  --LOCO=FALSE \
  --chrom=6 \
  --minMAF=0 \
  --minMAC=3 \
  --is_fastTest=TRUE \
  --GMMATmodelFile=${outpath}out_s1_0.125_2000.rda \
  --varianceRatioFile=${outpath}out_s1_0.125_2000.varianceRatio.txt \
  --SAIGEOutputFile=${outpath}out_6:32651477:A:T.txt \
  --is_output_moreDetails=TRUE \
  --sparseGRMFile=${outpath}sparseGRM_relatednessCutoff_0.125_2000_randomMarkersUsed.sparseGRM.mtx \
  --sparseGRMSampleIDFile=${outpath}sparseGRM_relatednessCutoff_0.125_2000_randomMarkersUsed.sparseGRM.mtx.sampleIDs.txt \
  --relatednessCutoff=0.125 \
  --condition=6:32651477:A:T

Could you provide any feedback on this issue.