Question: update fam file VCF -> plink bfile
0
gravatar for biosigil.blog
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
0
gravatar for chrchang523
13 months ago by
chrchang5234.8k
United States
chrchang5234.8k wrote:

--update-parents (http://www.cog-genomics.org/plink/1.9/data#update_indiv ), --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.

Help
Access

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