Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Microarrays and missing data, --chip #32

Open
stephenturner opened this issue Dec 23, 2019 · 0 comments
Open

Microarrays and missing data, --chip #32

stephenturner opened this issue Dec 23, 2019 · 0 comments

Comments

@stephenturner
Copy link
Collaborator

As described in the README, the --chip flag will limit the range to array SNPs only. How are missing genotypes handled in a multisample VCF file? For example, my multisample VCF file contains positions 1, 2, 3, 4, 5, but one particular sample is missing data (genotype not called) for positions 2 and 4. It appears in the results the range is still listed as 1, 2, 3, 4, 5 for all samples, even the one that was missing data. Is this the expected behavior? Desired behavior? I could limit the range to 1, 3, 5 for this particular sample by writing a new VCF extracting just this sample and removing missing genotypes from the VCF altogether, but that wouldn't allow me to analyze multiple samples at once, each missing a particular subset of the 5 positions on my hypothetical array. cc @vpnagraj @cneal13

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

No branches or pull requests

1 participant