Stephan Schiffels

Yes, that would be super useful for modelling. I might propose something if I find the time...

How about `skip`?

This sounds all great. Regarding aDNA data, just my five cents: This is so high-dimensional (a million genetic markers aren't uncommon) that it wouldn't fit into the exact same framework...

Yes, good point. Haplogroups could go into this, for sure, and they are already quite interesting. And of course, if we could even automatically download the full data somehow through...

Same here, packunused marks lots and lots of packages as redundant, when they're really not. Just a random thought: Does it properly search also through the additional files imported by...

Are you sure you have a correct mask file here? It seems you're giving it a VCF instead of a BED file for the mask.

could you post the error message, not just the Traceback, please?

Well, if the error occurred with the command `chromosome = fields[0]`, it means that you have an empty line there that you're trying to parse. Your input files must be...

This means that one of your VCFs is not in the right shape. My program expects the genotypes (e.g. `0|1`) in the 10th columns (so index 9), and that doesn't...

As above, please check that your VCF file is the shape that my scripts expects it. See my previous comment.