This is a test version of Biostars. For the public version, visit https://www.biostars.org.
vcflib-popstats to calculate FIS of a subpopulation

Dear Biostars

I used the vcflib-popstats to calculate FIS of a subpopulation but I have same FIS number for positions with the different number of hets (4th column). also, which formula did it use to the calculation Fis?

would you have any idea?

Thanks 

 target_allele_frequency    expected_heterozygosity observed_heterozygosity number_of_hets  number_of_homozygous_ref    number_of_homozygous_alt    target_Fis
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   27  0   0   1.00E-05
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   29  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   30  0   0   1.00E-05
    0.5 0.5 1   25  0   0   1.00E-05
    0.5 0.5 1   19  0   0   1.00E-05
    0.5 0.5 1   28  0   0   1.00E-05
vcflib popstats
• 1,906 views • link
updated by Zev.Kronenberg Zev.Kronenberg • written by Sam Sam

1 answer

Hi Sam,

The genotype object, stored in var.cpp, contains the methods to calculate various statistics. I used the equation below to calculate FIS, where nhet is the count of the heterozygous genotypes, and ngeno is the number of called genotypes, af, is the allele frequency. 

  fis = ( 1 - ((nhet/ngeno) / (2*af*(1 - af))));
• 0 views • link
written by Zev.Kronenberg Zev.Kronenberg

Log in to answer this question.

More posts like this
  • Heatmap legend position change
    written by hellokwmin • hellokwmin
    Hello, This is code that I am using to generate a heatmap: ```r color_fun <- colorRamp2(c(-4, 0, 4), c("red", "black", "green")) df = read.csv("RAT7basedlog2version2.csv", header …
  • LInear Models with HLA-B alleles in R
    written by schagas • schagas
    Hello, I think this question might be more a statistics understanding than a code troubleshooting. Still, I am posting here because there are more R-users …
  • GEO Query error
    written by Ondina • Ondina
    Hello everyone, I want to use GEO Query but I'm getting an error and I'm stuck. This is the code: gse <- getGEO("GSE21653", GSEMatrix = …
  • Missing data per site
    written by anna • anna
    Hi, I want to calculate statistics of missing data per each site in my vcf file. Using `vcftools --missing-site` gives wrong stats for several sites. …
  • How to sum up value based on previous rows in R
    written by ishackm ishackm
    Hi all, I have this dataset. Value Date Type Difference Increase or Decrease 1 55.55556 2018-09-21 ESVI 0.0000000 0 2 51.11111 2018-10-03 ESVI -4.4444444 -1 …
  • Runs of homozygosity calculation
    written by BAGeno BAGeno
    Hi, I am using plink to calculate runs of homozygosity.I prepared plink format file with vcftools I ran following command plink --file Plink.input --homozyg --out …
  • Counting nucleotide frequency using perl script
    written by MAPK MAPK
    I have this perl script below to calculate sequence length and their frequency along with nucleotide frequency(`A`,`T`,`G` and `C`). This script works fine for a …
  • vcf-annotate not adding annotations
    written by wes3985 • wes3985
    I would like to annotate a vcf file with the output of a filter that detects false positives. I know a very similar question to …
  • Index Vcf File
    written by mad.cichlids mad.cichlids
    <p>Hi, I was trying to index my vcf file, i first sort it and then zip it before indexing as this post: <a href='/p/78444/'>Tabix -p …
  • Extracting Set Of Numbers From A File
    written by Cynthia • Cynthia
    <p>In a text file there are numbers. for eg:</p> <blockquote> <p>human 0 1 2 3 4 5 6 7 8 9 10 11 12 13 …