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.


read loc vcf simple.vcf.gz
sum get R2
keep where test >= 0.8
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)