Question: Calculate the genotype frequency
0
gravatar for SOHAIL
11 months ago by
SOHAIL260
Beijing Institute of Genomics, CAS.
SOHAIL260 wrote:

Hi everyone:

I have a VCF file having multiple population samples. Is there any way to calculate the observed counts of genotypes called and observed genotype frequency of a subset of individuals from a given population within VCF file?

For example (genotype freq) : 0|0 : 0.2, 0|1: 0.5, 1|1: 0.3

just an analogous to --freq option in VCFtools field that calculates the allele frequency of REF and ALT allelic forms at given site.

Thanks a lot!

Regards, sohail.

ngs vcf • 798 views
ADD COMMENTlink modified 11 months ago by Wietje190 • written 11 months ago by SOHAIL260
0
gravatar for Wietje
11 months ago by
Wietje190
Germany
Wietje190 wrote:

You might want to try PLINK which allows you to read in a VCF file (--vcf [INPUTFILE]), select only a few individuals (--keep [FILE]) and then get genotype frequencies with --freqx (with aditional option --nonfounders if you want to have it for all individuals in your list) or allele frequencis with --freq. Check the manual for more detailed information:

PLINK 1.9: https://www.cog-genomics.org/plink/1.9/

PLINK 2.0: https://www.cog-genomics.org/plink/2.0/

I hope this is of some help to you!

ADD COMMENTlink written 11 months ago by Wietje190
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: 1167 users visited in the last hour