Sib-pair Command: read vcf


ClassData Declaration command
Nameread 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)