Tag: bgen

How can I keep INFO value when convert bgen to VCF by using plink2?

How can I keep INFO value when convert bgen to VCF by using plink2? 1 I am working on file handling for GWAS. When I converted bgen to VCF by using plink2 with a commands below, all INFO (and also FILTER) columns became “.” in the output VCF files. A…

Continue Reading How can I keep INFO value when convert bgen to VCF by using plink2?

Using QCTOOL v2 to process UK Biobank .bgen files

Using QCTOOL v2 to process UK Biobank .bgen files – why so slow? 0 I’m currently using QCTOOL v2 to process imputed .bgen files from UK Biobank, however they seem to be processing very slowly. Is this normal? My command is pretty basic; I’m filtering out a list of SNPs…

Continue Reading Using QCTOOL v2 to process UK Biobank .bgen files

qctool to merge two bgen file fails with no clear reason to

Hi, I am trying to merge two bgen files using qctool as explained here. I am using qctool_v2.2.0. The command works but ends with an error: ❱ qctool -g bug/in2.bgen -s bug/in2.sample -merge-in bug/in1.bgen bug/in1.sample -og bla.bgen -os bla.sample Welcome to qctool (version: 2.2.0, revision: unknown) (C) 2009-2020 University of…

Continue Reading qctool to merge two bgen file fails with no clear reason to

P-values far too high for quantitative regenie phenotype

P-values far too high for quantitative regenie phenotype 0 Hi all, I’m having some trouble running regenie (v2.2.4) on a quantitative phenotype for a large cohort. I’m testing a standard height GWAS with heights rounded to the nearest integer. I’ve tried a few different tests to see where the issue…

Continue Reading P-values far too high for quantitative regenie phenotype

vcf to bgen conversion using qctool v2 yields 0 snps

Hi all, I have a vcf file that was extracted from UKB data using qctool (v2.0.6-Ubuntu16.04-x86_64) and contains data in the GP format. This contains a bunch of SNPs from a single chromosome. ❱ wc -l chromosome1.vcf 260 chromosome1.vcf Then I try to convert this file to .bgen again using…

Continue Reading vcf to bgen conversion using qctool v2 yields 0 snps

Why may BOLT-LMM and SAIGE (quantitative, linear-mixed model) yield different results when ran on the absolutely the same dataset?

As a validation experiment, I have run the same GWAS of a quantitative phenotype derived from the UKBiobank, alongside the genomic data from the UKBiobank, once using the program BOLT-LMM and once using SAIGE linear mixed model (with selected quantitative trait tag). I wanted to see if the results would…

Continue Reading Why may BOLT-LMM and SAIGE (quantitative, linear-mixed model) yield different results when ran on the absolutely the same dataset?

Calculate LD matrix from bgen file

formatting error: Calculate LD matrix from bgen file 1 Hello, I am new to plink and am learning as I go. I am trying to calculate an LD matrix for a list of variants while using a bgen file as my reference population. See the command below: ./plink2/plink2 –r2 bin…

Continue Reading Calculate LD matrix from bgen file

extract list of SNPs from multiple chr{1:22}.bgen files using plink2

extract list of SNPs from multiple chr{1:22}.bgen files using plink2 1 hello, I have extracted out list of snps based on the maf cutoff 0,,0.0001, 0.001,0.01,0.1,.55,1.0. I am running plink2 to extract this list from .bgen files for individual chromosomes using the following code plink2 –chr{1:22}.bgen –extract maf1_snps for imputed…

Continue Reading extract list of SNPs from multiple chr{1:22}.bgen files using plink2