Calculate Fst to draw Manhattan plot

Hi All,

I trying to calculate Fst. In the first step, I used the following script to calculate Fst with window 100 kb

./vcftools --vcf input_data.vcf --weir-fst-pop population_1.txt --weir-fst-pop population_2.txt --fst-window-size 100000 --fst-window-step 50000 --out pop1_vs_pop2

But in the next step, I do not know how to convert the output file to z-transform? Because before plotting the plot using the Manhattan package in R, z-transform must be performed.

Best Regard

Mostafa

snp r

• 5,457 views • link

updated 4.4 years ago by beausoleilmo

• written 7.9 years ago by mostafarafiepour

Maybe this link is relevant, http://membres-timc.imag.fr/Olivier.Francois/LEA/files/LEA_snmf.html

Next, we convert the FSTFST values into absolute values of zz-scores.

#convert fst values into z-scores (absolute values) 
n = dim(Q(obj.snmf, K = 2))[1]
fst.values[fst.values<0] = 0.000001
K = 2
z.scores = sqrt(fst.values*(n-K)/(1-fst.values))

• 1 views • link

written 7.9 years ago by zx8754

many thanks for your reply,

This is the output of the Fst stage:

enter image description here

I do not know exactly which of these columns should be z-transformed and how?

• 1 views • link

written 7.9 years ago by mostafarafiepour

Can someone explain why the Fst needs to be z-transformed?

• 0 views • link

written 4.4 years ago by beausoleilmo

0 answers

No answers yet.

More posts like this

Error:Could not open individual file
written 4.8 years ago by kk.mahsa

When I ran the following command to Fst calculation vcftools --gzvcf All_samples.vcf.gz --weir-fst-pop Pop1_list.txt --weir-fst-pop Pop2_list.txt --fst-window-size 50000 --fst-window-step 25000 --out pop1_pop2_50_25.windowed.weir.txt I get this …
manhattan plot
written 5.0 years ago by anithanagaraj93 •

I calculated fst values by using following command, "vcftools --vcf a.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --out ex" and got negative and -nan Fst values. if …
Negative Fst values from VCFtools result
written 5.9 years ago by nitinra •

Hello all, I am calculating Fst using VCFtools and here is the code I am running: vcftools --vcf input.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --out pop1_pop2 …
vcftools --fst-window- problems
written 6.0 years ago by florycamacho305 •

Hi all, I'm using vcftools to calculate fst : vcftools --vcf file.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --out pop1_pop2.FST which works. However, when I'm trying to …
vcftools not ouputting log file when run from perl
written 6.7 years ago by dec986

I am running 325 vcftools commands to generate Fst values, which obviously needs to be automated. An example: vcftools --vcf big.vcf --weir-fst-pop pop_lists/pop1.txt --weir-fst-pop pop_lists/pop2.txt …
Using BCFtools to obtain between population Fst values from vcf file
written 6.8 years ago by ishengomae

I'm at a point where I would want to compute between population Weir and Cockerham Fst values directly from a large vcf file. Apparently VCFtools …
How vcftools calculate Fst value?
written 7.5 years ago by seta

Hi all, I used vcftools for calculation the Fst value between two populations (pop1 and pop2) using the below command: ./vcftools --vcf file.vcf --weir-fst-pop population_1.txt …
Pairwise Fst between populations using vcftools
written 10.3 years ago by MC •

Am trying to calculate pairwise Fst between populations with my SNP data in a .vcf file using vcftools-0.1.14. The output I get is a list …
Calculating statistically significant outlier for Pairwise Fst obtained from VCFTools
written 11.8 years ago by Anurag •

Hi, I calculated pairwise Fst using VCFTools: vcftools --vcf input_data.vcf --weir-fst-pop population_1.txt --weir-fst-pop population_2.txt --out pop1_vs_pop2 What method I should use for statistical significance to …
Vcftools Fst Calculation For Single Snps
written 13.7 years ago by klimklam •

<p>Hi there, I would like to obtain the Weir & Cockerham Fst for single SNPs from my vcf files for wished populations. I thought to …