How does plink read bed files?
plink bed files uses 2 bits for each genotype data point, which is most efficient in terms of storage space, but not so convenient for numerical analysis. I am wondering how does plink read these bits into a matrix of int/double for linear algebra operations?
• 4,274 views
•
link
updated
by
Ram
•
written
by
kindlychung •
•
written
by
kindlychung •
1 answer
The short answer is that it usually doesn't. Whenever possible, PLINK 1.9 never unpacks the data; instead it uses bitwise operations and population count to perform computation directly on the 2-bit representation.
• 0 views
•
link
updated
by
Ram
•
written
by
chrchang523
•
written
by
chrchang523
Log in to answer this question.
More posts like this
-
Efficient way to merge multiple single sample vcf files with same variants / load into PLINK
written by dmitry.s.kolobkov •
-
Parallel in PLINK for linear association for SNP effects
written by landscape95
-
Using pysam to parse bamfile and fetch based on bedfile, looking to speed up code.
written by Daniel
-
How to generate vcf using samtools Version: 0.1.19-96b5f2294a?
written by lakhujanivijay
-
plink mergelist put many individuals into same binary file set
written by gulcek •
-
defining reference alleles in 23andMe file to create plink .ped file
written by gulcek •
-
Is There An Explaination For Each Of The Pieces Of Blast Output?
written by agatorano •
-
How Can I Apply A Pssm Efficiently?
written by glocke01
-
Representing Interaction Data Using Iranges/Granges For 5C/Hi-C Data Analysis
written by Naveen •
-
Metabolic Network Drawing
written by Kwet •
Have you gone through the source code?
I am trying to, but I am not that good at C. Could you tell what functions are involved so that I can focus on these?
I'd search the source code for
fopenandfread, since those will give you a clue where to start.