Skip to content

edit loadGenotypeRegion so it filter variants before reading dosages - #613

Open
jaempawi wants to merge 1 commit into
StatFunGen:mainfrom
jaempawi:load-only-kept-variants
Open

jaempawi wants to merge 1 commit into
StatFunGen:mainfrom
jaempawi:load-only-kept-variants

Conversation

@jaempawi

Copy link
Copy Markdown

keepVariantsPath and keepIndel were applied after the dosage matrix for the whole region was read, so asking for a few variants still loaded every genotype in the region (300 chr22 variants over the full chromosome tried to allocate 14.8 GB). The filters are now applied to the variant index before the read, so only the kept variants are loaded. Output is unchanged; a new test checks that the kept dosages equal a full load followed by subsetting.

…nt index before reading dosages, so only the requested variants are loaded

This branch has not been deployed

No deployments
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant