Class | Data Declaration command |
Name | read vcf |
Arguments | <data file name> ped_id |
Reads only the IDs from a VCF file header and sets up the corresponding pedigree in memory. If the ped_id modifier is present, the ID is split to give a pedigree and individual ID. Neither the marker names nor genotype data are read from the VCF file.
Example:
clear read loc vcf simple.vcf.gz sum get R2 keep where test >= 0.8 pack read vcf simple.vcf.gz merge vcf simple.vcf.gz
See also:
merge vcf | merge VCF genotype data |
<< (read cases) | Up to index | >> (read hapmap) |