EvalAdmix: Difference between revisions
No edit summary |
|||
Line 14: | Line 14: | ||
evalAdmix can be installed from [https://github.com/GenisGE/evalAdmix github] | evalAdmix can be installed from [https://github.com/GenisGE/evalAdmix github] | ||
< | <pre> | ||
git clone https://github.com/GenisGE/evalAdmix.git | git clone https://github.com/GenisGE/evalAdmix.git | ||
cd evalAdmix | cd evalAdmix | ||
make | make | ||
</ | </pre> | ||
==Quick start== | ==Quick start== | ||
Line 52: | Line 52: | ||
===Plink=== | ===Plink=== | ||
Genotype data files in binary PLINK format (.bed .fam .bim). | |||
===Beagle genotype likelhoods=== | ===Beagle genotype likelhoods=== | ||
The input file contains genotype likelihoods in a .beagle file format [http://faculty.washington.edu/browning/beagle/beagle.html]. | The input file contains genotype likelihoods in a .beagle file format [http://faculty.washington.edu/browning/beagle/beagle.html]. | ||
Line 72: | Line 71: | ||
==Output Files== | ==Output Files== | ||
The analysis performed by | The analysis performed by evalAdmix produces 1 correlation file: | ||
==Run command example== | ==Run command example== | ||
=== low depth sequencing data === | |||
Download the input file containing genotype in binary plink format | |||
<pre> | |||
wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.bed | |||
wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.bim | |||
wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.fam | |||
</pre> | |||
Run ADMIXTURE to obtain admixture proprotions | |||
<pre>admixture admixTjeck2.bed 3</pre> | |||
Run evalAdmix | |||
<pre>./evalAdmix -plink admixTjeck2 -fname admixTjeck2.3.P -qname admixTjeck2.3.Q -P 20 </pre> | |||
=== low depth sequencing data === | === low depth sequencing data === | ||
Download the input file containing genotype likelihoods in beagle format | Download the input file containing genotype likelihoods in beagle format |
Revision as of 14:04, 19 July 2019
evalAdmix allows to evaluate the results of an admixture analysis. It takes as input the genotype data (either called genotypes in plink files or genotype likelihoods beagle files) used in the admixture analysis and the frequency and admixture propotions (P and Q files) generated.
The output is a pairwise correlation of residuals matrix between individuals The correlation will be 0 in case of a good fit of the data to the admixture model. When something is wrong, individuals from the same population will be positively correlated; and individuals from different populationts but that share one or more ancestral populations as admixture sources will have a negative correlation. Positive correlation between a pair of individuals might also be due to relatedness.
Download and Installation
evalAdmix can be installed from github
git clone https://github.com/GenisGE/evalAdmix.git cd evalAdmix make
Quick start
./evalAdmix -beagle inputBeagleFile.gz -fname myoutfiles.fopt.gz -qname myoutfiles.qopt -P 10
./evalAdmix -plink inputPlinkPrefix -fname myoutfiles.fopt.gz -qname myoutfiles.qopt -P 10
- -bealge beagle file of genotype likelihoods
- -fname number of clusters
- -o prefix of output file names
- -P Number of threads used
Parameters
./evalAdmix
Arguments: -plink path to binary plink file (excluding the .bed) -beagle path to beagle file containing genotype likelihoods (alternative to -plink) -fname path to ancestral population frequencies file -qname path to admixture proportions file -o name of the output file Setup: -P 1 number of threads -autosomeMax 23 autosome ends with this chromsome -nIts 5 number of iterations to do for frequency correction; if set to 0 calculates correlation without correction (fast but biased) -useSites 1.0 proportion of sites to use to calculate correlation of residuals -misTol 0.05 tolerance for considering site as missing when using genotype likelihoods. Use same value as used in NGSadmix to keep compatibility when using genotype likelihoods (-beagle) -minMaf 0.05 minimum minor allele frequency to keep site. Use same value as used in NGSadmix to keep compatibility when using genotype likelihoods (-beagle)
Input File
Plink
Genotype data files in binary PLINK format (.bed .fam .bim).
Beagle genotype likelhoods
The input file contains genotype likelihoods in a .beagle file format [1]. and can be compressed with gzip.
BAM files
If you have BAM files you can use ANGSD to produce genotype likelihoods in .beagle format. Please see Creation of Beagle files with ANGSD
VCF files
If you already have made a VCF file that contains genotype likehood information then it should be possible to convert .vcf files with genotype likelihoods to .beagle file via vcftools [2]
vcftools --vcf input.vcf --out test --BEAGLE-GL --chr 1,2
Chromosome has to be specified.
You can also use bcftools' [3] 'query' option for generating a .beagle file from a .vcf file.
Output Files
The analysis performed by evalAdmix produces 1 correlation file:
Run command example
low depth sequencing data
Download the input file containing genotype in binary plink format
wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.bed wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.bim wget http://pontus.popgen.dk/albrecht/open/admixTjeck/plink/admixTjeck2.fam
Run ADMIXTURE to obtain admixture proprotions
admixture admixTjeck2.bed 3
Run evalAdmix
./evalAdmix -plink admixTjeck2 -fname admixTjeck2.3.P -qname admixTjeck2.3.Q -P 20
low depth sequencing data
Download the input file containing genotype likelihoods in beagle format
wget popgen.dk/software/download/NGSadmix/data/Demo2input.gz wget popgen.dk/software/download/NGSadmix/data/Demo2pop.info
Execute NGSadmix to obtain admixture proportions
./NGSadmix -likes Demo2input.gz -K 3 -P 20 -o myoutfiles -minMaf 0.05
- Input file = Demo2input.gz
- Ancestral Populations K=3
- Computer cores = 20 (-P 20).
- Output prefix = myoutfiles (-o myoutfiles)
- SNPs with MAF > 5% (-minMaf 0.05)
Run evalAdmix
./evalAdmix -beagle Demo2input.gz -fname myoutfiles.fopt.gz -qname myoutfiles.qopt -P 20
- Input file = input.gz
- Ancestral Populations frequency file myoutfiles.fopt.gz
- Computer cores = 20 (-P 20).
Plot results in R
# read population labels and estimated admixture proportions pop<-read.table("Demo2pop.info",as.is=T) q<-read.table("myoutfiles.qopt") # order according to population and plot the NGSadmix reults ord<-order(pop[,1]) barplot(t(q)[,ord],col=2:10,space=0,border=NA,xlab="Individuals",ylab="Demo2 Admixture proportions for K=3") text(tapply(1:nrow(pop),pop[ord,1],mean),-0.05,unique(pop[ord,1]),xpd=T) abline(v=cumsum(sapply(unique(pop[ord,1]),function(x){sum(pop[ord,1]==x)})),col=1,lwd=1.2) r<-read.table("output.corres.txt") image(as.matrix(r))