Question: update fam file VCF -> plink bfile
gravatar for
13 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 • 838 views
ADD COMMENTlink modified 13 months ago by chrchang5234.8k • written 13 months ago by biosigil.blog0
gravatar for chrchang523
13 months ago by
United States
chrchang5234.8k wrote:

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

ADD COMMENTlink written 13 months ago by chrchang5234.8k
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: 759 users visited in the last hour