Segmentation fault when reading BGEN file

[jielab]

Dear Xiaowei:

When I use the latest version of rvtest and the BGEN input file, why the log says "Loaded [ 0 ] samples from VCF files"? There is NO VCF file at all. Also, the BGEN file has two part, one is .bgen and the other part is .sample. So, how do I specify the .sample file?

In the phenotype file, does my header line must have the following five as the header names: "fid iid motid fatid sex"? It seems that the command options in RVTEST is very similar to PLINK, such as --covar-name --pheno-name. But in PLINK, the phenotype file only needs the first two columns "FID IID" and then the rest could be anything.

Thank you very much & best regards, jie

Below is the LOG file:

Thank you for using rvtests (version: 20170818, git: 72ab2236f3be72afc5329dcea9a770747aa13330) .... .... Effective Options --inBgen ukb_imp_chr20_v2.bgen --out chr20.rvtest --pheno bmi.EUR.ped --pheno-name bmi --single score

[INFO] Program version: 20170818 [INFO] Analysis started at: Wed Sep 13 21:46:39 2017 [INFO] Loaded [ 0 ] samples from VCF files [INFO] Loaded [ 334348 ] sample pheontypes Segmentation fault

[dajiangliu]

Did you try --meta score option? --single score is less actively maintained for now.

[jielab]

still got the same error. I think the segmentation fault error came from reading the input file, RVTESTS has not gone to --meta or --single part yet. I still don't understand why the log file says "from VCF files" even when i specify it to read a .bgen file. And i don't know where to specify the .sample file.

BTW, I used PLINK2a to read this .bgen file, and there is no problem.

best regards, Jie

[zhanxw]

Hi Jie,

I will look into this issue.

Our initial goal is to use bgen files as genotype inputs and ped files as phenotype inputs. I will see how to incorporate .sample file.

Xiaowei

On Sep 13, 2017, at 10:01 PM, jiehuang001 [email protected] wrote:

still got the same error. I think the segmentation fault error came from reading the input file, RVTESTS has not gone to --meta or --single part yet. I still don't understand why the log file says "from VCF files" even when i specify it to read a .bgen file. And i don't know where to specify the .sample file.

BTW, I used PLINK2a to read this .bgen file, and there is no problem.

best regards, Jie

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or mute the thread.

[zhanxw]

The current version of RVTESTS requires that .bgen file includes sample names. In that case a separate .sample file is not be needed. I just learned that many other .bgen files do not include sample names. In that case, I need to make RVTESTS support loading .sample file. I will work that, and update this issue.

[jielab]

I see. i will not use RVTESTS for .bgen files based association analyses, for now then. Hopefully you will release a new version soon that address the issues that I raised.

Best regards, Jie

[zhanxw]

@jiehuang001 This version (https://github.com/zhanxw/rvtests/releases/tag/v2.0.6) should support BGEN format. Please let me know if you have additional questions. Thanks.