How to run BOLT-LMM against a Plink simulated dataset?
0
0
Entering edit mode
4.0 years ago
b.ambrozio ▴ 30

I'm trying run BOLT-LMM against a simulated dataset, which were generated with plink --simulate, however no success:

./BOLT-LMM_v2.3.4/bolt --bfile=ds1_1gt10k --phenoUseFam --statsFile ds1_1gt10k.bolt --LDscoresUseChip --lmm
...
ERROR: LOCO mixed model analysis requires >= 2 chromosomes or chunks
       To run this analysis, set --numLeaveOutChunks

The documentation don't say nothing about numLeaveOutChunks parameter. I've tried some arbitrary numbers as well, and things started working, however, after some time running, it ended up with another error:

./BOLT-LMM_v2.3.4/bolt --bfile=ds1_1gt10k --phenoUseFam --statsFile ds1_1gt10k.bolt --LDscoresUseChip --lmm --numLeaveOutChunks=5
...
Intel MKL ERROR: Parameter 6 was incorrect on entry to DGELS .

Intel MKL ERROR: Parameter 6 was incorrect on entry to DGELS .

Intel MKL ERROR: Parameter 6 was incorrect on entry to DGELS .
Intercept of LD Score regression for ref stats:   0.000 (0.000)
Estimated attenuation: -nan (nan)
Intercept of LD Score regression for cur stats: 0.000 (0.000)
Calibration factor (ref/cur) to multiply by:      -nan (nan)

Time for computing Bayesian mixed model assoc stats = 0.802051 sec

Calibration stats: mean and lambdaGC (over SNPs used in GRM)
  (note that both should be >1 because of polygenicity)
Segmentation fault
BOLT-LMM PLINK • 1.8k views
ADD COMMENT

Login before adding your answer.

Traffic: 2276 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6