Question: update fam file VCF -> plink bfile
gravatar for
24 months ago by
biosigil.blog0 wrote:

I have a large VCF (WGS >7000 samples) that I converted to plink bfile with this command

plink --vcf in.vcf --make-bed --out plink

A fam file is produced, but does not reflect the relationships, phenotypes, or genders

I've replaced the fam file with a correct one containing these information. Then I calculate allele frequency with this command

plink --bfile plink --freq counts --out plink

I've noticed that the allele counts are incorrect, like for variant A plink reports an allele count of 1 and when I manually parse the VCF, I find the allele count is 3 (for founders only)

So something went wrong.

How do I convert a VCF into a plink bfile while replacing/updating the fam file, so I don't get these errors.

bfile plink vcf • 1.4k views
ADD COMMENTlink modified 24 months ago by chrchang5236.3k • written 24 months ago by biosigil.blog0
gravatar for chrchang523
24 months ago by
United States
chrchang5236.3k wrote:

--update-parents ( ), --update-sex, and --pheno can be used for this purpose.

ADD COMMENTlink written 24 months ago by chrchang5236.3k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1024 users visited in the last hour